如何在光栅堆栈数据中做连续的差异？

Question

假设我有大RasterBrick,我想做连续的差异,这意味着layer2-layer1,layer3-layer2,layer4-layer3,......

我正在考虑使用Raster包,但我无法做到这一点.可以有人建议我使用(猜测)代码吗？我没有提供可重现的数据,认为问题很明显.提前致谢.

Answer 1

以下似乎有效:

library(raster)
r1 = raster()
r2 = raster()
r3 = raster()
dim(r2)
r1[] = runif(180*360)
r2[] = runif(180*360)
r3[] = runif(180*360)
r = stack(r1,r2,r3)
dim(r)
rd = calc(r, fun = diff)
dim(rd)