Le portail francophone de la géomatique

Salut Jérôme C,

Je confirme OpenStreetMap utilise bien la projection de Mercator, cependant etant un vrai newbie en SIG, je m excuse d'avance si mes questions semblent evidente

Apres avoir fait une recherche sur cette projection (wikipedia pour pas le citer) je suis tombé sur deux equations permettant d'obtenir les coordonnées x et y a d'un point a partir de sa lat/lon ma question est la suivante :
si j applique cette convertion pour mes 3 points est ce que je peux par la suite utiliser une simple formule de calcul de distance d'un point a une droite ? ou est ce que j ai rien compri lol (dans ce cas comment s y prendre svp )

Bonjour,
Je suis d'accord avec Jérôme. Une fois qu'on travaille dans la projection, on calcule suivant la géométrie cartésienne en plan, et il n'y a juste qu'à appliquer l'éventuelle correction des distances, due à la projection, exactement de la même façon que lorsqu'on travaille en lambert (LT93 par exemple). Si les distances sont petites, cette correction est négligeable.
Cordialement

Bonjour,
Les formules indiquées ont pour paramètres la lat et la lon. Je croyais que vous travaillez dans la projection Mercator, c'est à dire avec un X et un Y cartésien.
L'utilisation des coordonnées géographiques est possible, mais fait appel à la géométrie sphérique.
Je vais essayer de détailler non propos. Votre chemin est caractérisé un segment de droite A-B. Ce trajet est-il le plus court chemin de A à B tel que le prendrait un avion (l'orthodromie) ou est-ce le segment de droite traçé avec une règle sur la plan entre les points A et B? (la loxodromie).
N'hésitez pas à me détailler les hypothèses, j'essayerai d'être plus précis. Je ne connais pas OpenStreeMap, et je ne pourrai vous répondre que sur le calcul lui-même.

Si cette opération doit être répétitive, ce qui semble être le cas, ça vaut le coup de mettre au point une petite formule.
Cordialement.

wouw eh ben ca y est je suis perdu

en lisant la formule sur wikipedia ( je cite : Les équations suivantes déterminent les coordonnées x et y d'un point sur une carte de Mercator à partir de sa latitude φ et de sa longitude λ (avec λ0 au centre de la carte) )
je pensais justement que la formule me retournait x et y par rapport a une Lat / Lon passé en parametre

Je vais essayer de detailler la situation :

J'ai une map (apparement en projection Mercator)
Lorsque j exporte les données provenant de la region selectionnée, ces dernieres sont caracterisées de la facon suivante :
- une liste de chemin
- chaque chemin est composé d'une liste de noeud ayant pour parametres une latitude et une longitude

J'ai une trace (provenant de la region selectionnée) composé de la meme facon
- un chemin
- une liste de noeud qui compose se chemin

l'objectif c'est de prendre 2 noeuds de ma premiere source (les données provenant directement de la map) de former un segment, de prendre par exemple le 1e noeud de ma trace et de calculer la distance entre ce dernier et le segment formé

les distances sont assez courtes (quelques km pas plus)


ps : j ai trouvé ce petit bout de code pour faire la conversion

Private Function CalcTileXY(ByVal lat As Single, ByVal lon As Single, ByVal zoom As Long) As Point

        Dim xf As Single = (lon + 180) / 360 * 2 ^ zoom
        Dim yf As Single = (1 - Math.Log(Math.Tan(lat * Math.PI / 180) + 1 / Math.Cos(lat * Math.PI / 180)) / Math.PI) / 2 * 2 ^ zoom
        CalcTileXY.X = CLng(Math.Floor(xf))
        CalcTileXY.Y = CLng(Math.Floor(yf))

    End Function

Dernière modification par Dada97.2 (Wed 10 December 2008 17:21)

C'est tout a fait ca !
Je vous remercie pour l'attention que vous portez a mon probleme

re
Juste quelques precisions : le zoom qui apparait dans ma fonction correspond au niveau de zoom de ma map je suppose ?
apparement d'apres le lien suivant : http://wiki.openstreetmap.org/index.php … _tilenames
j'obtiens tilex et tiley c'est donc a partir de ca que je fais mes calculs ?
l'interpretation le resultat obtenu dependra de l'echelle de ce zoom ?

Cordialement.

Dernière modification par Dada97.2 (Wed 17 December 2008 02:59)

Bonjour,
Quelle joie de retrouver son forum préféré
Concernant le sujet en cours, je résume la question
Soit un réseau R composé de noeuds reliés par des segments. et des chemins C constitués de noeuds reliés par des segments. Ces deux ensembles sont connus en coordonnées géographiques Lon/lat.
Le problème est de calculer la distance de chaque noeud de C par rapport au segment le plus proche du réseau R.

Si telle est la question, comme prévu j'ai établi une formule de calcul de la distance en mètres d'un point à un segment.
Ja vais rédiger la méthode, et je mettrai un lien très vite.
Pour mémoire, l'utilisation d'une projection conduit à des écarts non négligeables
Cordialement

Bonjour,
Comme promis, voila l'explication détaillée
http://www.dlzlogic.com/Geodesie.pdf
Il est bien évident que toute remarque ou observation sera la bien venue.
bon courage.

Il est vrai que je n'ai pas fait de fonction qui calcule la distance d'un point à un segmant, mais j'ai expliqué la formule. Ou plutôt j'ai explique comment, lorsqu'on connait dans un triangle sphérique la  longueur AM et l'angle A, si ce triangle satisfait aux hypothèses de Legendre (excès sphérique) on peut calculer en géométrie plane.

Le source qui génère la fenêtre affiché est très grand, mais sans intérêt, puisqu'il ne s'agit que de saisie, transformation d'unités, affichage etc. Les fonctions permettant de transformer du LT93 en WGS sont dans une librairie dont j'ai cité l'auteur, et qui est libre.

Voila les lignes qui calculent la longueur de la projection

  double Zap=CalcAzimut(LA,MA,LP,MP);
  double ArcAP=CalcArc(LA,MA,LP,MP); // radians
fprintf(espion,"Azimut    Zap =%f g\n",Zap);
  double Zbp=CalcAzimut(LB,MB,LP,MP);
  double ArcBP=CalcArc(LB,MB,LP,MP); // radians
fprintf(espion," azimut  Zbp=%f g \n", Zbp);
// Calcul de¨PH , plus courte distance de P à AB
// l'angle en A est la différence des azimuts AB AP
   double AngA=Zap-ZAB;
   if (AngA > 100) AngA=200.0-AngA;
   double eAP=(double)0.5*ArcAP*ArcAP*sin(AngA*M_PI/200.0) * 200 / M_PI; // ATTENTION
   double DisAP=CalcLongArc(LA,MA,LP,MP);
   double DisPHa=fabs(DisAP*sin((AngA-eAP/3)*M_PI/200.0));
fprintf(espion,"A partir de A AngA=%f g DisAP=%0.2f m. DisPHa=%0.2f m. \n",AngA,DisAP,DisPHa);

   double AngB=Zbp-ZBA;
   if (AngB > 100) AngB=200.0-AngB;
   double eBP=(double)0.5*ArcBP*ArcBP*sin(AngB*M_PI/200.0) * 200 / M_PI;  // ATTENION
   double DisBP=CalcLongArc(LB,MB,LP,MP);
   double DisPHb=fabs(DisBP*sin((AngB-eBP/3)*M_PI/200.0));
fprintf(espion,"A partir de B AngB=%f g DisBP=%0.2f m. DisPHb=%0.2f m. \n",AngB,DisBP,DisPHb);
   eProj->Text=Ftoa((DisPHa+DisPHb)/2,2);


Seules les fonction concernant les deux problèmes de base me paraissaient dignes d'être jointes, le reste n'est que du petit calcul informatique. Mais bien sûr, je peux le rajouter. Il y a tout de même une autre application qui mériterait un module, c'est le calcul d'un point N à partir d'un point A à une distance donnée D et suivant un azimut Z. C'est un peu plus compliqué. Mais si ça intéresse quelqu'un ? à suivre.
Bonne soirée.

Dernière modification par Pierre Dolez (Fri 19 December 2008 21:41)