<div dir="ltr"><div>Thanks Even and Greg.</div><div><br></div><div>My problem is more oriented to <b>projected</b> CRS. Now I see that I have not mentioned this term before.</div><br><div>I agree with you that the ellipsoidal height should refer to the proper geographical CRS. In this case, the "base one" for the projected CRS in the compound CRS. Something similar applies to other compound CRS combinations, where the geographical CRS reference should match the projected and the vertical, right?<br></div><div><br></div><div>My problem is locating points from a point cloud initially with GPS coordinates. It is possible to convert from GPS coordinates in WGS84 - 3D (EPSG:4979) to any State Plane (projected) + NAVD88. Is it possible also to convert from WGS84 - 3D to the State Plane without any vertical CRS.</div><div><br></div><div><span style="font-family:monospace">This is contrary to a gravity-based height, where if someone gives you 
just the elevation at 2 points, you can for example know from which 
point towards which one water would flow.</span></div><div><br></div><div>I need the elevation to note the point in a 3D space, not to show the gravity-based height. In that case, of course, I need a gravity-based model.<br></div><div><br></div><div>I would like to store those points in the projected CRS in a LAS file.</div><div>If I just define the projected CRS (without vertical component) in the LAS, what should I expect to be the Z values?</div><div><br></div><div><span style="font-family:monospace">For projected CRS, yes, since there is currently no Projected 3D CRS in EPSG</span></div><div><span style="font-family:monospace"><br></span></div><div><span style="font-family:monospace"><font face="arial,sans-serif">Having Projected 3D CRS in EPSG would be a solution... is there any plan for that? There were more than 7000 projected CRSs in EPSG. I hardly imagine that they duplicate that number.</font><br></span></div><div><br></div><div>Related to this topic (I guess), if I run this "experiment", <br></div><div><span style="font-family:monospace">echo 40 -90     0 | cs2cs EPSG:4979 EPSG:8734 -f "%.6f" ==> 2343275.132900 1213912.862240 107.186796</span></div><div><span style="font-family:monospace">echo 40 -90 10000 | cs2cs EPSG:4979 EPSG:8734 -f "%.6f" ==> 2343275.132900 1213912.862240 32915.520130</span></div><div><span style="font-family:monospace">(EPSG:8734 = <a title="View Details" href="https://beta.epsg.org/crs_8734/NAD83-Illinois-West-ftUS-NAVD88-height-ftUS.html?sessionkey=23506x2d03">NAD83 / Illinois West (ftUS) + NAVD88 height (ftUS)</a>)</span></div><div>I see there is no change in the x and y coordinates due to the change in z. If the altitude is not changing the horizontal projection, can I assume that the ellipsoidal height in EPSG:4979 is perpendicular to the projected CRS? (I know it is very unlikely).<br></div><div><br></div><div><span style="font-family:monospace"><span style="font-family:arial,sans-serif">Thanks</span><br></span></div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature">.___ ._ ..._ .. . ._.  .___ .. __ . _. . __..  ... .... ._ .__<br>Entre dos pensamientos racionales <br>hay infinitos pensamientos irracionales.<br><br></div></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, 1 Sep 2020 at 15:25, 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"><u></u>
<div style="font-family:"monospace";font-size:9pt;font-weight:400;font-style:normal">
<p style="margin:0px;text-indent:0px">Javier,</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">> Does exist any EPSG or WKT for a Vertical CRS that means "ellipsoidal</p>
<p style="margin:0px;text-indent:0px">> height"?</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">No, and this is by design. Ellipsoidal height doesn't make sense, from a referencing point of view, if the longitude,latitude of the point where you measure it is not specified as well, hence the need for a Geographic3D CRS. This is contrary to a gravity-based height, where if someone gives you just the the elevation at 2 points, you can for example know from which point towards which one water would flow.</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">See <a href="http://docs.opengeospatial.org/is/18-010r7/18-010r7.html" target="_blank">http://docs.opengeospatial.org/is/18-010r7/18-010r7.html</a></p>
<p style="margin:0px;text-indent:0px">"""</p>
<p style="margin:0px;text-indent:0px">3.1.19 ellipsoidal height</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">    distance of a point from the reference ellipsoid along the perpendicular from the reference ellipsoid to this point, positive if upwards or outside of the reference ellipsoid</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">    Note 1 to entry: Only used as part of a three-dimensional ellipsoidal coordinate system or as part of a three-dimensional Cartesian coordinate system in a three-dimensional projected coordinate reference system, but never on its own.</p>
<p style="margin:0px;text-indent:0px">"""</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">Very early versions of the EPSG dataset (like 10 years ago or more) add codes for vertical CRS with ellipsoidal height, but this was removed in later versions. Some remains of that are still found in the GeoTIFF 1.0 specification that hard-coded such dataset, but they have disappeared since, and are no longer allowed in GeoTIFF 1.1. I've seen recently also some WKT-1 string for LAS files that were using a ellipsoidal vertical CRS as part of a compound CRS, but this is no longer considered as geodetic-savy practice.</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">> I don't want to "promoteTo3D" the 2D CRS  adding the third axis, because it</p>
<p style="margin:0px;text-indent:0px">> removes the authority.</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">For projected CRS, yes, since there is currently no Projected 3D CRS in EPSG</p>
<p style="margin:0px;text-indent:0px">For geographic CRS, to the extent where there's a 3D version of the one you promote to 2D registered, promoteTo3D should attach the code of the 3D CRS:</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">$ projinfo EPSG:4326 --3d</p>
<p style="margin:0px;text-indent:0px">[....snip ...]</p>
<p style="margin:0px;text-indent:0px">    ID["EPSG",4979]]</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">Even</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">-- </p>
<p style="margin:0px;text-indent:0px">Spatialys - Geospatial professional services</p>
<p style="margin:0px;text-indent:0px"><a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a></p></div></blockquote></div>