Рассчитать тренд в большом растровом стеке в R

Question

Я хочу вычислить тренд растрового стека в R, используя ежедневные данные за 50 лет для всей Европы.

time <- 1:nlayers(gimms.sum) 
fun=function(x) { if (is.na(x[1])){ NA } else { m = lm(x ~ time); summary(m)$coefficients[2] }}
gimms.slope=calc(gimms.sum, fun)

Код выше был взят отсюда [1]. Тем не менее, я получаю сообщение, что вектор 8,6 ГБ и не может выделить память

[1] https://matinbrandt.wordpress.com/2013/11/15/pixel-wise-time-series-trend-anaylsis-with-ndvi-gimms-and-r/

Как я могу решить эту проблему?

Answer 1

Я предполагаю, что вы используете растровую версию 2.7-15. В этой версии есть ошибка в настройке управления памятью. Вы можете исправить это с помощью rasterOptions(maxmem=1e09); или используйте версию 2.8-4 (достигла CRAN сегодня, но еще не скомпилирована для win и mac)