RasterBrick各時間間隔 (不規則) でバイオマスを表す時系列のラスター レイヤーを含む大きなオブジェクトがあります。私が必要としているのは、2 つの連続した期間の間のバイオマスの差の時系列です ( difference in biomass = current biomass layer - last weeks biomass layer)。私の考えでは、ループまたは適用関数の 1 つを使用して、時系列の前の関数と減算の関数RasterBrickをそれぞれ適用します。RasterLayerはRasterBrick順序付けられているため、実際のタイム スタンプは必ずしも重要ではありません。例を見つけようとしましたが、進歩に失敗しました。任意のポインターは非常に高く評価されます。
以下に私の状況の簡単な例を示します。
library(raster)
RasterStack への 10 個のラスタ レイヤのランダム セット
r <- raster(ncol=10, nrow=10)
;b <- brick( sapply(1:10, function(i) setValues(r, rnorm(ncell(r), i, 3))))
今RasterBrick、結果を含む9つのレイヤーが必要です
バイオマスの差 = 現在のバイオマス層 - 先週のバイオマス層。