<div dir="auto">You can calculate and pass in your own "--writers.gdal.bounds" and get aligned pixels.<div dir="auto"><br></div><div dir="auto">(By doing this, I am able to round trip a DEM into a point cloud in PDAL and back into a DEM losslessly.)<br><div dir="auto"><br></div><div dir="auto"><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, May 7, 2020, 16:16 Håkon Åmdal <<a href="mailto:hakon@aamdal.com">hakon@aamdal.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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>
_______________________________________________<br>
pdal mailing list<br>
<a href="mailto:pdal@lists.osgeo.org" target="_blank" rel="noreferrer">pdal@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/pdal" rel="noreferrer noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/pdal</a></blockquote></div>