#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, #是否顯示點(diǎn)
     fillcol = hcl.colors(60, "Sunset"))  #填充顏色

#改變圖形的形狀
#install.packages("polyclip")
library(polyclip)
# Circle環(huán)狀
s <- seq(0, 2 * pi, length.out = 3000)
circle <- list(x = 0.5 * (1 + cos(s)),
               y = 0.5 * (1 + sin(s)))
plot(tiles, pch = 19,
     col.pts = "white",
     border = "black",
     fillcol = hcl.colors(60, "Sunset""),
     clipp = circle)

library(ggvoronoi)
library(ggplot2)

set.seed(1)
x <- sample(1:600, size = 100)
y <- sample(1:600, size = 100)
dist <- sqrt((x - 200) ^ 2 + (y - 200) ^ 2)

df <- data.frame(x, y, dist = dist)

ggplot(df, aes(x, y)) +
  stat_voronoi(geom = "path",
               color = 6,    
               lwd = 0.7,   
               linetype = 1) 
  geom_point()

#添加熱圖
ggplot(df, aes(x, y, fill = dist)) +
  geom_voronoi() +
  geom_point() +
  scale_fill_gradient(low = "#20B2AA",
                      high = "#9370DB") #顏色

#添加邊界線
ggplot(df, aes(x, y, fill = dist)) +
  geom_voronoi() +
  stat_voronoi(geom = "path") +
  geom_point()+
  scale_fill_gradient(low = "#20B2AA",
                      high = "#9370DB")
 #theme(legend.position = "none") # 去掉右側(cè)圖例標(biāo)簽