Ticket #7 (closed defect: fixed)

Opened 5 years ago

Last modified 3 years ago

Dutch Coordinate System Not Working

Reported by: lxnyce Owned by: warmerdam
Priority: major Milestone:
Component: default Version: unspecified
Keywords: double stereographic Cc:

Description

Please see this bug for an overview :  http://trac.osgeo.org/gdal/ticket/2487

We are contacting customer about getting the prj generated by ESRI. I'll repost the basic problem though. From a customer :

Just so we know what we are talking about: the Dutch coordinate system is ‘Rijksdriehoeksstelsel’. Luckily that is commonly abbreviated to ‘RD’.

Most of that is just parameters, and we can fill them out basically. But what is fundamentally lacking is the actual projection formula used in that coordinate system. It’s the ‘Double Stereographic’ projection, used worldwide only in NL and on Nova Scotia. It is different from the regular stereographic projection, one that is supported….

Please find enclosed two point shapefiles, depicting the same points. One is in RD, the other one is in geographic coordinates on the WGS84 datum. Please also find enclosed a PDF (published by the Dutch Cadastre, the ‘owners’ of the coordinate system) describing the specs of the coordinate system (chapter 4.2.1 for latlong à X, Y, and 4.2.2 for X, Y à latlong)

Attachments

Sampledataset.zip Download (25.6 KB) - added by lxnyce 5 years ago.
rd_wgs84_proj_deviations.gif Download (16.3 KB) - added by janh 4 years ago.
rd_wgs84_proj_deviations.2.gif Download (16.3 KB) - added by janh 4 years ago.
rd_wgs84_errors.gif.gif Download (25.7 KB) - added by janh 4 years ago.

Change History

Changed 5 years ago by lxnyce

in reply to: ↑ description   Changed 5 years ago by bart

In my opinion, PROJ.4 does support the RDNAP projection. The parameters i use are:

proj=sterea lat_0=52.15616055555555 lon_0=5.38763888888889 k=0.9999079 x_0=155000 y_0=463000 towgs84=565.040,49.910,465.840,-0.40939,0.35971,-1.86849,4.0772 ellps=bessel

First time i discussed this ith Frank the 7 parameter datum transformation was not working correctly but i think starting version 4.6 this works pretty well.

Last check i did was today, compared it with the results from MapInfo? 9.0 difference was less the 1 meter !

Replying to lxnyce:

Please see this bug for an overview :  http://trac.osgeo.org/gdal/ticket/2487 We are contacting customer about getting the prj generated by ESRI. I'll repost the basic problem though. From a customer : Just so we know what we are talking about: the Dutch coordinate system is ‘Rijksdriehoeksstelsel’. Luckily that is commonly abbreviated to ‘RD’. Most of that is just parameters, and we can fill them out basically. But what is fundamentally lacking is the actual projection formula used in that coordinate system. It’s the ‘Double Stereographic’ projection, used worldwide only in NL and on Nova Scotia. It is different from the regular stereographic projection, one that is supported…. Please find enclosed two point shapefiles, depicting the same points. One is in RD, the other one is in geographic coordinates on the WGS84 datum. Please also find enclosed a PDF (published by the Dutch Cadastre, the ‘owners’ of the coordinate system) describing the specs of the coordinate system (chapter 4.2.1 for latlong à X, Y, and 4.2.2 for X, Y à latlong)

Changed 4 years ago by janh

Changed 4 years ago by janh

Changed 4 years ago by janh

  Changed 4 years ago by janh

This ticket can be closed. The parameters as given by Bart give a correct projection within 20 cm of the official transformation ( http://www.rdnap.nl/download/setup%20cc41.exe). Attached a map with deviations between this proj transformation and the official one for 425 core measurement points ( https://rdinfo.kadaster.nl/rd/kernnet_meter.xls), and the official "error map", showing the measurement errors of the Dutch grid ( http://www.rdnap.nl/download/coor-en-proj.zip, p. 45, in Dutch). Both maps are nearly identical.

  Changed 3 years ago by warmerdam

  • status changed from new to closed
  • resolution set to fixed

Closing based on Jan's advice.

Note: See TracTickets for help on using tickets.