目前,我必须从大量的栅格中创建大量的瓦片(由一个形状良好的多多边形定义)。
每个瓦片都是一个多边形,而栅格只是我磁盘上的tiff文件。
现在,我首先找到与瓦片相交的所有栅格的文件名,并将它们放在一个列表中。
然后,我应用以下过程
images = []
for file in files:
images.append(rasterio.open( file)
mosaic = merge(images)剩下的步骤是用定义瓷砖的形状多边形来裁剪“马赛克”栅格。
但是,通过执行合并步骤,马赛克变成了一个numpy数组,因此我无法对其应用掩码方法。
发布于 2018-09-30 23:16:10
可以对特定平铺的形状多边形使用rasterio.merge.merge()的bounds参数。
array_clipped_mosaic, out_transform = rasterio.merge.merge(datasets, bounds = *polygon.bounds)https://stackoverflow.com/questions/52569641
复制相似问题