[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