CMIP6多.nc文件在R语言中的合并方法求助

一、优先推荐terra包（简洁稳定，自动处理地理维度） #

terra对NetCDF格式的空间数据支持更友好，能自动对齐lon/lat和时间维度，避免维度不匹配与对象缺失问题：

1. 批量读取所有.nc文件

   # 安装并加载包
   install.packages("terra")
   library(terra)
   
   # 获取当前目录下所有.nc文件路径
   nc_files <- list.files(pattern = "\\.nc$", full.names = TRUE)
   
   # 批量读取为SpatRaster对象
   rasters <- lapply(nc_files, rast)
   

2. 合并并输出为单个.nc文件
   由于你的文件来自同一模型、空间范围一致，直接合并时间维度即可：

   # 合并所有栅格对象（自动对齐维度）
   merged_raster <- c(rasters, recursive = TRUE)
   
   # 输出为新的NetCDF文件
   writeCDF(merged_raster, filename = "merged_monthly_temp.nc", varname = "temp", overwrite = TRUE)
   

二、ncdf4+abind手动处理（适合精细控制场景） #

针对你用abind时遇到的问题，通过显式处理维度来解决：

1. 读取模板文件并提取核心维度

   install.packages("ncdf4")
   library(ncdf4)
   library(abind)
   
   # 读取第一个文件作为维度模板
   nc_template <- nc_open(nc_files[1])
   temp_template <- ncvar_get(nc_template, "temp") # 假设温度变量名为"temp"
   lon <- ncvar_get(nc_template, "lon")
   lat <- ncvar_get(nc_template, "lat")
   time_template <- ncvar_get(nc_template, "time")
   nc_close(nc_template)
   

2. 批量读取并合并数据

   # 初始化合并后的温度数组
   merged_temp <- temp_template
   time_combined <- time_template
   
   # 循环处理剩余文件
   for (file in nc_files[-1]) {
     nc <- nc_open(file)
     temp_current <- ncvar_get(nc, "temp")
     time_current <- ncvar_get(nc, "time")
     nc_close(nc)
   
     # 校验维度匹配（确保lon/lat数量一致）
     if (all(dim(temp_current)[1:2] == dim(merged_temp)[1:2])) {
       merged_temp <- abind(merged_temp, temp_current, along = 3) # 按时间维度（第三维）合并
       time_combined <- c(time_combined, time_current)
     } else {
       stop(paste("文件", file, "空间维度不匹配，请检查"))
     }
   }
   

3. 创建新的NetCDF文件并写入数据

   # 定义维度对象
   lon_dim <- ncdim_def(name = "lon", units = nc_template$dim$lon$units, vals = lon)
   lat_dim <- ncdim_def(name = "lat", units = nc_template$dim$lat$units, vals = lat)
   time_dim <- ncdim_def(name = "time", units = nc_template$dim$time$units, vals = time_combined)
   
   # 定义温度变量
   temp_var <- ncvar_def(name = "temp", units = nc_template$var$temp$units, dim = list(lon_dim, lat_dim, time_dim), missval = nc_template$var$temp$missval)
   
   # 写入并保存文件
   nc_out <- nc_create("merged_monthly_temp.nc", vars = list(temp_var))
   ncvar_put(nc_out, temp_var, merged_temp)
   nc_close(nc_out)
   

常见问题排查 #

• 维度错误：用dim()检查每个文件的温度变量维度，确保lon/lat的数量完全一致，仅时间维度为追加项
• lon对象缺失：读取文件时需显式提取lon变量，确认所有文件的经度变量名统一（比如都是"lon"而非"longitude"）

内容的提问来源于stack exchange，提问作者Ethan D.