I have a rasterstack (5 raster layers) that actually is a time series raster.
r <- raster(nrow=20, ncol=200) s <- stack( sapply(1:5, function(i) setVa
Thanks a lot @SimonO101 !
This code works.
r <- raster(nrow=10, ncol=10) s <- stack( sapply(1:5, function(i) setValues(r, rnorm(ncell(r), i, 3) )) ) s[1:3]<-NA vals <- values(s) coord <- xyFromCell(s,1:ncell(s)) combine <- cbind(coord,vals) write.table(combine,"xyvalues.txt")