R kullanarak belirli raster değerlerine renk atama

Ben değerleri sıfıra yakın bir harita oluşturmak levelplotiçin harika rasterVispaket işlevini kullanıyorum . Kırmızıdan maviye palet kullanarak bir tarama dosyası çiziyorum, ancak haritadaki sıfıra değerlere gri atamaya çalışıyorum.

Özellikle, bu rakamın renklerini çoğaltmak istiyorum:

Kırmızıdan maviye bir ölçek kullanıldığında, ancak sıfır değerlerinin gri ile renklendirildiğine dikkat edin.

Şu anda haritalarım şöyle görünüyor:

Ve bu kodu yeniden oluşturmak için kod (dosya https://www.dropbox.com/s/cypfdu1eaz2fuok/r.annual.tif?dl=0 adresinde bulunabilir ):

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

Yukarıdaki haritamdaki tüm sıfır değerlere nasıl gri atayabilirim?

Önceden tanımlanmış iki renk paletini ve istediğiniz gri rengi şu şekilde birleştirerek kendi renk paletinizi tanımlayabilirsiniz: [kırmızı turuncu sarı palet] + [sıfır gri renk] + [blues paleti] . Sıfır değerin gri renkle hizalanmasını nsağlamak için önceden tanımlanmış renk paletleri için aynı sayıda renk ( ) kullanmanız gerekir.

Aşağıdaki yorum kodunu deneyin:

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

Grafiğinizin yayınladığınız şekil örneğine benzemediğine dikkat edin, çünkü rasterinizde birçok Nas ile daha seyrek veriler var