Re: [R-sig-Geo] reading tiff with jump?

2008-08-05 Thread Agustin Lobo
Thanks to your input, this is what I'm doing and I think it's ok: require(rgdal) fname <- "/media/SWISNIFE1/GEODATABOLIVIAmini/geotif/p233_r070_ETM_20010825.tif" GDALinfo(fname=fname) output of GDALinfo: rows3736 columns 5108 bands 6 ll.x694431 ll.y8410014 res

Re: [R-sig-Geo] reading tiff with jump?

2008-08-05 Thread Ashton Shortridge
On Tuesday 05 August 2008, Roger Bivand wrote: > On Tue, 5 Aug 2008, Agustin Lobo wrote: > > Hi! > > > > is it possible to read a tiff file with readGDAL but just keeping > > one pixel every p columns and q lines? That is, > > a systematic sampling at reading. This is good to save memory > > and st

Re: [R-sig-Geo] reading tiff with jump?

2008-08-05 Thread Roger Bivand
On Tue, 5 Aug 2008, Agustin Lobo wrote: Hi! is it possible to read a tiff file with readGDAL but just keeping one pixel every p columns and q lines? That is, a systematic sampling at reading. This is good to save memory and still be able of displaying and performing statistics. As images are la

[R-sig-Geo] reading tiff with jump?

2008-08-05 Thread Agustin Lobo
Hi! is it possible to read a tiff file with readGDAL but just keeping one pixel every p columns and q lines? That is, a systematic sampling at reading. This is good to save memory and still be able of displaying and performing statistics. As images are large, reading in the entire image and then