RDNAPTRANS in Excel

Heeft of weet iemand wellicht een Excel bestand om RD coordinaten naar WGS84 om te zetten en vice versa? Dat moet volgens mij toch mogelijk zijn?

RDNAPTRANS in Excel zul je niet vinden. Je kunt bij het Kadaster als ontwikkelaar de omrekenformules opvragen inclusief de correctiegrids. Maar dat heb je niet zomaar in Excel ingebouwd. Defensie levert een programma dat gebaseerd is op RDNAPTRANS maar die mag je niet commercieel gebruiken.

Er zijn wel Excelbestanden (ook websites) waar je WGS84 naar RD om kan rekenen maar die zijn niet heel nauwkeurig. Je kunt van alles vinden als je binnen Google zoekt op ā€œExcel WGS RDā€.

Een andere mogelijkheid is QGIS gebruiken waarbij je je project instelt op RD en dan een CSV inlaadt vanuit WGS84. Er is volgens mij ook een projectie die wel gebruik maakt van een correctiegrid maar weet niet of dit op RDNAPTRANS is gebaseerd.

De nauwkeurige transformatie van RD naar WGS84 bestaat uit de volgende stappen:

  1. RD&NAP (EPSG:7415) naar ETRS89 (EPSG:4937), dit is RDNAPTRANS
    Aanvragen download RDNAPTRANS2018
  2. ETRF2000 (EPSG:7931) naar ITRF2008 (EPSG:7911), dit is een tijdsafhankelijke transformatie naar ITRS waarbij je dus een kalenderdatum op moet geven voor het juiste resultaat.
    http://etrs89.ensg.ign.fr/pub/EUREF-TN-1.pdf
  3. ITRF2008 (EPSG:7911) gelijkstellen aan WGS84-G1762 (EPSG:7665) voor 1 juli van het betreffende jaar.
    ftp://ftp.nga.mil/pub2/gandg/website/wgs84/NGA.STND.0036_1.0.0_WGS84.pdf

Dit in Excel programmeren kan wel, maar dat is niet eenvoudig, dus ik verwacht niet dat iemand dat gedaan heeft. Het is handiger om bestaande software te gebruiken om een Excel-bestand in te lezen, te transformeren en weer naar Excel te exporteren.

Met de meeste software gaat de transformatie van RD naar WGS84 op twee punten fout:

  1. Het correctiegrid van RDNAPTRANS wordt vaak nog overgeslagen, dat geeft een fout tot 0,25 m;
    Voor controle van het juist toepassen van RDNAPTRANS is er een validatieservice:
    Validatieservice - NSGI
  2. ETRS89 wordt vaak gelijkgesteld aan ITRS en WGS84, dat geeft een fout van inmiddels 0,8 m.
    Voor controle hiervan kun je kijken of software de coƶrdinaten wel verandert bij deze transformatie en of dit ook daadwerkelijk andere coƶrdinaten oplevert als je een andere kalenderdatum opgeeft, die gelijk zijn aan die van deze website:
    EUREF Permanent GNSS Network

Met de open-source software PROJ kan de transformatie correct uitgevoerd worden. Het commando voor RDNAPTRANS staat in de RDNAPTRANS-download. Het commando voor ETRF2000 naar ITRF2014 is:

cct -o output.txt -d 9 -I +proj=pipeline
  +step +proj=axisswap +order=2,1,3,4
  +step +proj=cart +ellps=GRS80
  +step +proj=helmert +x=0.0547 +y=0.0522 +z=-0.0741
    +rx=-0.001701 +ry=-0.010290 +rz=0.016632 +s=0.00212
    +dx=0.0001 +dy=0.0001 +dz=-0.0019
    +drx=-0.000081 +dry=-0.000490 +drz=0.000792 +ds=0.00011
    +convention=coordinate_frame +exact
  +step +proj=cart +inv +ellps=GRS80
  +step +proj=axisswap +order=2,1,3,4	 
  input.txt

Testpunt:
input.txt: 52.000000000 5.000000000 43.7000 10.750
(kalenderdatum 2020-10-01 => 2020.75 - referentie epoche 2010.0 = 10.75)

output.txt: 52.000005186 5.000007667 43.7179 10.7500
Ā 

Transformatie van ITRS naar WGS84 raad ik af, omdat ik het gebruik van WGS84 af raad. Ik adviseer ITRF2014 (ā‰ˆWGS84-G1762) te gebruiken.

4 likes

Edo Schrama destijd TU delft, heeft ongeveer 10 jaar geleden een conversieprogje in excel gemaakt en via internet gedeeld. Volgens mij werkt dat uitstekend.

Ik ken dat Excel-programma van E.J.O. Schrama niet, maar ik vermoed sterk dat het een benadering gebruikt. Als het 10 jaar geleden gemaakt is gebruikt het in ieder geval niet de procedure van de nieuwste RDNAPTRANSā„¢ voor de transformatie van RD naar ETRS89 en ook niet de actuele parameters voor de transformatie van ETRS89, via ITRS naar WGS84.

