R Trama paquete de recorrer los pixeles

Question

Quiero realizar un análisis de tendencia para cada píxel sobre la dimensión de tiempo. Por lo tanto, escribí un bucle en cada píxel de hacer el cálculo, pero esto es extremadamente lento. Estoy seguro de que debe haber una manera más rápida de realizar el cálculo.

También las dimensiones del archivo se 1440,720,480. Es este un factible tamaño de R? El problema es que no puedo encontrar la función estadística implementada en IDL que de lo contrario utilizar.

gridfile = 'test.nc'
data = brick(gridfile)

require(fume)

mk = matrix(numeric(0), 720,1440) 

for(i in 1:72) {
  for(j in 1:144) {
  temp = mkTrend(as.vector(data[i,j,]))
  mk[i,j] = temp$`Sen's Slope`
  }
}

Answer 1

Paquete raster proporciona una función calc que se ejecuta una función sobre toda la trama de píxeles, volviendo a la trama con los resultados; se puede tratar con grandes rásteres (de memoria). raster también proporciona cluster/multi-núcleo de la funcionalidad, aspecto en función de beginCluster; tiene la intención de trabajar con calc.