[Proj] Netherlands RD (28992) to UTM Zone 32 (3044) produces shifted coordinates

[Jan Hartmann]

Hi Jelmer,

The correct parameters for the grid shift from epsg:28992 to  WGS84 are 
described here:

http://trac.osgeo.org/geotiff/ticket/22#comment:3

There are several different formulas in the EPSG database, but the most 
recent one is your number 1 (= 4). It has been put in the PROJ4-database 
by Frank Warmerdam, as you can read in the posting above. Numbers 2 and 
3 are older definitions, although they deliver about the same results.  
They should be discarded, but will probably hang around yet for a long 
time in older software.

Groeten,

Jan

On 30-7-2013 10:49, Jelmer Baas wrote:
>
> Hello,
>
> For a new project we need to convert coordinates from EPSG 28992 to 
> EPSG 3044 with MapInfo, Oracle AND a few Proj.4 programs (PostGIS, a 
> custom DLL, Proj4js). The problem we are having is that we're seeing 
> some shifted coordinates.
>
> First test point (in the Netherlands)
>
> Origin               135745.924           517131.216
>
> Proj4js              236345.39474983042   5839505.550959114
>
> Postgis(1)           236345.390264227     5839505.55017783
>
> Postgis(2)           236345.394817119     5839505.55055913
>
> MapInfo              236345.47            5839504.058
>
> CS2CSN               236345.39            5839505.55
>
> DLL old(4)            236345.396465174     5839505.54496813
>
> DLL new(3)           236345.390259674     5839505.55009311
>
> Oracle               236345.49107424      5839504.12246656
>
> Test point (in Germany)
>
> Origin              254870.45            415123.56
>
> MapInfo             350265.960           5731670.259
>
> Postgis(1)          350263.4072329       5731667.880
>
> Postgis(2)           350263.4132521       5731667.880
>
> DLL new(3)           350263.407243136     5731667.88019826
>
> DLL old(4)           350263.403514925     5731667.88565864
>
> Oracle(5)           350265.967728329     5731670.33312477
>
> CS2CS                350263.41            5731667.88 43.92     cs2cs 
> +init=EPSG:28992 +to +init=EPSG:3044
>
> It seems MapInfo and Oracle more-or-less agree, as do all Proj.4-based 
> products. I don't know which are "right", meaning I also don't know 
> which ones are wrong -- let alone find the fault.
>
> Hopefully someone can help me with this problem, or at least point me 
> in the right direction.
>
> (1): "+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 
> +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel 
> +towgs84=565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725 
> +units=m +no_defs "
>
> (2) Another server where I updated the spatial_ref_system: 
> "+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 
> +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m 
> +towgs84=565.2369,50.0087,465.658,-0.406857330322398,0.350732676542563,-1.8703473836068,4.0812 
> +no_defs no_defs"
>
> (3) A custom built DLL for usage in VB6. "+proj=sterea 
> +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.999908 
> +x_0=155000 +y_0=463000 +ellps=bessel +units=m 
> +towgs84=565.2369,50.0087,465.658,-0.406857330322398,0.350732676542563,-1.8703473836068,4.0812 
> +no_defs <>"
>
> (4) Same DLL. "+proj=sterea +lat_0=52.15616055555555 
> +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 
> +ellps=bessel 
> +towgs84=565.417,50.3319,465.552,-0.398957,0.343988,-1.8774,4.0725 
> +units=m +no_defs <>"
>
> (5)select sdo_cs.transform(SDO_GEOMETRY('POINT(254870.45 415123.56)', 
> 90112),25832)  from dual;
>
> (5b) The Oracle definition contains some  tags I don't see anywhere in 
> Proj.4 (like the Bessel 1841 part):
>
> PROJCS["Netherlands National System", GEOGCS [ "Netherlands Bessel", 
> DATUM ["Netherlands Bessel", SPHEROID ["Bessel 1841", 6377397.155, 
> 299.1528128]], PRIMEM [ "Greenwich", 0.000000 ], UNIT ["Decimal 
> Degree", 0.01745329251994330]], PROJECTION ["Stereographic"], 
> PARAMETER ["Scale_Factor", 0.999908], PARAMETER ["Central_Meridian", 
> 5.387639], PARAMETER ["Latitude_Of_Origin", 52.156161], PARAMETER 
> ["False_Easting", 155000.000000], PARAMETER ["False_Northing", 
> 463000.000000], UNIT ["Meter", 1.000000000000]]
>
> A side note is that transformation via my DLL doesn't work properly 
> (i.e., coordinates are shifted tens of meters), unless I **first** 
> transform to WGS84 and then to 3044. I assume this is caused by an 
> easy-to-fix error, but I don't know what that is.
>
> P.S.: PostGIS, Proj4JS and my epsg file are all up-to-date.
>
> Thanks for any help,
>
> With kind regards,
>
> Jelmer Baas
>
>
> _______________________________________________
> Proj mailing list
> Proj at lists.maptools.org
> http://lists.maptools.org/mailman/listinfo/proj
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.maptools.org/pipermail/proj/attachments/20130730/8edf42dd/attachment-0001.htm