kijk eens op Regiolab-Delft data services

een Excel formule die hetzelfde resultaat geeft als de website van TU Delft Regiolab Data Services:

(waarbij cel G3 de x in RD bevat:)
=ALS(G3>0,(52.15517+(SOM((((((3235.65389*((H3-463000)10^-5))+(-32.58297((G3-155000)10^-5)^2)+(-0.2475((H3-463000)10^-5)^2)+(-0.84978((G3-155000)10^-5)^2((H3-463000)10^-5))+(-0.0655((H3-463000)10^-5)^3)+(-0.01709((G3-155000)10^-5)^2((H3-463000)10^-5)^2)+(-0.00738((G3-155000)10^-5))+(0.0053((G3-155000)10^-5)^4)+(-0.00039((G3-155000)10^-5)^2((H3-463000)10^-5)^3)+(0.00033((G3-155000)10^-5)^4((H3-463000)10^-5))+(-0.00012((G3-155000)10^-5)((H3-463000)*10^-5)))))))/3600)),"")

(waarbij cel H3 de y in RD bevat:)
=ALS(H3>0,(5.387206+(SOM((((((5260.52916*((G3-155000)10^-5))+(105.94684((G3-155000)10^-5)((H3-463000)10^-5))+(2.45656((G3-155000)10^-5)((H3-463000)10^-5)^2)+(-0.81885((G3-155000)10^-5)^3)+(0.05594((G3-155000)10^-5)((H3-463000)10^-5)^3)+(-0.05607((G3-155000)10^-5)^3((H3-463000)10^-5))+(0.01199((H3-463000)10^-5))+(-0.00256((G3-155000)10^-5)^3((H3-463000)10^-5)^2)+(0.00128((G3-155000)10^-5)((H3-463000)10^-5)^4)+(0.00022((H3-463000)10^-5)^2)+(-0.00022((G3-155000)10^-5)^2)+(0.00026((G3-155000)*10^-5)^5))))))/3600)),"")

1 like

Dit is een zogenaamde reeksontwikkeling die een benadering geeft van de transformatie tussen RD en ETRS89 / WGS84 met een nauwkeurigheid van ongeveer een meter (controleer zelf). Dit mag geen RDNAPTRANSā„¢ genoemd worden, want toestemming voor het gebruik van deze merknaam wordt alleen verleend bij implementatie met een nauwkeurigheid van 1 mm.

4 likes

Je hebt gelijk. Maar voor Geocaching oid van wellicht gemakkelijk in gebruik.

En de omgekeerde versie ervan, van WGS naar RD:

(waarbij C11 Latitude en D11 Longitude bevat:)
=155000+(190094.945 * (0.36 * (D11 - 5.38720621))) + (-11832.228 * (0.36 * (C11 - 52.1551744)) * (0.36 * (D11 - 5.38720621))) + (-144.221 * (0.36 * (C11 - 52.1551744))^2 * (0.36 * (D11 - 5.38720621))) + (-32.391 * (0.36 * (D11 - 5.38720621))^3) + (-0.705 * (0.36 * (C11 - 52.1551744))) + (-2.34 * (0.36 * (C11 - 52.1551744))^3 * (0.36 * (D11 - 5.38720621))) + (-0.608 * (0.36 * (C11 - 52.1551744)) * (0.36 * (D11 - 5.38720621))^3) + (-0.008 * (0.36 * (D11 - 5.38720621))^2) + (0.148 * (0.36 * (C11 - 52.1551744))^2 * (0.36 * (D11 - 5.38720621))^3)

(waarbij C11 Latitude en D11 Longitude bevat:)
=463000+(309056.544 * (0.36 * (C11 - 52.1551744))) + (3638.893 * (0.36 * (D11 - 5.38720621))^2) + (73.077 * (0.36 * (C11 - 52.1551744))^2 ) + (-157.984 * (0.36 * (C11 - 52.1551744)) * (0.36 * (D11 - 5.38720621))^2) + (59.788 * (0.36 * (C11 - 52.1551744))^3 ) + (0.433 * (0.36 * (D11 - 5.38720621))) + (-6.439 * (0.36 * (C11 - 52.1551744))^2 * (0.36 * (D11 - 5.38720621))^2) + (-0.032 * (0.36 * (C11 - 52.1551744)) * (0.36 * (D11 - 5.38720621))) + (0.092 * (0.36 * (D11 - 5.38720621))^4) + (-0.054 * (0.36 * (C11 - 52.1551744)) * (0.36 * (D11 - 5.38720621))^4)

1 like

Oei! De transformatie van ETRS89 / WGS84 naar RD is een stuk slechter, fouten tot ca.:

  • 70 meter binnen de rechthoek ruim om Nederland (50-56Ā°NB, 2-8Ā°OL)
  • 35 meter binnen de Nederlandse EEZ van de Noordzee
  • 4 meter op land binnen de Nederlandse rijksgrens
2 likes