在R中通过OPeNDAP获取THREDDS的ncdf文件后如何正确写入本地？

你的代码错误点 #

• 仅创建了空的NetCDF变量框架，但没有读取OPeNDAP链接中的实际数据并写入本地文件，生成的file.nc只是一个无数据的空壳
• 未复制原文件的维度定义、全局属性、变量属性，缺少这些核心元数据的NetCDF文件无法被正常解析

R中可行的实现方式 #

方法1：用ncdf4包完整读写 #

需要手动处理维度、数据、属性的复制，步骤如下：

library(ncdf4)

# 1. 打开OPeNDAP链接
url <- "https://thredds.met.no/thredds/dodsC/metpparchivev3/2020/01/01/met_analysis_1_0km_nordic_20200101T00Z.nc?x[663:1:674],y[732:1:750],latitude[732:1:750][663:1:674],longitude[732:1:750][663:1:674],altitude[732:1:750][663:1:674],air_temperature_2m[0:1:0][732:1:750][663:1:674]"
ncin <- nc_open(url)

# 2. 提取原文件的维度信息
dims <- ncin$dim

# 3. 创建本地NetCDF文件，包含维度和变量结构
filepath <- "file.nc"
ncout <- nc_create(filepath, vars = ncin$var, dims = dims)

# 4. 复制全局属性
for (att_name in names(ncin$atts)) {
  ncatt_put(ncout, 0, att_name, ncin$atts[[att_name]])
}

# 5. 逐个读取变量数据并写入新文件，同时复制变量属性
for (var_name in names(ncin$var)) {
  # 读取变量数据
  var_data <- ncvar_get(ncin, var_name)
  ncvar_put(ncout, var_name, var_data)
  # 复制变量属性
  for (att_name in names(ncin$var[[var_name]]$atts)) {
    ncatt_put(ncout, var_name, att_name, ncin$var[[var_name]]$atts[[att_name]])
  }
}

# 6. 关闭文件连接
nc_close(ncin)
nc_close(ncout)

方法2：用terra包简化操作（更推荐） #

terra原生支持OPeNDAP协议，处理流程更简洁，还会自动保留元数据：

library(terra)

# 直接读取OPeNDAP子集为SpatRaster对象
r <- rast(url)

# 写入本地NetCDF文件
writeRaster(r, "file_terra.nc", format = "CDF", overwrite = TRUE)

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