ggmap中带有自定义颜色代码的Voronoi镶嵌？

Posted 2023-02-16

【中文标题】ggmap中带有自定义颜色代码的Voronoi镶嵌？

【英文标题】：Voronoi tesselation in ggmap with custom color codes?

【发布时间】：2021-05-07 04:59:42

【问题描述】：

我一直在尝试在 ggmap 中绘制 voronoi 镶嵌，其中每个块的颜色将被赋予十六进制代码，例如 #FFCC00。到目前为止我想出的代码如下：

library(ggmap)
library(ggforce)

b <- get_map(c(2.09174, 50.52550, 7.36819, 53.68320),
             maptype = "toner",
             source = "stamen",
             zoom = 8)

lon <- c(3.76779, 5.31313, 3.48031, 3.90727, 4.15682)
lat <- c(51.2219, 52.0808, 50.7684, 51.2684, 50.9502)
hex_col <- c("#5A586E", "#47967F", "#4EB22E", "#9E82C5", "#ADCFAD")
to_plot <- data.frame(lon, lat, hex_col)


ggmap(b, base_layer = ggplot(data = to_plot,
                             aes(x = lon,
                                 y = lat))) +
  geom_voronoi_tile(aes(fill = hex_col)) +
  scale_fill_identity() +
  geom_voronoi_segment()

但是，当我添加fill = hex_col 参数时，会出现错误警告：

Warning message:
Computation failed in `stat_voronoi_tile()`:
There is at most one point, data or dummy, inside
the given rectangular window. Thus there are
insufficiently many points to triangulate/tessellate. 

我不确定如何修复，因为在添加参数之前，地图显示没有错误。因此我的问题是：如何将自定义颜色编码的 voronoi 镶嵌添加到 ggmap 上？

提前致谢！

【参考方案1】：

问题可能是 geom_voronoi_tile 期望 voronoi 多边形是闭合的，而您的数据集缺少外部边界。 一个快速的替代方法是回退到 ggvoronoi::geom_voronoi()。

library(ggmap)
library(ggforce)
library(ggvoronoi)
library(tidyverse)

b <- get_map(c(2.09174, 50.52550, 7.36819, 53.68320),
             maptype = "toner",
             source = "stamen",
             zoom = 8)

lon <- c(3.76779, 5.31313, 3.48031, 3.90727, 4.15682)
lat <- c(51.2219, 52.0808, 50.7684, 51.2684, 50.9502)
hex_col <- c("#5A586E", "#47967F", "#4EB22E", "#9E82C5", "#ADCFAD")
to_plot <- data.frame(lon, lat, hex_col)

ggmap(b, base_layer = ggplot(data = to_plot,
                             aes(x = lon,
                                 y = lat))) +
  geom_voronoi(mapping = aes(fill = hex_col), alpha = 0.5) +
  scale_fill_identity() +
  geom_voronoi_segment()