[Proj] Central Meridian (lon_0) sign
Richard didier
didier.richard at ign.fr
Tue Dec 14 05:28:31 EST 2004
Hi all,
When looking at the epsg file ($PROJ_HOME/share/proj) for the french Lambert I
noticed that lon_0 is negative whilst the longitude is at east of Greenwich.
The following tests show that something is wrong with the negative value :
Transformation between 2 lamberts :
cs2cs -f "%.2f" +init=epsg:27561 +to +init=epsg:27572 \
<<EOF
750000.00 300000.00
EOF
gives:
750283.12 2600360.77 0.00
which is correct.
Transformation between a lambert and its base geographic system :
cs2cs +init=epsg:27561 +to +init=epsg:4807 \
<<EOF
1029705.083 272723.849
EOF
gives:
3d39'45.975"E 50dN 0.000
which is wrong, it should be :
8d20'14.025"E 50dN 0.000
Transformation between a geographic system and its lambert :
cs2cs -f "%.3f" +init=epsg:4807 +to +init=epsg:27562 \
<<EOF
2d45'43.846"E 46d37'35.430"N
EOF
gives:
990233.155 193379.657 0.000
which is wrong, it should be :
632542.063 180804.162 0.000
I obtain the expected correct values by changing the sign of lon_0 (from
-2.33722917 to 2.33722917).
Any idea ?
--
Didier RICHARD, chef de projet "Diffusion des données numériques"
Institut Geographique National tel: +33/0 1.43.98.83.23
2/4 Avenue Pasteur fax: +33/0 1.43.98.81.71
F-94165 St-Mande Cedex e-mail: Didier.Richard at ign.fr
More information about the Proj
mailing list