如何在不生成断点的情况下修改栅格图的连续图例？

这是一张年均温地图，单位为摄氏度/10，数据取自CHELSA V1数据集。

我尝试了以下代码绘制地图：

min_val <- -15000
max_val <- 0

plot(r, col = rev(terrain.colors(10)), at = seq(min_val, max_val, length.out = 11))

但操作后仅生成了断点，地图显示无变化，图例仅出现9个断点，颜色渐变效果丢失。

尝试以下代码也无法正常工作：

plot(r, range=c(-20000, 0))

运行后得到空白图。

运行以下代码同样得到空白图：

min_val <- -10
max_val <- 0
plot(r, range=c(min_val, max_val))

先确认栅格数据r的实际数值范围，运行以下代码：

summary(r)

这能帮你明确数据的真实最小值和最大值，避免设置的参数和实际数据完全不匹配导致异常。

修复颜色渐变丢失问题 #

你设置的at参数是11个断点（对应10种颜色），但图例显示9个断点，大概率是数据里没有落在部分区间的数值。建议换用更多颜色数来提升渐变自然度，同时确保断点覆盖数据真实范围，也可以试试更适合温度数据的配色：

# 先安装viridis包（未安装时执行）
# install.packages("viridis")
library(viridis)

# 用summary(r)得到的真实最小值替换下面的real_min
real_min <- ... 
real_max <- 0

plot(r, col = viridis(100), at = seq(real_min, real_max, length.out = 101))

修复空白图问题 #

空白图说明你设置的range完全不在数据实际数值范围内。比如你的数据单位是摄氏度/10，真实数值应该是类似-150到0（对应-15℃到0℃），你之前设置的-20000或-10都和真实范围不匹配，导致没有像素落在指定区间。

用summary(r)拿到真实范围后，再调整参数：

# 假设summary得到的最小值是-120（对应-12℃）
plot(r, range = c(-120, 0), col = rev(terrain.colors(20)))

如果用raster包，也可以试试rasterVis包的levelplot，对连续栅格可视化更灵活：

# install.packages("rasterVis")
library(rasterVis)
levelplot(r, col.regions = rev(terrain.colors(100)), at = seq(real_min, real_max, length.out = 101))

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