Ben ik toch nog een keer… Misschien een nieuw topic waard, maar ik ga voor nu toch de vraag eerst hier stellen.
Vanaf de ETRS89-coördinaten wil ik nog 1 stap verder naar ITRF2014. ook weer ene hoop ingelezen, maar ik kom er nog niet uit.
Volgens de discussie hier: RDNAPTRANS in Excel - #3 door Jochem - RDNAPTRANS - Geoforum
Bestaat de transformatie uit 3 stappen:
- Eerst weer naar een geocentrische coördinatenop basis van GRS80.
(Moet dat met uitgerekende hoogte of ook weer met de fictieve hoogte van 43m) - Daarna het tijdsverschil uitrekenen in jaren tussen de huidige datum en 1-1-2010 (Voor de voorbeeldberekening wordt er 2020-10-01 genomen als huidige datum)
- Dan de juiste helmert parameters uitrekenen, met als basis de volgende waarden:
array(‘tx’=>0.0547,‘ty’=>0.0522,‘tz’=>-0.0741,‘rx’=>deg_to_rad1(+0.001701/3600),‘ry’=>deg_to_rad1(+0.010290/3600),‘rz’=>deg_to_rad1((-0.016632)/3600),‘s’=>+0.00000000212);
En als snelheid/jaar: array(‘tx’=>0.00010,‘ty’=>0.00010,‘tz’=>-0.00190,‘rx’=>deg_to_rad1(+0.0000810/3600),‘ry’=>deg_to_rad1(+0.0004900/3600),‘rz’=>deg_to_rad1((-0.0007920)/3600),‘s’=>+0.000000000110);
Waar kan ik de bron vinden van deze waarden? Ik heb ze nu uit de excel-discussie en uit de sas-implementatie van wgs84toitrf2014
En na de helmert-transformatie, weer terug naar ITRF2014, ook weer m.b.v. GRS80.
En dan zou ik op het juiste antwoord moeten uitkomen (volgens de SAS-implementatie).
ETRS 51.999994778222, 4.9999922685, 42.98187019
WGS84 ≈ ITRS2014: 52, 5, 43
Ik heb allerlei manieren geprobeerd, maar ik kom op andere antwoorden
Uitgangspunt ETRS89
‘lat’ => float 51.999994778222
‘lon’ => float 4.9999922685
‘h’ => float 42.98187019
Geocentric:
‘x’ => float 3920013.6292309
‘y’ => float 342956.22040882
‘z’ => float 5002836.8721211
Helmert-parameters na toepassen epoch-correctie (10.75):
‘tx’ => float -0.145225
‘ty’ => float -0.147725
‘tz’ => float 3.724475
‘rx’ => float -7.768569583705E-7
‘ry’ => float -4.6995050568092E-6
‘rz’ => float 7.5959347040671E-6
‘s’ => float -2.177975E-7
Na Helmert:
‘x’ => float 3919988.5121045
‘y’ => float 342990.10530378
‘z’ => float 5002852.3928946
Ellipsiodal (radian):
‘lat’ => float 0.90757534651209
‘lon’ => float 0.087275462441856
‘h’ => float 41.6441
Zou ITRS2014 moeten zijn:
‘lat’ => float 52.000236945266
‘lon’ => float 5.0005156529709
‘h’ => float 41.6441
Is er weer iemand die me op weg kan helpen?