绘制具有多个图层的栅格数据

Question

我需要创建一个带有多个栅格图层的图形。这些层是二进制数据（0和1），而某些层只有一个值（0或1）。对于只有一个值的图层，图形无法正确显示。

例如，在下图中。第3层应为绿色。

＃

r1 <- raster(nrow=10, ncol=10)
r1[] <- sample(c(0,1),100,replace=TRUE)

tst <- ratify(r1[[1]])
rat <- levels(tst)[[1]]
rat$legend <- c("Absence","Presence")
levels(tst) <- rat

r1 <- ratify(r1)
levels(r1) <- rat

r2 <- raster(nrow=10, ncol=10)
r2[] <- 0
r2 <- ratify(r2)
levels(r2) <- rat

r3 <- raster(nrow=10, ncol=10)
r3[] <- 1
r3 <- ratify(r3)
levels(r3) <- rat

r4 <- raster(nrow=10, ncol=10)
r4[] <- sample(c(0,1),100,replace=TRUE)
r4 <- ratify(r4)
levels(r4) <- rat

r <- stack(r1,r2,r3,r4)
library(rasterVis)
levelplot(r,col.regions=c("grey","darkgreen"))

Answer 1

我将R版本从R 3.0.2更新为R 3.1，并将软件包“raster”和“rasterVis”更新。它现在有效。