连接空间数据时如何解决球面几何故障

Question

我有一个 shapefile（有几个多边形）和一个带坐标的数据框。我想将数据框中的每个坐标分配给 shapefile 中的多边形。因此，要在数据框中添加一个具有多边形名称或 ID 的列这是 link to the data

library(sf)
library(readr)
shape <- read_sf("data/mesopelagic_regions/GlasgowMesopelagicProvinces_v1_2017.shp")
data<- read_csv("data/data.csv")

但是当我尝试加入他们时，我总是得到错误

pts = st_as_sf(data, coords = c("dec_lon", "dec_lat"), crs= 4326)

st_join(pts, shape)

我尝试了 over() 函数和其他技巧，例如 st_make_valid()，但我总是遇到此错误： Error in s2_geography_from_wkb(x, oriented = oriented, check = check) : Evaluation error: Found 30 features with invalid spherical geometry.

这是最近的一个问题（在我的代码工作之前），但现在我无法使用 sf 包来完成这个任务，我总是以这个错误结束。我更新了库以查看它是否有帮助，但我无法使其工作。

非常感谢你在这件事上的帮助

Answer 1

您有两个选择：

通过脚本中的 sf::sf_use_s2(FALSE) 关闭 s2 处理；理论上，行为应该恢复到 1.0
修复多边形对象的球形几何体；这将取决于您的错误的实际性质。

我无法访问您的文件并确定，但这段代码在过去对我有帮助：

yer_object$geometry <- yer_object$geometry %>%
  s2::s2_rebuild() %>%
  sf::st_as_sfc()

连接空间数据时如何解决球面几何故障

How to resolve spherical geometry failures when joining spatial data

r

sf

s2