计算R中大型栅格堆栈中的趋势

Question

我想使用整个欧洲50年的每日数据来计算R中的栅格堆栈的趋势。

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 Gb，无法分配内存

[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进行编译）