RStudio 4.0.3如何读取GLODAP压缩包、提取变量实现环境变量空间绘图

依赖环境准备 #

• 需提前安装2个核心功能包，直接运行以下代码即可完成安装，无需额外跳转第三方资源站：

install.packages(c("terra", "ncdf4"))

说明：terra包负责空间栅格读写与制图，ncdf4包支持NetCDF格式文件的元数据解析，两个包均为CRAN官方维护版本，适配所有R版本环境。

核心操作步骤 #

1. 定义资源路径与目标变量清单 #

直接复制运行，无需修改基础配置：

# 数据集固定直连地址
data_url <- "https://www.nodc.noaa.gov/archive/arc0107/0162565/2.2/data/0-data/mapped/GLODAPv2_Mapped_Climatology.tar.gz"
# 需提取的4个目标变量
target_vars <- c(
  "GLODAPv2.OmegaAinsitu",
  "GLODAPv2.OmegaCinsitu",
  "GLODAPv2.talk",
  "GLODAPv2.tco2"
)

2. 流式读取压缩包，不落地全量大文件 #

• 操作逻辑：通过R的临时目录机制建立网络流式通道，仅从压缩包中提取核心NetCDF数据文件，不将3GB以上的完整tar.gz包存储到本地磁盘，解决大文件下载、磁盘占满导致的读取失败问题。

# 创建R系统临时工作目录，操作完成后自动清理
temp_dir <- tempdir()
# 建立网络流式连接，设置300秒超时适配大文件读取
con <- url(data_url, open = "rb", timeout = 300)
# 从压缩包中定向提取核心NetCDF文件，跳过说明文档、冗余格式文件
untar(
  tarfile = con,
  exdir = temp_dir,
  files = "GLODAPv2_Mapped_Climatology.nc"
)
close(con)
# 定位本地临时路径下的NetCDF文件
nc_path <- file.path(temp_dir, "GLODAPv2_Mapped_Climatology.nc")

注意：如果运行后提示文件不存在，可临时去掉untar()函数中的files参数，先列出压缩包内全部文件路径，匹配到真实的NetCDF文件名后再调整提取规则即可。

3. 定向加载目标变量，避免全量内存占用 #

• 此前读取无有效返回数据的核心原因：全量加载整个NetCDF文件（包含所有深度层、全部数十个变量）导致内存溢出，触发R的静默截断。该步骤仅定向读取指定的4个目标变量，不加载冗余数据。

# 仅加载目标变量为空间栅格对象
data_rast <- rast(nc_path, subds = target_vars)
# 验证加载结果：打印后可查看每个变量的数值范围、空间分辨率、投影信息
print(data_rast)

如果打印结果显示变量层存在有效数值区间（非全NA），则说明加载成功。

4. 空间可视化制图 #

直接运行以下代码即可生成2行2列的多变量全球空间分布图，自动匹配WGS84投影的边界叠加：

# 设置绘图面板布局
par(mfrow = c(2, 2), mar = c(4, 4, 2, 1))
# 逐变量绘制空间图
for (var_idx in 1:nlyr(data_rast)) {
  # 绘制变量栅格
  plot(data_rast[[var_idx]], main = names(data_rast)[var_idx])
  # 叠加全球海岸线边界，提升可读性
  world_map <- terra::vect("world")
  terra::plot(world_map, add = TRUE, lwd = 0.8)
}

如果运行时提示无法加载世界边界矢量，可直接删除world_map和叠加绘图的两行代码，不影响核心栅格的绘制结果。

常见问题排查 #

• 网络连接中断：可适当延长url()函数中的timeout参数值，最高设置为600，适配低带宽读取场景。
• 加载后变量全为NA：运行nc <- nc_open(nc_path); print(nc$var)列出NetCDF内部的真实变量名，替换原有target_vars向量即可，多数情况是压缩包内变量的层级路径存在细微差异。
• 内存不足：读取栅格时可通过lyrs参数指定深度层范围，例如rast(nc_path, subds = target_vars, lyrs = 1:15)，仅加载上层1000米深度的数据，降低内存占用。

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