sf::st_intersection与OSM数据相交操作失败问题排查

问题根源 #

你的bbox多边形顶点顺序错误，导致生成的多边形内部区域实际被识别为外部空间（sf遵循OGC标准，多边形顶点需按顺时针/逆时针的正确环方向定义，反向环会让空间逻辑完全颠倒）。虽然可视化看起来两个图形重叠，但空间交集检测会判定它们没有真正的内部重叠，因此返回空结果。

验证方法 #

先检查几何有效性，确认问题：

# 检查bbox多边形是否有效
sf::st_is_valid(bbox$geometry)
# 查看具体错误原因
sf::st_is_valid(bbox$geometry, reason = TRUE)

大概率会返回类似“Invalid ring orientation”（环方向无效）的错误提示。

修复步骤 #

1. 修正bbox顶点顺序 #

将顶点调整为正确的闭合环顺序（比如顺时针：左下→右下→右上→左上）：

bbox = data.frame(
  x = c(-8664207, -8663609, -8663609, -8664207), 
  y = c(4807913, 4807913, 4808799, 4808799), 
  name = 'bbox'
) |>
  sf::st_as_sf(coords = c('x', 'y'), crs = 3857) |>
  dplyr::group_by(name) |>
  dplyr::summarize(do_union = FALSE) |>
  sf::st_cast('POLYGON')

2. 强制修复几何有效性（可选） #

如果调整顺序后仍存在问题，用st_make_valid修复所有无效几何：

# 修复两个几何的有效性
co_multipoly_valid = sf::st_make_valid(co_multipoly$geometry)
bbox_valid = sf::st_make_valid(bbox$geometry)

# 执行交集操作
sf::st_intersection(co_multipoly_valid, bbox_valid)

测试结果 #

调整顶点顺序后，重新运行st_intersection就能得到正确的交集结果，与可视化显示的重叠区域一致。

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