nc文件格式数据批量提取求均值（R语言）

nc文件格式数据批量提取求均值

install.packages("netcdf4")#如果已经下载可略过

library(ncdf4)

ncfile <- nc_open("cru_ts4.05.1951.1960.pre.dat.nc")#打开nc文件

pre<- ncvar_get(ncfile,"pre")#提取向量pre

dim(pre)#pre的维度最后一列用来确定循环的次数
c <- c()#这一步很重要，先确定空的cell,后面才能批量叠加，这个在好多地方都可以运用到。
for(i in 1:120){bb<- pre[ , ,i]

b <- bb[553:598,245:264]

c <- rbind(c,b)} 
#for循环代码i为时间即为维度的最后一列，而 bb[553:598,245:264]表示自己需要提取的nc的部分。
d <- rep(1:120,each=46)#分组的序号each为598-553+1。即为提取的一个时间维度的行数。
e <- data.frame(d,c)#把序号加载进入数据框内
e[is.na(e)] <- 0
#对数据中的NA赋值为0，也可以进行插值都可，特殊情况才会出现NA，绝大部分需要研究的区域的nc格式的数据都有数值。
f <- split(e,e$d)#根据序号对数据框进行分类

g <-lapply(f,sum)#对每一类（类=时间）进行分类其中加进去序号的数值了

h <- seq(46,5520,46)#序号数值计算

j <- data.frame(g)

k <- t(j)#对算出来的结果进行转置

l <- k-h#减去序号的数值得到每一时间序列的和
m <- l/920#除以（598-553+1）*（264-245+1）#得到均值。即每一时间维度得均值。

下节为多个在同一目录下的nc文件数据提取（有时下载nc文件特别多尤其是天尺度的，需要进行目录顺序读取）