Cara menetapkan warna ke nilai raster tertentu menggunakan R

Saya menggunakan levelplotfungsi rasterVispaket luar biasa untuk membuat peta yang nilainya berbeda sekitar nol. Saya merencanakan file raster menggunakan palet merah ke biru, tetapi saya mencoba untuk menetapkan nilai abu-abu ke nol di peta.

Secara khusus, saya ingin mereproduksi warna gambar ini:

Perhatikan bahwa di mana skala merah ke biru digunakan, tetapi nilai nol telah diwarnai dengan abu-abu.

Saat ini, peta saya terlihat seperti ini:

Dan ini adalah kode untuk memperbanyaknya (file tersedia di 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)

Bagaimana saya bisa menetapkan abu-abu ke semua nilai nol di peta saya di atas?

Anda dapat menentukan palet warna Anda sendiri dengan menggabungkan dua palet warna yang telah ditentukan dan warna abu-abu yang Anda inginkan sebagai: [palet kuning oranye merah] + [nol warna abu-abu] + [palet blues] . Untuk mendapatkan nilai nol selaras dengan warna abu-abu Anda harus menggunakan jumlah warna yang sama ( n) untuk palet warna yang telah ditentukan.

Coba kode yang dikomentari di bawah ini:

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

Perhatikan bahwa plot Anda tidak akan terlihat seperti contoh gambar yang Anda poskan karena raster Anda memiliki lebih banyak data yang jarang dengan banyak Nas