Kleine afwijking bij RDNAPTRANS2018 in PHP

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?