如何对比两个LineString类型GeoDataFrame并找出几何差异？

前提准备：统一坐标系 #

首先必须确保两个GeoDataFrame的**坐标系（CRS）**一致，否则空间对比会完全失效：

import geopandas as gpd

# 以post_df的坐标系为基准，转换new_df的CRS
new_df = new_df.to_crs(post_df.crs)

方法1：用空间连接（sjoin）筛选单方向差异要素 #

适合单独提取某一个GeoDataFrame中独有的要素（比如你蓝色圈出的new_df差异部分）：

# 空间连接，找出new_df中与post_df有空间交集的要素
matched_records = gpd.sjoin(new_df, post_df, how="inner", predicate="intersects")
# 从new_df中过滤掉已匹配的记录，剩下的就是差异部分
new_only_diff = new_df[~new_df.index.isin(matched_records.index)]

方法2：用几何叠加（overlay）生成双向对称差异 #

如果需要同时获取两个GeoDataFrame互相独有的几何部分（对应你预期的效果），直接用overlay的对称差模式：

# 生成两个GeoDataFrame的几何对称差，即双方独有的所有几何区域
symmetric_diff = gpd.overlay(new_df, post_df, how="symmetric_difference")

处理几何精度问题 #

如果因为坐标精度误差导致匹配失败，可以给几何添加微小缓冲来消除精度影响（缓冲值根据数据坐标单位调整，比如经纬度用0.00001，米单位用0.01）：

new_df['geometry'] = new_df.geometry.buffer(0.00001)
post_df['geometry'] = post_df.geometry.buffer(0.00001)

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