iT邦幫忙

2022 iThome 鐵人賽

DAY 20
0
AI & Data

氣象食材系列 第 20

[Day 20]進階預報資料-數值預報模式-區域預報模式(WRF-3公里)(grib2)

  • 分享至 

  • xImage
  •  

GRIB2資料,詳細介紹可以看這

可以解GRIB2檔案格式的套件列幾個如下,
pygrib, grib2io, gdal(沒錯,又是gdal), xarray, pynio

如果還記得昨天的內容,應該可以猜得到上述套件編譯過程中皆有連結grib2的函式庫。

這邊介紹之前常用的pygrib,安裝放式如下
conda install -c conda-forge pygrib

pip3 isntall pygrib

安裝完成後就開始解資料吧!

到氣象局開放資料官網下載xml,並打開xml,可以看到實際下載位置之連結(source標記)。
就下載吧,一份檔案大約接近200MB。
https://ithelp.ithome.com.tw/upload/images/20221004/20150923ODbnn0qRCu.jpg
https://ithelp.ithome.com.tw/upload/images/20221004/20150923djOQttbVJO.jpg

因為之前已經有下載過颱風預報的資料,故這邊就用這一份資料呈現

import pygrib as pb
grbdata = pb.open("day20/M-A0064-21091006-048.grb2")
print(grbdata.select()) #可以看有哪些變數,然後找到對應的index。
print(grbdata.select()[1].keys()) #可以看到選擇變數的key值

上述程式中的print(grbdata.select()),可以得到下圖(部分),黃色螢光筆的部分即為資料對應的index。
https://ithelp.ithome.com.tw/upload/images/20221004/20150923Cxb9iGM8cv.jpg
print(grbdata.select()[1].keys()),可以得到下圖(部分),有非常多的key可以選。
https://ithelp.ithome.com.tw/upload/images/20221004/20150923jhdsNFyNPG.jpg

經過這些key的選取之後,就可以得到不同的資訊值

print(grbdata.select()[1]["values"])

這樣就可以取得index為1的變數數值。至於這些數值式格點上的值,要視覺化的話,還需要格點的經緯度。
這些就待由明天介紹吧。


上一篇
[ Day 19] 簡介相同資料格式可用不同套件讀取
下一篇
[Day 21]進階預報資料-數值預報模式-區域預報模式(WRF-3公里)視覺化(grib2)
系列文
氣象食材30
圖片
  直播研討會
圖片
{{ item.channelVendor }} {{ item.webinarstarted }} |
{{ formatDate(item.duration) }}
直播中

尚未有邦友留言

立即登入留言