使用R将多年月度温度栅格转为NetCDF时保留变量名

刚好做过类似的需求，要把这种按月份命名波段的栅格转成NetCDF还不丢变量名，其实关键就是在导出的时候指定好变量名参数就行，下面是完整的可运行代码：

1. 加载所需包 #

首先得把处理栅格和NetCDF的包装上并加载：

# 安装包（如果还没装的话）
install.packages(c("raster", "ncdf4"))
# 加载包
library(raster)
library(ncdf4)

2. 读取并合并多年栅格 #

读取你的年度栅格文件，然后合并成一个包含所有月份的栅格对象，这里原波段的变量名会自动保留：

# 读取单年度温度栅格
r1 <- brick('TavgM_1981.grd')
r2 <- brick('TavgM_1982.grd')

# 合并所有年度栅格成一个对象
merged_temp <- stack(r1, r2)

# 可以验证一下变量名是否保留
print(names(merged_temp))

3. 导出为NetCDF并保留变量名 #

这一步是核心！用writeRaster导出时，通过varname参数指定为栅格的原有变量名，这样NetCDF里每个变量的名字就会和原栅格波段名完全一致：

# 导出为NetCDF文件，保留原变量名
writeRaster(
  merged_temp,
  filename = "Monthly_Temperature_1981-1982.nc",
  format = "CDF",  # 指定格式为NetCDF
  varname = names(merged_temp),  # 关键：用原栅格的变量名作为NetCDF的变量名
  overwrite = TRUE,  # 如果文件已存在则覆盖
  NAflag = -9999  # 设置缺失值标记，按需调整
)

小提示 #

如果你的年度栅格文件很多，可以用list.files()批量读取，不用一个个手动定义对象：

# 批量读取所有.grd格式的温度文件
temp_files <- list.files(pattern = "TavgM_\\d{4}\\.grd$", full.names = TRUE)
merged_temp <- stack(temp_files)

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