You could get the "bbox" list from mercantile.tiles and then get the bounds for each tile using mercantile.bounds and then check the intersection with your irregular polygon, this would be fairly slow depending on the bounds / zoom. I'm not aware of any other libraries.
_______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
