<div dir="ltr"><div dir="ltr">Hi Even,<br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Oct 22, 2020 at 10:41 AM Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Matthias,<br>
<br>
(would be more a topic for proj mailing list)<br></blockquote><div><br></div><div>(noticed just after sending, but immediately after figured that the audience will likely be similar)<br></div><div> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
> I am trying to apply a vertical grid shift to a coordinate. I am<br>
> extrapolating from the example in the documentation at<br>
> <a href="https://proj.org/operations/transformations/vgridshift.html" rel="noreferrer" target="_blank">https://proj.org/operations/transformations/vgridshift.html</a><br>
> <br>
> The command I am entering is<br>
> <br>
> echo 46.9480 7.4474 590 2020 | cct +proj=vgridshift<br>
> +grids=./geoid+HTRANS_ETRS2020.gtx<br>
> <br>
> the result is:<br>
> <br>
>  46.9480000000    7.4474000000          -inf     2020.0000<br>
> <br>
> I was hoping for the altitude being corrected by roughly 50 meters (instead<br>
> of being -inf).<br>
<br>
If the coordinate is located at Kornhausplatz in Bern, Switzerland, rather <br>
than at Dollo in Ethiopia, you'd better switch your 2 first figures ;-)<br></blockquote><div><br></div><div>Fair point, I am indeed more interested in Bern results than in Ethiopian results.</div><div>Unfortunately The result doesn't change significantly<br></div><div><br></div><div>echo 7.4474 46.9480 590 2020 | cct +proj=vgridshift +grids=./geoid+HTRANS_ETRS2020.gtx       <br>  7.4474000000   46.9480000000          -inf     2020.0000<br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
When using PROJ pipelines, the coordinate order is longitude,latitude / <br>
easting,northing , unless there are +proj=axisswap operations changing that. <br>
When using higher level PROJ interfaces like cs2cs or projinfo with EPSG code, <br>
then you'll have to comply with EPSG axis order.<br></blockquote><div><br></div><div>Thanks for the context, good to know !<br></div><div>Matthias<br></div><div> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
Even<br>
<br>
-- <br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
</blockquote></div></div>