<div dir="ltr"><div>Hi</div><div><br></div><div>I am testing a conversion from WGS72 (EPSG:4322) to Cape (EPSG:4222), and I expect to get a noticeable difference in Z, due to the significant datum shift.</div><div>However I do not get it.</div><div><br></div><div>To avoid any misunderstanding, I am converting everything to 3D with this command:</div><div><span style="font-family:monospace">as3d () { projinfo EPSG:$1 --single-line -o WKT2:2019 --3d -q } </span></div><div><br></div><div>Using projinfo with "--spatial-test intersects" returns a reasonable transformation (I think), that composes two Helmert transformations.</div><div></div><div></div><div>What am I doing wrong?</div><div><br></div><div><div>Pivoting manually via EPSG:4326 gives reasonable values in Z<br></div></div><div><br></div><div>Run with PROJ 7.2.1 compiled in Ubuntu 18.04, without any grid file.<br></div><div><br></div><div>Thanks<br></div><div><br></div><div><span style="font-family:monospace">$ echo -33.9 18.4 0 | cs2cs "$(as3d 4322)" "$(as3d 4222)" -d 8<br>-33.89987620   18.40079770 0.00000000</span></div><div><span style="font-family:monospace"><br>$ echo -33.9 18.4 0 | cs2cs "$(as3d 4322)" "$(as3d 4326)" -d 8 | cs2cs "$(as3d 4326)" "$(as3d 4222)" -d 8<br>-33.89987619 18.40079770 -34.18419224<br>$ echo -33.9 18.4 0 | cs2cs "$(as3d 4322)" "$(as3d 4326)" -d 8<br>-33.89996800      18.40015389 -3.05061427<br>$ echo -33.9 18.4 0 | cs2cs "$(as3d 4326)" "$(as3d 4222)" -d 8<br>-33.89990820       18.40064382 -31.13371642<br></span></div><div><span style="font-family:monospace"><br></span></div><div><span style="font-family:monospace">$ projinfo -s "$(as3d 4322)" -t "$(as3d 4222)" --spatial-test intersects<br>Candidate operations found: 5<br>-------------------------------------<br>Operation No. 1:<br><br>unknown id, WGS 72 to WGS 84 (2) + Inverse of Cape to WGS 84 (1), 11 m, South Africa - mainland onshore.<br><br>PROJ string:<br>+proj=pipeline<br>  +step +proj=axisswap +order=2,1<br>  +step +proj=unitconvert +xy_in=deg +z_in=m +xy_out=rad +z_out=m<br>  +step +proj=push +v_3<br>  +step +proj=cart +ellps=WGS72<br>  +step +proj=helmert +x=0 +y=0 +z=4.5 +rx=0 +ry=0 +rz=0.554 +s=0.219<br>        +convention=position_vector<br>  +step +proj=helmert +x=136 +y=108 +z=292<br>  +step +inv +proj=cart +a=6378249.145 +rf=293.4663077<br>  +step +proj=pop +v_3<br>  +step +proj=longlat +a=6378249.145 +rf=293.4663077<br>  +step +proj=unitconvert +xy_in=rad +z_in=m +xy_out=deg +z_out=m<br>  +step +proj=axisswap +order=2,1<br></span>...<br><br></div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature">.___ ._ ..._ .. . ._.  .___ .. __ . _. . __..  ... .... ._ .__<br>Entre dos pensamientos racionales <br>hay infinitos pensamientos irracionales.<br><br></div></div></div>