如何在柵格堆棧數據中進行連續差異？

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)