如何在python中首先合並然后裁剪柵格

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)