[pdal] upsample

Andrew Bell andrew.bell.ia at gmail.com
Fri Aug 14 15:36:20 PDT 2020


Alison,

This would be simple to implement from scratch. I could write you a custom
filter if that's what you want. Contact me off list if you'd like to
discuss.

On Fri, Aug 14, 2020, 5:22 PM Bradley Chambers <brad.chambers at gmail.com>
wrote:

> Good call with GDAL reader/writer, but this only gives you 1 point per
> cell. Another partial solution would be to use filters.locate, searching
> for the max in the Z dimension (see the docs on how to specify this). You'd
> have to split using filters.splitter first, and would have the same issue
> in the end of having a single point per cell. I'm just not sure we have a
> filter that could take that value and apply it to the rest of the points in
> the cell. You could write a Python script and call that with filters.python.
>
> On Fri, Aug 14, 2020 at 4:05 PM Connor Manning <connor at hobu.co> wrote:
>
>> Sounds like you could do this in raster land since you are interested in
>> 1x1 meter cells.  Perhaps the GDAL writer (
>> https://pdal.io/stages/writers.gdal.html) with a `resolution` of 1 and
>> an `output_type` of "max".  And then to reinterpret your output as a point
>> cloud: the GDAL reader (https://pdal.io/stages/readers.gdal.html).
>>
>> On Fri, Aug 14, 2020 at 3:12 PM Alison Cramer <
>> alisoncramer at nevada.unr.edu> wrote:
>>
>>> I have a point cloud of an abandoned mine land with a waste rock pile.
>>> The waste rock pile has rills (soil removed by erosion) caused by large
>>> rainfall events.  This causes differences in elevation (z values) within 1
>>> square meter.  I would like to perform a moving window of 1m along the
>>> waste rock pile.  This would hopefully give me "cells" of 1m by 1m.  My
>>> goal is to find the maximum z value in every 1m cell.  Once I find the max
>>> z value in each 1m cell, I would like to replace the current z values with
>>> the specific maximum z value calculated for each cell.  Then export to a
>>> new point cloud.  Does anyone have an idea as to how I would go about this
>>> in pdal?
>>> _______________________________________________
>>> pdal mailing list
>>> pdal at lists.osgeo.org
>>> https://lists.osgeo.org/mailman/listinfo/pdal
>>
>> _______________________________________________
>> pdal mailing list
>> pdal at lists.osgeo.org
>> https://lists.osgeo.org/mailman/listinfo/pdal
>
> _______________________________________________
> pdal mailing list
> pdal at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/pdal
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20200814/4b4ba0c5/attachment-0001.html>


More information about the pdal mailing list