Geopandas中从大多边形移除内部小多边形及pygeos报错解决

报错原因 #

你虽已安装PyGEOS，但当前Python环境未正确识别到它，或是GeoPandas未关联上PyGEOS。另外，你用的symmetric_difference并非实现需求的最佳操作——要得到带孔洞的大多边形，应该用差集操作而非对称差。

修复方案 #

1. 确保PyGEOS正常可用 #

先验证PyGEOS是否能被导入：

import pygeos
print(pygeos.__version__)

若仍报NameError，按以下方式修复环境：

• 卸载重装PyGEOS：

  pip uninstall pygeos -y
  pip install pygeos
  

• 若用conda环境，优先用conda安装：

  conda install -c conda-forge pygeos
  

• 确认GeoPandas版本与PyGEOS兼容（GeoPandas 0.10+建议搭配PyGEOS 0.10+）

2. 替换为正确的几何操作 #

要从摩洛哥多边形中移除卡萨布兰卡多边形，使用overlay的how='difference'才是正确选择。修改后的代码：

import geopandas as gpd

# 加载摩洛哥数据
world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
geomap_world = world[world['name'] == 'Morocco']

# 构造GeoDataFrame（假设geomap是你已加载的卡萨布兰卡数据）
df_casablanca = gpd.GeoDataFrame({'geometry': geomap.geometry, 'label': ['casablanca']})
df_morocco = gpd.GeoDataFrame({'geometry': geomap_world.geometry, 'label': ['morocco']})

# 执行差集：从摩洛哥多边形中减去卡萨布兰卡
result = df_morocco.overlay(df_casablanca, how='difference')

# 查看结果
print(result)
result.plot()

3. 更直接的替代方案 #

如果overlay仍有问题，可直接对几何对象执行差集操作：

# 获取单个几何对象（确保数据集只有一个目标多边形）
morocco_geom = geomap_world.geometry.iloc[0]
casablanca_geom = geomap.geometry.iloc[0]

# 若卡萨布兰卡是多个多边形，先合并为单个几何
# casablanca_geom = geomap.geometry.unary_union

# 生成带孔洞的多边形
holed_polygon = morocco_geom.difference(casablanca_geom)

# 转为GeoDataFrame
result_gdf = gpd.GeoDataFrame({'geometry': [holed_polygon]}, crs='EPSG:4326')
result_gdf.plot()

注意事项 #

• 确保两个数据集的CRS完全一致（均为EPSG:4326），若不一致需先调用to_crs转换
• 若卡萨布兰卡包含多个离散多边形，需先通过unary_union合并为单个几何，再执行差集操作

内容的提问来源于stack exchange，提问作者bravopapa