<div dir="ltr"><div>Instead of "+init=EPSG:32611" or "+proj=utm +zone=11 +type=crs", you might try the <a href="https://spatialreference.org/ref/epsg/32611/proj4/">fully expanded PROJ string</a>:</div><div><br></div><div><blockquote style="margin:0 0 0 40px;border:none;padding:0px"><div>+proj=utm +zone=11 +ellps=WGS84 +datum=WGS84 +units=m +no_defs<br></div></blockquote></div><div><br></div><div>So your pipeline would be:</div><div><br></div><div><blockquote style="margin:0px 0px 0px 40px;border:none;padding:0px"><div><i>{"pipeline": [{"type" : "readers.las",</i></div><div><i> "filename": "input_Ellipsoidal.laz"},</i></div><div><i> {"type":"filters.reprojection",</i></div><div><i> "in_srs":"EPSG:32611",</i></div><div><i> "out_srs":"</i><span style="background-color:rgb(255,255,0)">+proj=utm +zone=11 +ellps=WGS84 +datum=WGS84 +units=m +no_defs</span><i> +geoidgrids=g2003conus.gtx"},</i></div><div><i> {"type" : "writers.las",</i></div><div><i> "filename": "output_Geoid.laz",</i></div><div><i> "compression": "laszip",</i></div><div><i> "a_srs": "EPSG:32611+5703"</i></div><div><i> }]}</i></div></blockquote><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Apr 28, 2023 at 10:09 AM Matt Beckley <<a href="mailto:matt.beckley@earthscope.org">matt.beckley@earthscope.org</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">Hello,<div><br></div><div>In older versions of PDAL I used to be able to run a pipeline like:</div><div><br></div><blockquote style="margin:0px 0px 0px 40px;border:none;padding:0px"><div><i>{"pipeline": [{"type" : "readers.las",</i></div><div><i> "filename": "input_Ellipsoidal.laz"},</i></div><div><i> {"type":"filters.reprojection",</i></div><div><i> "in_srs":"EPSG:32611",</i></div><div><i> "out_srs":"+init=EPSG:32611 +geoidgrids=g2003conus.gtx"},</i></div><div><i> {"type" : "writers.las",</i></div><div><i> "filename": "output_Geoid.laz",</i></div><div><i> "compression": "laszip",</i></div><div><i> "a_srs": "EPSG:32611+5703"</i></div><div><i> }]}</i></div></blockquote><div><br></div><div>This would run fine, and the output would have the geoid properly applied.</div><div><br></div><div>However, newer versions of PDAL (e.g. 2.5.3) uses newer version of PROJ, which no longer likes its string in the form, "+init=epsg:XXXX". It will give the warning:</div>"Warning 1: +init=epsg:XXXX syntax is deprecated. It might return a CRS with a non-EPSG compliant axis order."
<div><br></div><div>I adjusted the PROJ string to remove this convention:</div><div><br></div><div><div><i>{"pipeline": [{"type" : "readers.las",</i></div><div><i> "filename": "input_Ellipsoidal.laz"},</i></div><div><i> {"type":"filters.reprojection",</i></div><div><i> "in_srs":"EPSG:32611",</i></div><div><i> "out_srs":"</i>+proj=utm +zone=11 +type=crs<i> +geoidgrids=g2003conus.gtx"},</i></div><div><i> {"type" : "writers.las",</i></div><div><i> "filename": "output_Geoid.laz",</i></div><div><i> "compression": "laszip",</i></div><div><i> "a_srs": "EPSG:32611+5703"</i></div><div><i> }]}</i></div></div><div><br></div><div>This runs without warning or errors, but the output grid <b>does not</b> have the geoid applied. Most of the documentation on this type of operation I've seen still refers to the older PROJ syntax of, "+init=epsg:XXXX".</div><div><br></div><div>Is there a way to apply a geoid with the new PROJ syntax, and hopefully I am just missing something with the proper syntax to do this operation on the newer versions of PDAL?</div><div><br></div><div>Thanks!</div></div>
_______________________________________________<br>
pdal mailing list<br>
<a href="mailto:pdal@lists.osgeo.org" target="_blank">pdal@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/pdal" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/pdal</a><br>
</blockquote></div></div>