Advertisement
Guest User

Untitled

a guest
Jun 29th, 2016
63
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.64 KB | None | 0 0
  1. r <- as.data.frame(cellStats(x,mean))
  2.  
  3. library(raster) #load raster package
  4.  
  5. rasterlist <- list("path 1","path 2",...) #create list of raster file paths
  6. outlist <- list() #create empty list to store outputs from loop
  7.  
  8. for (i in 1:length(rasterlist)) { #for each raster in rasterlist
  9. r <- raster(rasterlist[[i]]) #read element i of rasterlist into R
  10. val <- getValues(r) #get raster values
  11. m <- mean(val,na.rm=T) #remove NAs and compute mean
  12. outlist[[i]] <- c(rasterlist[[i]],m) #store raster path with mean
  13. return("complete")
  14. }
  15.  
  16. df <- data.frame(do.call(rbind,outlist)) #convert list to data frame
  17. colnames(df) <- c("raster path","mean")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement