Geopandas合并多边形遇TopologyException：无法为外壳分配自由孔洞

这个错误源于你的多边形数据存在拓扑无效问题（比如自相交、孔洞未正确关联到外壳、几何方向错误等），导致unary_union无法正常处理。以下是具体解决步骤：

1. 修复所有多边形的拓扑有效性 #

优先使用shapely的make_valid工具，它能处理绝大多数拓扑错误，覆盖范围比buffer(0)更广：

import geopandas as gpd
from shapely.validation import make_valid

# 将你的FeatureCollection转换为GeoDataFrame
gdf = gpd.GeoDataFrame.from_features(your_feature_collection['features'])

# 批量修复无效几何
gdf['geometry'] = gdf.geometry.apply(lambda geom: make_valid(geom) if not geom.is_valid else geom)

如果make_valid无法解决个别自相交问题，再尝试用buffer(0)补漏：

gdf['geometry'] = gdf.geometry.apply(lambda geom: geom.buffer(0) if not geom.is_valid else geom)

2. 按road_name分组合并 #

修复后直接用dissolve方法完成分组合并：

merged_gdf = gdf.dissolve(by='road_name', aggfunc='first').reset_index()

若dissolve仍报错，手动分组并再次校验每组几何：

merged_data = []
for road_name, group in gdf.groupby('road_name'):
    # 二次确保组内几何有效
    group['geometry'] = group.geometry.apply(make_valid)
    # 执行合并操作
    union_geom = group.geometry.unary_union
    merged_data.append({'road_name': road_name, 'geometry': union_geom})

merged_gdf = gpd.GeoDataFrame(merged_data, crs=gdf.crs)

3. 针对性排查（若上述方法失效） #

报错提示的坐标138.6334, -34.8695附近可能存在严重拓扑问题，提取该区域数据单独检查：

# 提取报错坐标周边的多边形
problem_gdf = gdf.cx[138.632:138.635, -34.8705:-34.8685]
# 可视化查看问题（需matplotlib支持）
problem_gdf.plot()

重点排查该区域是否存在：

• 孔洞与外壳交叉重叠
• 多个多边形边界完全重合导致冲突
• 极端坐标精度误差

针对这类问题，可尝试先简化几何以减少坐标点干扰：

gdf['geometry'] = gdf.geometry.simplify(tolerance=0.0001)  # tolerance根据数据精度调整

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