使用polygon裁剪栅格，rasterio.mask报错：ValueError: Input shapes do not overlap raster.

背景：栅格文件tif 和 一个保证被栅格完全覆盖的polygon(shp格式) ，使用rasterio.mask报错：ValueError: Input shapes do not overlap raster.（栅格与形状不重叠）

我用的fiona打开的polygon,代码：

def clipRasterByShapefile(src, shpdatafile, dst, nodata=0):

    with fiona.open(shpdatafile, "r", crs=from_epsg(32718)) as shapefile:
        features = [feature["geometry"] for feature in shapefile]
    src = rasterio.open(src)
    

    print(shapefile.crs, src.crs)
    
    out_image, out_transform = rasterio.mask.mask(src, features,
                                             all_touched=True,
                                             crop=True,
                                             nodata=nodata)

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

    output_file = rasterio.open(dst, "w", compress="LZW", **out_meta)
    output_file.write(out_image)
    output_file.close()

这个错偶尔会有，原因是polygon的crs为None，可以打印出来检查一下。

解决办法：把polygon的crs转换为与tif相同的crs

from fiona.crs import from_epsg
from rasterio.warp import transform_geom

features = [transform_geom(from_epsg(4326), src.crs, geom) for geom in features]

哎，以后全面chatgpt时代这种小问题记录不知道还有没有意义