相关疑难解决方法(0)

我正在寻找一个简单的(如果存在的)算法来找到球体表面上一组点的Voronoi图.源代码会很棒.我是德尔福人(是的,我知道......),但我也吃C代码.

我试图在中央子午线不同于0的Robinson投影中投影世界地图。根据此StackOverFlow线程，这应该很容易（尽管示例使用sp）。

这是我的可复制代码：

library(sf)
library(ggplot2)
library(rnaturalearth)

world <- ne_countries(scale = 'small', returnclass = 'sf')

# Notice +lon_0=180 instead of 0
world_robinson <- st_transform(world, crs = '+proj=robin +lon_0=180 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs')

ggplot() +
geom_sf(data = world_robinson)

这就是结果。多边形从投影的一侧向另一侧闭合。

尝试具有sp相同的效果。我还尝试了一个shapefile，该文件仅包含来自http://www.naturalearthdata.com/的海岸线多边形（无政治边界），效果相似。

我试图在Mac OS X和Ubuntu 18.04上的两个独立的R安装上运行我的代码段。