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

Jose Gonçalves jagoncal at gmail.com
Tue Jul 30 03:53:10 PDT 2013


Hello

There are a few possible reasons for the differences you find. In the
Oracle definition:

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]]

you have the central point (5.387639 E, 52.156161 N) rounded with 6 decimal
places. The sixth decimal place is 11 cm in latitude, so it may justify
differences of that order. The ellipsoid definition is correct, so it is no
reason for differences.
The projection you are using in PROJ.4 is "Oblique Stereographic
Alternative" (+proj=sterea). PROJECTION ["Stereographic"] probably
corresponds to +proj=stere, which is different.

Another possible reason are the datum transformation parameters. Check if
the parameters used by Mapinfo and Oracle in the Bursa-Wolf transformation
are the same you are using with cs2cs. Different sets of parameters justify
differences of a few meters.

Regards

José Gonçalves





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/7e9e7e73/attachment.html>


More information about the Proj mailing list