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

Question

假設我有大RasterBrick，我想連續做差，這意味着二層，layer1的，三層，二層，layer4-三層，...... ..

我想用光柵包本，但我無法做到。請問任何人可以給我一個工作（猜測）的代碼嗎？我沒有提供可重複的數據，認爲問題是清楚的。 在此先感謝。

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)