如何使用 R 将 Grib1 转换为 Netcdf？

Question

我刚刚下载了一些 grib 格式的气候数据。 我想使用“R”将其转换为 NetCDF 格式。

此外，由于文件由不同的变量组成，我想一次将一个变量提取到单个文件中。

Answer 1

如果没有您的特定文件，很难回答这个问题。 您应该考虑制作可重现的示例，尤其是当您要发布到 R 板时。

对于 R，请查看library(raster)和library(ncdf4) 。 我只是抓取了我看到的第一个 grib1 文件，并整理了一个快速示例。

library(raster)
library(ncdf4)

download.file(url = 'ftp://ftp.hpc.ncep.noaa.gov/grib/20130815/p06m_2013081500f030.grb', destfile = 'test.grb')

(r <- raster('test.grb'))

n <- writeRaster(r, filename = 'netcdf_in_youR_comp.nc', overwrite = TRUE)

Answer 2

1. RNOMADS

Rnomads 包有一个函数 readgrib 为外部库提供包装器，允许读取 grib 文件

2.转换为netcdf

如果 GRIB 数据位于常规经纬度网格上，那么可能更简单的方法是转换为 netcdf 作为更发达的阅读支持（您可能已经习惯使用它）

您可以通过多种方式转换 grib，最简单的两种方式是

首席数据官：

cdo -f nc copy test.grb test.nc 

如果需要 netcdf4 约定，请使用“-f nc4”。

ECCODES（在 mac 上使用 brew install eccodes 安装）

grib_to_netcdf -o test.nc test.grb

Answer 3

您可以使用计算机上安装的 ncl

library(ncdf)
 system(ncl_convert2nc  xxxx.grb, internal = TRUE)
 my.nc <- open.ncdf("result.nc")
 print(my.nc)