如何为 R 中的点坐标构建 Voronoi 多边形？

编程问答 2022-05-15

问题描述

我在阿尔卑斯山的观测站有多个点（2000+）。我想用它们来代表最近的地理区域，即不靠近另一个观测站。我做了一些研究，并认为使用 Varanoi 多边形可能是最好的方法。

尝试在 R 中构建这些后，多边形图与我在 R 中的图形不太匹配。

我附上了我正在试验的样本数据点，以及显示这些点不同图形的两个图像。

我需要做些什么不同的事情来确保点对齐？

积分：

Longitude:
15.976667 12.846389 14.457222 13.795556  9.849167 16.055278 13.950833 15.666111  9.654722 15.596389 13.226667 15.106667 13.760000 12.226111  9.612222 17.025278  9.877500 15.368056 13.423056 12.571111 16.842222 13.711667 14.003056 12.308056 13.536389

Latitude:
48.40167 48.14889 47.56778 46.72750 47.45833 48.04472 47.82389 47.49472 47.35917 48.64917 48.25000 48.87139 47.87444 47.42806 47.20833 47.77556 47.40389 47.87583 47.53750 46.77694 47.74250 46.55000 48.37611 47.38333 47.91833

图片： Leaflet中25个样本点的地图：

enter image description here

Voronoi 情节：

enter image description here

显然这两个不是同一个图像，所以我一定是做错了什么。这是我用来生成 Voronoi 图和传单地图的代码。

meta25%>% 
  st_as_sf(coords = c("Longitude","Latitude"),crs = sp::CRS("+proj=longlat +datum=WGS84")) %>% 
  mapview()

m1 = matrix(meta25$Longitude,meta25$Latitude,ncol=2,nrow=25) %>% st_multipoint()
voronoi_grid <- st_voronoi(m1)
plot(voronoi_grid,col = NA)
plot(m1,add = TRUE,col = "blue",pch = 16)

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

gisleafletrrvoronoi