Geopandas/Rasterio掩膜报Input shapes do not overlap raster错误如何解决

DevOps 小满

2026-9-25

问题排查与解决

首先明确：该报错不是因为未给全球人口栅格明确指定WGS 84坐标系导致。Rasterio读取tif格式栅格时会自动解析文件内置的坐标系统参数，你已经确认栅格本身CRS为EPSG:4326，无需额外指定。

以下是可能的成因及对应解决方法：

成因1：未筛选到正确的福特波塔尔矢量要素
你是通过遍历属性表所有字段匹配FORT PORTAL字符串完成筛选，可能出现匹配为空、匹配到其他同名要素的情况。
解决方法：执行print(len(gdf))检查筛选结果的要素数量，若输出为0则说明未匹配到有效要素，需要调整字段匹配逻辑，指定行政名称对应的字段做精确匹配，避免遍历所有字段的误判。
成因2：矢量与栅格实际范围无交集
即使二者CRS一致，也可能因为栅格范围本身不覆盖目标区域导致报错。
解决方法：分别打印二者的边界对比：
```
# 打印福特波塔尔矢量四至：(最小经度, 最小纬度, 最大经度, 最大纬度)
print(gdf.total_bounds)
# 打印人口栅格四至
print(gpw.bounds)
```
福特波塔尔的实际经纬度范围大致为东经30.1°_{30.4°、北纬0.5°}0.8°，如果栅格范围不在该区间、或者矢量坐标经纬度顺序颠倒（如纬度在前经度在后，数值变成0.5°_{0.8°经度、30.1°}30.4°纬度），都会出现无交集报错。若为坐标顺序颠倒，可通过shapely.ops.transform交换x/y坐标顺序后重新赋值给gdf的geometry列。
成因3：栅格读取对象或参数异常
你在代码中先读取了一次栅格的窗口数据gpw_region，如果后续对gpw原始数据集对象做过修改、或者你打开的gpw本身是已经裁剪过的乌干达局部栅格而非全球栅格，也可能出现范围不匹配。另外rasterio的mask方法传入GeoSeries类的geometry序列偶发识别异常。
解决方法：重新执行gpw = rio.open('你的全球栅格路径')读取原始栅格，调用mask方法时将geometry转为列表：shapes=list(gdf.geometry)；也可将矢量和栅格导入QGIS可视化，直观确认二者是否存在空间重叠。