Section courante

A propos

Section administrative du site

Une des fonctions les plus communes de la géographie et des systèmes modernes, c'est le calcul de la distance géographique entre deux coordonnées de Longitude et de Latitude. Il n'y a aucune nécessité de grande connaissance en trigonométrie pour arriver à se genre de calcul dans le format qu'on le souhaite, Km, Miles ou Miles Nautiques. Ainsi, si vous savez les coordonnées suivantes :

Ville Latitude Longitude
Montréal 45 31N 73 34O
Paris 48 50N 2 20E

A l'aide du code source Visual C++ .NET 2003 suivant, vous trouverez la réponse que vous souhaitez :

  1. #include "stdafx.h"
  2. #include "iostream"
  3. #include "NxMath.h"
  4.  
  5. double CoordToDeltaKm(
  6.     double Q1Latitude,double Q1LatiDeg, char Q1LatiDirection,
  7.     double Q1Longitude,double Q1LongDeg, char Q1LongDirection,
  8.     double Q2Latitude,double Q2LatiDeg, char Q2LatiDirection,
  9.     double Q2Longitude,double Q2LongDeg, char Q2LongDirection
  10. ) {
  11.     double a1,b1,a2,b2,RawDelta;
  12.     a1=(Q1Latitude+(Q1LatiDeg/60)) * NxPiF64/180;
  13.     if(Q1LatiDirection=='N') a1=-a1;
  14.     b1=(Q1Longitude+(Q1LongDeg/60)) * NxPiF64/180;
  15.     if(Q1LongDirection=='O') b1=-b1;
  16.     a2=(Q2Latitude+(Q2LatiDeg/60)) * NxPiF64/180;
  17.     if(Q2LatiDirection=='N') a2=-a2;
  18.     b2=(Q2Longitude+(Q2LongDeg/60)) * NxPiF64/180;
  19.     if(Q2LongDirection=='O') b2=-b2;
  20.      RawDelta = NxMath::acos(NxMath::cos(a1) * NxMath::cos(b1) * NxMath::cos(a2) * NxMath::cos(b2) + NxMath::cos(a1)*NxMath::sin(b1)*NxMath::cos(a2)*NxMath::sin(b2) + NxMath::sin(a1)*NxMath::sin(a2));
  21.     return RawDelta * 6378.0;
  22. }
  23.     
  24. double CoordToDeltaStatuteMiles(
  25.     double Q1Latitude,double Q1LatiDeg, char Q1LatiDirection,
  26.     double Q1Longitude,double Q1LongDeg, char Q1LongDirection,
  27.     double Q2Latitude,double Q2LatiDeg, char Q2LatiDirection,
  28.     double Q2Longitude,double Q2LongDeg, char Q2LongDirection
  29. ) {
  30.     double a1,b1,a2,b2,RawDelta;
  31.     a1=(Q1Latitude+(Q1LatiDeg/60)) * NxPiF64/180;
  32.     if(Q1LatiDirection=='N') a1=-a1;
  33.     b1=(Q1Longitude+(Q1LongDeg/60)) * NxPiF64/180;
  34.     if(Q1LongDirection=='O') b1=-b1;
  35.     a2=(Q2Latitude+(Q2LatiDeg/60)) * NxPiF64/180;
  36.     if(Q2LatiDirection=='N') a2=-a2;
  37.     b2=(Q2Longitude+(Q2LongDeg/60)) * NxPiF64/180;
  38.     if(Q2LongDirection=='O') b2=-b2;
  39.      RawDelta = NxMath::acos(NxMath::cos(a1)*NxMath::cos(b1)*NxMath::cos(a2)*NxMath::cos(b2) + NxMath::cos(a1)*NxMath::sin(b1)*NxMath::cos(a2)*NxMath::sin(b2) + NxMath::sin(a1)*NxMath::sin(a2));
  40.     return RawDelta * 3963.1;
  41. }
  42.     
  43. double  CoordToDeltaNauticalMiles(
  44.     double Q1Latitude,double Q1LatiDeg, char Q1LatiDirection,
  45.     double Q1Longitude,double Q1LongDeg, char Q1LongDirection,
  46.     double Q2Latitude,double Q2LatiDeg, char Q2LatiDirection,
  47.     double Q2Longitude,double Q2LongDeg, char Q2LongDirection
  48. ) {
  49.     double a1,b1,a2,b2,RawDelta;
  50.     a1=(Q1Latitude+(Q1LatiDeg/60)) * NxPiF64/180;
  51.     if(Q1LatiDirection=='N') a1=-a1;
  52.     b1=(Q1Longitude+(Q1LongDeg/60)) * NxPiF64/180;
  53.     if(Q1LongDirection=='O') b1=-b1;
  54.     a2=(Q2Latitude+(Q2LatiDeg/60)) * NxPiF64/180;
  55.     if(Q2LatiDirection=='N') a2=-a2;
  56.     b2=(Q2Longitude+(Q2LongDeg/60)) * NxPiF64/180;
  57.     if(Q2LongDirection=='O') b2=-b2;
  58.      RawDelta = NxMath::acos(NxMath::cos(a1)*NxMath::cos(b1)*NxMath::cos(a2)*NxMath::cos(b2) + NxMath::cos(a1)*NxMath::sin(b1)*NxMath::cos(a2)*NxMath::sin(b2) + NxMath::sin(a1)*NxMath::sin(a2));
  59.     return RawDelta * 3443.9;
  60. }
  61.      
  62. int _tmain(int argc, _TCHAR* argv[])
  63. {
  64.      std::cout << "Distance entre Montréal et Paris en Km: " <<
  65.                     CoordToDeltaKm(45, 31,'N', 73, 34,'O',    48, 50,'N', 2,  20,'E') << std::endl;
  66.      std::cout << "Distance entre Montréal et Paris en Miles: " <<
  67.                     CoordToDeltaStatuteMiles(45, 31,'N', 73, 34,'O',    48, 50,'N', 2,  20,'E') << std::endl;
  68.      std::cout << "Distance entre Montréal et Paris en Miles Nautique: " <<
  69.                     CoordToDeltaNauticalMiles(45, 31,'N', 73, 34,'O',    48, 50,'N', 2,  20,'E') << std::endl;
  70.      return 0;
  71. }

on obtiendra le résultat suivant :

Distance entre MontrÚal et Paris en Km: 5510.17
Distance entre MontrÚal et Paris en Miles: 3423.85
Distance entre MontrÚal et Paris en Miles Nautique: 2975.3

Dernière mise à jour : Vendredi, le 2 août 2013