いくつかの気候データを grib 形式でダウンロードしました。「R」を使用してNetCDF形式に変換したい。
さらに、ファイルはさまざまな変数で構成されているため、一度に 1 つの変数を個別のファイルに抽出したいと考えています。
いくつかの気候データを grib 形式でダウンロードしました。「R」を使用してNetCDF形式に変換したい。
さらに、ファイルはさまざまな変数で構成されているため、一度に 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)
1.RNOMADS
パッケージ Rnomads には関数 readgrib があり、外部ライブラリにラッパーを提供して、grib ファイルを読み取ることができます
2. netcdf への変換
GRIB データが通常の緯度経度グリッド上にある場合、おそらくより簡単な方法は、より開発された読み取りのサポートとして netcdf に変換することです (そして、おそらく既に使用に慣れているでしょう)。
いくつかの方法でグリブを変換できます。最も簡単な 2 つは次のとおりです。
CDO:
cdo -f nc copy test.grb test.nc
netcdf4 規則が必要な場合は、「-f nc4」を使用します。
ECCODES (brew install eccodes を使用した Mac インストールの場合)
grib_to_netcdf -o test.nc test.grb