Używamextract funkcja zraster pakiet do wyodrębniania danych z plików rastrowych przy użyciu obszaru zdefiniowanego przez pliki kształtu. Mam jednak problemy z ilością pamięci wymaganą przez ten proces. Mam dużą liczbę plików shapefiles (~ 1000). Pliki rastrowe są duże (~ 1.6 gb)

Mój proces to:

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)

Czy są jakieś zmiany, które mogę wprowadzić, aby zmniejszyć obciążenie pamięci? Próbowałem załadować mniej plików shapefiles, które działały przez około 5 minut, zanim pamięć ponownie się powiększyła. Jego komputer quad core 2,4 ghz z ramieniem 8 gb