如何使用spplot绘制其他栅格？

Question

我想将SpatialPolygonsDataFrame绘制为一个半透明的主要对象（右侧有图例），但我想将其他栅格（山体阴影）作为背景绘制 - 只是为了制作更好的地图。我需要类似的东西：

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)