如何在python中先合并然后裁剪光栅

如何在python中先合并然后裁剪光栅,python,rasterio,Python,Rasterio,目前,我必须从大量的光栅创建大量的平铺(由一个形状多边形定义) 每个磁贴都是一个多边形,光栅只是我磁盘上的tiff文件 现在,我首先找到与平铺相交的光栅的所有文件名,并将它们放入列表中 然后,我应用以下程序 images = [] for file in files: images.append(rasterio.open( file) mosaic = merge(images) 剩下的步骤是裁剪“马赛克”光栅和定义平铺的形状多边形 但是,通过执行合并步骤,马赛克变成了一个nump

目前,我必须从大量的光栅创建大量的平铺(由一个形状多边形定义)

每个磁贴都是一个多边形,光栅只是我磁盘上的tiff文件

现在,我首先找到与平铺相交的光栅的所有文件名,并将它们放入列表中

然后,我应用以下程序

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

mosaic =  merge(images)
剩下的步骤是裁剪“马赛克”光栅和定义平铺的形状多边形


但是,通过执行合并步骤,马赛克变成了一个numpy数组,因此我无法对其应用mask方法。

您可以将rasterio.merge.merge()的边界参数用于特定瓷砖的形状多边形

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

可以将rasterio.merge.merge()的bounds参数用于特定平铺的形状多边形

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

谢谢你,这是一种优雅的方式!谢谢你,这是一种优雅的方式!