[pdal] Trouble reading projection from some 3DEP tiles

Jim Klassen klassen.js at gmail.com
Wed Sep 11 11:36:16 PDT 2024


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
-------------- next part --------------
A non-text attachment was scrubbed...
Name: USGS_LPC_MN_UpperMissRiver_B22_334_5271.laz.json
Type: application/json
Size: 26476 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20240911/49e54b99/attachment-0001.json>


More information about the pdal mailing list