如何在R中使用colorRampPalette实现非对称正负值色阶？

你的核心问题是非对称数据范围下，直接用均匀渐变的色阶导致负值区域颜色区分度低，且0点未对齐到红-蓝的分界（白色）。下面先解释colorRampPalette的用法，再给出修正代码。

关于colorRampPalette的用法 #

colorRampPalette()是用来生成颜色渐变函数的工具：

• 输入一组参考颜色（比如brewer.pal(11, "RdBu")的11种颜色），返回一个新函数；
• 这个新函数接受整数n，会生成n种颜色，按输入参考颜色的顺序做线性插值，形成平滑渐变。

你之前的代码问题：
brewer.pal(11, "RdBu")本身是红→白→蓝的渐变序列，rev()后变成蓝→白→红，但直接生成20种均匀渐变的颜色，完全没考虑你的数据区间分布（17个负值区间，3个正值区间），导致蓝色系颜色被压缩，负值区域难以区分，同时0点没有准确对应到白色分界。

修正后的代码 #

library(RColorBrewer)
library(fields)
library(maps)

# 定义断点
brks <- seq(-85, 15, 5)
# 统计负值、正值区间的数量
neg_interval_num <- sum(brks < 0)  # 17个负值区间（从-85~-80到-5~0）
pos_interval_num <- sum(brks > 0)  # 3个正值区间（从0~5到10~15）

# 构建蓝→白→红的参考调色板（反转RdBu得到正确的蓝到红顺序）
rdbu_blue_to_red <- rev(brewer.pal(11, "RdBu"))

# 生成负值区间的颜色：从深蓝到浅白，共17种
neg_colors <- colorRampPalette(rdbu_blue_to_red[1:6])(neg_interval_num)
# 生成正值区间的颜色：从浅白到深红，共3种
pos_colors <- colorRampPalette(rdbu_blue_to_red[6:11])(pos_interval_num)

# 合并颜色，去掉重复的白色（避免色阶出现两次白）
colorscale <- c(neg_colors[-length(neg_colors)], pos_colors)

# 绘图
windows()
image.plot(x=lon, y=lat, inputvar,
           xlim=c(12,25), ylim=c(44,50),
           col=colorscale, breaks=brks, lab.breaks=brks,
           xlab="Longitude (deg)", ylab="Latitude (deg)")
title(main="Annual Bias Difference Summer Days")
mtext("Actual Difference (days)", side=4, adj=0.5, line=0.4)
map(database="world", add=TRUE)
abline(h=seq(44,50,1), v=seq(12,24,2), lty=2, col="grey", lwd=0.5)

代码说明 #

1. 区间统计：先算出负值和正值各自的区间数量，确保颜色分配和数据区间匹配；
2. 调色板拆分：把RdBu反转成蓝→白→红的序列，拆分出「蓝到白」「白到红」两段，分别对应负值和正值区域；
3. 渐变生成：针对负值的17个区间生成足够的蓝色渐变，针对正值的3个区间生成红色渐变，保证每个区间的颜色区分度；
4. 去重合并：去掉两段渐变重复的白色，避免色阶出现冗余分界。

内容的提问来源于stack exchange，提问作者kettaren ray