<div dir="ltr"><div dir="ltr">Hi Jochem,<div><br></div><div>So if I'm not entirely mistaken, you'll need something like the Forth "swap"</div><div>operator, or the PostScript "exch" operator to allow something like this:</div><div><br></div><div>###  The inverse part (read from bottom up, and read pop as push, push as pop) ###<br><br></div><div>step omit_fwd proj=push v_3              # restore the geophysical height<br><br># The bread and butter<br>step omit_fwd proj=tmerc inv lat_0=12.180658675 lon_0=-68.251802281 k_0=1 x_0=23209.5600 y_0=21423.9900 ellps=intl<br>step omit_fwd proj=cart ellps=intl<br>step omit_fwd proj=helmert x=-366.1939 y=-115.0688 z=-776.7039 rx=20.963080 ry=16.462749 rz=-14.276379 s=-12.80900 convention=coordinate_frame exact<br>step omit_fwd proj=cart inv ellps=GRS80<br><br># The great swap show...<br>step omit_fwd proj=push v_3              # restore the geometrical height<br>step omit_fwd proj=exch                  # swap the two heights, to get the geometrical to the top-of-stack<br>step omit_fwd proj=pop v_3               # save the geophysical height<br>step omit_fwd proj=vgridshift grids=egm2008_0_45.tif inv<br>step omit_fwd proj=pop v_3                 # save the geometrical height</div><div># End of inverse part - read from here and upwards</div><div><br></div><div>In that case, I think the start of a solution would be to open an issue</div><div>over at <a href="https://github.com/OSGeo/PROJ/issues">https://github.com/OSGeo/PROJ/issues</a>, requesting an</div><div>exch/swap operator, and summarizing this thread</div><div><br></div><div><br></div><div><br></div></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg-4092377013765372164"><div lang="NL" style="overflow-wrap: break-word;">
</div>

</div></blockquote></div></div>