Hello,

I came across an interesting phenomenon when dealing with a large number of
rasters sets (~1300 with each one being 250x250).  What happened was that I
was unable to create the brick using a list of rasters due to time
constraints whereas when creating a brick from an 3D array the
solution comparatively was lightning fast.  Here is an example:

library(raster)

library(raster)
V <- rnorm(5000000)
A <- array(V, c(100,100,500))

L <- unlist(apply(A, 3, list), recursive = FALSE)
L <- lapply(L, raster)

system.time(A <- brick(A))
#   user  system elapsed
#   0.48    0.10    0.58

system.time(L <- brick(L))
#   user  system elapsed
#  44.80    0.22   46.21

Is this huge speed difference the expected result or am I missing something?

Thank you,
-Jonathan

	[[alternative HTML version deleted]]

