<div dir="ltr"><div>Javier, useful link, I think it did provide me with some higher level of clarity :)</div><div><br></div><div>Thanks,</div><div><br></div><div>Mathieu<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Nov 10, 2020 at 5:10 PM Javier Jimenez Shaw <<a href="mailto:j1@jimenezshaw.com">j1@jimenezshaw.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"><div dir="ltr"><div>Hi Mathieu</div><div><br></div><div>I do not know if I understand your question correctly. In Switzerland the geoid is above the ellipsoid (about 50 meters). That is a positive geoid undulation. So the orthometric height (over the geoid) is smaller than the ellipsoidal height.<br></div><div><br></div><div>This may clarify it:<br><a href="https://www.ngs.noaa.gov/GEOID/GSVS/" target="_blank">https://www.ngs.noaa.gov/GEOID/GSVS/</a></div><div><br></div><div>Cheers<br></div><div><div><div dir="ltr">.___ ._ ..._ .. . ._.  .___ .. __ . _. . __..  ... .... ._ .__<br>Entre dos pensamientos racionales <br>hay infinitos pensamientos irracionales.<br><br></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, 10 Nov 2020 at 10:39, Mathieu Pellerin <<a href="mailto:nirvn.asia@gmail.com" target="_blank">nirvn.asia@gmail.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"><div dir="ltr">Greetings,<br><br>I’m hoping the mailing list can enlighten me regarding handling of vertical shift in PROJ.<br><br>Using this Switzerland vertical shift grid (<a href="https://www.dropbox.com/s/6uzpg3gr47q1p1k/CHGeo04.gtx?dl=0" target="_blank">https://www.dropbox.com/s/6uzpg3gr47q1p1k/CHGeo04.gtx?dl=0</a>), if I type in the following command to adjust vertical height (100 meters here):<br><br>echo 7.9995 46.7949 100 | cct +proj=vgridshift +grids=CHGeo04.gtx<br><br>PROJ returns the following result:<br><br>7.9995000000   46.7949000000       49.8406           inf<br><br>The adjusted height is 49.8 meters. This puzzles me as the grid’s value at the provided longitude and latitude (i.e.  7.9995 46.7949) is 50.179. I would have assumed that this meant a vertical adjustment of +50.179, not -50.179 as observed by the PROJ returned values.<br><br>Is PROJ’s result the correct one based on the grid? Or should an inverse transformation (i.e. cct -t) be used here? If so, what’s the rationale of whether to pick forward or inverse transformation?<br><br>Regards,<br><br>Mathieu Pellerin</div>
_______________________________________________<br>
PROJ mailing list<br>
<a href="mailto:PROJ@lists.osgeo.org" target="_blank">PROJ@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/proj" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/proj</a><br>
</blockquote></div>
</blockquote></div>