<html><head><meta http-equiv="Content-Type" content="text/html; charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On Jun 29, 2020, at 10:18 AM, Stephen V. Mather <<a href="mailto:svm@clevelandmetroparks.com" class="">svm@clevelandmetroparks.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class="">Hi All,</div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class=""><br class=""></div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class="">I am processing a dataset for which I have a digitized building footprints. It'd be even cooler if I had vegetation, but beggers != choosers.</div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class=""><br class=""></div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class="">In principle, I could use filters.overlay, but this is a massive dataset, so I am curious if it makes more sense to rasterize, tile things up, and use filters.colorization.</div></div></blockquote><div><br class=""></div><div>Define massive ;)</div><div><br class=""></div><div>I think your instinct to rasterize the footprints is the right one here. filters.overlay is going to do point-in-polygon for every point/poly combo. Most are quickly thrown out, but you're going to be checking every one. A raster mask is going to be much quicker, and you can control the resolution with gdal_rasterize.</div><br class=""><blockquote type="cite" class=""><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class=""><span style="font-size: 12pt;" class="">Happy to test and report back on what is fastest, but in case someone has a hunch, or an alternative recommended approach, I thought I would ask.</span></div></blockquote><br class=""></div><div>I'm curious in how different things might be, but I suspect the approaches will separate quite quickly for significantly sized polygon and point cloud sets.</div><div><br class=""></div><div>Howard</div><br class=""></body></html>