如何使用Python的Geopandas移除合并区域后的残留线条？

我通过以下代码读取并筛选两个Shapefile的目标区域：

import geopandas as gpd

sg = gpd.read_file('shapefile1.shp')
sd = gpd.read_file('shapefile2.shp')

kw = sg[sg['region_name'] == 'region1']
tg = sd[sd['region_name'] == 'region2']

region1和region2的初始形态：

执行合并并绘图：

mr = kw.geometry.reset_index(drop=True).union(tg.geometry.reset_index(drop=True))
mr.plot()

得到合并后的区域：

但查看边界时发现内部残留多余线条：

mr.boundary.plot()

已尝试simplify(tolerance=0.01)和buffer(0)方法，均未解决。需要得到仅含外部轮廓、内部无多余线条的结果。

试试以下几种针对性处理方法：

1. 优先合并单图层内的几何
   如果kw或tg包含多个几何对象，直接使用union可能因拓扑问题残留线条。先通过unary_union合并单个图层内的所有几何，再做跨图层合并：

kw_union = kw.geometry.unary_union
tg_union = tg.geometry.unary_union
mr = kw_union.union(tg_union)

2. 双向Buffer清理拓扑
   单纯的buffer(0)无法处理细微的重叠或缝隙时，用微小的正向Buffer再反向Buffer，强制清理内部冗余：

# 微小膨胀后收缩，消除拓扑异常
mr_clean = mr.buffer(0.001).buffer(-0.001)
# 按需简化边界
mr_clean = mr_clean.simplify(0.01)
mr_clean.boundary.plot()

3. 修复几何有效性
   原始几何或合并后的几何可能存在拓扑无效（比如自相交），用Shapely的make_valid修复后再提取外部轮廓：

from shapely.validation import make_valid

# 修复几何有效性
mr_valid = make_valid(mr)
# 若修复后为MultiPolygon，提取面积最大的外部轮廓（适配单一大区域场景）
if mr_valid.geom_type == 'MultiPolygon':
    mr_clean = max(mr_valid.geoms, key=lambda g: g.area)
else:
    mr_clean = mr_valid
mr_clean.boundary.plot()

4. 检查并修复原始几何
   先确认原始筛选出的几何是否有效，无效几何会导致合并后异常：

# 检查几何有效性
print("kw几何是否全部有效：", kw.geometry.is_valid.all())
print("tg几何是否全部有效：", tg.geometry.is_valid.all())

# 修复无效几何
kw['geometry'] = kw.geometry.apply(make_valid)
tg['geometry'] = tg.geometry.apply(make_valid)

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