[pdal] Specifying SRS of output GeoTIFF

[Howard Butler]

The output units of the nc_ground.tiff 'resolution' parameter is now going to be in degrees. A value of '1.0' is probably not going to be what you want.

Write your raster in UTM and then use gdalwarp to put it into whatever coordinate system you want.

> On Apr 12, 2021, at 5:35 PM, Andreas Yankopolus <andreas at yank.to> wrote:
> 
> 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
> _______________________________________________
> pdal mailing list
> pdal at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/pdal