RDNAPTRANS in Excel

mluijben · 30 oktober 2020 om 12:09

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?

Anton · 30 oktober 2020 om 12:16

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.

Jochem · 30 oktober 2020 om 14:12

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

RD&NAP (EPSG:7415) naar ETRS89 (EPSG:4937), dit is RDNAPTRANS
Aanvragen download RDNAPTRANS2018
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
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:

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

leo · 29 januari 2021 om 14:05

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

Jochem · 29 januari 2021 om 14:53

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.

Thibaud · 23 februari 2022 om 10:24

kijk eens op Regiolab-Delft data services

mathieu-kool · 14 juni 2022 om 12:23

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)),"")

Jochem · 15 juni 2022 om 11:27

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.

mathieu-kool · 15 juni 2022 om 11:49

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)

Jochem · 17 juni 2022 om 15:18

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

Jochem · 7 oktober 2022 om 09:25

Er staat een tikfout in een van de coëfficienten 144.221 moet 114.221 zijn. Zonder die tikfout is de transformatie van ETRS89 / WGS84 naar RD een stuk beter: fouten tot een halve meter.