m пытается найти ограничивающие полигоны всех целых в двумерной сетке с большим значением отсутствия данных (1e6). Я'у нас есть список дыр, работающих с помощью Scipy 'с меткой. Не окунаясь в гдальs полигонализировать, есть ли простой способ генерировать ограничивающие полигоны? Я вижу, что есть matplotlib.pylab.contour, но он пытается нарисовать сюжет, который я действительно не делаюне хочу Любая рекомендация о том, как получить ограничивающие полигоны для каждой метки (желательно с возможностью немного упростить полигоны, если это возможно)? Я'Я уверен, что могу написать что-то, что будет проходить границы каждой помеченной дыры, но есть ли что-то, что уже существует?

from osgeo import gdal
from scipy import ndimage

dem_file = gdal.Open('dem.tif')
dem = dem.file.GetRasterBand(1).ReadAsArray()

# Get a binary image of the no-data regions.  The no-data value is large
bin = dem > 9e5

# Find all the wholes.  Anything with a label > 0.
labels, num_labels = ndimage.measurements.label(bin)
num_labels
1063

# The hole's label and size. Skip 0 as that label has all the valid data.
holes = [(label, sum(labels==label)) for label in range(1, num_labels)]
holes[:3]
[(1, 7520492),
 (2, 1),
 (3, 1),]

например вместо того, чтобы вести счетя ищу границы всех этих белых областей, как показано в qgis, что было сделано с помощью gdal_polygonalize.py.I '