如何获取API返回的NetCDF文件？R语言调用相关问题求助

问题分析 #

API返回的是二进制NetCDF文件，直接用httr::content()默认解析会得到raw对象，无法直接使用；download.file()报错是因为默认以文本模式保存二进制文件，导致文件损坏，和API的filename参数无关。

方法一：用httr直接保存文件（推荐） #

使用httr::write_disk()在请求时直接将二进制响应写入文件，避免中间处理出错：

library(httr)

# API请求链接
api_url <- "https://dataset.api.hub.geosphere.at/v1/grid/historical/inca-v1-1h-1km?parameters=RR&parameters=T2M&parameters=UU&parameters=VV&parameters=RH2M&start=2014-01-01&end=2014-02-01&bbox=47.25392576%2C11.34458542%2C47.278834%2C11.44352138&output_format=netcdf&filename=file_try"

# 发送请求并直接保存到本地文件
response <- GET(
  url = api_url,
  accept("application/x-netcdf"), # 明确指定接受NetCDF格式，可选但更稳妥
  write_disk("output_file.nc", overwrite = TRUE) # 保存路径，允许覆盖已有文件
)

# 检查请求是否成功
stop_for_status(response)

# 用ncdf4读取文件
library(ncdf4)
nc_data <- nc_open("output_file.nc")

方法二：先获取raw内容再写入文件 #

如果需要先处理响应内容再保存，可指定content()的as = "raw"参数，再用writeBin写入文件：

library(httr)

api_url <- "https://dataset.api.hub.geosphere.at/v1/grid/historical/inca-v1-1h-1km?parameters=RR&parameters=T2M&parameters=UU&parameters=VV&parameters=RH2M&start=2014-01-01&end=2014-02-01&bbox=47.25392576%2C11.34458542%2C47.278834%2C11.44352138&output_format=netcdf&filename=file_try"

response <- GET(api_url, accept("application/x-netcdf"))
stop_for_status(response)

# 获取二进制raw内容
raw_content <- content(response, as = "raw")

# 写入NetCDF文件
writeBin(raw_content, "output_file.nc")

# 读取验证
library(ncdf4)
nc_data <- nc_open("output_file.nc")

方法三：修复download.file()的使用方式 #

核心是设置mode = "wb"以二进制模式保存文件，否则二进制文件会被损坏：

api_url <- "https://dataset.api.hub.geosphere.at/v1/grid/historical/inca-v1-1h-1km?parameters=RR&parameters=T2M&parameters=UU&parameters=VV&parameters=RH2M&start=2014-01-01&end=2014-02-01&bbox=47.25392576%2C11.34458542%2C47.278834%2C11.44352138&output_format=netcdf&filename=file_try"

# 二进制模式保存
download.file(api_url, destfile = "output_file.nc", mode = "wb")

# 读取验证
library(ncdf4)
nc_data <- nc_open("output_file.nc")

关键说明 #

• API中的filename参数只是服务器返回的文件名建议，不会影响本地保存的文件名，无需担心冲突。
• 必须以二进制模式处理NetCDF文件，文本模式会破坏文件结构，导致ncdf4无法识别。
• 使用stop_for_status()可以在请求失败时直接抛出错误，便于调试。

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