使用Geopandas实现匹配ArcGIS Dissolve的多段线去重叠效果

之前尝试方法的效果偏差原因 #

你之前使用的gpd.dissolve()、shapelyunary_union操作无法达到预期，核心是两类方法的设计目标和你的需求完全错位：

• 带统一填充列的gpd.dissolve()默认会将同组要素合并为单一多部件几何，哪怕后续调用explode()拆分多部件，拆分过程会在所有重叠、交叉位置打断原线节点，破坏原始线段的连续形态
• 全量线的unary_union操作本质是拓扑融合工具，会自动在所有线的交叉、重叠位置切割节点，把原本连续的原线段拆分为离散拓扑边，完全不符合“不拆分原有线段”的约束

匹配需求的实现逻辑（严格对齐你的约束） #

整个流程全程不对原始线几何做节点修改、切割操作，仅通过空间重叠判定筛选冗余要素：

• 为所有原始线添加唯一标识，从流程起点锁定原几何的完整坐标
• 逐要素检测与已保留要素的长段重叠关系，将重叠长度占比超过阈值（适配pysheds的随机冗余特征，通常取0.9，即重叠长度占较短线总长度的90%以上判定为冗余）的要素筛除
• 端点正常相接、仅短段重叠的合法溪流拓扑关系全部保留
• 最终输出结果全部为独立单线要素，无多部件聚合几何，无切割产生的新线段

可直接复用的代码实现 #

import geopandas as gpd
from shapely.ops import unary_union

# 为原始数据添加临时唯一标识，锁定原始几何关系
stream_0["original_id"] = stream_0.index
keep_id_list = []
# 重叠判定阈值：重叠长度占较短线长度的90%以上判定为冗余
overlap_ratio_threshold = 0.9

for idx, current_row in stream_0.iterrows():
    current_geom = current_row.geometry
    # 第一个要素直接纳入保留列表
    if not keep_id_list:
        keep_id_list.append(idx)
        continue
    # 快速计算当前线与已保留要素集合的重叠长度
    kept_geom_set = unary_union(stream_0.loc[keep_id_list, "geometry"].values)
    overlap_geom = current_geom.intersection(kept_geom_set)
    overlap_len = overlap_geom.length
    # 计算重叠占比，低于阈值即判定为非冗余要素
    shorter_len = min(current_geom.length, overlap_len)
    current_overlap_ratio = overlap_len / shorter_len if shorter_len !=0 else 0
    if current_overlap_ratio < overlap_ratio_threshold:
        keep_id_list.append(idx)

# 基于保留id过滤结果，全程不修改任何原始几何
deduped_stream = stream_0.loc[keep_id_list].drop(columns=["original_id"]).reset_index(drop=True)

实现效果对齐说明：这个逻辑和你在ArcGIS Pro中使用Dissolve工具、取消“创建多部件要素”选项后的处理效果完全一致——仅移除重叠的冗余重复线段，不切割原线、不聚合多部件，100%保留原线段的坐标形态。

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