如何使用R将6小时分辨率NetCDF数据集转换为日尺度NetCDF数据集？

以下是基于terra包（R中处理栅格/NetCDF数据的高效工具）的完整实现方案，适配ERA5 6小时雪深数据的日尺度聚合需求：

步骤1：安装并加载依赖包 #

# 首次运行时执行安装
install.packages(c("terra", "ncdf4"))

# 加载工具包
library(terra)
library(ncdf4)

步骤2：读取原始6小时分辨率NetCDF数据 #

# 替换为你的数据集路径
nc_file <- "1980-2022_era5_snow_depth_6h.nc"
# 读取为栅格时间序列对象
snow_rast <- rast(nc_file)

# 验证时间维度，确认是6小时间隔
print(time(snow_rast))

步骤3：按日聚合数据 #

ERA5的6小时数据时间戳通常为00:00、06:00、12:00、18:00 UTC，雪深的日尺度统计一般取日最大值（可根据需求替换为均值/最小值）：

# 提取日期部分（忽略小时），保持UTC时区避免偏移
daily_groups <- as.Date(time(snow_rast), tz = "UTC")

# 按日期分组聚合，自动忽略缺失值
daily_snow <- tapp(snow_rast, daily_groups, fun = function(x) max(x, na.rm = TRUE))

# 查看聚合后的时间维度（应为每日一个时间戳）
print(time(daily_snow))

步骤4：保存为日尺度NetCDF文件 #

# 替换为输出文件路径
output_file <- "1980-2022_era5_snow_depth_daily.nc"

# 保存时保留原数据的属性信息（需匹配你的变量名/单位）
writeCDF(daily_snow, 
         filename = output_file,
         varname = "sd",  # ERA5雪深默认变量名为sd
         unit = "m",      # ERA5雪深单位为米
         longname = "Daily Maximum Snow Depth (Aggregated from ERA5 6-hour data)",
         overwrite = TRUE)

关键说明 #

• tapp()是terra专为时间序列聚合设计的函数，效率远高于手动循环分组
• 若需其他聚合逻辑，直接替换fun参数即可（如mean取日均值）
• 必须指定tz = "UTC"处理时间，避免时区转换导致的日期错位
• 若原数据包含多变量，可先用subset(snow_rast, "sd")筛选出雪深变量再处理

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