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

[Jan Hartmann]

I don't think that that is the problem here. The differences are only a 
few centimeters, while errsors resulting from missing +towgs84 
parameters are in the order of 100 meters.

Jan

On 30-7-2013 11:57, Milo van der Linden (gmail) wrote:
> Hello Jelmer,
>
> 1. I would strongly advice to transform everything to wgs84 first.
> 2. wgs84 and etrs89 are based upon different spheroids/datums. 
> Transforming from wgs84 TO etrs89 will not be correct if your etrs89 
> definition does not have +towgs84 defined. And by default I believe it 
> isn't. Take a look at this ticket: 
> http://trac.osgeo.org/proj/ticket/15 and see if it is of any help.
>
> So, what I would do is:
>
> epsg:28992 -> wgs84 (visual check) -> etrs89
>
> Good luck!
>
>
>
>
> 2013/7/30 Jelmer Baas <baas at speerit.nl <mailto:baas at speerit.nl>>
>
>     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 <mailto:Proj at lists.maptools.org>
>     http://lists.maptools.org/mailman/listinfo/proj
>
>
>
> _______________________________________________
> 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.osgeo.org/pipermail/proj/attachments/20130730/025e6d7a/attachment.html>