r - 栅格堆栈的密度图

Question

我想使用rasterVisR-CRAN 中的包来实现栅格堆栈的密度图。我想更改表示栅格值分布的颜色，并添加一个图例。使用这个包可以做这些事情吗？我没有找到任何例子来做到这一点。这是我的示例代码

require(raster)
require(rasterVis)
slogo <- stack(system.file("external/rlogo.grd", package="raster")) 
densityplot(slogo,col=c("orange","skyblue1","cyan"))

Answer 1

大多数rasterVis功能都是基于lattice包的。因此，您应该使用par.settings提供您首选主题的参数。在这里，您将找到有关lattice主题的更多信息。

关于您的问题，请尝试以下代码：

library(rasterVis)
f <- system.file("external/test.grd", package="raster")
r <- raster(f)
s <- stack(r, r+500, r-500)
densityplot(s, par.settings = theEconomist.theme)