<div dir="ltr">Thank you so much for your help it is much clearer now and a perfect timing with the new release.<div>That is good to know about WGS84 and the resulting uncertainty, it is slightly worrying, I will investigate that further.<br>Thank again<br><br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div style="font-size:12.8px;color:rgb(80,0,80)"><p><b><span style="color:rgb(102,102,102)">Julien Schroder | Spatial Analyst | Photographer</span></b><u></u></p></div><div><a style="color:rgb(34,34,34);font-size:12.8px">julien.schroder@gmail.com</a><font color="#500050"><span style="font-size:12.8px"> | +1 514 797 1464<br></span></font></div><div><font color="#500050"><span style="font-size:12.8px"><br></span></font></div><div style="font-size:12.8px"><div style="font-size:12.8px;color:rgb(80,0,80)"><a href="http://arctic-mood.com/" style="color:rgb(17,85,204)" target="_blank">Website</a> | <a href="https://www.instagram.com/arctic.mood/?hl=en" style="color:rgb(17,85,204)" target="_blank">Instagram</a> | <a href="https://www.linkedin.com/in/julien-schroder-2225a578" style="color:rgb(17,85,204)" target="_blank">Linkedin</a> | <a href="https://www.facebook.com/julien.schroder.1/" style="color:rgb(17,85,204)" target="_blank">Facebook</a> | <a href="http://twitter.com/Julienschroder" style="color:rgb(17,85,204)" target="_blank">Twitter</a></div></div></div></div></div></div></div></div></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Le lun. 4 nov. 2019 à 07:36, Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> a écrit :<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Brian,<br>
<br>
I've been working on that very subject those last days. Pending work is in<br>
<a href="https://github.com/OSGeo/PROJ/pull/1710" rel="noreferrer" target="_blank">https://github.com/OSGeo/PROJ/pull/1710</a><br>
and<br>
<a href="https://github.com/OSGeo/proj-datumgrid/pull/60" rel="noreferrer" target="_blank">https://github.com/OSGeo/proj-datumgrid/pull/60</a><br>
<br>
to add support for GEOID18 (as well as older GEOID versions for specialized <br>
purposes)<br>
For getting a value in the grids, PROJ will use the Datum of the GeographicCRS <br>
indicated in the EPSG database for the transformation record between NAVD88 <br>
and the NAD83 realizations you mention below, possibly doing a horizontal <br>
adjustment to go back and forth between the initial and final CRS and this <br>
interpolation CRS. It will not go through WGS84.<br>
<br>
So for example (with the 2 above pull requests in place), transforming between <br>
NAD83(2011) + NAVD88 height and NA83(2011) (3D) shows that no horizontal <br>
adjustment is attempted:<br>
<br>
projinfo -s EPSG:6318+5703 -t EPSG:6319 --spatial-test intersects -o PROJ<br>
Candidate operations found: 2<br>
-------------------------------------<br>
Operation n°1:<br>
<br>
unknown id, Inverse of NAD83(2011) to NAVD88 height (3), 0.015 m, USA - CONUS <br>
- onshore<br>
<br>
PROJ string:<br>
+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert <br>
+xy_in=deg +xy_out=rad +step +proj=vgridshift +grids=g2018u0.gtx +multiplier=1 <br>
+step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1<br>
<br>
-------------------------------------<br>
Operation n°2:<br>
<br>
unknown id, Inverse of NAD83(2011) to NAVD88 height (2), 0.02 m, USA - Alaska<br>
<br>
PROJ string:<br>
+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert <br>
+xy_in=deg +xy_out=rad +step +proj=vgridshift +grids=g2012ba0.gtx <br>
+multiplier=1 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step <br>
+proj=axisswap +order=2,1<br>
<br>
Even<br>
<br>
> GEOID18 supersedes GEOID12B in CONUS and Puerto Rico.  GEOID12B should<br>
> be used in Alaska, Hawaii, Guam, and the Commonwealth of the Northern<br>
> Mariana Islands.  American Samoa should not be used for anything as<br>
> ASVD02 is no longer valid after earthquakes.<br>
> <br>
> A note on Geoids, Horizontal (Geometric) Datums and Vertical Datums.  To<br>
> get an orthometric height (aka elevation) in a particular vertical datum<br>
> you should use coordinates in a particular geometric datum with lat,lon,<br>
> and ellipsoid height along with a geoid model to get an orthometric<br>
> height.  For the US here are the appropriate datums and geoid models. <br>
> Note that the corresponding geoid models should be used based on the<br>
> datum of your coordinates. Also note the the original NAD 83 (1986) was<br>
> a 2D datum while all other realizations were 3D, aka geometric.<br>
> <br>
> Also note that if you are concerned about accuracy, WGS84 should be<br>
> avoided if at all possible.   WGS84 has a 3-5 meter uncertainty so any<br>
> time you use coordinates or use it as a pivot in transformations you<br>
> introduce this uncertainty into your coordinates.<br>
> <br>
> Table on geoids and NAD 83 realizations<br>
> <br>
> NAD 83 (realization)  Geoid Model     Regions<br>
> NAD 83 (1986)         N/A<br>
> NAD 83 (HARN)         GEOID96/GEOID99<br>
> NAD 83 (FBN)  GEOID03<br>
> NAD 83 (CORS96)       GEOID03<br>
> NAD 83 (NSRS2007)     GEOID09<br>
> NAD 83 (2011)         GEOID12B        HI, AK, Guam, CNMI<br>
> NAD 83 (2011)         GEOID18         CONUS, PRVI<br>
> <br>
<br>
-- <br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
_______________________________________________<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>