clusterR возвращает "аргумент нулевой длины"

Question

Я использовал R raster Пакет долгое время, но сейчас я действительно не могу разобраться с этой проблемой clusterR. Я должен вычислить индекс SPI для растра netCDF. Это делается для каждой ячейки с использованием временных рядов ячейки и возвращением временных рядов индекса SPI для этой ячейки.

Пример входного файла (около 4 МБ) можно найти здесь.

Смотрите код ниже:

library(raster)
library(SPEI)

calcspi <- function(pr) { #this function calculates the SPI index for each timeseries of values
    pr <- as.numeric(pr)
    if (all(is.na(pr[1:20]))) { #Check that this is not an NA cell
        outspi <- rep(NA, length(pr))
    } else {
        outspi <- fitted(spi(pr, 12, na.rm=TRUE))
    }
    return(outspi)
}

b <- brick("input_crop.nc", varname="pr")
readAll(b) #As requested in the comments

###THIS WORKS BUT IS SLOW:
bc <- calc(b, calcspi)

###THIS DOES NOT:
beginCluster(n=4)

bc <- clusterR(b, calc, args=list(fun="calcspi"))
#[1] "argument is of length zero"
#attr(,"class")
#[1] "snow-try-error" "try-error"
#Error in clusterR(b, calc, args = list(fun = "calcspi")) : cluster error

endCluster()

###THIS DOESN'T EITHER:
beginCluster(n=4)

f <- function(x) calc(x, calcspi)
bc <- clusterR(b, f)
#[1] "argument is of length zero"
#attr(,"class")
#[1] "snow-try-error" "try-error"
#Error in clusterR(b, f) : cluster error

endCluster()

traceback() совершенно бесполезен в этом случае. Что случилось?

Answer 1

Это сработало для меня:

b <- mybrick
#readAll(b) #As requested in the comments
#parallel processing 
ff <- function(x) calc(x, calcspi)
beginCluster(8)
bc <- clusterR(b, fun = ff,export='calcspi')
endCluster()