R sf中与Esri Identity工具等价的操作实现方法咨询

Identity工具的核心能力正好匹配需求：保留输入图层X的全部边界范围，与叠加图层Y相交的区域继承Y的属性与细分边界，未相交区域保留X原有属性，目前有多种可落地的实现方案：

1 桌面端工具实现（QGIS） #

QGIS原生内置同名Identity工具，无需额外开发即可直接使用，操作路径：

• 打开QGIS后在顶部菜单栏选择处理 > 工具箱
• 在工具箱搜索框输入Identity，找到矢量叠加分类下的Identity工具
• 输入图层选择待保留全部边界的X数据，叠加图层选择Y数据，设置输出路径后运行即可得到完全匹配需求的结果

2 代码实现（Python + GeoPandas） #

如果需要批量处理或者嵌入自动化流程，可以用GeoPandas的叠加分析接口组合实现，示例代码如下：

import geopandas as gpd
import pandas as pd

# 读取空间数据，支持shp、geojson、GPKG等常见格式
x_layer = gpd.read_file("your_x_data.geojson")
y_layer = gpd.read_file("your_y_data.geojson")

# 校验坐标系一致，不一致则转换
if x_layer.crs != y_layer.crs:
    y_layer = y_layer.to_crs(x_layer.crs)

# 计算相交部分，继承Y的属性
intersection_part = gpd.overlay(x_layer, y_layer, how="intersection")
# 计算X中不与Y相交的部分，保留X原有属性
difference_part = gpd.overlay(x_layer, y_layer, how="difference")
# 合并两部分得到最终结果
identity_result = gpd.GeoDataFrame(
    pd.concat([intersection_part, difference_part], ignore_index=True),
    crs=x_layer.crs
)

# 导出结果
identity_result.to_file("identity_result.shp", encoding="utf-8")

如果X和Y存在重名字段，可在读取后通过y_layer = y_layer.rename(columns={"old_name": "y_old_name"})的方式重命名Y的字段，避免属性覆盖。

3 其他GIS平台实现 #

• 国产商业GIS平台（SuperMap、MapGIS等）均内置同名Identity叠加分析工具，操作逻辑和ArcGIS完全一致，直接调用即可
• 空间数据库（PostGIS）可通过ST_Intersection+ST_Difference+UNION ALL的SQL逻辑组合实现等价效果

注意事项 #

• 叠加运算前需确认两个图层坐标系完全一致，否则会出现计算错误
• 图层存在拓扑错误（如自相交、重叠面）时，建议先做拓扑修复再做叠加，减少无效结果
• 大体积数据运算时建议添加空间索引，提升计算效率

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