deldir包绘制Voronoi图
#install.packages("deldir")
library(deldir)
# data
set.seed(1)
x <- runif(60)
y <- runif(60)
# Calculate Voronoi Tesselation and tiles
tesselation <- deldir(x, y)
tiles <- tile.list(tesselation)
plot(tiles, pch = 19,
border = "black", #边界颜色
showpoints = TRUE, #是否显示点
fillcol = hcl.colors(60, "Sunset")) #填充颜色