[pdal] Trouble reading projection from some 3DEP tiles

Andrew Bell andrew.bell.ia at gmail.com
Wed Sep 11 14:12:59 PDT 2024


Is this a different issue than the one you reported earlier in the day?

On Wed, Sep 11, 2024 at 3:49 PM Jim Klassen via pdal
<pdal at lists.osgeo.org> wrote:
>
> I've noticed that for some of the 3DEP tiles that PDAL doesn't calculate
> a 4326 BBOX.  This seems to be related to Proj not liking the how the
> projection is specified in the LAZ file.  If I force the projection as
> EPSG:6344+5703 in the reader then "pdal info" works as expected.  My
> main question is how can I tell if the issue is with the CRS in the LAZ
> file or the problem is with PDAL/Proj not being flexible enough in
> parsing the CRS?
>
> Steps to reproduce:
>      wget
> https://rockyweb.usgs.gov/vdelivery/Datasets/Staged/Elevation/LPC/Projects/MN_UpperMissRiver_B22/MN_UpperMissRiver_2_B22/LAZ/USGS_LPC_MN_UpperMissRiver_B22_334_5271.laz
>
>      pdal info --all USGS_LPC_MN_UpperMissRiver_B22_334_5271.laz
>
> Observations:
>      Under "stac" there is a error with message "Failed to create STAC
> Feature with missing key. 'EPSG:4326'".
>
>      Under stats.bbox there is no 4326 bbox calculated.
>
> And with --verbose 8 I get:
>
> (PDAL Debug) Debugging...
> (pdal info Debug) Executing pipeline in stream mode.
> (pdal info filters.stats Debug) GDAL debug:
> OGRSpatialReference::Validate: Coordinate system of GeographicCRS in the
> WKT definition is different from the one of the authority. Unsetting the
> identifier to avoid confusion
> (pdal info filters.stats Debug) GDAL debug:
> OGRSpatialReference::Validate: Parsing error : syntax error, unexpected
> AUTHORITY, expecting number. Error occurred around:
> UM["North American Vertical Datum
> 1988",AUTHORITY["EPSG","5103"]],UNIT["Meter",1
>                                          ^
> (pdal info filters.hexbin Debug) GDAL debug:
> OGRSpatialReference::Validate: Coordinate system of GeographicCRS in the
> WKT definition is different from the one of the authority. Unsetting the
> identifier to avoid confusion
> (pdal info filters.hexbin Debug) GDAL debug:
> OGRSpatialReference::Validate: Parsing error : syntax error, unexpected
> AUTHORITY, expecting number. Error occurred around:
> UM["North American Vertical Datum
> 1988",AUTHORITY["EPSG","5103"]],UNIT["Meter",1
> _______________________________________________
> pdal mailing list
> pdal at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/pdal



-- 
Andrew Bell
andrew.bell.ia at gmail.com


More information about the pdal mailing list