[pdal] [EXTERNAL] Re: PDAL Python bindings

[Smith, Michael ERDC-RDE-CRREL-NH]

Stefan,

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.

Mike


----
Michael Smith
US Army Corps
Remote Sensing GIS/Center
michael.smith at usace.army.mil


From: <pdal-bounces at lists.osgeo.org<mailto:pdal-bounces at lists.osgeo.org>> on behalf of Stefan Ziegler <stefan.ziegler.de at gmail.com<mailto:stefan.ziegler.de at gmail.com>>
Date: Wednesday, August 12, 2015 at 3:02 PM
To: Andrew Bell <andrew.bell.ia at gmail.com<mailto:andrew.bell.ia at gmail.com>>
Cc: "pdal at lists.osgeo.org<mailto:pdal at lists.osgeo.org>" <pdal at lists.osgeo.org<mailto:pdal at lists.osgeo.org>>
Subject: [EXTERNAL] Re: [pdal] PDAL Python bindings
Resent-From: Michael Smith <michael.smith at usace.army.mil<mailto:michael.smith at usace.army.mil>>

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.

Stefan

On Wed, Aug 12, 2015 at 8:36 PM, Andrew Bell <andrew.bell.ia at gmail.com<mailto:andrew.bell.ia at gmail.com>> wrote:
On Wed, Aug 12, 2015 at 12:04 PM, Stefan Ziegler <stefan.ziegler.de at gmail.com<mailto:stefan.ziegler.de at gmail.com>> wrote:
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).

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.

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().

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.

--
Andrew Bell
andrew.bell.ia at gmail.com<mailto:andrew.bell.ia at gmail.com>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20150812/494ec77f/attachment.html>