Je réalise actuellement un site web aéronautique, et dans ce cadre je souhaite proposer des nouvelles fonctionnalités à mes utilisateurs.
L'une de ces nouvelles fonctions sera de pouvoir cliquer sur une carte et d'afficher les informations à l'emplacement géographique du clic.
J'ai donc trouvé sur le site de l'IGN des formules pour convertir un X,Y Lambert vers coordonnées WGS84 mais je n'arrive pas à les appliquer correctement, je dois faire une erreur quelque part.
Lien IGN :
http://professionnels.ign.fr/DISPLAY/000/526/701/5267019/NTG_71.pdf (page 12)
http://www.ign.fr/DISPLAY/000/526/700/5267002/transformaton.pdf (page 9)
Ma projection est défini comme suit :
a=6378137; //demi grand axe de l'ellipsoide (m)
e=1/298.257221976; //première excentricité de l'ellipsoide
//paramètres de projections
l0=deg2rad(2); //longitude de référence
phi0=deg2rad(47); //latitude d'origine en radian
phi1=deg2rad(45); //1er parallele automécoique
phi2=deg2rad(49); //2eme parallele automécoique
x0=0; //coordonnées à l'origine
y0=0; //coordonnées à l'origine
Je calcul "n" et "c" grace à ces valeurs comme suit :
n=0,73150235178492
c=11747978,099289
Mon point d'essai est l'aérodrom de Lognes Emerainville situé :
Latitude : 48,822
Longitude : 2,623
X : 45656.295287915
Y : 202915.47936158
Ces valeurs sont issu du calcul GPS -> Lambert et son correctes puisque le résultat attendu est correct :
"n" et "c" sont normalement correctes pour la même raison, car utilisé à l'identique dans la formule WGS->Lambert
Je cherche déjà à calculer la longitude qui est le plus simple :
g = arctan(X-Xs/Ys-Y)
l = phi0 +g/n
avec Xs=0 et Ys=0 soit :
g=arctan(X/-Y)
X/-Y = -0,224921428
y = Atan(-0,224921428)=-0,221239655
Long=phi0 + y/n
y/n=-0,302445582
Long=0,820304748437335-0,302445582
Long=0,517859166140255
Au lieu de 2,623...
Je n'arrive pas à comprends mon erreur !
Merci d'avance pour votre aide précieuse !