Re: [R-sig-Geo] Problem with NDVI difference with subset image

2013-06-24 Thread Sarah Goslee
Hi, On my system, it's a writeRASTER() options issue. I don't use raster very often, but I played with the options enough to come up with a working solution. ## this doesn't work in the right way - something about the format ## used by writeRaster isn't compatible with lssub() #Create geotiff w

Re: [R-sig-Geo] Problem with NDVI difference with subset image

2013-06-23 Thread ASANTOS
Hi Dra. Goslee, I work in OS linux and I like very much your lssub function, despites your notice. I think that the problems are in linux system, because: > landc<- stack(c("stackIm1.sample.tif","tackIm2.sample.tif")) Error in .local(.Object, ...) : `/home/asantos/Documentos/Sen

Re: [R-sig-Geo] Problem with NDVI difference with subset image

2013-06-23 Thread Sarah Goslee
Hi, What happens? Do you get an error? Where? What is your sessionInfo()? As it says in the lssub() help, this function was written for a particular purpose, is only known to work on linux, and may not be widely applicable. It's a "use at your own risk" kind of function, and you may well be be

[R-sig-Geo] Problem with NDVI difference with subset image

2013-06-23 Thread Alexandre Santos
Dear Members,     I'm having trouble calculating NDVI difference, first the images Geotiff can not view using Ubuntu 4.12 64-bit and therefore can not get to the NDVI, follow an example: require(raster) require(sp) require(rgdal) require(landsat) # #Create raster r <- raster(nc=1000, nr=10