[Proj] Proj command usage

Herve Baille hb at meteorage.com
Tue Sep 4 09:21:11 PDT 2007


Hello

New in this list.
We are changing our own projection functions for a standard tool, 
libproj4, on Sun Sparc Solaris 10.
We project (latitude=45, longitude=0) to a Lambert conic conform 
projection called "extended Lamber II" used in France.
Our program gives (x=415728.737817, y=2002668.173285), which is very 
close to the result of the National Geograpgic Institute.
The "proj" command (with correct projection parameters):
echo 0 45 | proj +proj=lcc +lon_0=2.3372917 +lat_0=46.8 +lat_2=45.8989 
+lat_1=47.696 +x_0=600000 +y_0=2200000 +ellps=clrk80
gives:
415663.74       2002665.67
We are not very far on Y axis but 65 meter away on X axis.
If I specify the XYZ offset between WGS84 and NTF (dx=-168, dy=-60, 
dz=320) with option +towgs84
the command:
echo 0 45 | proj +proj=lcc +lon_0=2.3372917 +lat_0=46.8 +lat_2=45.8989 
+lat_1=47.696 +x_0=600000 +y_0=2200000 +ellps=clrk80 
+towgs84=-168,-60,320,0,0,0,0
gives exatly the same result
415663.74       2002665.67
as if the option +towgs84 is not used.

Where could be the mistake?
Is there an other way to specify the translation?
Thanks
Herve




More information about the Proj mailing list