Discussion :

[lhbilmarocain]

Bonjour tout le mode je suis bloqué sur une partie,
C'est que je veux convertir la distance exprimée en Latitude & Longitude en Km ou en mètre !

J'ai un point de repère que je fais entrer (lat1, lon1)

J'ai un fichier que je lis sur matlab comprenant des centaines de coordonnées exprimées en (lat2, lon2) aussi.

Alors à partir du premier point je veux faire une mesure de distance entre les différents points sur le fichier puis convertir le résultat en KM ou en M


Merci de votre aide.

[christophe_halgand]

Bonjour,

Une soustraction devrait répondre à ta première question.
Tu devrais trouver ta seconde réponse ici.

Cordialement,
Christophe Halgand

[lhbilmarocain]

Merci pour votre réponse,

Sinon, je reviens pour proposer ma solution, comme ça tout le monde en profitera :-)

Code MATLAB :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
function [a,c,dlat,dlon]=haversine11(lat1,lon1,lat2,lon2)
    dlat = radians(lat2-lat1);
    dlon = radians(lon2-lon1);
    lat1 = radians(lat1);
    lat2 = radians(lat1);
    a = (sin(dlat/2)).^2 + cos(lat1) * cos(lat2) * (sin(dlon/2)).^2;
    c = 2 * atan2(sqrt(a),sqrt(1-a));
    c = c*6372.8;
   fprintf('distance: %4f km\n',6372.8 * c);
end
 
function rad = radians(degree) 
    rad = degree * pi / 180;
end

Bonne journée.

[souviron34]

Pour convertir un couple (latitude/longitude) en mètres, il faut connaître des choses essentielles :

• Te demande-t-on la distance à la surface du globe (courbe) ou la distance linéaire ?
  
  
• Quel que soit ce qu'on te demande, il te faudra :
  
  
  
  • le rayon de la Terre utilisé
  • l'ellipticité retenue
  • la position du centre de la Terre
  
  
  C'est ce que'on appelle le "datum"... (par exemple WGS84)
  
  
• Un complément du point 1 est : si l'on te demande la distance linéaire, avec quelle projection et quels paramètres de projection ? (Lambert, stéréographique, ..., quelle latitude de référence, etc)

La transformation de coordonnées Lat/Lon en coordonnées rectangulaires n'est malheureusement pas simple du tout...

[lhbilmarocain]

Oui merci pour les explications mais je l'ai déjà réalisé.
Merci

[souviron34]

Oui merci pour les explications mais je l'ai déjà réalisé.
Merci

Oui mais elle n'est pas générale (valeurs codées en dur), et fausse

Faire des dlat et dlon et prendre un actangente ne donne pas la solution mathématique...

Mais puisque tu es fier de toi, c'est l'essentiel, hein ?

J'espère juste que ça ne servira pas à guider le prochain avion que je prend...

[prgasp77]

Merci pour votre réponse,

Sinon, je reviens pour proposer ma solution, comme ça tout le monde en profitera :-)

Code MATLAB :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
function [a,c,dlat,dlon]=haversine11(lat1,lon1,lat2,lon2)
    dlat = radians(lat2-lat1);
    dlon = radians(lon2-lon1);
    lat1 = radians(lat1);
    lat2 = radians(lat1);
    a = (sin(dlat/2)).^2 + cos(lat1) * cos(lat2) * (sin(dlon/2)).^2;
    c = 2 * atan2(sqrt(a),sqrt(1-a));
    c = c*6372.8;
   fprintf('distance: %4f km\n',6372.8 * c);
end
 
function rad = radians(degree) 
    rad = degree * pi / 180;
end

Bonne journée.

Salut, on tente tous de faire ça la première fois qu'on a affaire à des coordonnées gps . Mais ce n'est pas la bonne solution, elle induira des erreurs importantes.
Pour moi, la solution la plus simple est de convertir les données WGS84 en UTM. Mais cette technique ne fonctionne que pour des points espacé d'une distance de l'ordre de 1000km (et quand je dis « de l'ordre », j'entends bien « ordre de grandeur »). L'avantage de l'UTM, c'est que tu travailles directement avec des distances (voir l'article Wikipedia).

On trouve des fonctions de conversions écrites en C, Fortran, Java, ... un peu partout sur le web.

Cordialement,

[lhbilmarocain]

Non ça marche parfaitement ! car y a une GROSSE fonction dessous qui fait un grand travail !
J'ai comparé les résultats avec Google Maps et ça donne le même résultat

[prgasp77]

J'ai comparé les résultats avec Google Maps et ça donne le même résultat

Pour l'avoir déjà fait je me répète : non ça ne marchera pas. Tente à une autre latitude pour te convaincre Il faut savoir que la Terre est loin d'être sphérique, c'est plutôt un genre de patatoïde.