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