<div dir="ltr"><br><div>There are many changes that are coming shortly in the EPT reader for version 2.1 that should reduce the number of points retained.  In the more recent code, we're trying to store the minimum number of points necessary to meet the query.  That said, without knowing exactly what you're doing and what you were doing before, it's hard to comment.  I think I would wait and see if things are improved in version 2.1 before worrying too much about taking other steps.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Mar 16, 2020 at 8:27 AM Balázs Dukai <<a href="mailto:B.Dukai@tudelft.nl" target="_blank">B.Dukai@tudelft.nl</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">Hello,<br>
<br>
I'm using PDAL's C++ API (v2.0.1) to query an EPT and get the point <br>
cloud within a polygon. I do this for millions of polygons that are <br>
often adjacent (building footprints), and the EPT covers a whole <br>
country. The EPT is stored as LAZ.<br>
<br>
I noticed that there was a relatively large increase in memory use <br>
(cannot tell exactly, sorry ) after switching to PDAL+EPT. This made me <br>
wonder what exactly is read into memory when doing a polygon query from <br>
EPT? I made a drawing of what I think is happening , could someone <br>
confirm or reject this? See <br>
<a href="https://www.dropbox.com/s/ppa24o666ulunxc/ept_octree.jpg?dl=0" rel="noreferrer" target="_blank">https://www.dropbox.com/s/ppa24o666ulunxc/ept_octree.jpg?dl=0</a>.<br>
<br>
In the image you'll see an octree of an EPT, with the root cell on the <br>
top and the leaf that contains the polygon on the bottom. The red <br>
splotch is the polygon. The darker blue cubes are the octree cells of <br>
the EPT that I think are kept completely in memory in order to get all <br>
the points within that red polygon. Is that correct?<br>
<br>
Then if this is the case, how could I reduce the memory use in PDAL? And <br>
speed up extracting the polygon-extents from the EPT?<br>
<br>
Many thanks,<br>
<br>
Balázs Dukai<br>
<br>
_______________________________________________<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 clear="all"><div><br></div>-- <br><div dir="ltr">Andrew Bell<br><a href="mailto:andrew.bell.ia@gmail.com" target="_blank">andrew.bell.ia@gmail.com</a></div>