[pdal] filters.reprojection with WKT and no EPSG

[Kirk Waters - NOAA Federal]

Dear PDAL (and maybe GDAL) gurus,
I'm trying to figure out a reprojection issue where the output projection
lacks an EPSG code. I think it should work with a WKT, but I get this error:
(pdal pipeline filters.reprojection Error) GDAL failure (1) No inverse
operation

The projection is from geographic NAD83(CORS96) to UTM zone 10 in the same
system.
Here's a simplification of the pipeline I'm trying to run:
[
  {
    "type": "readers.las",
    "filename": "
https://coast.noaa.gov/htdata/lidar1_z/tmp/20180316_WA_Olympic_Peninsula_C1_2017_46123H5107_1.copc.laz
"
  },
  {
    "type": "filters.reprojection",
    "in_srs": "EPSG:6783+5703",
    "out_srs": "COMPD_CS[\"NAD83(CORS96) / UTM Zone 10N + NAVD88
height\",PROJCS[\"NAD83(CORS96) / UTM Zone
10N\",GEOGCS[\"NAD83(CORS96)\",DATUM[\"NAD83_Continuously_Operating_Reference_Station_1996\",SPHEROID[\"GRS
1980\",6378137,298.257222101,AUTHORITY[\"EPSG\",\"7019\"]],AUTHORITY[\"EPSG\",\"1133\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"6782\"]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"latitude_of_origin\",0],PARAMETER[\"central_meridian\",-123],PARAMETER[\"scale_factor\",0.9996],PARAMETER[\"false_easting\",500000],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"Easting\",EAST],AXIS[\"Northing\",NORTH]],VERT_CS[\"NAVD88
height\",VERT_DATUM[\"North American Vertical Datum
1988\",2005],UNIT[\"metre\",1],AXIS[\"Gravity-related
height\",UP],AUTHORITY[\"EPSG\",\"5703\"]]]"
  },
  {
    "type": "writers.copc",
    "pipeline": true,
    "forward": "header, vlr",
    "offset_x": "auto",
    "offset_y": "auto",
    "offset_z": "auto",
    "extra_dims": "all",
    "filename": "pipeline10.copc.laz"
  }
]

If I switch the output srs to be "EPSG:6339+5703", I get no errors (6339 is
UTM zone 10 for NAD83(2011)), though I don't think it actually shifted from
CORS96 to 2011. GDAL osr.SpatialReference doesn't dislike the CORS96 utm 10
wkt, but it throws an exception if I try to AutoIndentifyEPSG().

Is this the intended behavior or am I doing something wrong? I'm guessing
it's related to the lack of an EPSG code, but perhaps I've gone down the
wrong path.

Note that the input file in the example uses EPSG:6782, but I simplified
the pipeline by removing the geoid application part.

Thanks for any pointers.

Kirk Waters, PhD
NOAA Office for Coastal Management
Applied Sciences Program
coast.noaa.gov/digitalcoast
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20260324/437faded/attachment.htm>