郭曉東　王世濤　趙斐苗

[摘要]該應(yīng)用軟件采用Borland公司C++Builder6和OpenGL開(kāi)發(fā)軟件,利用Delaunay三角形插入算法,完成了三維雨量面的模擬生成,具備三維雨量面的等雨量線和點(diǎn)雨量輸出的基本功能。

[關(guān)鍵詞]三維雨量面 Delaunay GIS OpenGL Bitmap

我國(guó)是一個(gè)農(nóng)業(yè)大國(guó),對(duì)自然氣候依賴比較重,容易受到農(nóng)業(yè)災(zāi)害的影響。農(nóng)業(yè)災(zāi)害最常規(guī)的表現(xiàn)形式就是干旱或洪澇災(zāi)害,而干旱、洪澇災(zāi)害的重要指標(biāo)之一就是面雨量。如果能及時(shí)掌握了比較準(zhǔn)確的面雨量,對(duì)于農(nóng)業(yè)的抵抗災(zāi)害、江河的抗洪能力以及水庫(kù)的蓄洪規(guī)模都具有重要的意義。

隨著非常規(guī)探測(cè)技術(shù)的迅猛發(fā)展,以自動(dòng)氣象站降水觀測(cè)資料為主,結(jié)合雷達(dá)、衛(wèi)星估測(cè)降水資料輔助修正模式正在形成,通過(guò)資料融合、四維同化來(lái)進(jìn)行估算已成為流域面雨量估算技術(shù)的發(fā)展趨勢(shì)。而以GIS為平臺(tái),實(shí)現(xiàn)數(shù)據(jù)分析和顯示一體化管理則代表了面雨量業(yè)務(wù)系統(tǒng)的發(fā)展方向。通過(guò)水文模擬與水文觀測(cè)的比較,可以對(duì)面雨量估算結(jié)果進(jìn)行間接檢驗(yàn)。如何將雨量和上述資料融合在一起,提高降水量精度和降水?dāng)?shù)量的預(yù)測(cè),是一個(gè)亟待探討和研究的重要問(wèn)題。

三維雨量面,顧名思義就是固定時(shí)間段固定區(qū)域內(nèi)任意一點(diǎn)上的降水不流失,整個(gè)區(qū)域所有點(diǎn)降水將形成一個(gè)立體的降水體。在不考慮地貌特征,認(rèn)為地表是均勻平整的情況下,水體的上表面就是該區(qū)域某時(shí)段內(nèi)的三維雨量面。構(gòu)建三維雨量面要求有大量的數(shù)據(jù)支撐。

三維雨量面可以更直觀的反映面雨量的分布特征,善于表現(xiàn)降水空間分布的復(fù)雜性,同時(shí)可以非常方便的加入或融合到GIS系統(tǒng)中去。如果結(jié)合GIS地理信息提供的地貌資料,風(fēng)場(chǎng)、溫度場(chǎng)、氣壓場(chǎng)等綜合資料,可以進(jìn)一步揭示空間信息與降雨之間的關(guān)系,大大提高降水量預(yù)測(cè)的精度,還可以準(zhǔn)確估算流域的平均降水量,從而更好的預(yù)測(cè)、預(yù)警災(zāi)情,對(duì)于農(nóng)業(yè)的抵抗災(zāi)害、江河的抗洪能力以及水庫(kù)的蓄洪規(guī)模都具有重要的意義。

本文采用利用Borland公司的C++Builder6和OpenGL三維開(kāi)發(fā)軟件(如圖1所示),開(kāi)發(fā)設(shè)計(jì)了河南省雨量面模擬、分析軟件,再現(xiàn)了河南省某時(shí)段的三維雨量面,現(xiàn)簡(jiǎn)介如下。

一、雨量面的數(shù)據(jù)支持

三維雨量面的構(gòu)建必須有大量的數(shù)據(jù)作為支撐,必須選取合適的插入法。經(jīng)比較,本文選取了Delaunay三角形插值法。它是氣象部門和GIS研究都常用的一種插值方法,其優(yōu)點(diǎn)如下:

1.三角形插值法可以較好地表現(xiàn)雨量、地形、地貌和物體表面形狀特征,準(zhǔn)確性可以滿足要求;作為主要的數(shù)字地形模型表示法,它的生成算法已趨于成熟。

2.三角網(wǎng)保證了相同的觀測(cè)站點(diǎn)生成的三角形插值圖像唯一,可重復(fù)性好。易于更新,可適應(yīng)各種分布密度的數(shù)據(jù),減小運(yùn)算量。

3.結(jié)構(gòu)良好,數(shù)據(jù)結(jié)構(gòu)簡(jiǎn)單,數(shù)據(jù)冗余度小,存儲(chǔ)效率高,可以表示線性特征和迭加任意形狀的區(qū)域邊界。

