使用GDAL.Warp基于Shapefile多边形裁剪栅格时无法正确裁剪的问题

从你提供的截图能看到，用Shapefile多边形裁剪栅格后，裁剪线外的区域仍残留着值为1的像素——哪怕你尝试设置不同的dstNodata值（0或者np.nan），结果都没有改善。我之前处理遥感影像裁剪时也碰到过类似的情况，给你几个实用的排查和解决方向：

1. 先确认坐标系是否完全匹配 #

这是最容易踩的坑：如果你的Shapefile和待裁剪栅格的坐标系（投影）不一致，GDAL自动重投影的过程可能会出现错位，导致裁剪线和栅格像素边缘不贴合，最终留下多余的像素。你可以用gdalinfo命令查看栅格的坐标系，用ogrinfo查看Shapefile的坐标系，确保两者的EPSG代码完全一致。

2. 保证Nodata值和栅格数据类型兼容 #

你提到设置0和np.nan都无效，这里要注意数据类型的匹配：

• 如果你的栅格是整数类型（比如uint8、int16这类），np.nan是浮点数，根本无法被写入到整数栅格中，这时候必须选择一个整数类型的nodata值（比如0、255这类不会和有效数据冲突的值）；
• 如果是浮点数栅格，np.nan才会生效，此时设置整数nodata可能会被自动转换为浮点数，但效果可能不如预期。

3. 修改代码，添加关键参数强化裁剪效果 #

试试给gdal.Warp加上cutlineBlend=0参数，这个参数可以让GDAL严格按照裁剪线边界进行掩膜，避免边缘的羽化效果；另外建议指定输出数据类型和输入一致，避免类型转换导致的nodata失效。修改后的代码示例如下：

from osgeo import gdal

# 根据栅格数据类型选择适配的nodata值，这里以整数栅格用0为例
nodata_val = 0

gdal.Warp(
    "rasterOut.jp2",
    "rasterIn.jp2",
    cutlineDSName="mask.shp",
    cropToCutline=True,
    dstNodata=nodata_val,
    cutlineBlend=0,  # 严格按照裁剪线裁剪，无羽化过渡
    outputType=gdal.GDT_Byte,  # 这里要和输入栅格的数据类型匹配，比如输入是uint8就用GDT_Byte
)

4. 检查Shapefile的拓扑有效性 #

如果你的Shapefile多边形存在自相交、重叠、几何断裂等拓扑错误，GDAL的裁剪逻辑也会出现异常。你可以用QGIS的“检查几何有效性”工具，或者ogr2ogr的修复命令来处理Shapefile的拓扑问题，确保多边形是闭合且有效的。

如果以上方法都没能解决问题，你可以试试另一种稳妥的替代方案：先用gdal.Rasterize把Shapefile转成和待裁剪栅格分辨率、范围一致的掩膜栅格，再用gdal.Translate结合掩膜来完成裁剪。

备注：内容来源于stack exchange，提问作者ilFonta