我在使用我的数据填充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, 创