<div dir="ltr"><div>If you want retile a mosaic from 1 sqkm tiles to e.g. 750 m * 750 m tiles, I would use python for calculating the new bounding boxes and invoke pdal. Just the same as with GDAL. I can also make system calls but I prefer gdal/ogr python bindings.<br><br></div>Stefan<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Aug 12, 2015 at 8:36 PM, Andrew Bell <span dir="ltr"><<a href="mailto:andrew.bell.ia@gmail.com" target="_blank">andrew.bell.ia@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span class="">On Wed, Aug 12, 2015 at 12:04 PM, Stefan Ziegler <span dir="ltr"><<a href="mailto:stefan.ziegler.de@gmail.com" target="_blank">stefan.ziegler.de@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>Use case 1: We have approx. 1000 tiles (one sqkm each) and need to reproject them. Afterwards we nee to adjust the bounding boxes to "nice" coordinates again. Let's say the source bbox is something like 620'000/240'000, 621'000/241'000. After reprojecting it will be 2'620'000.55/1'240'000.34, 2'621'000.58/1'241'000.37. But we want a bbox like this: 2'620'000/1'240'000, 2'621'000/1'241'000. I can achieve this with available pdal features (tindex, translate etc.) and some python script for the loop with system calls (for pdal).<br></div></div></div></div></blockquote><div><br></div></span><div>Are you just saying that you want to write python instead of some shell to invoke pdal?  Not sure how a python API would help with that.  This is probably best done with a pipeline (XML) that you could invoke with `pdal pipeline`, but maybe I'm not understanding.</div><span class=""><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>Use case 2: Some small (wsgi) webservice like requesting lidar data with a bounding/polygon as parameter value. Again I could do this wie os.system().<br></div></div></div></div></blockquote><div><br></div></span><div>There are lots of tools out there to do this in various ways.  Hobu Inc. makes greyhound (not sure on its release status), which can serve a single file of data based on a query, but we'll have something much fancier and more efficient shortly.  I still don't understand what kind of an API would be beneficial for this.  You may have to explain a bit more.</div><span class="HOEnZb"><font color="#888888"><div><br></div></font></span></div><span class="HOEnZb"><font color="#888888">-- <br><div>Andrew Bell<br><a href="mailto:andrew.bell.ia@gmail.com" target="_blank">andrew.bell.ia@gmail.com</a></div>
</font></span></div></div>
</blockquote></div><br></div>