黔西南州區域自動站氣象資料繪圖系統的研制

王文勇，袁洪釗，周 雪，陳 林，何汝凡

(1．貴州省黔西南自治州氣象局，貴州 興義 562400；2．貴州省安龍縣氣象局，貴州 安龍 552400)

?

黔西南州區域自動站氣象資料繪圖系統的研制

王文勇1，袁洪釗1，周 雪2，陳 林2，何汝凡2

(1．貴州省黔西南自治州氣象局，貴州 興義 562400；2．貴州省安龍縣氣象局，貴州 安龍 552400)

隨著氣象服務的不斷推進，用戶對氣象服務產品的要求越來越高，提高氣象服務產品的質量對于基層氣象臺站來說顯得至關重要。目前，區域氣象站的建設已覆蓋到各鄉鎮，將這些區域氣象站資料用直觀的色斑圖來顯示，可有效提高氣象服務質量和服務水平。系統綜合使用VB、Fortran、GrADS等技術，實現對區域自動站氣象資料的查詢、處理與繪圖等功能，并可分縣制作繪圖邊界和色斑圖，為決策氣象服務提供精美素材。

VB語言；Fortran語言；GrADS繪圖；區域自動站；氣象資料

1 引言

近年來，區域自動氣象站建設已覆蓋到各鄉鎮，中小尺度的區域氣象觀測系統已經形成，如何有效地利用這些豐富的氣象資料開展氣象服務，制作出內容豐富、直觀的氣象服務產品，已成為當前縣級氣象部門亟需解決的事情。

在眾多的繪圖軟件中，GrADS(Grid Analysis and Display System)是當今氣象界廣泛使用的一種數據處理和顯示軟件，其具有功能強大、顯示快速、出圖類型多樣化、圖形美觀等特點。雖然GrADS的操作使用命令行方式，但GrADS以其強大的數據分析能力、靈活的環境設置、豐富的出圖類型等特點，已成為國內外氣象界通用的標準圖形環境之一。使用GrADS繪圖時，由于需要人工輸入命令，操作起來比較繁瑣，但結合VB語言進行軟件開發，可以很好地解決這個問題[1]。

Visual Basic是一種由微軟公司開發的包含協助開發環境的事件驅動編程語言，從任何標準來說，VB都是世界上使用人數最多的語言。VB擁有圖形用戶界面(GUI)和快速應用程序開發(RAD)系統，程序員可以輕易地使用DAO、RDO、ADO連接數據庫，輕松地創建ActiveX控件，或輕松地使用VB提供的組件快速建立一個應用程序[2]。GrADS以命令行方式運行，VB中的shell函數為執行GrADS程序提供了可能，免去開發人員不熟悉GrADS的麻煩[3]。

2 實現方法

2.1.1 混合編程 綜合使用VB語言、Fortran語言、GrADS繪圖等技術，實現對區域自動站氣象資料的查詢、處理與繪圖等功能。在VB語言開發環境中，主要使用shell函數分別調用Fortran程序和GrADS程序，如調用GrADS程序的命令格式為：Shell("gradsnc.exe -blc rain.gs", vbHide)，調用Fortran程序的命令格式為：Shell("rain.exe", vbHide) ，其中rain.exe為Fortran語言編譯程序。

2.1.2 進程等待的實現 VB使用shell函數時，采取多進程方式執行。VB程序執行完shell函數后，會立刻轉入下一語句，而不等待外部進程執行結束，導致數據查詢未結束就開始繪圖，或者圖形還未繪制成功就顯示圖形。系統采用進程等待方法，即等待調用的外部進程運行結束后，再繼續執行后面的語句。關鍵代碼如下：

hWnd=OpenProcess(&H100000,0,pid)

Do

DoEvents

Loop Until WaitForSingleObject(hWnd,0)=0

2.2 資料查詢方法

我又打量了一下屋里的擺設，心里更慌了，因為這顯然不是在我家里。我不在自家床上睡覺又會在哪里呢，我努力回憶了幾次，可腦子里始終是一截黏糊糊的空白。

