我有同样的问题连接空间数据时如何解决球面几何故障
我以前可以使用的代码现在不能用于从 GEOS 到 s2 的转换。
我不确定如何创建此错误的可重现示例。 sf 包附带的示例数据集可以很好地使用此代码,因为我确信它们也已更新为可以与 s2 包一起使用。最好的
ranges <- st_read("Data/Range maps/PRISM_shorebird_ranges_breeding/PRISM_shorebird_ranges_breeding.shp")
ranges2 <- ranges %>%
group_by(binomial) %>%
select(-seasonal, -Shape_Area) %>%
summarize()
Error in s2_geography_from_wkb(x, oriented = oriented, check = check) :
Evaluation error: Found 2 features with invalid spherical geometry.
[1] Loop 0 is not valid: Edge 319 has duplicate vertex with edge 322
[2] Loop 0 is not valid: Edge 53052 has duplicate vertex with edge 53055.
我尝试了以下解决方案连接空间数据时如何解决球面几何故障
您有两个选择:
1. 通过脚本中的 sf::sf_use_s2(FALSE) 关闭 s2 处理;理论上,行为应该恢复到 1.0 版本之前的行为
这确实有效,但似乎不是一个好的长期解决方案。我的其他代码会继续遇到这个问题吗?如何防止这种情况发生并使用更新的 sf 包的完整功能?
2. 修复多边形对象的球面几何形状;这将取决于你错误的实际性质。”
ranges$geometry <- ranges$geometry %>%
s2::s2_rebuild() %>%
sf::st_as_sfc()
Error in s2_geography_from_wkb(x, oriented = oriented, check = check) :
Evaluation error: Found 22 features with invalid spherical geometry.
[1] Loop 2 is not valid: Edge 846 has duplicate vertex with edge 853
[2] Loop 10 is not valid: Edge 0 has duplicate vertex with edge 4
... and 20 more.
如您所见,该解决方案对我不起作用。如何解决这个球面几何错误?