RDNAPTRANS in Excel

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