如何识别缓冲区图层重叠（相交）区域并解决st_intersection运行报错问题

报错原因 #

1. 逻辑不匹配：你需要计算250米、500米两个缓冲区图层的相交区域，但现有代码仅向st_intersection()传入了单个b500图层，此时函数会默认计算该图层内部所有要素的n元交集，不仅不符合业务需求，也更容易触发拓扑计算错误。
2. 底层触发原因：CPL_nary_intersection(x)报错通常来自几何无效、精度不匹配、拓扑异常三类问题，单独设置精度未生效通常是因为未将精度设置后的结果覆盖原图层，或者精度阈值设置不合理。

解决方案 #

第一步：修正交集计算逻辑（优先执行） #

你需要将两个缓冲区图层都传入st_intersection()，计算二元交集即可拿到两类缓冲区的重叠区域，示例代码如下：

# 假设250米缓冲区图层名为b250，替换为你实际的图层名即可
b_overlap <- st_intersection(b500, b250)

# 管道写法
b_overlap <- b500 %>% 
  st_intersection(b250)

第二步：处理几何有效性问题（如第一步仍报错） #

st_make_valid()比调整精度对拓扑错误的修复成功率更高，执行以下代码修复两个图层的无效几何后再计算交集：

# 检查是否存在无效几何，返回值大于0即存在
sum(!st_is_valid(b500))
sum(!st_is_valid(b250))

# 修复无效几何
b500_valid <- st_make_valid(b500)
b250_valid <- st_make_valid(b250)

# 用修复后的图层计算交集
b_overlap <- st_intersection(b500_valid, b250_valid)

第三步：正确配置精度参数（如仍有报错） #

st_set_precision()需要将修改后的结果赋值覆盖原对象才会生效，如果你使用的是米制投影，可参考以下配置：

# 精度设置为1米，可根据你的数据精度调整阈值
b500_valid <- b500_valid %>% st_set_precision(1e3)
b250_valid <- b250_valid %>% st_set_precision(1e3)

b_overlap <- st_intersection(b500_valid, b250_valid)

特殊场景补充（仅针对需要计算单图层内部要素重叠的需求） #

如果你确实需要计算500米缓冲区内部多个要素的重叠区域，可先筛选出存在重叠的要素再做交集，避免全量n元计算触发报错：

# 查找所有存在重叠的要素索引
overlap_idx <- st_overlaps(b500_valid) %>% unlist() %>% unique()
# 仅对存在重叠的要素做交集计算
b500_internal_overlap <- b500_valid[overlap_idx, ] %>% st_intersection()

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