|
Lat,long -> Lambert93
|
 Informations rapides |
|
|
|
Informations détaillées
|
Description :
Projection de coordonnées géographiques en Lambert 93.
Ce code permet de placer des coordonnées sur une carte.
Portion de 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
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
|
/*
références :
Méthode de calcul pour une projection de type lambert conique conforme sécante (
NTG_71.pdf):
http://www.ign.fr/affiche_rubrique.asp?rbr_id=1700&lng_id=FR
*/
//variables:
//système WGS84
$a=6378137; //demi grand axe de l'ellipsoide (m)
$e=0.08181919106; //première excentricité de l'ellipsoide
//paramètres de projections
$l0=$lc=deg2rad(3); //longitude de référence
$phi0=deg2rad(46.5); //latitude d'origine en radian
$phi1=deg2rad(44); //1er parallele automécoïque
$phi2=deg2rad(49); //2eme parallele automécoïque
$x0=700000; //coordonnées à l'origine
$y0=6600000; //coordonnées à l'origine
//coordonnées du point à traduire
$phi=deg2rad($latitude);
$l=deg2rad($longitude);
//calcul des grandes normales
$gN1=$a/sqrt(1-$e*$e*sin($phi1)*sin($phi1));
$gN2=$a/sqrt(1-$e*$e*sin($phi2)*sin($phi2));
//calculs de slatitudes isométriques
$gl1=log(tan(pi()/4+$phi1/2)*pow((1-$e*sin($phi1))/(1+$e*sin($phi1)),$e/2));
$gl2=log(tan(pi()/4+$phi2/2)*pow((1-$e*sin($phi2))/(1+$e*sin($phi2)),$e/2));
$gl0=log(tan(pi()/4+$phi0/2)*pow((1-$e*sin($phi0))/(1+$e*sin($phi0)),$e/2));
$gl=log(tan(pi()/4+$phi/2)*pow((1-$e*sin($phi))/(1+$e*sin($phi)),$e/2));
//calcul de l'exposant de la projection
$n=(log(($gN2*cos($phi2))/($gN1*cos($phi1))))/($gl1-$gl2);//ok
//calcul de la constante de projection
$c=(($gN1*cos($phi1))/$n)*exp($n*$gl1);//ok
//calcul des coordonnées
$ys=$y0+$c*exp(-1*$n*$gl0);
//calcul des coordonnées lambert
$x93=$x0+$c*exp(-1*$n*$gl)*sin($n*($l-$lc));
$y93=$ys-$c*exp(-1*$n*$gl)*cos($n*($l-$lc));
|
|
| Informations & Services |
 |
|
|
|
 Newsletter
Autre Sources
Tips
| Convertisseurs |
Aucune Catégorie
|
|
