[OSRS-PROJ] how use "proj"????

[Frank Warmerdam]

Gracia-Lorca Gregory wrote:
> I would use "proj" to convert points in lambert 2 to lambert 3,
> I have done :
>  
> proj +init=esri:27572 +proj=lcc +lat_1=44.10000000000001 
> +lat_0=44.10000000000001 +lon_0=2.33722917 +k_0=0.999877499 +x_0=600000 
> +y_0=3200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0  <<EOF
> 807000  1842000
> EOF
>  
> but the result is
>  
> *  *
>  
> Any idea?
> Can you tell me how use "proj" or "cs2cs"?

Greg,

With proj you could accomplish this as a two stage process if there is no datum
change required.

warmerda at gdal2200[3]% proj -I +init=esri:27572 | proj +proj=lcc +lat_1=44.10000000000001 +lat_0=44.10000000000001 
+lon_0=2.33722917 +k_0=0.999877499 +x_0=600000 +y_0=3200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0
807000 1842000

result:
806696.79       3142093.82

With the cs2cs commmand it can be done in one program:

cs2cs +init=esri:27572 +to +proj=lcc +lat_1=44.10000000000001 +lat_0=44.10000000000001 +lon_0=2.33722917 
+k_0=0.999877499 +x_0=600000 +y_0=3200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0
807000 1842000

result:
806707.71       3141745.34 -165.43

The difference in results in presumably because cs2cs is applying the
datum shift, though I haven't verified this.

Best regards,
-- 
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Programmer for Rent


----------------------------------------
PROJ.4 Discussion List
See http://www.remotesensing.org/proj for subscription, unsubscription
and other information.