Nous utilisons des cookies pour vous garantir la meilleure expérience sur notre site. Si vous continuez à utiliser ce dernier, nous considèrerons que vous acceptez l'utilisation des cookies. J'ai compris ! ou En savoir plus !.
banniere

Le portail francophone de la géomatique


Toujours pas inscrit ? Mot de passe oublié ?
Nom d'utilisateur    Mot de passe              Toujours pas inscrit ?   Mot de passe oublié ?

Annonce

Suite à un problème technique intervenu entre le 22 et le 23 mars, nous avons du procéder dans la soirée du 25 mars, à la restauration de la base de données du 24 mars (matinée).

En clair, nous avons perdu vos contributions et inscriptions du dimanche 24 et du lundi 25 mars.
Nous vous prions de nous excuser.

#1 Fri 19 May 2000 12:32

Sylvain Théry
Invité

2°20' 14,025 -> metres

Bonjour,

Ca ressemble a un probleme d'origine. Le decallage quand on est dans une
projections geographique classique est de 2° 20' 14,025, mais dans un
repere kilometrique, savez-vous combien ca fait ?


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

 

Pied de page des forums

Powered by FluxBB