[pdal] Specifying SRS of output GeoTIFF

Andreas Yankopolus andreas at yank.to
Mon Apr 12 15:35:15 PDT 2021


I’m trying to perform PDAL operations and write the resulting GeoTIFF in WGS84 (EPSG:4326

Using this sample data for example: http://fatra.cnr.ncsu.edu/foss4g2017/nc_tile_0793_016_spm.zip

This is my pipeline:

{
  "pipeline": [
    "nc_tile_0793_016_spm.las",
    {
      "type": "filters.range",
      "limits": "Classification[2:2]"
    },
    {
      "type": "filters.reprojection",
      "out_srs": "EPSG:4326"
    },
    {
      "type": "writers.gdal"
      "filename": "nc_ground.tiff",
      "gdaldriver": "GTiff",
      "dimension": "Z",
      "output_type": "idw",
      "resolution": "1.0",
    }
  ]
}

I don’t get a sensible TIFF with the reprojection filter in place. Without it, output is as expected in the point cloud’s SRS of ESPG:32119. Am I thinking about the reprojection filter incorrectly? Maybe it needs to be integrated into the writers.gdal block?

Thanks,

Andreas
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20210412/19e4247c/attachment.html>


More information about the pdal mailing list