【R】worldmet

1. はじめに

worldmetは、NOAA Integrated Surface Databaseにアクセスして容易にデータを取得できるパッケージです。世界35,000 地点の情報が載っています。所在地は map　で確認できます。

2. インストール

CRANからインストールできます。

install.packages("worldmet")

3. つかってみる

成田空港のメタデータを取得します。

library(worldmet)

getMeta(site = "narita")

> getMeta(site = "narita")
# A tibble: 1 x 13
  usaf   wban  station     ctry  st    call  latitude longitude `elev(m)` begin     
  <chr>  <chr> <chr>       <chr> <chr> <chr>    <dbl>     <dbl>     <dbl> <date>    
1 476860 99999 NARITA INTL JA    NA    RJAA      35.8      140.        43 1978-05-19
# ... with 3 more variables: end <date>, code <chr>, dist <lgl>

通常は、緯度経度で指定します。ドイツ・ブラウンシュバイクの情報を取得してみます。

info <- getMeta(lat = 52.2637014559553, lon = 10.527587681798869)
info

> info
# A tibble: 10 x 15
   usaf   wban  station     ctry  st    call  latitude longitude `elev(m)` begin     
   <chr>  <chr> <chr>       <chr> <chr> <chr>    <dbl>     <dbl>     <dbl> <date>    
 1 103480 99999 BRAUNSCHWE~ GM    NA    EDVE      52.3     10.6       89.9 1931-01-02
 2 103430 99999 CELLE       GM    NA    ETHC      52.6     10.0       39.3 1936-01-02
 3 103380 99999 HANNOVER    GM    NA    EDDV      52.5      9.68      55.8 1926-01-01
 4 102460 99999 FASSBERG    GM    NA    ETHS      52.9     10.2       74.7 1936-01-02
 5 103340 99999 WUNSTORF    GM    NA    ETNW      52.5      9.43      57   1937-05-07
 6 103350 99999 BUCKEBURG   GM    NA    ETHB      52.3      9.08      70.1 1973-01-02
 7 104360 99999 KASSEL CAL~ GM    NA    EDVK      51.4      9.38     277.  1988-12-23
 8 105540 99999 ERFURT      GM    NA    EDDE      51.0     11.0      316.  1931-01-01
 9 102240 99999 BREMEN      GM    NA    EDDW      53.0      8.79       4.3 1929-08-01
10 101490 99999 HAMBURG FI~ GM    NA    EDHI      53.5      9.84       6.7 1988-10-28
# ... with 5 more variables: end <date>, code <chr>, longr <dbl>, latr <dbl>,
#   dist <dbl>

地図とともに表示されるのがいいですね。

風向情報を取得してみます。場所は、コードで指定します。

met_braunschweig <- importNOAA(code = "103480-99999", year = 2020)
head(met_braunschweig)

> head(met_braunschweig)
# A tibble: 6 x 21
  code    station   date                latitude longitude  elev    ws    wd air_temp
  <fct>   <fct>     <dttm>                 <dbl>     <dbl> <dbl> <dbl> <dbl>    <dbl>
1 103480~ BRAUNSCH~ 2020-01-02 05:00:00     52.3      10.6  89.9  1.5    100     -4  
2 103480~ BRAUNSCH~ 2020-01-02 06:00:00     52.3      10.6  89.9  1.05   110     -4  
3 103480~ BRAUNSCH~ 2020-01-02 07:00:00     52.3      10.6  89.9  0.75   230     -3  
4 103480~ BRAUNSCH~ 2020-01-02 08:00:00     52.3      10.6  89.9  1.05   100     -2.5
5 103480~ BRAUNSCH~ 2020-01-02 09:00:00     52.3      10.6  89.9  1.05   230      1  
6 103480~ BRAUNSCH~ 2020-01-02 10:00:00     52.3      10.6  89.9  3.35   260      3.5
# ... with 12 more variables: atmos_pres <dbl>, visibility <dbl>, dew_point <dbl>,
#   RH <dbl>, ceil_hgt <dbl>, cl_1 <dbl>, cl_2 <dbl>, cl_3 <dbl>, cl <dbl>,
#   cl_1_height <dbl>, cl_2_height <dbl>, cl_3_height <dbl>

openaitパッケージを使うと容易に表現できます。

library(openair)
windRose(met_braunschweig)

4. さいごに

世界の気象データがこんなに簡単に取得できるなんて驚きです。