如何在ggplot中用RGB颜色绘制地理点并添加连续RGB图例？

用户希望通过RGB颜色值表示每个地理点的三个变量占比，编写了如下代码：

data <- data.frame(
  Variable = c("Point1", "Point2", "Point3"),
  R = c(0.4, 0.6, 0.3),
  G = c(0.3, 0.2, 0.4),
  B = c(0.3, 0.2, 0.3),
  x = c(1, 2, 4),
  y = c(5, 6, 7)
)
color_plot <- rgb(data$R, data$G, data$B)

ggplot() +
  geom_point(data = data, aes(x = x, y = y, color = color_plot))  

但生成的图例为离散型，无实际意义，需要添加带有R、G、B名称的连续颜色刻度。

方法1：使用ggnewscale添加独立连续图例 #

该方法可为R、G、B三个分量分别创建连续颜色图例，点的实际颜色由三个分量共同决定：

先安装并加载依赖包：

install.packages(c("ggplot2", "ggnewscale"))
library(ggplot2)
library(ggnewscale)

修改后的代码：

# 绘制基础散点图，直接指定点的RGB颜色
p <- ggplot(data, aes(x = x, y = y)) +
  geom_point(aes(color = I(rgb(R, G, B))), size = 5)  # I()用于跳过映射，直接设置颜色

# 添加R分量的连续图例
p <- p +
  new_scale_color() +
  geom_point(aes(color = R), size = 0) +  # 空点仅用于生成图例
  scale_color_gradientn(name = "R分量", colors = c("black", "red"), limits = c(0,1))

# 添加G分量的连续图例
p <- p +
  new_scale_color() +
  geom_point(aes(color = G), size = 0) +
  scale_color_gradientn(name = "G分量", colors = c("black", "green"), limits = c(0,1))

# 添加B分量的连续图例
p <- p +
  new_scale_color() +
  geom_point(aes(color = B), size = 0) +
  scale_color_gradientn(name = "B分量", colors = c("black", "blue"), limits = c(0,1))

print(p)

最终会生成三个独立的连续图例，每个图例从黑色（0值）对应到原色（1值），清晰展示各分量的数值变化。

方法2：三元图例（适用于三变量总和为1的占比场景） #

如果R、G、B为总和为1的占比数据，可使用ggtern包创建三元坐标轴图例，同时叠加地理坐标：

install.packages("ggtern")
library(ggtern)

ggtern(data, aes(R, G, B)) +
  geom_point(aes(x = x, y = y), inherit.aes = FALSE, size = 5, color = rgb(data$R, data$G, data$B)) +
  labs(x = "R分量", y = "G分量", z = "B分量") +
  theme_bw()

该图的三元坐标轴展示三个分量的比例关系，散点在右侧x-y轴上显示地理位置，颜色对应RGB组合，适合直观展示三变量的比例分布。

方法3：自定义RGB组合图例（进阶） #

若需要用单个图例展示RGB组合变化，可生成参考图例并与主图组合：

# 生成RGB颜色参考数据
rgb_ref <- expand.grid(R = seq(0,1,0.1), G = seq(0,1,0.1), B = seq(0,1,0.1))
rgb_ref$color <- rgb(rgb_ref$R, rgb_ref$G, rgb_ref$B)

# 绘制主地理散点图，预留右侧空间放图例
main_plot <- ggplot(data, aes(x, y)) +
  geom_point(color = rgb(data$R, data$G, data$B), size = 5) +
  theme(plot.margin = margin(1, 4, 1, 1, "cm"))

# 绘制RGB参考图例
legend_plot <- ggplot(rgb_ref, aes(R, G)) +
  geom_tile(aes(fill = color)) +
  scale_fill_identity() +
  labs(x = "R", y = "G") +
  annotate("text", x = 1.1, y = seq(0,1,0.2), label = seq(0,1,0.2)) +  # 标注B值
  theme_minimal() +
  theme(axis.text.y = element_blank(), axis.ticks.y = element_blank())

# 组合主图和图例
library(gridExtra)
grid.arrange(main_plot, legend_plot, ncol = 2, widths = c(3,1))

该方法通过额外的图例图展示R-G平面的颜色变化，同时标注B值，适合直观展示RGB组合的关系。

内容的提问来源于stack exchange，提问作者El toto