[R-sig-Geo] Mask from raster with less RAM

[Frederico Faleiro]

Dear all,

I would like to generate a mask from a raster, but my workflow needs a lot
of RAM to process big rasters. I need this polygon mask to use in another's
rasters.
Do you know another approach that needs less RAM?

# reproducible example
library(raster)
# read data to create mask
r <- raster(system.file("external/test.grd", package="raster"))
r[!is.na(r)] <- 1
pol <- rasterToPolygons(r, dissolve = T) #  a lot of RAM to process the data
# apply the mask in another raster
 r2 <- raster(extent(r), res(r))
r2[ ] <- 1
r2.mask <- mask(r2, pol)

Cheers!

-- 
Frederico Faleiro
Postdoctoral Researcher in the INCT-EECBio (https://www.eecbio.ufg.br/)
Department of Ecology | Institute of Biological Science | Federal
University of Goiás | Brazil
RG: https://www.researchgate.net/profile/Frederico_Faleiro
CV: http://lattes.cnpq.br/4926404840659003

	[[alternative HTML version deleted]]