r - 将 animate() 与 R 栅格中的一系列水平图一起使用

Question

我有 25 个年度土地覆盖栅格的时间序列。由于这是分类数据，我使用levelplot(inputRaster)（rasterVis 库的一部分）来绘制单个栅格。但是，我想按栅格库的功能顺序绘制年度栅格。animate当我使用

rasStack <- stack(listOfRasters) animate(rasStack)

结果没有分类图例。简而言之：我如何结合 和 的levelplot功能animate？

Answer 1

函数animate只接受栅格对象作为输入。您可以尝试saveGIF为级别图设置动画：

library(raster)
library(rasterVis)
library(animation)
library(classInt)

r <- raster(ncol=40, nrow=20)
r[] <- rnorm(n=ncell(r))
s <- stack(x=c(r, r*r, r*r*r, r*r*r*r))

classes <- classIntervals(values(r), n=5, style="fisher", precision = 3)
brks <- classes$brks
brks <- round(brks, 2)

saveGIF({
  for(i in c(1:nlayers(s))){
    l <- levelplot(s[[i]], colorkey=list(at=brks, labels=c(as.character(brks))), margin=FALSE)
    plot(l)
  }
}, interval=0.2, movie.name="animation.gif")