Как сделать последовательную разницу в данных растрового стека?

Предположим, что у меня большой RasterBrick, и я хотел бы сделать последовательное различие, что означает слой2-слой1, слой3-слой2, слой4-слой3, …… ..

Я думал использовать для этого пакет Raster, но не могу этого сделать. Может ли кто-нибудь предложить мне рабочий (угадывающий) код, пожалуйста? Я не предоставляю воспроизводимые данные, полагая, что проблема ясна. Заранее спасибо.


person WAW    schedule 06.03.2016    source источник


Ответы (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)
person Edzer Pebesma    schedule 06.03.2016