<div><div dir="auto">Let me try that tomorrow. You'll hear from me. Thanks for quick feedback!</div></div><div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 7 May 2020 at 17:24, Andrew Bell <<a href="mailto:andrew.bell.ia@gmail.com">andrew.bell.ia@gmail.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 <span style="color:rgb(32,33,36);font-size:0.875rem;letter-spacing:0.2px;font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif;white-space:nowrap">Håkon,</span><div><span style="outline:none"><span name="Håkon Åmdal" style="color:rgb(32,33,36);font-size:0.875rem;font-weight:bold;display:inline;vertical-align:top;letter-spacing:0.2px;line-height:20px"><br></span></span></div>I know nothing about gdalwarp, but I would think you could get perfect alignment by using the bounds option in `writers.gdal`.  Are you having a problem despite using this option?</div><br><div class="gmail_quote"></div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, May 7, 2020 at 5:16 PM Håkon Åmdal <<a href="mailto:hakon@aamdal.com" target="_blank">hakon@aamdal.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);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">pdal@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/pdal" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/pdal</a></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr">Andrew Bell<br><a href="mailto:andrew.bell.ia@gmail.com" target="_blank">andrew.bell.ia@gmail.com</a></div>
</blockquote></div></div>