[Proj] Datum shift not working

Even Rouault even.rouault at spatialys.com
Sat Aug 30 03:40:22 EST 2014


Paul,

as far as I know, the proj utility doesn't do datum shifts. You should rather 
use cs2cs (and as per 
http://trac.osgeo.org/proj/wiki/FAQ#WhydoIgetdifferentresultswith4.5.0and4.6.0, 
"PROJ.4 will now only apply a datum shift if both the source and destination 
coordinate system have valid datum shift information") :

$ cs2cs +proj=longlat +ellps=clrk80 +towgs84=0,0,0 +no_defs +to +proj=utm 
+zone=40 +ellps=clrk80 
+towgs84=-180.624,-225.516,173.919,-0.81,-1.898,8.336,16.7101 +units=m 
+no_defs -f "%.6f"
59.28632619 24.27334106
731801.730465	2685882.517256 82.812593

$ cs2cs +proj=longlat +ellps=clrk80 +towgs84=0,0,0 +no_defs +to +proj=utm 
+zone=40 +ellps=clrk80 +units=m +no_defs -f "%.6f"
59.28632619 24.27334106
732098.167370	2686171.824604 0.000000

Even


Le samedi 30 août 2014 10:26:28, Paul Boettinger a écrit :
> Hello,
> i want to do a transformation with the EPSG Code 3440
> In my epsg file 3440 references the following parameters:
> +proj=utm +zone=40 +ellps=clrk80
> +towgs84=-180.624,-225.516,173.919,-0.81,-1.898,8.336,16.7101 +units=m
> +no_defs
> 
> Executing proj 4.8.0_0 from the MacPorts Repository on OSX 10.9.4 makes no
> differences using the transformation with or without the datum shift. Any
> help welcome,
> 
> Paul
> 
> proj +proj=utm +zone=40 +ellps=clrk80
> +towgs84=-180.624,-225.516,173.919,-0.81,-1.898,8.336,16.7101 +units=m
> +no_defs -f "%.6f" 59.28632619 24.27334106
> 732098.167370	2686171.824604
> 
> proj +proj=utm +zone=40 +ellps=clrk80 +units=m +no_defs -f "%.6f"
> 59.28632619 24.27334106
> 732098.167370	2686171.824604
> _______________________________________________
> Proj mailing list
> Proj at lists.maptools.org
> http://lists.maptools.org/mailman/listinfo/proj

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com


More information about the Proj mailing list