[pdal] Trouble reading projection from some 3DEP tiles

[Jim Klassen]

Yes, separate issues, so separate email threads.

On 9/11/24 16:12, Andrew Bell wrote:
> 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
>
>