栅格堆栈的密度图

Question

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

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 参数来提供您喜欢的主题。 Here 您将找到有关 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)

栅格堆栈的密度图

Density plot of raster stack

r

r-raster

rastervis