R中光柵圖中黑色的基圖NA值

Question

我有一個光柵化的 S4 對象，它是通過在數據集上使用 rasterize 函數創建的。 放入其中的數據集有一些 NA 值。 在 R 中使用基本繪圖函數時，我想用黑色描繪這些 NA 單元格。

嘗試了以下代碼：

library(RColorBrewer)

my.colors = colorRampPalette(rev(brewer.pal(9, "RdBu")))

plot(rasterized object, col= ifelse(is.na(rasterized object@data@values),"black", my.colors(255)), box=F, add=F)

# col(ifelse) 的靈感來自使用 NA 數據繪制條件顏色

這不起作用，許多具有真實值的網格單元被用黑色描繪。

沒有錯誤消息。

鑒於我想使用R 的基本繪圖語法而不是 ggplot ，可以改變什么？

Answer 1

我猜你正在使用raster包並且光柵化的 S4 對象是一個RasterLayer

如果是這樣，光柵包（版本 3.4.5）的plot函數有一個colNA參數。

plot(r, col = my.colors(255), colNA = "black")

這是一個可重現的示例：

library(raster)
library(RColorBrewer)
r <- raster(nrows=10, ncols=10)
r <- setValues(r, c(rep(NA, 10), 1:90)) 

my.colors = colorRampPalette(rev(brewer.pal(9, "RdBu")))

plot(r, col = my.colors(255), colNA = "black")