R中的栅格图例：如何为特定值上色？

我在R中绘制了一个高程模型

r <- raster("example.dem")
plot(r, col = topo.colors(20))

其中一些值低于海平面（0），我希望这些值显示为红色。如何在plot（）中为特定颜色分配特定范围？

这是一个非常简单的例子

library(raster)
data(volcano)
volcanoR <- raster(volcano)

#making colors below 100 red and above 180 blue in this example

breakpoints <- c(94,100,120,140,160,180,195)
colors <- c("red","white","white","white","white","blue")
plot(volcanoR,breaks=breakpoints,col=colors)

您只需要向图传递一个断点矢量和一个颜色矢量以匹配断点。请检查RColorbrewer包装，以获取一些非常好的内置色带。还要检查classInt包装以制作断点。