So weisen Sie bestimmten Rasterwerten mit R Farben zu

Ich benutze die levelplotFunktion des fantastischen rasterVisPakets, um eine Karte zu erstellen, deren Werte um Null herum abweichen. Ich zeichne eine Rasterdatei mit einer Rot-Blau-Palette, versuche jedoch, Nullwerte in der Karte Grau zuzuweisen.

Insbesondere möchte ich die Farben dieser Figur reproduzieren:

Beachten Sie, dass dort, wo eine Rot-Blau-Skala verwendet wird, die Nullwerte jedoch grau gefärbt wurden.

Derzeit sehen meine Karten folgendermaßen aus:

Und dies ist der Code, um ihn zu reproduzieren (Datei verfügbar unter https://www.dropbox.com/s/cypfdu1eaz2fuok/r.annual.tif?dl=0 ):

# Load required packages
library(rasterVis)

# open file
r.annual <- raster("Downloads/r.annual.tif")

# Set color palette
myTheme=rasterTheme(region=brewer.pal('RdBu', n=11))

# Plot
levelplot(annual.mask, par.settings=myTheme, margin=F)

Wie kann ich allen Nullwerten in meiner obigen Karte Grau zuweisen?

Sie können Ihre eigene Farbpalette definieren, indem Sie zwei vordefinierte Farbpaletten und Ihre gewünschte Graufarbe wie folgt verketten: [Rot-Orange-Gelb-Palette] + [Null-Grau-Farbe] + [Blau-Palette] . Um den Nullwert an der grauen Farbe auszurichten, müssen Sie nfür die vordefinierten Farbpaletten die gleiche Anzahl von Farben ( ) verwenden.

Versuchen Sie den kommentierten Code unten:

# Load libraries
library('rasterVis')

# open raster file
r.annual <- raster("r.annual.tif")

# Set color palette
zeroCol <-"#B3B3B3" # (gray color, same as your figure example)
reds <- rev(brewer.pal('YlOrRd', n = 7))
blues <- brewer.pal('Blues', n = 7)

myTheme <- rasterTheme(region = c(reds, zeroCol, blues))

# Plot
levelplot(r.annual, par.settings = myTheme, margin = FALSE, main = expression("Precipitation" ~ (mm ~ year^{-1})))

Beachten Sie, dass Ihr Plot nicht wie das von Ihnen veröffentlichte Beispiel aussieht, da Ihr Raster mit vielen Nas spärlichere Daten enthält