R-如何使用36个文件中的数据填充数组

时间:2015-07-19 19:26:55

标签: arrays r multidimensional-array field fill

我在使用我的数据填充3D数组时遇到问题。

我想做什么: 我有36个带有变量的“.nc”文件,其中我将值> x计算为24个时间步长。那36个.nc应该显示为矩阵(bycol或row)x 24小时。 因此我的阵列应该暗淡(6,6,24)。我的问题是到目前为止我的代码只为每个6x6填充一个值。但至少我有24种不同的价值观。

我尝试了不同的方法,但这是我最新的代码:

Data<- array(NA,dim=c(6,6,24))



for (j in 0:36) {
nc=open.ncdf(paste("/xxxx/Data#",j,".nc",sep=""))

Var=get.var.ncdf (nc,"Var")
close(nc)

for (k in 1:6) { 
for( l in 1:6)      {
for(i in 1:24)          {

Data[k,l,i]<- (sum(Var[,,i]>x)) 
            }
        }       
}       

我也尝试过:

Data[,,i]<- (sum(Var[,,i]>x)) 

没有k,l循环,但这也不起作用。

提前谢谢!

BRGDS, 创

0 个答案:

没有答案