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

Tue Jul 30 02:57:28 PDT 2013

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>

>  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.osgeo.org/pipermail/proj/attachments/20130730/0b4c008b/attachment.html>