R/leafletで風を可視化

3413 ワード
R wgrib2 leaflet R テキストリンク

wgrib2をインストール。

wget https://www.ftp.cpc.ncep.noaa.gov/wd51we/wgrib2/wgrib2.tgz.v3.0.2
tar xvfz wgrib2.tgz.v3.0.2
cd grib2
vim makefile # 125行目（gcc）と126行目（gfortran）をコメントアウト
make # wgrib2/wgrib2が実行ファイル

grib2jsonをインストール。

git clone git clone https://github.com/cambecc/grib2json.git
cd grib2json
mvn package
cd target
tar xvfz grib2json-0.8.0-SNAPSHOT.tar.gz # grib2json-0.8.0-SNAPSHOT/bin/grib2jsonが実行ファイル

GFS 1.0 Degreeをダウンロード＆前処理。

wget https://nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/gfs.20210307/00/gfs.t00z.pgrb2.1p00.f000
wgrib2 ./gfs.t00z.pgrb2.1p00.f000
wgrib2 ./gfs.t00z.pgrb2.1p00.f000 -match ":*GRD:10 m above" -grib tmp.grib2
wgrib2 -V tmp.grib2
grib2json --names --data --o tmp.json tmp.grib2

R/leafletで可視化（QiitaでHTMLを貼り付けることはできないんだろうか…）。

library(leaflet)
library(leaflet.extras2)

content <- "tmp.json"
leaflet() %>%
  addTiles() %>%
  addVelocity(content=content) %>%
  setView(lng=135, lat=35, zoom=2)

wgrib2とかgfsとか久々すぎて...