Question

目前，我必须从大量栅格中创建很多图块（由有形的多边形定义）。

每个图块都是一个多边形，而栅格只是我磁盘上的tiff文件。

现在，我首先找到与图块相交的栅格的所有文件名，并将它们放在列表中。

然后我执行以下步骤

images = []
for file in files:
   images.append(rasterio.open( file)

mosaic =  merge(images)

剩下的步骤是使用定义图块的形状多边形来裁剪“马赛克”栅格。

但是，通过执行合并步骤，镶嵌图变成了一个numpy数组，因此我无法对其应用mask方法。

Answer 1

您可以将rasterio.merge.merge（）的bounds参数与特定图块的匀称多边形一起使用。

array_clipped_mosaic, out_transform = rasterio.merge.merge(datasets, bounds = *polygon.bounds)