
<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Connor—<br class=""><div><br class=""><blockquote type="cite" class=""><div class="">You can intersect your AOI with the GeoJSON boundaries to determine your overlaps: <a href="https://raw.githubusercontent.com/hobu/usgs-lidar/master/boundaries/resources.geojson" class="">https://raw.githubusercontent.com/hobu/usgs-lidar/master/boundaries/resources.geojson</a></div><div class=""><div dir="ltr" class=""><div class=""><br class=""></div><div class="">This is updated automatically daily here: <a href="https://github.com/hobuinc/usgs-lidar/tree/master/boundaries" class="">https://github.com/hobuinc/usgs-lidar/tree/master/boundaries</a></div></div></div></blockquote><div><br class=""></div><div>This is perfect. It will be a straightforward task to write a Python script that uses GDAL to perform the overlap calculations to determine the relevant collects.</div></div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">Kirk—</div><div class=""><br class=""></div><div class=""><blockquote type="cite" class="">I could have this wrong, but I think the json file that has the boundaries for the EPT datasets can be found at <a href="https://github.com/hobuinc/usgs-lidar/blob/master/boundaries/boundaries.topojson" class="">https://github.com/hobuinc/usgs-lidar/blob/master/boundaries/boundaries.topojson</a>. That would allow you to figure out which datasets are within an AOI. Standing up a service that would let you specify an AOI and then spit back all the points, merging the different collects as it goes, has some inherent problems. It's technically feasible; the issue is that some things shouldn't be mixed. For example, pre and post-event collections for hurricanes.</blockquote><br class=""></div><div class="">Point taken on being careful with the merges. I’m using the LPC data to get a rough idea of average building/tree heights for radio propagation calculations, converting the LPC data to surface-height raster files with resolutions between 3m to 30m.</div><div class=""><br class=""></div><div class="">To date I’ve been converting LPC tiles to a raster tiles in parallel, then merging the raster tiles with GDAL as a sequential operation. I’m often processing several hundred square km of data at once, and telling PDAL to load all that data into memory at once doesn’t work out. Maybe there’s a trick I’m missing?</div><div class=""><br class=""></div><div class="">Thanks,</div><div class=""><br class=""></div><div class="">Andreas</div><div class=""><br class=""></div></body></html>