[pdal] [EXTERNAL] Re: PDAL Python bindings

[Stefan Ziegler]

Thanks, I "knew" that there is some magic :) Gonna try that.

Stefan

On Wed, Aug 12, 2015 at 9:22 PM, Smith, Michael ERDC-RDE-CRREL-NH <
Michael.Smith at erdc.dren.mil> wrote:

> 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> on behalf of Stefan Ziegler <
> stefan.ziegler.de at gmail.com>
> Date: Wednesday, August 12, 2015 at 3:02 PM
> To: Andrew Bell <andrew.bell.ia at gmail.com>
> Cc: "pdal at lists.osgeo.org" <pdal at lists.osgeo.org>
> Subject: [EXTERNAL] Re: [pdal] PDAL Python bindings
> Resent-From: Michael Smith <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>
> wrote:
>
>> On Wed, Aug 12, 2015 at 12:04 PM, Stefan Ziegler <
>> 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
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20150812/e372355a/attachment.html>