在VB中，事先通過調用OpenCn模塊連接數據庫，然后調用OpenRs模塊打開數據庫記錄集。具體調用方法如下：

OpenCn(Cip, pubs, users, pw)，其中Cip是服務器IP地址，pubs是數據庫名稱，users是用戶名，pw是密碼。

OpenRs(sqlstr)，其中sqlstr是查詢語句，可根據實際需要編寫。系統主要從貴州省氣象局服務器上查詢數據，SQL語句格式如下：

Sqlstr="select StationID,sum(R1H) as sums from dbo." & tabb & " where ( StationID in ( " & aa & ") and (ObservTime between′" & (t1) & "' and ′" & (t2) & "′ )) group by StationID"

其中tabb代表數據庫表tabHourData或tabRainMinute,StationID、R1H分別代表站名和小時雨量，aa代表臺站名稱,t1、t2分別代表查詢開始時間和結束時間。數據查詢輸出格式為純文本，命名為rain.txt，其存儲格式如圖1所示。

圖1 rain.txt文件格式Fig.1 File format of rain.txt

2.3 站點資料處理

根據GrADS繪制站點資料圖流程，事先使用Fortran語言將rain.txt文件處理成GrADS能讀取的二進制數據文件rain.dat，然后將Fortran處理代碼編譯成可供VB調用的可執行程序rain.exe。

數據處理完成后，為rain.dat編寫數據描述文件rain.ctl。由于繪制站點圖需要一個地圖映射文件rain.map，此文件通過執行“stnmap.exe -i rain.ctl”命令生成，在輸入rain.ctl時需要輸入該文件的完整路徑。

再準備供站點繪圖插值用的格點數據文件grid.dat及grid.ctl。根據用戶圖形顯示需要，編寫站點繪圖的批處理文件rain.gs。完成以上準備工作后，即可調用GrADS軟件繪制站點圖。

3 各功能模塊的實現

3.1 系統主界面

系統運行主界面簡單整潔，左側為設置和文本顯示區，右側為圖形顯示區，如圖2所示。

圖2 系統主界面Fig.2 Main interface of system

3.2 繪圖文件準備

3.2.1 地圖文件 利用各縣市邊界的經緯度數據制作GrADS地圖文件，分別命名為：qxn.map、56792.map、57900.map、57902.map、57905.map、57906.map、57907.map、57908.map、57909.map。根據選擇的區域，繪制出相應的地圖。

3.2.2 格點數據文件 對不同區域進行繪圖時，所需要的格點數據文件不同，插值對格點的數量及分布有嚴格的要求。系統為各縣市分別制作了相應的格點數據文件，分別命名為：qxngrid.dat、56792grid.dat、57900grid.dat、57902grid.dat、57905grid.dat 57906

grid.dat、57907grid.dat、57908grid.dat、57909grid.dat。

3.2.3 站點信息文件 繪制站點圖需要站點信息，將黔西南州所有區域站按縣市進行分類，按要求制作為站點信息文件，分別命名為：qxn.txt、56792.txt、57900.txt、57902.txt、57905.txt、57906.txt、57907.txt、57908.txt、57909.txt。軟件在查詢數據時，也是通過讀取這些文件確定繪圖站點。

3.2.4 地圖屏蔽文件 繪制美觀的色斑圖，自然離不開對圖形的美化，利用動力論壇中網友提供的程序province-basemap.gs對地圖以外的區域進行屏蔽，只保留地圖內的圖像，這樣畫出來的圖既美觀又整潔。province-basemap.gs程序需要將地圖的邊界經緯度按要求制作成***_out.txt文件，系統分別為各縣市制作了此文件，分別命名為：qxn_out.txt、56792_out.txt、57900_out.txt、57902_out.txt、57905_out.txt、57906_out.txt、57907_out.txt、57908_out.txt、57909_out.txt。

3.3 主要實現功能

3.3.1 多種繪圖類型 不同的繪圖類型有不同的色斑，繪圖類型主要有：1 h降水量、3 h降水量、6 h降水量、12 h降水量、24 h降水量、累計小時降水量、降水距平百分率、氣溫等共8個類型，用戶可根據需要選擇不同的繪圖類型。

