我有一个有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
正如您所看到的,所有列都在“年,小时和日”中。有没有办法在一栏中写下“半小时”的时间步骤,另一栏中的年份和当天的相同?
答案 0 :(得分: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
,而{{ 1}}还有2个维度,其中Tair
为维度1,land
为维度2.因此,要将它们一起输出,首先需要time_counter
它们,因为它们的共享维度是列({ {1}}),然后进行转置。