一、環境准備
(1)、python3環境
(2)、wgirb工具(用於讀取grib1文件),下載地址: ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib
(3)、wgirb2工具(用於讀取grib2文件),下載地址:ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib2
二、環境測試
工具下好了,windows下無需安裝,怎么看是否能用呢?很簡答,以wgrib為例,下載好了之后在我的文件夾下是這樣的:

可以看到里面有一個應用程序,那么從命令行進入該文件夾下,運行命令:wgrib.exe就可以查看一些使用幫助了,如下圖所示

可以看到很多提示信息,也就是教我們如何使用這個工具。wgrib2也是同樣的查看。
接下來就可以使用python來讀取grib1和grib2數據了(其實就是在python里面調用cmd命令而已)
三、讀取數據
grib1和grib2數據我了解不多,說錯的地方還請大家指正,它們是兩種氣象數據,在我看來,它相當於一個大倉庫,里面是一個個“小房間”,每個房間都有編號,每個“小房間”存放的是一個氣象變量,例如溫度、濕度、壓力等等,那么如何讀取我們需要的氣候變量呢?下面分別以grib1和grib2數據為例,來讀取溫度變量的值。二者讀取方式相似,只是使用的工具不同。
剛剛說了,每一個“小房間”存放的是一個氣象變量且“小房間”有編號,那么我們有2中方法進行讀取,一種是根據“小房間”編號進行讀取,另一種是根據變量名稱進行匹配讀取,一般來說后者使用的較多,畢竟使用名稱是肯定錯不了的。下面進行分別介紹。
1、按編號讀取
既然要根據編號讀取,那么首先要知道編號才行,還好wgrib工具中也提供了這和功能。查看“小房間”編號的代碼如下:
import os
# 將當前目錄轉到wgrib工具目錄
os.chdir("G:\mygrib\wgrib1")
# 執行cmd命令,查看各個“小房間”具體信息
os.system(r"wgrib G:\mete_data\W_NAFP_C_ECMF_20170713190330_P_C3E07131200071312001 -v")
運行得到的結果如下圖所示:
如上圖所示,我們可以得到很多有用的信息,例如:2m高的溫度變量所在的“小房間”編號是82,接下來就可以根據這個編號82來進行讀取2m溫度變量,讀取代碼如下:
import os
# 將當前目錄轉到wgrib工具目錄
os.chdir("G:\mygrib\wgrib1")
# 執行cmd命令,查看各個“小房間”具體信息
# os.system(r"wgrib G:\mete_data\W_NAFP_C_ECMF_20170713190330_P_C3E07131200071312001 -v")
# 知道了編號82,那么就可以根據這個編號把數據讀到一個txt里面,注意:wgirb工具不支持讀到csv文件,而wgrib2工具支持
os.system("wgrib G:\mete_data\W_NAFP_C_ECMF_20170713190330_P_C3E07131200071312001 -d 82 -h -text -o 2m_tmp.txt")
讀取結果如下圖所示:

除此之外,wgrib工具還支持“變量名稱匹配讀取”,在這里就不介紹了,可以看下一部分的利用wgrib2工具進行“變量名稱匹配讀取”二者是類似的。
2、按名稱匹配讀取
一般來說,用名稱匹配讀取這個方法的占多數,因為雖然grib文件可以看做“倉庫”,一個個氣象變量可以看做“小房間”,但是不同的grib文件內部的“小房間”分布可能是不一樣的,這樣的話同一氣候變量的編號也就不確定,無法單靠一個編號批量讀取很多文件。而利用名稱匹配讀取就不存在這個問題,因為氣象變量的名稱是不會變的。以讀取grib2文件里的2m處的溫度為例,這里使用wgrib2工具,wgrib2工具功能更強大,可以直接讀出為csv文件,代碼如下:
import os
# 將當前目錄轉到wgrib2工具目錄
os.chdir("D:\mydownload\wgrib2")
# 進行氣象變量名稱匹配,讀取2m處溫度這個氣象變量到2m_tmp.csv文件中
os.system("wgrib2 G:\mete_data\Z_NAFP_C_BABJ_20170709000000_P_CNPC-T639-GMFS-HNEHE-03900.grib2 -match ':TMP:2 m' -csv 2m_tmp.csv")
讀取結果如下圖所示:

csv格式的數據更便於我們處理,倒數第一列就是溫度值。
同理,其他的氣候變量同樣可以使用匹配的方式讀取出來,只是改變下代碼-match后面的氣象變量即可,大家可能會問這些氣象變量的簡寫如何得知,這些都是固定的,目前為止我只用了溫度和降水,其他的我也不知道
,有需要的話可以去查一查。
結束
這兩個工具的功能還有很多,目前我用到的只是一丟丟,有其他需要的童鞋可以使用wgrib.exe和wgrib2.exe命令來進行查看幫助,發掘更多功能。在linux下同樣可以使用這兩個工具進行讀取,只是在安裝和命令上有一點小小的差異,大家感興趣的可以一起討論下
,水平有限,肯定有寫的不對的地方,歡迎指正!
