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")