В настоящее время мне нужно создать множество тайлов (определяемых фигурным мультиполигоном) из большого количества растров.
Каждая плитка представляет собой многоугольник, а растры - это просто файлы TIFF на моем диске.
Прямо сейчас я сначала нахожу все имена файлов растров, пересекающих плитку, и помещаю их в список.
Затем я применяю следующую процедуру
images = []
for file in files:
images.append(rasterio.open( file)
mosaic = merge(images)
Остающийся шаг - обрезать растр мозаики с помощью фигурного многоугольника, определяющего плитку.
Однако, выполнив этап слияния, мозаика стала массивом numpy, поэтому я не могу применить к нему метод маски.