根据空间点数据创建边界多边形以在传单中进行绘制

Question

你好（在这里映射新手！）

我四处游历，无法解决似乎棘手的问题。

我有基本的XY（坐标）数据：

我要做的是基于不重叠且具有一定大小限制的坐标数据创建相邻的多边形（这样它们就不会永远延伸到海洋中）。

请原谅我糟糕的MS Paint技术，但所需的结果将是这样的：

我有一个标记为陆/海界面的多边形，所以多边形也不能重叠。

我正在使用Leaflet使这些地图具有交互性，它不是用于任何统计分析，而是用于提供概述。

最终目的是使每个多边形都以变量（例如温度）上色，并覆盖生态数据。

一些示例数据：

    > data[1:10,]
   Station  Lat_dec  Long_dec Surface_T
1      247 50.33445 -2.240283     15.19
2      245 50.58483 -2.535217     14.11
3      239 50.16883 -2.509250     15.41
4      225 50.32848 -2.765967     15.34
5      229 50.63900 -2.964800     14.09
6      227 50.33757 -3.303217     15.12
7      217 50.16657 -3.563817     15.13
8      207 49.66683 -3.556550     15.04
9      213 50.16512 -3.824667     14.97
10     219 49.83707 -3.815483     14.78

产生图1的代码是基本的传单脚本：

leaflet() %>% 
  addProviderTiles('Esri.OceanBasemap'
  ) %>% 
  addCircleMarkers(data = data,
                   lng = ~Long_dec,
                   lat = ~Lat_dec,
                   radius = 2
  ) %>%
  addPolygons(data = Land,
              weight = 1,
              color = 'black')

整天让我感到沮丧的是，大多数示例都使用下载的多边形（例如，看起来像是经典的美国，而不是制造它们）

任何帮助深表感谢！ （或者我问得太多！）

吉姆

Answer 1

这是开始的事情：

library(sf)
library(dplyr)

#create sf object with points
stations <- st_as_sf( df, coords = c( "Long_dec", "Lat_dec" ) ) 

#create voronoi/thiessen polygons
v <- stations %>% 
  st_union() %>%
  st_voronoi() %>%
  st_collection_extract()

library(leaflet)
leaflet() %>% 
  addTiles() %>% 
  addCircleMarkers( data = stations ) %>%
  addPolygons( data = v )