我有ncdf文件,可以下载here。该数据集包含两个变量(TA和DMI)。 TA值表示自1800-01-01 00:00:00以来单位MONTHS的日期。如何将TA值转换为日期(as.Date)?
这是我的代码:
x<-open.ncdf('dmi.cdf')x$
dmi<-get.var.ncdf(x, varid = 'DMI')
ta<-get.var.ncdf(x, varid = 'TA')
date<-as.Date(ta, origin=c('1800-01-01'))
这里我从我的代码中给出了我的数据集的第一行的第12行的结果。我应该按月而不是白天。
> ta[1:12]
[1] 853.0192 854.0198 855.0205 856.0211
[5] 857.0218 858.0225 859.0231 860.0238
[9] 861.0244 862.0251 863.0257 864.0264
> as.Date(ta[1:12], origin=c('1800-01-01'))
[1] "1802-05-04" "1802-05-05" "1802-05-06" "1802-05-07"
[5] "1802-05-08" "1802-05-09" "1802-05-10" "1802-05-11"
[9] "1802-05-12" "1802-05-13" "1802-05-14" "1802-05-15"
我预期的结果(即年份,我不确定它是正确的1802)
[1] "1802-01-01" "1802-02-01" "1802-03-01" "1802-04-01"
[5] "1802-05-01" "1802-06-01" "1802-07-01" "1802-08-01"
[9] "1802-09-01" "1802-10-01" "1802-11-01" "1802-12-01"