<div dir="ltr"><div>Dear PDAL (and maybe GDAL) gurus,</div><div>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:</div><div>(pdal pipeline filters.reprojection Error) GDAL failure (1) No inverse operation</div><div><br></div><div>The projection is from geographic NAD83(CORS96) to UTM zone 10 in the same system.</div><div>Here's a simplification of the pipeline I'm trying to run:</div><div>[<br>  {<br>    "type": "readers.las",<br>    "filename": "<a href="https://coast.noaa.gov/htdata/lidar1_z/tmp/20180316_WA_Olympic_Peninsula_C1_2017_46123H5107_1.copc.laz">https://coast.noaa.gov/htdata/lidar1_z/tmp/20180316_WA_Olympic_Peninsula_C1_2017_46123H5107_1.copc.laz</a>"<br>  },<br>  {<br>    "type": "filters.reprojection",<br>    "in_srs": "EPSG:6783+5703",<br>    "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\"]]]"<br>  },<br>  {<br>    "type": "writers.copc",<br>    "pipeline": true,<br>    "forward": "header, vlr",<br>    "offset_x": "auto",<br>    "offset_y": "auto",<br>    "offset_z": "auto",<br>    "extra_dims": "all",<br>    "filename": "pipeline10.copc.laz"<br>  }<br>]</div><div><br></div><div>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(). </div><div><br></div><div>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.</div><div><br></div><div>Note that the input file in the example uses EPSG:6782, but I simplified the pipeline by removing the geoid application part.</div><div><br></div><div>Thanks for any pointers.</div><div><br></div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><font face="arial, helvetica, sans-serif">Kirk Waters, PhD </font></div><div><font face="arial, helvetica, sans-serif">NOAA Office for Coastal Management<br></font></div><div><font face="arial, helvetica, sans-serif">Applied Sciences Program      </font></div><div><div><font face="arial, helvetica, sans-serif"><a href="http://coast.noaa.gov/digitalcoast" target="_blank">coast.noaa.gov/digitalcoast</a></font></div></div><div><br></div><div><br></div></div></div></div></div>