<div dir="ltr"><div dir="ltr">Thanks all for the help!</div><div dir="ltr"><br></div><div dir="ltr">Dinking around a little bit more, I found through trial and error that eliminating the +vunits seemed to solve it:</div><div dir="ltr"><br></div><div dir="ltr">echo 565890.947 4192011.754 <b>55</b> | cs2cs +proj=utm +zone=10 +datum=NAD83 +units=m +no_defs +geoidgrids=us_noaa_g2018u0.tif +to +proj=utm +zone=10 +datum=WGS84 +units=m +no_defs<br>565890.95 4192011.75 <b>22.79</b><br></div><div dir="ltr"><br></div><div>Somehow, the addition of +vunits=m on the target angers proj:</div><div><br></div><div>echo 565890.947 4192011.754 <b>55</b> | cs2cs +proj=utm +zone=10 +datum=NAD83 +units=m +no_defs +geoidgrids=us_noaa_g2018u0.tif +to +proj=utm +zone=10 +datum=WGS84 +units=m <b>+vunits=m </b>+no_defs<br>565890.91 4192011.91 <b>55.00</b><br></div><div dir="ltr"><br></div><div dir="ltr"><div>Adding +vunits=m to the source seems not to matter. Wish I could say more on why this is the case; is it a bug or is it a malformed request?</div><div><br></div><div>Thanks!</div><div>Patrick</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Jul 1, 2020 at 7:11 PM Brian Shaw <<a href="mailto:brian.shaw@noaa.gov">brian.shaw@noaa.gov</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>
Hi Patrick<br>
<br>
NOAA does maintain the NAD 83, NAVD 88 and geoid models, looks like
the typo might be in the readme from anothers email since NGA
maintains WGS84 and EGM08. <br>
<br>
I apologize as Im not completely sure on the format of your command
but I think I have it you are passing UTM coordinates in NAD83 to
geoid18 to get an NAVD88 height except that instead of passing an
ellipsoid height you are passing an orthometric height. <br>
<br>
To calculate the ellipsoid height you would use the formula h=H+N
where h (ellipsoid height, in this case NAD 83), H (orthometric
height, in this case NAVD 88) and N (geoid height, in this case
GEOID18).<br>
<br>
You should be providing an ellipsoid height to the geoid to get an
orthometric height. Maybe its because its the end of my work day
that I cant figure out how to do this in one step but I believe if
you provide 0 as the height this will then give you the geoid height
for that location (in geoid18) and you can just add this to your
orthometric height to get the ellipsoid height. You should get a
negative number since the geoid is all negative for CONUS.<br>
<br>
more information about geoid18 can be found here:<br>
<a href="https://geodesy.noaa.gov/GEOID/GEOID18/" target="_blank">https://geodesy.noaa.gov/GEOID/GEOID18/</a><br>
<br>
Hopefully this helps if I think of a way to do it in one command Ill
respond again if someone else doesnt provide a solution.<br>
<br>
Cheers<br>
Brian<br>
<br>
<div>On 7/1/2020 5:49 PM, Patrick Young
wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">Hello,
<div><br>
</div>
<div>I am trying to figure out how to convert some elevations I
have into heights above ellipsoid via the cs2cs tool. I am
using proj 7.0.1, and I've run projsync --all to snag all the
grids.</div>
<div><br>
</div>
<div>Here is an example of what isn't working for me:</div>
<div><br>
</div>
<div>echo 565890.947 4192011.754 55 | cs2cs +proj=utm +zone=10
+datum=NAD83 +units=m +no_defs +geoidgrids=us_noaa_g2018u0.tif
+to +proj=utm +zone=10 +datum=NAD83 +units=m +vunits=m
+no_defs<br>
565890.95 4192011.75 55.00<br>
</div>
<div><br>
</div>
<div>I am confused as to why the elevation value isn't changing
for me. My (extremely limited) understanding is that I am
asking for 55m to be converted from a height referenced
against the NADV88 geoid to one referenced against the NAD83
ellipsoid. Am I using the wrong geoid grid? From the grid
readme, it mentions "us_nga_g2018u0.tif : Conterminous US,
referenced to NAD83(2011). Vertical datum: NAVD88"
(presumably the nga <-> noaa is a typo in the readme) so
I think I'm using the right one.</div>
<div><br>
</div>
<div>Thanks for any help,</div>
<div>Patrick</div>
</div>
<br>
<fieldset></fieldset>
<pre>_______________________________________________
PROJ mailing list
<a href="mailto:PROJ@lists.osgeo.org" target="_blank">PROJ@lists.osgeo.org</a>
<a href="https://lists.osgeo.org/mailman/listinfo/proj" target="_blank">https://lists.osgeo.org/mailman/listinfo/proj</a>
</pre>
</blockquote>
<br>
<pre cols="72">--
*************************************
Brian Shaw
Rocky Mountain Regional Advisor (CO, MT, WY)
NOAA's National Geodetic Survey (NGS)
Phone # 240-533-9522</pre>
</div>
</blockquote></div></div>