如何在R中构建循环下载REASv3.2 CO2数据并生成raster对象

前置依赖 #

首先安装并加载所需R包，优先使用terra替代传统raster包，运算效率更高，也可直接导出兼容raster格式的结果：

# 安装包（首次运行执行）
install.packages(c("terra", "ncdf4"))
# 加载包
library(terra)
library(ncdf4)

操作步骤 #

• 提前将1950-2015年的REASv3.2 CO2数据netCDF（.nc）文件下载到本地同一个文件夹中，确保文件命名包含对应年份标识（官方默认命名规则通常已满足该要求）
• 按以下逻辑构建循环，批量提取指定变量并输出栅格结果

完整代码示例 #

# 1. 配置基础参数
# 本地nc文件存储文件夹路径，替换为你自己的路径
nc_dir <- "C:/your/path/to/REAS_data/"
# 目标变量名，和文件内变量名完全一致
target_var <- "REASv3.2_CO2-BF_TOTAL"
# 目标年份范围
year_range <- 1950:2015
# 输出栅格文件存储路径
output_dir <- "C:/your/path/to/output/"
dir.create(output_dir, recursive = TRUE, showWarnings = FALSE)

# 2. 初始化存储多波段栅格的空对象
co2_stack <- rast()

# 3. 循环处理每年数据
for (y in year_range) {
  # 匹配对应年份的nc文件，适配官方默认命名规则，如果你的文件名不同需要调整匹配规则
  nc_file <- list.files(nc_dir, pattern = paste0("REASv3.2.*", y, ".*\\.nc$"), full.names = TRUE)
  if (length(nc_file) == 0) {
    message(y, "年对应数据文件不存在，跳过")
    next
  }
  # 读取nc文件中的指定变量
  r <- rast(nc_file, subds = target_var)
  # 给单波段栅格命名为对应年份
  names(r) <- paste0("CO2_", y)
  # 加入多波段栈
  co2_stack <- c(co2_stack, r)
  
  # 如果需要单独输出每年的单波段raster文件，取消注释下方代码
  # writeRaster(r, filename = paste0(output_dir, "REAS_CO2_", y, ".tif"), overwrite = TRUE)
}

# 4. 输出整体多波段raster文件（tif格式兼容所有GIS和R的raster包读取）
writeRaster(co2_stack, filename = paste0(output_dir, "REAS_CO2_1950_2015.tif"), overwrite = TRUE)

# 如果需要转为传统raster包的RasterStack对象，执行下方代码即可
# co2_raster_stack <- raster::stack(co2_stack)

注意事项 #

• 读取文件报错时首先检查变量名是否和nc文件内的变量完全一致，可执行ncdf4::nc_open(nc_file)查看文件内所有变量名
• 自动批量下载可配合download.file函数，根据官方文件下载链接的规律构造每年的下载地址后放入循环内执行即可

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