<div dir="ltr">One trick I like regarding accessing overview levels is to use the open option to snag a handle direct to the overview level (see Even's post <a href="http://erouault.blogspot.com/2014/10/warping-overviews-and-warped-overviews.html">http://erouault.blogspot.com/2014/10/warping-overviews-and-warped-overviews.html</a>); then you can run polygonize against a much smaller dataset.  If it's a COG presumably it has internal overviews.</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Jan 13, 2020 at 12:22 PM Brian <<a href="mailto:mulcahy.brianp@gmail.com">mulcahy.brianp@gmail.com</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">So I am not sure if this is the right place or not. Feel free to remove if it is not. <div><br></div><div>The problem:  I am needing to see if an extent (vector polygon, 4 vertices) has any interaction with a raster stored in the cog format.</div><div><br></div><div>I thought about using raster_polygonize and running a intersect but that is going to be way too slow and I don't need the detail at which the raster is stored in. </div><div>So I was thinking of some alternatives.</div><div>Is there anyway to output just the polygon of the raster at a certain overview level? </div><div>Is there any way to run intersects on a raster without converting it? </div><div>I know I can retrieve the extent of the raster and quickly determine if an extent has any interaction with that raster extent. I also know that I can send points(x,y) and retrieve the value of the raster at that point is there any way to manipulate /game this? Maybe sending the extent as individual points with added points in between the vertices and just running multiple value at location operations? For my need, some false positives would be acceptable. If it helps raster data is of rivers and inundations.</div><div><br></div><div>Has anyone else come up with some type of process for this? </div></div>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a></blockquote></div>