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 http://erouault.blogspot.com/2014/10/warping-overviews-and-warped-overviews.html); then you can run polygonize against a much smaller dataset. If it's a COG presumably it has internal overviews.
On Mon, Jan 13, 2020 at 12:22 PM Brian <[email protected]> wrote: > So I am not sure if this is the right place or not. Feel free to remove if > it is not. > > 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. > > 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. > So I was thinking of some alternatives. > Is there anyway to output just the polygon of the raster at a certain > overview level? > Is there any way to run intersects on a raster without converting it? > 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. > > Has anyone else come up with some type of process for this? > _______________________________________________ > gdal-dev mailing list > [email protected] > https://lists.osgeo.org/mailman/listinfo/gdal-dev
_______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
