r - 如何使用 spplot 绘制额外的栅格？

Question

我想将 SpatialPolygonsDataFrame 绘制为半透明的主对象(右侧有图例)，但我想绘制额外的栅格(山体阴影)作为背景 - 只是为了制作更好的 map 。我需要这样的东西:

spplot(polygons, sp.layout = list(list("raster", myRaster)))

但查看?spplot，似乎无法在sp.layout 中指定栅格。我不能将栅格指定为主要对象，因为主要对象是 polygons dataFrame，我想在右侧的图例中绘制其值比例。

如何在 spplot 中绘制额外的栅格？

Answer 1

这是一种方法。不过，可能有一种更简洁的方法可以在不绘制多边形对象两次的情况下实现它……

library(sp)
library(rasterVis)
r <- raster(nrow=18, ncol=36)
r[] <- runif(ncell(r)) * 10
r[r > 8] <- NA
pol <- rasterToPolygons(r, function(x) x > 6)
spplot(pol) + levelplot(r) + spplot(pol)

或者:

library(latticeExtra)
spplot(pol) + spplot(r) + spplot(pol)

---

编辑

根据@OscarPerpiñán 的评论，更好的方法是:

spplot(pol) + as.layer(spplot(r), under = TRUE)