用Rasterio和Fiona裁剪光栅

Question

我想用分水岭作为边界来剪辑光栅数据。我用Fiona和Rasterio试过了。

这里有我的代码：

import fiona
import rasterio
from rasterio.rio.clip import clip
with fiona.open("oreto_bacino2.shp", "r") as shapefile:
    geoms = [feature["geometry"] for feature in shapefile]

with rasterio.open("cn.asc") as src:
    out_image, out_transform = clip (geoms, src, crop=True)
    out_meta = src.meta.copy()

out_meta.update({
    "driver": "GTiff",
    "height": out_image.shape[1],
    "width": out_image.shape[2],
    "transform": out_transform
})

with rasterio.open("cn_masked.tif", "w", **out_meta) as dest:
    dest.write(out_image)

，这是我所犯的错误：

complete_var = '_%s_COMPLETE' % (prog_name.replace('-', '_')).upper()

这会引发此错误：

AttributeError：'rasterio._io.RasterReader‘对象没有属性’替换‘

有人知道我做手术的正确方法吗？

Answer 1

rio函数(如rasterio.rio.clip.clip函数)不能作为rasterio.rio.clip.clip函数调用，因为它们已经附加了一个命令行接口(所有这些装饰器)，并且为从那里获取参数进行了优化。

您可以使用查看剪辑函数的代码来查看它是如何实现剪裁和模仿的。它基本上是从几何图形的边界创建一个窗口，并使用该窗口从源栅格读取(窗口读取)。

或者，您可以遵循使用文档中的示例的rasterio.mask.mask(..., crop=True)，不仅要剪辑到边界，还可以屏蔽几何之外的像素，这可能更接近您想要做的事情。