使用sf包提取多边形内的点时报无效球面几何错误如何解决

你的gefahren多边形数据存在16个拓扑无效的几何要素，常见问题包括重复顶点、边重叠、自相交、多边形环未正常闭合等。sf v1.0及以上版本默认使用s2引擎处理WGS84这类地理坐标系数据，对几何有效性的校验规则远严于之前的平面计算模式，因此触发报错。

第一步：修复无效几何 #

首先安装依赖包lwgeom，使用st_make_valid()函数自动修复多边形的拓扑错误：

install.packages("lwgeom") # 未安装时执行
library(sf)
library(lwgeom)

# 修复多边形几何
gefahren_valid <- st_make_valid(gefahren)

第二步：执行点筛选 #

你需要提取多边形范围内的点，使用st_intersects()做索引筛选比st_intersection()效率更高，不需要额外做几何裁剪，仅返回匹配的点要素：

# 筛选多边形范围内的点，sparse=FALSE返回布尔矩阵方便直接索引
point_in_poly <- geb_pers_bunk_coord_filt[st_intersects(geb_pers_bunk_coord_filt, gefahren_valid, sparse = FALSE), ]

备选方案：临时关闭s2引擎 #

如果修复几何后仍有报错，且你的研究范围较小（你提供的坐标范围仅几百米，平面计算误差可忽略），可以临时关闭s2引擎使用平面计算模式：

# 关闭s2引擎
sf_use_s2(FALSE)

# 执行筛选
point_in_poly <- geb_pers_bunk_coord_filt[st_intersects(geb_pers_bunk_coord_filt, gefahren_valid, sparse = FALSE), ]

# 操作完成后改回默认设置，避免影响其他空间计算
sf_use_s2(TRUE)

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