Heyo,
I have a script that uses the SentinelSat API to download images that are then processes using Snappy. Now I have found an edge case where the footprint I am searching in overlaps (just barely) with a Sentinel-1 image. So that image is downloaded. However it is such a small overlap that when I try to subset in Snappy it gives me an empty image.
I have found an awkward solution I convert the image footprint and the subsetting footprint to shapely and only download if the area over intersection is large enough. I am wondering if there is a cleaner method of doing this?
p_subset = shapely.wkt.loads(subset)
for product in products:
for i, r in product.iterrows():
P_img = shapely.wkt.loads(r.footprint)
if P_img.intersection(p_subset).area > 0.001:
api.download(i)
Thanks!
Z