注意

この記事は、2年ほど前にQiitaに投稿した記事の再投稿になります。
若干内容が古くなっている可能性があります。

インストールをして、簡単な動作確認をしてみます。
踏み込んだ内容はありませんので、あくまで初心者の方用の記事です。

wgrib2 インストール

インストールはこちら。

いっぱいファイルがあって何を落とせばいいかわからない方は、最初にwgrib2.exeだけを落としてください。
wgrib2.exeだけで実行すれば、実行に必要な(不足している)ファイルがダイアログで出力されるので、それらを同じディレクトリに入れて、パスを通せばOKです。

wgrib2 動かしてみる

今回は、NOMADSから風速データを落としてきます。
GDAS 0.25(grib filter) => gdas.2018xx と進みます。

今回は取り合えず**地上10m地点のU方向, V方向の風**を取得します。
ですので、チェックボックスは以下をチェック。
- 10m above ground
- UGRD
- VGRD

そして、落としたファイルをwgrib2で見てみます。
Shellで下記を実行。


$ wgrib2 [ファイル名]
1:0:d=2018041200:UGRD:10 m above ground:anl:
2:778510:d=2018041200:VGRD:10 m above ground:anl:

出力の内容はつまり、
1：地上10mでU方向のデータ
2：地上10mでV方向のデータ
ということです。うまく取れてるようです。

ちなみに、Optionで*-V*を指定すると、


$ wgrib2 -V [ファイル名]
1:0:vt=2018041200:10 m above ground:anl:UGRD U-Component of Wind [m/s]:
ndata=1038240:undef=0:mean=-0.0738295:min=-19.6813:max=28.6987
grid_template=0:winds(N/S):
  lat-lon grid:(1440 x 721) units 1e-06 input WE:NS output WE:SN res 48
  lat 90.000000 to -90.000000 by 0.250000
  lon 0.000000 to 359.750000 by 0.250000 #points=1038240

2:778510:vt=2018041200:10 m above ground:anl:VGRD V-Component of Wind [m/s]:
ndata=1038240:undef=0:mean=0.173403:min=-26.8313:max=23.2687
grid_template=0:winds(N/S):
  lat-lon grid:(1440 x 721) units 1e-06 input WE:NS output WE:SN res 48
  lat 90.000000 to -90.000000 by 0.250000
  lon 0.000000 to 359.750000 by 0.250000 #points=1038240

と、データの詳細な情報が見れます。

更に、-dオプションをつけて、抜き出すデータの番号を指定すると、


$ wgrib2 -d 1 [ファイル名]
1:0:d=2018041200:UGRD:10 m above ground:anl:

このように、指定された番号に対応するデータのみ取り出せます。
なので、例えばU方向のデータだけ取り出して別のファイルに保存したければ、


$ wgrib2 -d 1 [ファイル名] -order we:ns -ieee u.dat(好きな名前でいいです)

とすればOKです。

grib2json インストール

grib2json インストール

このリポジトリをローカルにもってきます。


$ git clone https://github.com/weacast/weacast-grib2json.git

として、srcにあるファイルだけどこか別のディレクトリに保存すればOKです。

ちなみに、そのままだと
Error: Unable to access jarfile だとか、何かしらエラーが出てしまいます。
ですので、私は grib2json.cmd を


  @echo off
  "java" -Xmx512M -jar "%~dp0grib2json.jar" %*

と書き換えたらうまいこと動きました。

grib2json 動かしてみる

先ほど落とした風のファイルを、json形式で出力してみます。


$ grib2json [ファイル名]
[
  {
  "header":{
    "discipline":0,
    "gribEdition":2,
    "gribLength":778510,
    "center":7,
    "subcenter":0,
    "refTime":"2018-04-12T00:00:00.000Z",
    "significanceOfRT":1,
    "productStatus":0,
    "productType":1,
    "productDefinitionTemplate":0,
    "parameterCategory":2,
    "parameterNumber":2,
    "parameterUnit":"m.s-1",
    "genProcessType":2,
    "forecastTime":0,
    "surface1Type":103,
    "surface1Value":10.0,
    "surface2Type":255,
    "surface2Value":0.0,
    "gridDefinitionTemplate":0,
    "numberPoints":1038240,
    "shape":6,
    "gridUnits":"degrees",
    "resolution":48,
    "winds":"true",
    "scanMode":0,
    "nx":1440,
    "ny":721,
    "basicAngle":0,
    "lo1":0.0,
    "la1":90.0,
    "lo2":359.75003,
    "la2":-90.0,
    "dx":0.25,
    "dy":0.25
  }
},
{
  "header":{
    "discipline":0,
    "gribEdition":2,
    "gribLength":748505,
    "center":7,
    "subcenter":0,
    "refTime":"2018-04-12T00:00:00.000Z",
    "significanceOfRT":1,
    "productStatus":0,
    "productType":1,
    "productDefinitionTemplate":0,
    "parameterCategory":2,
    "parameterNumber":3,
    "parameterUnit":"m.s-1",
    "genProcessType":2,
    "forecastTime":0,
    "surface1Type":103,
    "surface1Value":10.0,
    "surface2Type":255,
    "surface2Value":0.0,
    "gridDefinitionTemplate":0,
    "numberPoints":1038240,
    "shape":6,
    "gridUnits":"degrees",
    "resolution":48,
    "winds":"true",
    "scanMode":0,
    "nx":1440,
    "ny":721,
    "basicAngle":0,
    "lo1":0.0,
    "la1":90.0,
    "lo2":359.75003,
    "la2":-90.0,
    "dx":0.25,
    "dy":0.25
    }
  }
]

1, 2それぞれが別の連想配列になっています。
ではさらに、風速データのみをJSON形式で保存してみます。


$ grib2json --names --data -o file.json(好きなファイル名) [ファイル名]

すると、-oオプションの次に与えた任意のファイル名で、JSON形式のファイルが保存されます。

こうすると、何かと扱い易くて便利です。
(ちなみに、leaflet.jsのvelocityプラグインがなかなか楽しいので、是非grib2jsonを使って実行してみてください)

最後に

つい最近使うことになった上記二つのツールですが、便利なのに案外初心者向けの記事がなくて困ってしまいました。
ですので、この記事が少しでも次のプログラマーを目指す方の一助になればいいな、と思っています。

また、私自身ペーペーなので、訂正すべき個所があればコメントをお願いいたします。