如何用RColorBrewer实现以0为中心的矩阵可视化配色？

当然可以实现这个需求！你想要的是以0为中心的对称颜色映射——0对应白色，正负值分别对应红/蓝（或蓝/红），我们只需要调整颜色方案和色阶范围就能轻松搞定。

问题分析 #

你原来的代码用了brewer.pal(8,"RdBu")，但这个配色的色阶没有强制以0为中心，而且8种颜色的过渡不够细腻，导致0不一定刚好对应白色。我们需要做两个关键调整：

1. 创建蓝-白-红（或红-白-蓝）的对称渐变配色
2. 设置对称的数值范围（zlim），确保0正好落在色阶的中点，对应白色

修改后的完整代码 #

library(RColorBrewer)

# 生成数据（加set.seed让结果可重复）
set.seed(123)
vec = rbinom(10000, 1, 0.1) 
n = sum(vec) 
vec = ifelse(vec == 1, rnorm(n), 0) 
mat = matrix(vec, 100, 100)

# 自定义对称渐变配色：从深蓝到白再到深红，100种颜色过渡（更细腻）
my_palette <- colorRampPalette(c("darkblue", "white", "darkred"))(100)

# 计算数据的最大绝对值，设置对称的数值范围
max_abs_value <- max(abs(vec))

# 绘制图像，核心是zlim参数确保0在色阶中心
image(t(mat)[, nrow(mat):1], 
      col = my_palette, 
      xaxt = "n", 
      yaxt = "n", 
      frame.plot = TRUE, 
      useRaster = TRUE,
      zlim = c(-max_abs_value, max_abs_value)) # 关键：对称范围让0完美对应白色

关键细节说明 #

• 配色自定义：colorRampPalette能生成平滑的渐变颜色，如果你想要红-白-蓝的顺序，只需要把颜色参数改成c("darkred", "white", "darkblue")即可。
• 对称色阶：zlim = c(-max_abs_value, max_abs_value)让颜色映射的范围完全对称于0，这样数值0正好落在色阶的中点，精准对应白色。
• 可选方案：如果你想用RColorBrewer自带的对称配色，比如RdBu系列的11色版本（中间为白色），也可以这样写：

  pal <- brewer.pal(11, "RdBu")
  image(t(mat)[, nrow(mat):1], 
        col = pal, 
        xaxt = "n", 
        yaxt = "n", 
        frame.plot = TRUE, 
        useRaster = TRUE,
        zlim = c(-max_abs_value, max_abs_value))
  

  不过自定义渐变的颜色过渡会更细腻，适合展示连续数值的变化。

内容的提问来源于stack exchange，提问作者Robin Lindström