<div dir="ltr">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.<br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Aug 14, 2020 at 4:05 PM Connor Manning <<a href="mailto:connor@hobu.co" target="_blank">connor@hobu.co</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Sounds like you could do this in raster land since you are interested in 1x1 meter cells.  Perhaps the GDAL writer (<a href="https://pdal.io/stages/writers.gdal.html" target="_blank">https://pdal.io/stages/writers.gdal.html</a>) with a `resolution` of 1 and an `output_type` of "max".  And then to reinterpret your output as a point cloud: the GDAL reader (<a href="https://pdal.io/stages/readers.gdal.html" target="_blank">https://pdal.io/stages/readers.gdal.html</a>).<br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Aug 14, 2020 at 3:12 PM Alison Cramer <<a href="mailto:alisoncramer@nevada.unr.edu" target="_blank">alisoncramer@nevada.unr.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><span style="color:rgb(0,0,0);font-family:Verdana,Geneva,Helvetica,Arial,sans-serif;font-size:13.44px">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? </span><br></div></div>
_______________________________________________<br>
pdal mailing list<br>
<a href="mailto:pdal@lists.osgeo.org" target="_blank">pdal@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/pdal" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/pdal</a></blockquote></div>
_______________________________________________<br>
pdal mailing list<br>
<a href="mailto:pdal@lists.osgeo.org" target="_blank">pdal@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/pdal" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/pdal</a></blockquote></div>