我有一个netcdf
文件可以从中提取一些变量。提供者告诉我,时间变量(A)为以来的小时数1900-01-01 00:00:0.0
。我需要将此变量转换为正常日期并写入另一个txt文件。
f1 <- open.ncdf("C:\\Users\\data.nc")
# [1] "file C:\\Users\\data.nc has 3 dimensions:"
# [1] "longitude Size: 2"
# [1] "latitude Size: 2"
# [1] "time Size: 2920"
# [1] "file C:\\Users\\data.nc has 5 variables:"
# [1] "short stm[longitude,latitude,time] Longname:soil textre Missval:-32767"
时间:
A <- get.var.ncdf(nc=f1,varid="time")
head(A)
## [1] 990552 990558 990564 990570 990576 990582
另一个变量:
B1 <- get.var.ncdf(nc=f,varid="stm")
write.table(t(rbind(A,B1)),file="output.txt")
在写入文本输出文件之前,我需要您的帮助将时间变量转换为正常日期。
日期
如果A
是以来的小时数"1900-01-01"
,则将它们除以24,然后将此日期传递到origin
函数内的参数中as.Date
,如
A <- c(990552, 990558, 990564, 990570, 990576, 990582)
as.Date(A/24, origin = "1900-01-01")
## [1] "2013-01-01" "2013-01-01" "2013-01-01" "2013-01-01" "2013-01-02" "2013-01-02"
日期+时间
如果您还想要小时,请改用as.POSIXct
3600乘以(您可能还应该在tz
参数中指定时区)
as.POSIXct(A*3600, origin = "1900-01-01")
## [1] "2013-01-01 02:00:00 IST" "2013-01-01 08:00:00 IST" "2013-01-01 14:00:00 IST" "2013-01-01 20:00:00 IST" "2013-01-02 02:00:00 IST" "2013-01-02 08:00:00 IST"
时区
为了知道您的时区,请输入
Sys.timezone()
有关可能的时区的完整列表,请键入
OlsonNames()
添加/删除时间
无论哪种方式,您都可以通过将小时数乘以3600
并从结果中添加/删除来删除/添加小时,例如
as.POSIXct(A*3600, origin = "1900-01-01") - 2*3600
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句