Leyenda raster en R: ¿cómo colorear valores específicos?

Question

Tengo un modelo de elevación trazado en R

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

Algunos de los valores están por debajo del nivel del mar (0), y me gustaría que aparecieran en rojo. Cómo puedo asignar rangos específicos a colores específicos en plot()?

Answer 1

He aquí un ejemplo muy sencillo

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)

Sólo tiene que pasar al gráfico un vector de puntos de ruptura y un vector de colores que coincidan con los puntos de ruptura. Echa un vistazo a la RColorbrewer paquete para algunas rampas de color muy bonito construido. Compruebe también el classInt para crear los breakpoints.