Raster multi-band time series analysis

Hi Victor,

I don't have much experience with beginCluster; however, I have written 
a parallel version of raster::calc a little while ago that uses forking.
https://github.com/dutri001/bfastSpatial/blob/master/R/mc.calc.R

calc and overlay are analogue so that it shouldn't be too hard to extend 
this function to overlay, using an undefined number of rasterStack or 
Bricks as input. You will need a bit more checks and control flow at the 
beginning of the function though.
I'm not sure at all whether this is the most efficient way of 
parallelizing things, but it will for sure speed up your processing 
compared to single core processing.

Hope this helps,
Cheers,
Lo?c

Hi Lo?c,

Thank your very much for your answer. I am working on your suggestion.
But I still have a question.

The function raster::overlay works for multiple RasterBrick objects.
However, I noticed that raster::beginCluster works with raster::overlay
as long as a single RasterStack or RasterBrick is provided as the first
argument. Do you have any suggestion for parallel processing of multiple
RasterBrick objects?

Thank you!

Best,
Victor