我有一个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"
Run Code Online (Sandbox Code Playgroud)
时间:
A <- get.var.ncdf(nc=f1,varid="time")
head(A)
## [1] 990552 990558 990564 990570 990576 990582
Run Code Online (Sandbox Code Playgroud)
另一个变量:
B1 <- get.var.ncdf(nc=f,varid="stm")
write.table(t(rbind(A,B1)),file="output.txt")
Run Code Online (Sandbox Code Playgroud)
在写入文本输出文件之前,我需要您的帮助将时间变量转换为正常日期.
日期
如果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"
Run Code Online (Sandbox Code Playgroud)
日期+时间
如果你也想要小时,请改用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"
Run Code Online (Sandbox Code Playgroud)
时区
要了解您的时区,请键入
Sys.timezone()
Run Code Online (Sandbox Code Playgroud)
有关可能时区的完整列表,请键入
OlsonNames()
Run Code Online (Sandbox Code Playgroud)
添加/删除小时数
无论哪种方式,您都可以通过乘以小时数3600
和添加/删除结果来删除/添加小时数
as.POSIXct(A*3600, origin = "1900-01-01") - 2*3600
Run Code Online (Sandbox Code Playgroud)