我有一个每年25个土地覆盖栅格的时间序列。由于这是分类数据,我使用levelplot(inputRaster)
( rasterVis库的一部分)绘制一个栅格。但是,我想像光栅库的函数一样,按顺序绘制年度光栅。当我用
rasStack <- stack(listOfRasters) animate(rasStack)
结果没有一个明确的传说。简单地说:如何将levelplot
和animate
的功能结合起来
发布于 2017-12-07 11:23:03
函数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")
https://stackoverflow.com/questions/47698349
复制相似问题