Eu tenho usado oextract função doraster pacote para extrair dados de arquivos raster usando uma área definida por shapefiles. No entanto, estou tendo problemas com a quantidade de memória que esse processo está exigindo agora. Eu tenho um grande número de shapefiles (~ 1000). Os arquivos raster são grandes (~ 1.6gb)

Meu processo é:

shp <- mclapply(list.files(pattern="*.shp",full.names=TRUE), readShapePoly,mc.cores=6)
ndvi <- raster("NDVI.dat")
mc<- function(y) {
temp <- gUnionCascaded(y)
extract <- extract(ndvi,temp)
mean <- range(extract, na.rm=T )[1:2]
leng <- length(output)
}
output <- lapply(shp, mc)

Há alguma alteração que eu possa fazer para reduzir a carga de memória? Tentei carregar menos shapefiles que funcionaram por cerca de 5 min antes que a memória aumentasse novamente. É um computador quad core 2.4ghz com 8gb de ram