使用ggplot2绘制固定色阶填充等高图：超范围值显示问题

核心问题原因 #

ggplot2的连续填充色阶默认会将超出limits范围的数值标记为NA，导致这些区域无法正常显示（通常呈现空白或默认背景色）。要解决这个问题，需要通过oob（out of bounds）参数指定超出范围数值的处理方式。

解决方案：使用scales::squish压缩超出范围的值 #

scales包中的squish函数可以将超出色阶范围的数值映射到色阶的端点，这样超出区域会显示为色阶的最小/最大颜色，而非空白。

代码示例 #

假设你已经将栅格数据转换为包含x、y、value和分面变量的数据框（比如用raster::as.data.frame(your_raster, xy = TRUE)），修改后的代码如下：

基础连续色阶版本

library(ggplot2)
library(scales) # 需加载scales包

ggplot(data = raster_df, aes(x = x, y = y, fill = value)) +
  geom_raster() + # 或geom_tile，根据数据类型选择
  facet_wrap(~your_facet_var) + # 替换为你的分面变量
  scale_fill_continuous(
    limits = c(-40, 40), # 固定色阶范围
    oob = squish, # 关键：压缩超出范围的值到边界
    breaks = seq(-40, 40, 20), # 可选：自定义色阶断点
    labels = c("≤-40", "-20", "0", "20", "≥40") # 可选：明确标注超出范围的端点
  ) +
  labs(fill = "数值") # 可选：修改图例标题

自定义渐变颜色版本

如果需要使用多色渐变（比如蓝-白-红），可以用scale_fill_gradientn：

ggplot(data = raster_df, aes(x = x, y = y, fill = value)) +
  geom_raster() +
  facet_wrap(~your_facet_var) +
  scale_fill_gradientn(
    colors = c("darkblue", "white", "darkred"), # 自定义渐变颜色
    limits = c(-40, 40),
    oob = squish
  )

额外说明 #

• squish函数会将所有小于-40的数值统一映射到-40对应的颜色，大于40的数值映射到40对应的颜色，保证所有区域都有颜色填充。
• 如果需要更灵活的超出范围处理（比如用特定颜色标记超出区域），可以自定义oob函数，但squish是最常用的解决方案。

内容的提问来源于stack exchange，提问作者Prasad Thota