4

我有一个有 6 个维度的 netcdf 文件：

f=open.ncdf("C:\\BR_Ban.nc")

我读了变量date：

A = get.var.ncdf(nc=f,varid="date",verbose=TRUE)

然后我像读取变量Tair一样读取变量date。我将这两个变量写入两个文本文件。谁能帮我把它们都写到一个文本文件或excel中？这段代码只会写一个变量：

write.table(as.double(A),"C:\\folder\\shwon_Br_Ban_flux net.txt")

这个变量的问题date是我得到了一个看起来像这样的文本文件：

"x"
"1"  2004
"2"   1
"3"  0.5
"4"  2004
"5"    1
"6"    1

取决于

"157675" 2006
"157676" 365
"157677" 23.5
"157678" 2007
"157679" 366
"157680" 0

如您所见，全部在“年、时和日”列中。有没有办法在一个列中写“半小时一次”的时间步长，在另一列中写年份，并且在一天中写相同的？

4

1 回答 1

4 
f <- open.ncdf("C:\\BR_Ban.nc")
A <- get.var.ncdf(nc=f,varid="date")
B <- get.var.ncdf(nc=f,varid="Tair")
write.table(t(rbind(A,B)),file="output.txt")

正如 netCDF 文件的摘要中提到的那样，date有 2 个维度，维度 1（行）datedim和维度 2（列）time_counter，同时Tair也有 2 个维度，land维度 1 和time_counter维度 2。所以要将它们一起输出，您需要首先给rbind他们，因为他们的共享维度是列 ( time_counter)，然后转置。

于 2012-10-12T08:40:57.660 回答