Question

我有一个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")

我需要你的帮助才能在写入文本输出文件之前将时间变量转换为正常日期。

Answer 1

<强>日期

如果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

如何将数字转换为正常日期？

1 个答案: