Ich habe das benutztextract Funktion von derraster Paket zum Extrahieren von Daten aus Rasterdateien in einem durch Shapefiles definierten Bereich. Ich habe jedoch Probleme mit der Speichermenge, die dieser Prozess jetzt benötigt. Ich habe eine große Anzahl von Shapefiles (~ 1000). Die Rasterdateien sind groß (~ 1,6 GB)

Mein Prozess ist:

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)

Kann ich Änderungen vornehmen, um die Speicherauslastung zu verringern? Ich habe versucht, weniger Shapefiles zu laden, was ungefähr 5 Minuten gedauert hat, bevor sich der Speicher wieder auflöste. Es ist ein Quad-Core-Computer mit 2,4 GHz und 8 GB RAM