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

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)