Le portail francophone de la géomatique

Bonjour,

252.6395 Km (selon un calcul effectue par GpsUtility... dont je ne connais
pas la validite....)

Sinon, utilisez la formule diffusee sur la liste gps@texte-a-enlever.eurenet.com en novembre
dernier... Bon courage ! ;-)))

On pose :
- Gi : longitude Est du point i ;
- Li : latitude Nord du point i ;
- Hi : altitude du point i ;
- A : rayon moyen de la Terre ;

On supposera la Terre spherique : l'approximation est suffisante pour des
comparaisons de distance. Un calcul vraiment rigoureux supposerait de passer
par l'ellipsoide, voire le geoide, avec utilisation des integrales
elliptiques.

La direction geocentrique du point initial i est :
Zi = {cos(Gi)*cos(Li) ; sin(Gi*cos(Li) ; sin(Li) }

Le plan tangent a la sphere terrestre en i est defini par les deux vecteurs
:
Yi = { -sin(Gi) ; cos(Gi) ; 0 }
Xi = {-cos(Gi)*sin(Li) ; -sin(Gi)*sin(Li) ; cos(Li) }

Xi est oriente vers le Nord et Yi vers l'Est.

Le cosinus de la distance entre deux points U et V est donnee par le produit
scalaire Zu.Zv, soit :
D = cos(Gu-Gv)*cos(Lu)*cos(Lv)  sin(Lu)*sin(Lv) ;

Si l'on est familier du mille marin, on prend l'arc inverse de ce cosinus en
minutes de degres, sinon on prend cet arc en radians et on multiplie par le
rayon terrestre, eventuellement corrige de l'altitude moyenne.

Les vecteurs X et Y permettent d'obtenir le cap a suivre C pour aller de U
vers V, soit :
K*cos(C) = cos(Lu)*sin(Lv) - sin(Lu)*cos(Lv)*cos(Gv-Gu)
K*sin(C) = cos(Lv)*sin(Gv-Gu)

Avec : K*K D*D=1 et K >0, le point le plus distant etant aux antipodes.

Les signes de cos(C) et sin(C) permettent de trouver le bon quadrant (NE,
SE, SW ou NW), et leur rapport la valeur du cap vrai geographique.

Rahhhhhhhhhh ! ;-)))

Cordialement,

Sylvain THERY - CSN geographe IRD (ex-Orstom)

CERMES-OCCGE
Unite de Parasitologie et de Transmission des Schistosomes
B.P. 10887 Niamey - NIGER
http://www.ird.ne/cermes