stack many files without loading into memory
I think you are correct. s=stack(fn,quick=T) writeRaster(s,'localpath/2012data') would get the data local. I guess the trade off is that the file size is an order of magnitude bigger than if I saved them in an .RData file but much quicker to access. ds On Wed, Feb 4, 2015 at 12:51 PM, Dominik Schneider <
dominik.schneider at colorado.edu> wrote:
Wouldn't that keep the link to the server on which they are stored now? On Wed, Feb 4, 2015 at 12:50 PM, Michael Sumner <mdsumner at gmail.com> wrote:
Why not stack(fn) ? On Thu, 5 Feb 2015 06:41 Dominik Schneider < Dominik.Schneider at colorado.edu> wrote:
Hi -
I have some data on a server but would like to bring them local in a
somewhat compressed format that is still easy to access.
/Volumes/hD/2012 -> 100 geotiffs
~/project/data/ -> store those geotiffs here without needing server
access.
untested, I think I could do something like:
s=stack()
writeRaster(s,'2012stack')
fn=list.files('/Volumes/hD/2012',pattern='*.tif',full.names=T)
lapply(fn,function(f){
s=stack('2012stack')
r=raster(f)
names(r)=gsub(pattern='.tif',replacement='',basename(f))
s=addLayer(s,r)
writeRaster(s,'2012stack')
})
Or is it better to save to a .RData?
Is there a better way that doesn't require me to loop through each
geotiff
since I can't load it all into memory.
Thanks
Dominik
[[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo