<div dir="ltr"><div>Thanks, I "knew" that there is some magic :) Gonna try that.<br><br></div>Stefan<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Aug 12, 2015 at 9:22 PM, Smith, Michael ERDC-RDE-CRREL-NH <span dir="ltr"><<a href="mailto:Michael.Smith@erdc.dren.mil" target="_blank">Michael.Smith@erdc.dren.mil</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word;color:rgb(0,0,0);font-size:14px;font-family:Calibri,sans-serif">
<div>
<div>
<div>Stefan,</div>
<div><br>
</div>
<div>You can actually do that with the pdal split command. You can create a tile index and then feed that to the split command and set the length parameter on the split command.</div>
<div><br>
</div>
<div>Mike</div>
<div><br>
</div>
<div><br>
</div>
<div>
<div>
<div>----</div>
<div>Michael Smith</div>
</div>
<div>US Army Corps</div>
<div>Remote Sensing GIS/Center</div>
<div><a href="mailto:michael.smith@usace.army.mil" target="_blank">michael.smith@usace.army.mil</a></div>
<div><br>
</div>
</div>
</div>
</div>
<div><br>
</div>
<span>
<div style="font-family:Calibri;font-size:11pt;text-align:left;color:black;BORDER-BOTTOM:medium none;BORDER-LEFT:medium none;PADDING-BOTTOM:0in;PADDING-LEFT:0in;PADDING-RIGHT:0in;BORDER-TOP:#b5c4df 1pt solid;BORDER-RIGHT:medium none;PADDING-TOP:3pt">
<span style="font-weight:bold">From: </span><<a href="mailto:pdal-bounces@lists.osgeo.org" target="_blank">pdal-bounces@lists.osgeo.org</a>> on behalf of Stefan Ziegler <<a href="mailto:stefan.ziegler.de@gmail.com" target="_blank">stefan.ziegler.de@gmail.com</a>><br>
<span style="font-weight:bold">Date: </span>Wednesday, August 12, 2015 at 3:02 PM<br>
<span style="font-weight:bold">To: </span>Andrew Bell <<a href="mailto:andrew.bell.ia@gmail.com" target="_blank">andrew.bell.ia@gmail.com</a>><br>
<span style="font-weight:bold">Cc: </span>"<a href="mailto:pdal@lists.osgeo.org" target="_blank">pdal@lists.osgeo.org</a>" <<a href="mailto:pdal@lists.osgeo.org" target="_blank">pdal@lists.osgeo.org</a>><br>
<span style="font-weight:bold">Subject: </span>[EXTERNAL] Re: [pdal] PDAL Python bindings<br>
<span style="font-weight:bold">Resent-From: </span>Michael Smith <<a href="mailto:michael.smith@usace.army.mil" target="_blank">michael.smith@usace.army.mil</a>><br>
</div><div><div class="h5">
<div><br>
</div>
<blockquote style="BORDER-LEFT:#b5c4df 5 solid;PADDING:0 0 0 5;MARGIN:0 0 0 5">
<div>
<div>
<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>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>
<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><font color="#888888">
<div><br>
</div>
</font></span></div>
<span><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>
</div>
</div>
</blockquote>
</div></div></span>
</div>
</blockquote></div><br></div>