R 围绕点坐标创建非重叠六边形

我必须在坐标点周围创建非重叠、不规则的六边形（voronoi六边形）。我在某地和其他地方找到了一些类似问题的答案，但我没有设法适应这些解决方案

只使用sp、光栅和rgdal是可能的，因为我有一些部分数据（见下图），但我无法重现

点的坐标来自如下所示的shp：

class       : SpatialPointsDataFrame 
features    : 148392 
extent      : -179.9844, 179.6227, 53.00372, 79.90965  (xmin, xmax, ymin, ymax)
crs         : +proj=longlat +datum=WGS84 +no_defs 
variables   : 7
names       :   point, cell, grid,       lat,         lon,   cntr,         cont 
min values  : 2601723,  100, warp, 53.003719, -179.984421, Canada,       Arctic 
max values  : 3213512,  967, warp, 79.909653,  179.622711,    USA, NorthAmerica 

dput示例如下：

structure(list(point = c(3110786, 2861560, 2864634, 2864638, 
                         2864642, 2867700), cell = c("428", "282", "282", "282", "282", 
                                                     "282"), grid = c("warp", "warp", "warp", "warp", "warp", "warp"
                                                     ), lat = c(72.407028, 61.31004, 61.422215, 61.422215, 61.422215, 
                                                                61.534386), lon = c(-120.742577, -140.21254, -140.013031, -140.247162, 
                                                                                    -140.481308, -140.047607), cntr = c("Canada", "Canada", "Canada", 
                                                                                                                        "Canada", "Canada", "Canada"), cont = c("Arctic", "Arctic", "Arctic", 
                                                                                                                                                                "Arctic", "Arctic", "Arctic")), row.names = c(NA, 6L), class = "data.frame")

结果应该是一个shp，如下所示：

这些点来自shp，只需要六边形结构