[pdal] Writing multi-band rasters

[Andreas Yankopolus]

Doug,

> Geotiff supports multiple bands. If your workflow is set up to generate a single band at a time, you can use gdalbuildvrt https://gdal.org/programs/gdalbuildvrt.html <https://gdal.org/programs/gdalbuildvrt.html> to make a virtual multiband raster of the individual bands , then gdal_translate , https://gdal.org/programs/gdal_translate.html#gdal-translate <https://gdal.org/programs/gdal_translate.html#gdal-translate> to convert your virtual multiband raster into a multiband geotiff.

I’m seeing that I could create one PDAL pipeline per dimension, so that if the target GeoTIFF had six dimensions, I’d have six JSON files writing six GeoTIFFs, then use GDAL to perform the merge and output a six-band GeoTIFF.

Is there a more elegant way to do this in PDAL, creating a single pipeline that processes six dimensions of data and writes it to a multi-band raster? My end goal is converting LAS files into rasters with separate bands for terrain height, surface height, and surface characteristics (ground, foliage, structure, etc).

I can get around well in GDAL, particularly the C++ API, but am new to PDAL Hence I don’t have a good sense on how it’s typically employed for more complex processing tasks.

Thanks,

Andreas
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20210407/2e234e94/attachment.html>