<div dir="ltr"><div dir="ltr"><br><br>On Fri, May 3, 2019 at 8:58 PM Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> wrote:<br>><br>> On vendredi 3 mai 2019 13:01:01 CEST Alan Snow wrote:<br><div>[...]<br></div><div>></div>> The +init=ITRF2000:ITRF93 definition uses a +proj=helmert operation which<br>> operates in the cartesian/geocentric coordinate space. So you need to convert<br>> from geographic to geocentric coordinates first.<br>><br>> $ echo "-117.61748, 59.4967, 329.396" | src/cct -d 7 +proj=pipeline \<br>>   +step +proj=unitconvert +xy_in=deg +xy_out=rad \<br>>   +step +proj=cart +ellps=GRS80 \<br>>   +step +init=ITRF2000:ITRF93 \<br>>   +step +inv +proj=cart +ellps=GRS80 \<br>>   +step +proj=unitconvert +xy_in=rad +xy_out=deg<br><div><br></div><div></div><div>that means we (developers) need to check with proj_angular_input() if such a unit conversion is needed? AFAIU, units are now degrees as in EPSG:4326, and in this case a conversion to radians is no longer required, as in proj-4.</div><div><br></div><div>Markus M<br></div></div><div dir="ltr"><br></div></div>