Wolfram Function Repository

Instant-use add-on functions for the Wolfram Language

Function Repository Resource:

 FromRDCoordinates

Source Notebook

Convert Dutch RD coordinates to GeoPosition objects

Contributed by: Sander Huisman

ResourceFunction["FromRDCoordinates"][{x,y}]

converts the Dutch RD coordinates {x,y} to GeoPosition objects.

Details and Options

Converts the Dutch RD coordinates (Rijksdriehoekscoördinaten) {x,y} to a latitude-longitude pair in the WGS84 datum using a fourth-order conformal mapping obtained from 6,000 references points in the Netherlands.

Examples

Basic Examples (1) 

Convert an RD coordinate to a GeoPosition object:

In[1]:=
ResourceFunction["FromRDCoordinates"][{257688, 471277}]
Out[1]=

Applications (1) 

Plot the waterways of The Netherlands using data from www.pdok.nl given in RD coordinates:

In[2]:=
(* Evaluate this cell to get the example input *) CloudGet["https://www.wolframcloud.com/obj/80bc8baa-d2d1-4bc1-adc4-7b32ee148d9f"]
Out[3]=

Properties and Relations (2) 

Using FromRDCoordinates and the resource function ToRDCoordinates we can convert from and to RD coordinates:

In[4]:=
ResourceFunction["FromRDCoordinates"][ ResourceFunction["ToRDCoordinates"][{52.22, 6.89}]]
Out[5]=

The "origin" of the RD coordinates is somewhere in France:

In[6]:=
GeoGraphics[GeoMarker[ResourceFunction["FromRDCoordinates"][{0, 0}]], GeoRange -> 400000]
Out[6]=

The closest city is La Celle-Saint-Cyr:

In[7]:=
GeoNearest[Entity["City"], ResourceFunction["FromRDCoordinates"][{0, 0}]]
Out[7]=

Possible Issues (2) 

The coordinates are valid only within the border given by the following RD coordinates:

In[8]:=
bounds = {{141000, 629000}, {100000, 600000}, {80000, 500000}, {-7000, 392000}, {-7000, 336000}, {101000, 336000}, {161000, 289000}, {219000, 289000}, {300000, 451000}, {300000, 614000}, {259000, 629000}}; GeoGraphics[Polygon[ResourceFunction["FromRDCoordinates"] /@ bounds]]
Out[9]=

Asking for a coordinate outside this region—in this example, Brussels—spawns a message and extrapolates the fit:

In[10]:=
ResourceFunction["FromRDCoordinates"][{80515, 316108}]
Out[10]=

Publisher

SHuisman

Version History

  • 1.0.0 – 26 June 2019

Source Metadata

  • Citation:
    • de Kooker, E., "Het stelsel van de Rijksdriehoeksmeting en het European Terrestrial Reference System 1989–Gebruik en definitie.", Kadaster: Apeldoorn, 2000.
    • Schreutelkamp, F. H., Strang van Hees, G. L., "Benaderingsformules voor de transformatie tussen RD– en WGS84–kaartcoördinaten." Geodesia, no. 2, 64–69, 2001.

Related Resources

License Information