Right way to do this transformation/projection : ITRF2014 -> WGS84 UTM

--===============2062134248846790095==
Content-Type: multipart/alternative; boundary="000000000000e5b386059af15b0b"

--000000000000e5b386059af15b0b
Content-Type: text/plain; charset="UTF-8"

Hello,
First I want to apologize as this is a bit of a beginner question but after
hours of experimentation I am still confused.

I have coordinates in ITRF2014 (both in geocentric and geodetic
coordinates) with an associated EPOCH which I am trying to
transform/project to WGS84 UTM18N using pyproj/PROJ and EPSG codes. I am
working with surveying data, trying to maximize the accuracy. So far my
group made the approximation that coordinates in ITRF2014 were close enough
to WGS84 to not bother with the transformation but I would like to change
that.

I am looking for advice on how to do this properly :

   - Would you preferably use the geocentric coordinates for the
   transformation/projection?
   - Could I just do something like this :

#with EPSG 7789 +> ITRF2014 geocentric and EPSG 32618 WGS84 UTM18N

transformer = Transformer.from_crs(7789,32618)

X2,Y2, Z2, _ = transformer.transform(xx=X, yy=Y,zz=Z, tt=g.epoch)


During my experimentation I noticed that doing EPSG 7789 => EPSG 32618
holds different results (about 4cm diffrence) than  EPSG 7789  => EPSG 4919
=> EPSG 32614
Which surprised me as I thought that going from ITRF2014 to WGS84 involved
a step through ITRF2000.

In the same way converting between ITRF geocentric (EPSG 7789) and ITRF2014
lat long Height (EPSG 7912) gives me some 7cm differences which confused me
a lot and lead me to ask for help.

As a bonus question, do you know of a tool that support those
transformation and could help me validate my results?

Thanks a lot for your help!
Happy end of the year!
Regards

--000000000000e5b386059af15b0b
Content-Type: text/html; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable

<div dir=3D"ltr">Hello,<div dir=3D"ltr" class=3D"gmail_signature" data-smar=
tmail=3D"gmail_signature"><div dir=3D"ltr"><div dir=3D"ltr"><div dir=3D"ltr=
"><div dir=3D"ltr"><div style=3D"font-size:12.8px"></div></div></div></div>=
</div></div><div>First I want to apologize as this is a bit of a beginner q=
uestion but after hours of experimentation I am still confused.</div><div><=
br></div><div>I have coordinates in ITRF2014 (both in geocentric and geodet=
ic coordinates) with an associated EPOCH which I am trying to transform/pro=
ject to WGS84 UTM18N using=C2=A0pyproj/PROJ and EPSG codes. I am working wi=
th surveying data, trying to maximize the accuracy. So far my group made th=
e approximation that coordinates in ITRF2014 were close enough to WGS84 to =
not bother with the transformation=C2=A0but I would like to change that.</d=
iv><div><br></div><div>I am looking for advice on how to do this properly :=
</div><div><ul><li>Would you preferably use the geocentric coordinates for =
the transformation/projection?</li><li>Could I just do something like this =
:=C2=A0</li></ul></div><blockquote style=3D"margin:0 0 0 40px;border:none;p=
adding:0px"><blockquote style=3D"margin:0 0 0 40px;border:none;padding:0px"=
><blockquote style=3D"margin:0 0 0 40px;border:none;padding:0px"><div>#with=
 EPSG 7789=C2=A0+&gt; ITRF2014 geocentric and EPSG 32618 WGS84 UTM18N=C2=A0=
=C2=A0</div></blockquote></blockquote><blockquote style=3D"margin:0 0 0 40p=
x;border:none;padding:0px"><blockquote style=3D"margin:0 0 0 40px;border:no=
ne;padding:0px"><div>transformer =3D Transformer.from_crs(7789,32618)</div>=
</blockquote></blockquote><blockquote style=3D"margin:0 0 0 40px;border:non=
e;padding:0px"><blockquote style=3D"margin:0 0 0 40px;border:none;padding:0=
px"><div>X2,Y2, Z2, _ =3D transformer.transform(xx=3DX, yy=3DY,zz=3DZ, tt=
=3Dg.epoch)</div></blockquote></blockquote></blockquote><blockquote style=
=3D"margin:0 0 0 40px;border:none;padding:0px"><div><br></div></blockquote>=
During my experimentation I noticed that doing EPSG 7789 =3D&gt; EPSG 32618=
 holds different results (about 4cm diffrence) than=C2=A0

EPSG 7789=C2=A0 =3D&gt; EPSG 4919 =3D&gt; EPSG 32614

<br><div>Which surprised me as I thought that going from ITRF2014 to WGS84 =
involved a step through ITRF2000.</div><div><br></div><div>In the same way =
converting between ITRF geocentric (EPSG 7789) and ITRF2014 lat long Height=
 (EPSG 7912) gives me some 7cm differences which confused me a lot and lead=
 me to ask for help.</div><div><br></div><div>As a bonus question, do you k=
now of a tool that support those transformation and could help me validate =
my results?<br><br>Thanks a lot for your help!<br>Happy end of the year!<br=
>Regards<br><br></div></div>

--000000000000e5b386059af15b0b--

--===============2062134248846790095==
Content-Type: text/plain; charset="utf-8"
MIME-Version: 1.0
Content-Transfer-Encoding: base64
Content-Disposition: inline

X19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX18KUFJPSiBtYWls
aW5nIGxpc3QKUFJPSkBsaXN0cy5vc2dlby5vcmcKaHR0cHM6Ly9saXN0cy5vc2dlby5vcmcvbWFp
bG1hbi9saXN0aW5mby9wcm9qCg==

--===============2062134248846790095==--