I have to reclassify a large (> 8000x8000 pixels) geotif raster. I normally would do it through grass (or saga) but would like to try the raster package in R. The problem is that I'm a bit confused on how to do it. My understanding of the package is that it provides a way to operate in the rasters by pieces, not having to read the whole thing into memory. I've seen there's a function reclass() but how do I actually read the raster? If I use readAll() this is going to be as using readGDAL() directly, so this is not the way. Do I have to use a for() of readRows(), then reclass()? This seems a bit low-level. Is there a way of assigning the tif file to an R raster object and then just use reclass() and R would read the file by pieces as needed? That's what I would like to...
Agus
begin:vcard fn:Agustin Lobo n:Lobo;Agustin org:Institut de Ciencies de la Terra "Jaume Almera" CSIC adr:;;Lluis Sole Sabris s/n;Barcelona;;08028;Spain email;internet:agustin.l...@ija.csic.es url:http://www.ija.csic.es/gt/obster version:2.1 end:vcard
_______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo