如何通过Terra从含valid_time字段的ERA5 NetCDF提取时间信息

由于ERA5的NetCDF文件使用valid_time而非标准time作为时间变量名，terra::time()默认无法识别，返回全NA。你可以通过terra自带的工具直接提取并赋值时间，无需依赖ncdf4包，步骤如下：

1. 提取valid_time变量：使用terra::ncmeta()直接读取NetCDF文件中的valid_time字段值，该函数可以获取NC文件的元数据变量。
2. 转换时间格式：ERA5的valid_time通常是Unix时间戳（单位为秒），需要转换为POSIXct格式的时间对象。
3. 赋值给SpatRaster：通过terra::time<-将转换后的时间绑定到SpatRaster对象上，之后就能正常使用terra::time()获取时间信息。

完整代码示例：

library(terra)

ncfile <- 'filepath'
sample_rast <- terra::rast(ncfile)

# 提取valid_time原始值
valid_time_raw <- terra::ncmeta(ncfile, var = "valid_time")$value

# 转换为可读的时间格式（ERA5默认以1970-01-01为起始的秒级时间戳）
sample_dates <- as.POSIXct(valid_time_raw, origin = "1970-01-01", tz = "UTC")

# 将时间赋值给SpatRaster对象
terra::time(sample_rast) <- sample_dates

# 验证结果：查看SpatRaster信息，此时会显示时间字段
sample_rast

# 直接调用time()获取时间
terra::time(sample_rast)

注意事项 #

• 如果你的valid_time单位不是秒（比如小时），需要先对原始值进行转换，例如：valid_time_raw * 3600再传入as.POSIXct。
• 赋值后，SpatRaster的图层名称不会自动更新，但时间属性会被正确识别，后续时间相关的操作（如按时间筛选图层）都可以正常使用。

内容的提问来源于stack exchange，提问作者Jaken