<div dir="ltr">Hi there,<div><br></div><div>I'm currently working on a LIDAR data project where we classify "high vegetation" and buildings using PDAL. We convert them to rasters before we send them downstream.</div><div><br></div><div>Now, we are looking to combine the rasters generated with satellite imagery. For this project, we believe it's important to have an 1:1 mapping between the pixels in the different sources.</div><div><br></div><div>To align our satellite images with the underling projection, we use the "-tap" option in gdalwarp. However, no such feature exists for the GDAL-writer in PDAL. Which means the pixels are offset by up to half a meter (for UTM and 1.0 scale).</div><div><br></div><div>I think we can overcome the above issue by using PDAL to generate a much higher resolution TIF file (say 0.1) and then resample with gdalwarp and the "-tap" option. Or even just keep the resolution, just use gdalwarp to resample.</div><div><br></div><div>I was wondering if a similar "-tap" option belongs in PDAL, and how that would look like? I initially thought it was argument passed to the driver, but AFAIK the pixel alignment happens before it is passed to the driver. If I understand it correctly, it would then be an extra option to be passed to the "writers.gdal" component.</div><div><br></div><div>I'd love to contribute, but I'm no C++ developer.</div><div><br></div><div>Håkon</div></div>