Comment attribuer des couleurs à des valeurs raster spécifiques à l'aide de R

J'utilise la levelplotfonction du rasterVispackage génial pour créer une carte dont les valeurs divergent autour de zéro. Je trace un fichier raster à l'aide d'une palette rouge à bleu, mais j'essaie d'attribuer des valeurs de gris à zéro dans la carte.

Plus précisément, je veux reproduire les couleurs de cette figure:

Notez que lorsqu'une échelle de rouge à bleu est utilisée, mais que les valeurs nulles ont été colorées en gris.

Actuellement, mes cartes ressemblent à ceci:

Et voici le code pour le reproduire (fichier disponible sur 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)

Comment puis-je attribuer du gris à toutes les valeurs nulles dans ma carte ci-dessus?

Vous pouvez définir votre propre palette de couleurs en concaténant deux palettes de couleurs prédéfinies et la couleur grise de votre choix: [palette jaune orange rouge] + [couleur gris zéro] + [palette bleus] . Pour aligner la valeur zéro sur la couleur grise, vous devez utiliser le même nombre de couleurs ( n) pour les palettes de couleurs prédéfinies.

Essayez le code commenté ci-dessous:

# 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})))

Notez que votre tracé ne ressemblera pas à l'exemple de figure que vous avez publié, car votre raster contient plus de données clairsemées avec de nombreux Nas