Distance Point Segment !

Version imprimable

Bonjour tout le monde,
j' ai la problematique suivante :
j ai 3 coordonnées (latitude longitude) et je dois calculer la distance entre une des coordonnées vis a vis d'un segment formé avec les 2 autres
j'ai mis en place le code suivant qui me retourne une distance ! mais le probleme c'est que je ne sais pas dans quel unité exprimer cette distance (si c'est en mettre, kilometre etc) si le resultat ne correspond a aucun de ces types pouvez vous m'indiquer comment le convertir

Le code :
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 /* * Calcule la distance du TkPT par rapport au segment des 2pts .osm */ float distancePointSegment(void){ float a,b,A,B,xH,yH; float x1,x2,y1,y2,xM,yM; float distres=0; x1 = atof(coordonnees[1]); y1 = atof(coordonnees[0]); x2 = atof(coordonnees[3]); y2 = atof(coordonnees[2]); xM = atof(lon); yM = atof(lat); if (y2==y1){ a=(y2-y1)/(x2-x1); b=y1-a*x1; A=xM*(x2-x1)-(y2-y1)*(b-yM); B=(x2-x1)+a*(y2-y1); xH=A/B; distres=abs(xH-xM)*sqrt(1+pow((x2-x1)/(y2-y1),2)); }else{ xH=xM; yH=y1; distres=sqrt(pow((xH-xM),2)+pow((yH-yM),2)); } return distres; }
Merci pour votre aide

23/11/2008, 23h13
Zavonen

Citation:

J'ai 3 coordonnées (latitude longitude)

Chez moi 1+1 ça fait 2 !

Citation:

je dois calculer la distance entre une des coordonnées vis a vis d'un segment formé avec les 2 autres

Je sais calculer la distance d'un point à un autre, la distance d'un point à une droite, à un segment. Une coordonnée est un NOMBRE. Et entre le "entre" et le "vis à vis" on attrape le tournis.
Sérieusement il serait temps que les posts soient rédigés en un langage, sinon scientifiquement correct, au moins syntaxiquement correct.
NB: Si tu es étranger ne tiens pas compte de ma remarque.(mais les DOM c'est bien la France, non ?)
Pour finir comment veux tu qu'on réponde à ce genre de charabias.
Je pourrai reformuler ta question, mais je préférerais que tu le fasses.
24/11/2008, 01h10
snipes

:king: un simple : c'est mal formulé j'aurais compri !

Citation:

J'ai 3 coordonnées (latitude longitude)

Je voulais dire j'ai 3 points !
je te rassure chez moi aussi ca fait 2 : 1+1
24/11/2008, 08h34
Zavonen

Tu es ici en coordonnées 'sphériques' (rho, phi ,theta) en fait la valeur de rho est constante parce que tes points sont sur la surface d'une sphère. Donc sur une sphère (rho fixe) un point est effectivement caractérisée par deux angles et deux seulement. Mais ce que tu sembles ignorer c'est que sur la surface d'une sphère les 'segments' ne sont pas des portions de droite au sens usuel, ce sont des arcs de grands cercles. Les formules permettant de calculer la distance de deux points en fonction de leurs coordonnées sphériques sont données par les formules d'Haversine:http://en.wikipedia.org/wiki/Haversine_formula
Quand à la distance d'un point à un segment elle est donnée par:
d(P,[AB])=min(d(P,A),d(P,B),d(P,H)) où H est la projection orthogonale de P sur (AB), mais là encore il s'agit de géométrie sphérique, et les notions de 'segment', 'droite', 'disatnce' doivent être compris au sens de la géométrie sphérique.
24/11/2008, 11h52
Ulmo

Citation:

Envoyé par Zavonen

[...}la distance d'un point à un segment [...] est donnée par:
d(P,[AB])=min(d(P,A),d(P,B),d(P,H)) où H est la projection orthogonale de P sur (AB)[...]

C'est d(P,H) uniquement si H est entre A et B, sinon c'est bien min(d(P,A),d(P,B)).

Maintenant, la notion de "entre" sur une sphère...
24/11/2008, 14h10
FR119492

Salut!
En pratique, il faudrait connaître la taille du domaine dans lequel se trouvent tes points. S'ils sont très éloignés, c'est vrai, comme l'écrit Zavonen, que tu dois travailler en coordonnées sphériques et que la ligne reliant deux points est un arc de cercle, et non un segment de droite. En revanche, si les points sont proches, une projection (Mercator, Lambert ou autre) sera plus simple et ne produira qu'une erreur négligeable.
Jean-Marc Blanc
25/11/2008, 14h12
souviron34

encore un fois :roll: je répète : transformer des latitudes/longitudes en mètres dépend de la projection utilisée...

Il faut donc connaître la projection utilisée, transformer les lat/lon en coordonnées rectangulaires (mètres), et calculer la distance.

Sinon Jean-Marc a raison, ce sera des coordonnées sphériques (ou plutôt elliptiques), mais la vraie distance devra tenir compte de l'ellipsoide utilisé (WGS 84 ??) , avec son vrai rayon de la Terre, son vrai centre X,Y,Z, et sa vraie approximation pour cette partie du globe..... Très nettement plus compliqué et peut-être pas ce que tu veux...
26/11/2008, 11h13
Ulmo

Et je me répète aussi, mais sur la sphère, 2 points ne définissent pas un segment (il y en a deux).
26/11/2008, 17h19
pseudocode

:google: :arrow: Orthodromie

Bonjour

Je ravive ce sujet pour obtenir une réponse à la question suivante :

Comment connaitre de quel côté du segment se trouve le point ?

Voici ma source en Delphi :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
Function TMapForm.DistPntSegm(Pnt, Pnt1, Pnt2 : PointFloat) : TPntToSegm;
Var
 Len, ParamR : Double;
 PntR : PointFloat;
begin                              
  Result.x := 0; Result.y := 0; Result.Dist := RefMaxDist;
  Len := Dist2Pnts(Pnt1, Pnt2);
  ParamR := (((Pnt.X - Pnt1.X)*(Pnt2.X - Pnt1.X)) + ((Pnt.Y - Pnt1.Y)*(Pnt2.Y - Pnt1.Y)))/Sqr(Len);
  if (ParamR > 0) and (ParamR < 1) then
  begin  //  La projection du point est sur le segment aux coordonnées
    PntR.x := Pnt1.X + ParamR * (Pnt2.X - Pnt1.X);
    PntR.y := Pnt1.Y + ParamR * (Pnt2.Y - Pnt1.Y);
    Result.x := PntR.x; // X point projection de Pnt.X sur le segment
    Result.y := PntR.y; // Y point projection de Pnt.Y sur le segment
    Result.Dist := Dist2Pnts(Pnt, PntR); //  Longueur du segment
  end;
end;

Merci d'avance

24/05/2009, 14h52
souviron34

bonjour

une bonne source d'informations (entre autres) est la FAQ des UseNet Newsgroup, en particulier celle du newsgroup comp.graphics.algorithms :

http://www.faqs.org/faqs/graphics/algorithms-faq/

Voir les points 1.02 et 1.03

:D
05/01/2010, 14h57
nicomac3

Pb transformation de coord

Bonjour

je réalise une étude sur les radars utilisés en aéronautique et pour cela j'ai besoin de convertir des coordonnées de la manière suivante :
j'ai bien essayé de me débrouiller entre les solutions aux réponses apportées et le lien de newsgroup (super bien fait d'ailleurs!) mais je ne suis pas parvenu à m'en sortir.
voila mon pb :
je connais un point (mon radar) en latitude, longitude
pour chaque radar, je connais tout les plots l'entourant (les petits avions) en coordonnées polaires Rho et Theta

j'ai besoin de convertir l'emplacement défini en (rho, theta) avec comme origine ma latitude / longitude du radar en latitude / longitude du plot lui meme.

je ne sais pas si j'ai réussi à me faire comprendre. est ce que je dois utilisé une proj style WGS84 ou je ne connais strictement rien? je suis un peu perdu

merci pour votre aide
05/01/2010, 15h50
souviron34

je dirais qu'il faut connaître quel est le rapport de projection utilisé sur ton écran, ainsi que la projection utilisée.

Normalement, il doit y avoir quelque part cette spec, soit dans la doc soit dans le code du logiciel de radar, avec vraisemblablement une formule..

Ensuite, c'est simple...

Par contre, il faut soit le code soit la doc du radar et surtout de l'installation, avec l'écran..

PS: oui, vraisemblablement le système de coordonées sera WGS84..
05/01/2010, 17h26
nicomac3

Citation:

je dirais qu'il faut connaître quel est le rapport de projection utilisé sur ton écran, ainsi que la projection utilisée.

le radar utilise une projection orthogonale sur le sol. c'est cette distance en NM, que je récupère. de plus je ne sais pas s'il y a une intégration de la variation du NM (dépend un tt petit peu de lat/long) ou de la variation du rayon de la terre étant donné que la distance max des avions est d'environ 200 NM (360km)

pour l'écran, enfin la visualisation si j'ai bien compris j'utilise google earth

Citation:

Ensuite, c'est simple...

tu connais la formule magique??