3.3.2 在線實時繪圖 繪圖文件生成后，選擇時間范圍，點擊“在線繪圖”按鈕，系統即可進行數據查詢、處理與自動繪圖，實現一鍵繪圖功能，并將圖形顯示于軟件界面上，生成的圖形文件rain.gif和rain.gmf自動存放在軟件目錄下，其中rain.gif文件讀取方便，但硬拷貝效果差；rain.gmf則具有高質量、硬拷貝效果好等優點，用于在Word、PPT等軟件中進行服務材料制作時效果較好。

3.3.3 診斷繪圖 使用該系統時，如發現某個站點的數據是錯誤的，想對其進行修改后再繪圖，此時在線繪圖就很不方便。系統提供了數據修改功能，只要輸入需要修改數據的臺站號及正確值，點擊“數據修改”按鈕，即可將正確的數據保存在rain.txt文件中，再點診斷繪圖，就可以不重新查詢數據庫，直接讀取修改后的數據進行繪圖。若不知道正確值，可將其值改為缺測值-999.0。

3.3.4 按鄉鎮統計 很多鄉鎮都有兩個以上區域自動站，如只需按鄉鎮開展服務，系統設計了按鄉鎮最大(高)、鄉鎮最小(低)、鄉鎮平均等功能進行統計，默認為不選取。

3.3.5 文本顯示區域 當繪圖類型為雨量時，將所繪站點按從大到小進行排列。當繪圖類型為氣溫時，將所繪站點按從低到高進行排列。統計出來的文本經過簡單的編輯后，可用于編發氣象短信。

4 結語

綜合利用VB編程語言、GrADS繪圖軟件、Fortran語言等各自的特點，設計出一個實用的區域自動站氣象資料繪圖系統。系統實現分縣制作地圖文件，繪制色斑圖，并將地圖以外區域屏蔽掉，繪制生成的圖形產品真實美觀，在雨量實時監測方面更直觀、方便，使值班人員能夠集中精力關注本區域的雨量情況，決策者可從中直觀地了解氣象信息。同時，系統可生成圖形、表格、文字等3種產品，圖形和表格產品用于制作決策氣象服務材料，文字產品經過簡單修改后可用于發布手機短信。另外，系統還實現了區域自動站氣象資料統計功能。自其投入業務運行以來，大大減輕了值班人員的工作量，明顯提高了決策氣象服務現代化水平。

[1] 虞蘇青，鄧世友．淺析用VB調用Grads繪圖軟件的技巧——雷達回波拼圖的報警系統[J]．貴州氣象，2012，36(3)：48-49.

[2] 彭國倫．Fortran95程序設計[M]．北京：中國電力出版社，2002.

[3] 劉輝權，肖天貴，邱靜雅．氣象預報預警系統開發中VB與GrADS、Fortran混合編程技術[J]．成都信息工程學院學報，2011， 26(2)：208-214.A research of graphics system using regional weather station data in Qianxinan

WANG Wenyong1,YUAN Hongzhao1,ZHOU Xue2,CHEN Lin2,HE Rufan2

(1.Qianxinan Meteorological Bureau of Guizhou Province, Xingyi 562400, China;2.Anlong Meteorological Bureau of Guizhou Province, Anlong 552400, China)

With the advancement of meteorological service, user requirements for meteorological service products are getting higher and higher, To improve the quality of the meteorological service products is very important for the meteorological stations at basic level. At present, the construction of the regional meteorological station has been covered by the township, Drawing these regional weather station data into choroplethic graph, can effectively improve the quality and service level of meteorological service. This system can realize regional automatic weather data query, processing and mapping functions using VB, Fortran and GrADS technologies. And it can produce drawing boundary and choroplethic graph, provide fine material for decision-making meteorological service.

VB language; fortran language; GrADS language; regional automatic weather station; meteorological data

1003-6598(2016)04-0061-03

2016-03-08第

王文勇(1987—)，男(苗族)，助工，主要從事綜合氣象觀測及綜合管理工作，E-mail:568949994@qq.com。

P415.1+2

B