How to extract values from rasterstack with xy coordinates?

前端 未结 2 1109
终归单人心
终归单人心 2020-12-16 05:23

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         


        
2条回答
  •  太阳男子
    2020-12-16 06:12

    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") 
    

提交回复
热议问题