windows python读取grib2数据
一、環境準備
(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下同樣可以使用這兩個工具進行讀取,只是在安裝和命令上有一點小小的差異,大家感興趣的可以一起討論下,水平有限,肯定有寫的不對的地方,歡迎指正!
總結
以上是生活随笔為你收集整理的windows python读取grib2数据的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: mac-phpstorm 快捷键
- 下一篇: lvs+keepalived 集群