如何用Python修复自相交无效多边形并合并Shapefile多边形？

问题背景 #

尝试合并Shapefile中的多边形时，因部分多边形存在自相交拓扑错误，执行代码触发TopologyException（边位置冲突）和ValueError（无法从空值创建几何），合并操作直接失败。

核心问题分析 #

原代码先执行unary_union合并所有几何，但此时数据中已存在无效几何，直接合并会导致GeoPandas返回空几何对象，后续的修复步骤根本无法启动。正确逻辑是先修复每个几何的有效性，再执行合并操作。

修复代码 #

import geopandas as gpd
from shapely.validation import make_valid, explain_validity

# 读取目标Shapefile
pz32 = gpd.read_file("B://_Shp_robocze//test//test.shp")

# 1. 逐个检查并修复几何有效性
pz32['geometry'] = pz32['geometry'].apply(
    lambda geom: make_valid(geom) if not geom.is_valid else geom
)

# 2. 过滤修复后仍为空的无效几何（避免后续合并报错）
pz32 = pz32[pz32['geometry'].notna()]

# 3. 执行多边形合并
merged_geometry = pz32.unary_union

# 将合并结果转为GeoSeries输出
shp = gpd.GeoSeries([merged_geometry])
print(shp)

# 可选调试：查看单个几何的无效原因
# for idx, geom in enumerate(pz32['geometry']):
#     if not geom.is_valid:
#         print(f"几何{idx}无效原因：{explain_validity(geom)}")

关键说明 #

• make_valid()：Shapely的核心拓扑修复函数，可自动处理自相交、环方向错误等常见几何问题，返回符合OGC标准的有效几何对象。
• 先修复再合并：必须确保所有输入几何有效后，再执行unary_union，避免无效几何导致合并流程中断。
• 过滤空几何：极端无效的几何修复后可能仍为空值，提前过滤能避免后续操作触发空值错误。

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