如何在terra包中绘制栅格并将0值固定为白色

针对你的需求，核心是强制绘图时使用对称的数值范围（以0为中心），同时让调色板的中间色对应白色。以下是两种可行的方法：

方法一：指定zlim强制对称范围 #

直接在plot()中用zlim参数锁定你计算的对称范围，确保颜色映射基于这个范围，让0正好对应调色板的中间（白色）：

# 准备环境与数据（你的原有代码）
library(terra)
library(magrittr)

set.seed(1234)
rr <- rast(matrix(rnorm(400, 1.5, 1), nrow=20, ncol=20))

# 计算以0为中心的对称范围
the_range <- range(values(rr)) %>% abs %>% max %>% multiply_by(c(-1,1))

# 创建调色板（反转RdBu，让负值蓝、正值红）
the_palette_fc <- leaflet::colorNumeric(palette = "RdBu", domain = the_range, reverse = TRUE)

# 生成对应颜色序列
colors <- the_palette_fc(seq(the_range[1], the_range[2], length.out=31))

# 绘图时指定zlim，强制使用对称范围
plot(rr, col = colors, zlim = the_range, legend = TRUE)

方法二：自定义breaks精确控制颜色区间 #

如果需要更精细的区间控制，可以手动创建对称的断点，确保0作为分割点之一，让颜色严格对应正负区间：

# 沿用前面的rr、the_range对象
n_colors <- 31
# 创建对称断点（n个颜色对应n+1个断点）
breaks <- seq(the_range[1], the_range[2], length.out = n_colors + 1)

# 生成对应颜色
pal <- leaflet::colorNumeric(palette = "RdBu", domain = the_range, reverse = TRUE)
colors <- pal(seq(the_range[1], the_range[2], length.out = n_colors))

# 绘图时指定breaks和颜色
plot(rr, breaks = breaks, col = colors, legend = TRUE)

为什么之前的方法无效？ #

1. 步骤4未指定zlim：terra::plot()默认会使用栅格的实际数值范围（而非你计算的对称范围），导致颜色映射的中间点偏移，白色无法对应0。
2. 使用coltab错误：coltab仅适用于离散型数值栅格（整数类型），你的栅格是连续浮点型，因此无法正确映射。

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