同一目录下的多个文件顺序提取
前两部分讲解了如何提取一个nc文件的数据,以及如何算平均值。但是在很多情况下我们所收集到的资料为多个nc文件,且是大量的。这时候一个一个文件进行提取显然是不明智的,对于这种情况,根据网上大神的讲解。经由自己修改以后,写出了一个对于同一目录多个nc文件(不止限于nc文件各种文件通用)顺序读取且进行合成的R语言代码供大家学习。
fun <- function(x){
filenames <- list.files(x)
path <- paste(x,filenames,sep="")
len <- length(path)
c <- c()
for(i in 1:len){hdf5_data <- h5read(path[i],name = "burned_area")
#读取hdf5文件,这里可以换成读取各种文件格式的代码
for (m in 1:12){h5burn <- hdf5_data[[m]]
h5burned <- h5burn$burned_fraction
b <- h5burned[1165:1172,220:226]
c <- rbind(c,b)#这几步都与前两节读取一个nc文件形式相同
}
}
return(c)
}