4.三角網(wǎng)三角形數(shù)據(jù)可以直接用于三維物體的軟件生成。

本文將河南省121個(gè)氣象觀測(cè)臺(tái)站作為插入點(diǎn),根據(jù)經(jīng)緯度利用三角形插值法生成了河南省二維delaunay三角形網(wǎng)絡(luò)。程序設(shè)計(jì)過(guò)程如下圖所示。

如圖3所示,為生成的二維Delaunay三角形網(wǎng)絡(luò)。

二、三維開(kāi)發(fā)軟件的選擇及技術(shù)開(kāi)發(fā)

構(gòu)建三維的雨量面,常規(guī)的程序開(kāi)發(fā)軟件難以實(shí)現(xiàn),更難以將多個(gè)面、場(chǎng)疊加,最好采用專業(yè)的三維開(kāi)發(fā)軟件。當(dāng)前最主流的三維開(kāi)發(fā)軟件有DirectX和OpenGL,是兩款非常優(yōu)秀的軟件。由于氣象部門數(shù)據(jù)運(yùn)算需求,經(jīng)常要在UNIX、LIUNX等操作系統(tǒng)下操作,鑒于OpenGL(Open Graphic Library)具有良好的程序通用性和可移植性,本文選擇了采用OpenGL來(lái)解決面雨量的三維顯示。

三、雨量面的分析

本文對(duì)于三維雨量面的分析分析包含最基本的功能,主要有兩點(diǎn):點(diǎn)雨量的獲取和等雨量線的獲取。

1.點(diǎn)雨量的獲取。點(diǎn)雨量的獲取非常簡(jiǎn)單。我們只要知道了需要查找的點(diǎn)坐標(biāo)位置,即可獲得需要的坐標(biāo)點(diǎn)上像素顏色,顏色和雨量之間關(guān)系換算即可得到該點(diǎn)雨量。

2.等雨量線的獲取。

等雨量線是一種等值線。等值線是聯(lián)結(jié)具有相同數(shù)量指標(biāo)的空間分布現(xiàn)象與二維平面的圖形,不同數(shù)值的等值線組成的圖形可以精確表示自然現(xiàn)象的垂直變化和水平方向的強(qiáng)弱差異,是一項(xiàng)重要的分析指標(biāo)。由于空間的限制,等值線可以是閉合的也可以是不閉合的。等值線一般不相交。

本文提出一種新穎的等值線獲取方法,就是通過(guò)對(duì)比位圖色彩值來(lái)獲取等值線。上文構(gòu)建了三維的雨量面,同時(shí)生成輸出該雨量面的彩色正視圖。正視圖是三維雨量面在水平面上的正投影。由于雨量面是內(nèi)部連續(xù)的,雨量與色彩值是線性相關(guān)的,相同雨量值色彩值必定相同,在色彩值相同的地方,必定存在等雨量線。我們可以開(kāi)發(fā)編寫軟件,通過(guò)對(duì)雨量面正投影位圖的色彩分析,通過(guò)過(guò)濾生成等雨量線,設(shè)置多個(gè)雨量段對(duì)應(yīng)的顏色段閥,通過(guò)與非運(yùn)算,篩去不符合閥值的點(diǎn),即可輕易獲得任意等雨量線,僅生成一次雨量面,即可完成所有任務(wù),對(duì)于生成的等雨量線還可以方便保存為位圖。程序編寫過(guò)程如下。

四、結(jié)論

本文通過(guò)采用Delaunay三角形插值方法,保證了三維雨量面的數(shù)據(jù)需求;采用了OpenGL技術(shù),完成了三維雨量面的生成;采用位圖分析技術(shù)來(lái)獲取三維雨量面上的點(diǎn)雨量輸出和任意等雨量線的輸出,將格點(diǎn)插入轉(zhuǎn)換為對(duì)位圖的處理,完成了雨量面的基本數(shù)據(jù)輸出,不受格點(diǎn)限制,減輕了計(jì)算量,提高了工作效率。同時(shí)并為下一步引入GIS地面系統(tǒng)奠定基礎(chǔ)。

參考文獻(xiàn):

[1]R. K. L inslcg 等著. 劉光文等譯.工程水文學(xué).水利出版社, 1991.18～81.

[2]Hutchinson MF.Interpolating mean rainfall using thin plate smoothing splines [J ] . Int . J . Geogr. Inf . Sys. ,1995:385～403.

[3] 葉飛躍等,等值線線性追蹤算法及其在地圖重建中的應(yīng)用[J].計(jì)算機(jī)仿真,2003,20(6).