R 在 map 中用六边形替换点

发布于11月23日

在 map 上，我试图用六边形来替换点，但做不到.我试过用geom_hex，但没有得到同样的结果，有什么建议吗？

ggplot() + 
  geom_point(data=df, aes(x=x, y=y, color = value), linewidth = 1, na.rm = FALSE) +
  scale_colour_gradientn(colours = c('#80146E','#804AA4','#6771B8','#379BC2',
                                     '#3ABCBF','#79D4B8','white','#FCEE88','#F8C34B', 
                                     '#F7A72B', '#F4792D'),
                         values = scales::rescale(-10, -8, -6, -4, -2, -1, 0, 
                                                    0.3, 0.5, 0.7, 1))) + 
  theme_void() + 
  coord_map()

编辑

df <- head(df)
dput(df)
structure(list(x = c(-1.29180627034999, -1.19180627074999, -1.09180627114999, 
-0.991806271549991, -0.891806271949991, -0.791806272349991), 
    y = c(44.24152724085, 44.24152724085, 44.24152724085, 44.24152724085, 
    44.24152724085, 44.24152724085), value = c(-0.709841088006482, 
    -0.889560720544426, -0.998504169250709, -0.905970350930832, 
    -0.888156362836211, -0.960843878765758)), row.names = c(NA, 
6L), class = "data.frame")

推荐答案

我们可以通过获取带有海拔数据的法国栅格，对其进行下采样并转换为x，y，z格式来近似您的数据.

df <- raster::getData('alt', country = 'FRA', mask = TRUE) |>
  raster::aggregate(fact = 15) |>
  raster::rasterToPoints() |>
  as.data.frame() |>
  setNames(c('x', 'y', 'value'))

使用这些点的代码绘制的图类似于您自己的图像，但由于值的不同，其 colored颜色有所不同

library(ggplot2)

ggplot(df, aes(x, y, color = value)) + 
  geom_point(size = 0.8) +
  theme_void() + 
  scale_colour_gradientn(colours = c('#80146E','#804AA4','#6771B8','#379BC2',
                                     '#3ABCBF','#79D4B8','white','#FCEE88',
                                     '#F8C34B', '#F7A72B', '#F4792D'),
                         values = c(0, 0.01, 0.02, 0.1, 0.3, 0.6, 0.8, 0.85,
                                    0.9, 0.95, 1)) +
  coord_map()

如果你想换成六边形，你需要stat_hex_summary,fun = mean和aes(z = value).还记得把你的scale_color_gradientn号改为scale_fill_gradientn号

ggplot(df, aes(x, y, z = value)) +
  stat_summary_hex(fun = mean, color = 'white', bins = 60) +
  theme_void() + 
  scale_fill_gradientn(colours = c('#80146E','#804AA4','#6771B8','#379BC2',
                                     '#3ABCBF','#79D4B8','white','#FCEE88',
                                     '#F8C34B', '#F7A72B', '#F4792D'),
                         values = c(0, 0.01, 0.02, 0.1, 0.3, 0.6, 0.8, 0.85,
                                    0.9, 0.95, 1)) +
  coord_map()