等值線在線修改系統(tǒng)的設(shè)計(jì)與應(yīng)用

徐映雪，牛鑫艷，翟彩霞，薛　偉

(1．北京清流技術(shù)股份有限公司，北京 100073；2．北京清流地下水資源研究所，北京 100073)

等值線在線修改系統(tǒng)的設(shè)計(jì)與應(yīng)用

徐映雪1，2，牛鑫艷1，翟彩霞1，2，薛偉1

(1．北京清流技術(shù)股份有限公司，北京 100073；2．北京清流地下水資源研究所，北京 100073)

摘要：等值線面是水利信息化數(shù)據(jù)處理、成果展示及預(yù)測(cè)預(yù)警的重要方法。本文介紹了自主開(kāi)發(fā)的等值線在線修改系統(tǒng)Conmas，利用數(shù)據(jù)流技術(shù)通過(guò)對(duì)原始樣本數(shù)據(jù)的編輯實(shí)現(xiàn)了可視化的等值線自動(dòng)修改，極大提高了等值線繪制工作的便捷性；并結(jié)合水利信息化業(yè)務(wù)需要，改進(jìn)了裁剪功能，拓展了統(tǒng)計(jì)分析、剖面圖繪制、特殊區(qū)域識(shí)別等功能。等值線算法庫(kù)經(jīng)二次開(kāi)發(fā)后已成功應(yīng)用于水利部水情信息中心的地下水業(yè)務(wù)與分析系統(tǒng)、地質(zhì)災(zāi)害監(jiān)測(cè)預(yù)警平臺(tái)等項(xiàng)目中。

關(guān)鍵詞：插值；等值線；在線修改；數(shù)據(jù)流；Conmas

等值線面作為數(shù)據(jù)處理及成果展示的重要方法，在水文、氣象、地質(zhì)災(zāi)害監(jiān)測(cè)等許多方面有著深入的研究應(yīng)用，尤其是在地質(zhì)災(zāi)害預(yù)測(cè)領(lǐng)域中也有著廣泛的應(yīng)用[1]。它可以有效地表示雨量、地下水位、墑情、溫度等幾何場(chǎng)變量的大小、分布信息，等值線之間的距離反映了場(chǎng)變量的變化大小。傳統(tǒng)的手工繪圖方法計(jì)算量大、工作繁重。隨著20世紀(jì)中期電子計(jì)算機(jī)的出現(xiàn)，等值線的繪制也開(kāi)始進(jìn)入程序化和計(jì)算機(jī)化的新階段[2-4]，目前國(guó)內(nèi)外的商業(yè)繪圖軟件基本都有繪制等值線的功能。最常見(jiàn)的等值線繪圖軟件為美國(guó)Golden Software 公司編制的Surfer8.0[5]，其處理二維數(shù)據(jù)，主要面向等高線、三維表面圖的繪制，最大的不足為圖形與常用GIS軟件接口較差；另外美國(guó)MathWorks公司出品的商業(yè)數(shù)學(xué)軟件Matlab[6]，是用于算法開(kāi)發(fā)、數(shù)據(jù)可視化、數(shù)據(jù)分析以及數(shù)值計(jì)算的高級(jí)計(jì)算語(yǔ)言和交互式環(huán)境，Matlab對(duì)數(shù)據(jù)的處理功能強(qiáng)大，但對(duì)于繪制等值線來(lái)說(shuō)，需要編程實(shí)現(xiàn)，且對(duì)等值線的編輯功能較弱，對(duì)于初學(xué)者上手較難，該軟件主要面向氣象數(shù)據(jù)的處理和展示；此外還有美國(guó)Tecplot公司推出的數(shù)據(jù)處理軟件Tecplot系列軟件[7]、美國(guó)馬里蘭大學(xué)開(kāi)發(fā)的氣象數(shù)據(jù)分析與現(xiàn)實(shí)軟件GrADS、Esri開(kāi)發(fā)的ArcGIS[8]，國(guó)產(chǎn)軟件GeoStar 的GeoTIN 模塊、MapGIS[9]的DTM分析模塊等均能生成等值線，但它們?cè)诘戎稻€自動(dòng)生成、編輯、數(shù)據(jù)接口以及信息系統(tǒng)二次開(kāi)發(fā)集成等方面均有一定的不足之處。

本文主要介紹了基于GIS平臺(tái)開(kāi)發(fā)的一套等值線在線修改系統(tǒng)，利用數(shù)據(jù)流技術(shù)通過(guò)對(duì)原始樣本數(shù)據(jù)的編輯實(shí)現(xiàn)了可視化的等值線自動(dòng)修改，極大提高了等值線繪制工作的便捷性；結(jié)合水利信息化業(yè)務(wù)需要，改進(jìn)了裁剪功能，拓展了統(tǒng)計(jì)分析、剖面圖繪制、特殊區(qū)域識(shí)別等功能。等值線算法庫(kù)經(jīng)二次開(kāi)發(fā)后已成功應(yīng)用于水利部水情信息中心的地下水業(yè)務(wù)與分析系統(tǒng)、地質(zhì)災(zāi)害監(jiān)測(cè)預(yù)警平臺(tái)等項(xiàng)目中。

1等值線在線修改流程

1.1等值線生成基本步驟

等值線生成的步驟一般是遵循一定規(guī)律的、模式化的，這為實(shí)現(xiàn)等值線繪制過(guò)程的程序化提供了前提。首先利用原始的樣本散點(diǎn)數(shù)據(jù)創(chuàng)建矩形網(wǎng)格類型的數(shù)字模型(GRD模型)，進(jìn)而進(jìn)行等值點(diǎn)查找、等值線追蹤、等值面構(gòu)建等，等值線生成，基本步驟如下：

步驟一：獲取樣本散點(diǎn)數(shù)據(jù)文件，必須包括觀測(cè)點(diǎn)或取樣點(diǎn)的坐標(biāo)和觀測(cè)數(shù)據(jù)。

步驟二：生成矩形網(wǎng)格數(shù)據(jù)文件，即計(jì)算網(wǎng)格節(jié)點(diǎn)值，由散點(diǎn)數(shù)據(jù)經(jīng)插值(克里格、反距離加權(quán)等)計(jì)算而得。

步驟三：等值點(diǎn)計(jì)算，根據(jù)矩形網(wǎng)格數(shù)據(jù)線性內(nèi)插得到穿過(guò)網(wǎng)格線的等值點(diǎn)。

步驟四：等值線追蹤，即等值點(diǎn)相連，構(gòu)成閉合或與插值區(qū)域邊界線相交的曲線。

1.2等值線在線修改流程

樣本數(shù)據(jù)的多寡、空間分布、數(shù)據(jù)正確性、人工經(jīng)驗(yàn)等因素決定了等值線的形狀。因此實(shí)際工作中，往往是通過(guò)人工判斷等值線是否符合區(qū)域?qū)嶋H，符合區(qū)域?qū)嶋H則完成等值線繪制工作，否則需要編輯樣本散點(diǎn)數(shù)據(jù)，通過(guò)對(duì)樣本散點(diǎn)數(shù)據(jù)的反復(fù)編輯，最終生成符合區(qū)域?qū)嶋H情況的等值線圖。完成等值線的繪制工作后，往往還要根據(jù)實(shí)際的業(yè)務(wù)工作需要，進(jìn)行等值線的裁減、標(biāo)注、面積量算等。

因此本文提出利用數(shù)據(jù)流技術(shù)，建立樣本散點(diǎn)數(shù)據(jù)、矩形網(wǎng)格數(shù)據(jù)、等值線等結(jié)果數(shù)據(jù)之間的單向連通性和傳遞性，通過(guò)對(duì)原始樣本數(shù)據(jù)的編輯后自動(dòng)觸發(fā)該流程，最終實(shí)現(xiàn)等值線圖的自動(dòng)繪制和修改。

2等值線在線修改系統(tǒng)實(shí)現(xiàn)與應(yīng)用

2.1設(shè)計(jì)目標(biāo)

等值線在線修改系統(tǒng)的主要設(shè)計(jì)目標(biāo)可歸結(jié)如下：借助GIS技術(shù)實(shí)現(xiàn)散點(diǎn)圖可視化及在線編輯；建立散點(diǎn)數(shù)據(jù)和等值線圖之間數(shù)據(jù)流的連通性和單向傳遞性，以及自動(dòng)觸發(fā)機(jī)制，一旦樣本散點(diǎn)數(shù)據(jù)發(fā)生變化，則等值線圖以及基于此樣本數(shù)據(jù)生成的剖面圖、特殊區(qū)域等均實(shí)現(xiàn)同步更新計(jì)算功能。這樣便可直接在散點(diǎn)圖上對(duì)原始數(shù)據(jù)進(jìn)行可視化編輯，等值線圖實(shí)現(xiàn)自動(dòng)重繪和在線修改。

2.2系統(tǒng)結(jié)構(gòu)

等值線自動(dòng)繪制系統(tǒng)Conmas主要部分采用CS結(jié)構(gòu)設(shè)計(jì)，其具有表現(xiàn)元素多元化，交互性能好，安全性能易保證、客戶端響應(yīng)速度快等優(yōu)勢(shì)。程序總體框架由三個(gè)模塊構(gòu)成：數(shù)據(jù)管理模塊、數(shù)據(jù)處理模塊、圖形輸出模塊。

數(shù)據(jù)管理模塊主要實(shí)現(xiàn)創(chuàng)建數(shù)據(jù)庫(kù)并儲(chǔ)存數(shù)據(jù)，例如散點(diǎn)數(shù)據(jù)、剖線數(shù)據(jù)、裁剪邊界數(shù)據(jù)、網(wǎng)格數(shù)據(jù)等。數(shù)據(jù)按照點(diǎn)、線、面不同類型進(jìn)行分類儲(chǔ)存。該模塊為繪圖的基礎(chǔ)，提供數(shù)據(jù)支持和調(diào)用。

數(shù)據(jù)處理模塊主要實(shí)現(xiàn)對(duì)圖形的輔助計(jì)算功能。例如散點(diǎn)的插值計(jì)算、數(shù)據(jù)篩選、剖面數(shù)據(jù)的計(jì)算、圖形裁剪、基于網(wǎng)格的體積面積計(jì)算、等值線光滑算法、圖形的長(zhǎng)度面積計(jì)算等。該模塊是Conmas的核心技術(shù)部分，體現(xiàn)了等值線的高級(jí)應(yīng)用。

圖形輸出模塊主要實(shí)現(xiàn)圖形的繪制、編輯、保存及輸出功能。例如散點(diǎn)圖、等值線面圖、裁剪邊界圖、柵格底圖、剖面圖、Shp圖等圖形元件的繪制。該模塊輸出等值線圖，生成用戶所需要的最終效果圖。

2.3主要功能及應(yīng)用

除了等值線的在線修改基本功能之外，結(jié)合水利信息化業(yè)務(wù)要求開(kāi)發(fā)了分區(qū)統(tǒng)計(jì)、特殊區(qū)域繪制、剖面圖等功能。

(1)等值線在線修改

在實(shí)際工作中，監(jiān)測(cè)站點(diǎn)分布及其監(jiān)測(cè)數(shù)據(jù)的質(zhì)量往往會(huì)影響等值線的繪制，具體包括：①區(qū)域控制性監(jiān)測(cè)站點(diǎn)不足；②常規(guī)監(jiān)測(cè)站點(diǎn)變更，造成數(shù)據(jù)的不連續(xù)；③由于自動(dòng)監(jiān)測(cè)設(shè)備損壞等原因造成的漏測(cè)等，導(dǎo)致監(jiān)測(cè)數(shù)據(jù)不全；④監(jiān)測(cè)數(shù)據(jù)不準(zhǔn)確，比如地下水水位觀測(cè)井為生產(chǎn)結(jié)合井時(shí)，可能觀測(cè)的是動(dòng)水位數(shù)據(jù)；⑤統(tǒng)測(cè)類站點(diǎn)存在的站點(diǎn)變更或統(tǒng)測(cè)數(shù)據(jù)由于人為原因造成的不精確等。因此需要根據(jù)區(qū)域的地形、水文地質(zhì)情況等，通過(guò)站點(diǎn)的管理來(lái)繪制符合區(qū)域?qū)嶋H的等值線。

借助GIS平臺(tái)設(shè)計(jì)了針對(duì)監(jiān)測(cè)站散點(diǎn)數(shù)據(jù)的人機(jī)交互式的可視化操作界面(圖1)，包括可不使用一個(gè)或多個(gè)明顯不合理的監(jiān)測(cè)數(shù)據(jù)；可虛擬控制性監(jiān)測(cè)站點(diǎn)數(shù)據(jù)，在目標(biāo)區(qū)域點(diǎn)擊后輸入虛擬監(jiān)測(cè)數(shù)據(jù)即可；根據(jù)監(jiān)測(cè)數(shù)據(jù)的連續(xù)性可對(duì)時(shí)間序列不全的數(shù)據(jù)進(jìn)行插補(bǔ)；可直接修改不合理的監(jiān)測(cè)數(shù)據(jù)等。通過(guò)編輯操作散點(diǎn)數(shù)據(jù)即可實(shí)現(xiàn)等值線的自動(dòng)繪制和修改功能。

圖1　等值線在線修改系統(tǒng)界面

(2)等值線裁剪的改進(jìn)

目前大部分軟件對(duì)圖形的裁剪最常用的處理方法是基于矩形網(wǎng)格的裁剪方法，即以實(shí)際裁剪的邊界線所經(jīng)過(guò)的網(wǎng)格為邊界對(duì)圖形進(jìn)行裁剪[10]。該方法思路簡(jiǎn)單編程較易實(shí)現(xiàn)，但應(yīng)用此方法進(jìn)行圖形裁剪存在很大的缺陷，裁剪后的圖形邊界不光滑，存在明顯的鋸齒狀；網(wǎng)格單元的尺寸越大這樣的不光滑性越明顯。若基于網(wǎng)格無(wú)限細(xì)分思想進(jìn)行改進(jìn)則會(huì)造成計(jì)算效率低、占用空間大，而且在GIS平臺(tái)下，很難完美解決這一問(wèn)題。

為了避免上述缺點(diǎn)，我們采用了一種新算法——基于多邊形的裁剪方法，使繪圖邊界更加光滑更符合實(shí)際(圖2)。其中的多邊形都可以是一般多邊形，既可以是凹多邊形，也可以是有內(nèi)孔的多邊形。該算法使用單線性鏈表數(shù)據(jù)結(jié)構(gòu)，與其他使用雙鏈表或樹(shù)結(jié)構(gòu)的算法相比，具有最簡(jiǎn)單的結(jié)構(gòu)和最快的執(zhí)行速度。新算法的具體步驟描述如下：①首先指定裁剪框window和待裁剪圖形polygon,給定的裁剪框和待裁剪圖形均是具有方向的多邊形結(jié)構(gòu)，即多邊形的頂點(diǎn)是以一定順序(順或逆時(shí)針)輸入的；②對(duì)待裁剪的polygon進(jìn)行方向調(diào)整(與裁剪框的輸入方向一致)；③對(duì)待裁剪多邊形與裁剪框多邊形求交點(diǎn)集。對(duì)于交點(diǎn)集按照待裁剪多邊形的序號(hào)(pi)從小到大排序，然后記錄并保存由入邊到出邊之間的結(jié)點(diǎn)；同樣的交點(diǎn)集按照裁剪框Window邊的序號(hào)(wi)從小到大排序，然后記錄并保存由出邊到入邊之間的結(jié)點(diǎn)；④連接這些交點(diǎn)得到多個(gè)閉合的多邊形就是待裁剪多邊形polygon與裁剪框window 的交集。

圖2　基于網(wǎng)格法和多邊形法裁剪效果對(duì)比圖

(3)剖面圖繪制

剖面圖展示的是沿著某一剖線斷面的Z屬性值隨X軸坐標(biāo)的變化曲線(圖3)。根據(jù)不同走向的剖面圖，我們可以分析出沿某一剖線的數(shù)據(jù)分布特征、判斷異常的走向和分布范圍等信息。比如在水資源領(lǐng)域，我們可以用來(lái)反應(yīng)沿某一剖面方向地下水水位、埋深等監(jiān)測(cè)因素的變化情況。

剖面圖上的剖面數(shù)據(jù)點(diǎn)無(wú)外乎以下三類：一、剖線的首尾端點(diǎn)；二、剖線轉(zhuǎn)折點(diǎn)；三、剖線和網(wǎng)格線的交點(diǎn)。如圖4所示,做一條剖線和網(wǎng)格線交于E、G兩點(diǎn)；ABCD為4個(gè)網(wǎng)格節(jié)點(diǎn)(Z值已知)；F為剖線的轉(zhuǎn)折點(diǎn)；E、F、G均是剖面數(shù)據(jù)點(diǎn)。

圖3　基于等值線的剖面圖

圖4　剖線切割網(wǎng)格線圖

基于網(wǎng)格的剖面數(shù)據(jù)提取算法的基本原理是：首先由散點(diǎn)插值得到網(wǎng)格數(shù)據(jù)文件，然后對(duì)剖線附近的網(wǎng)格數(shù)據(jù)進(jìn)行線性內(nèi)插或雙線性內(nèi)插得出剖面數(shù)據(jù)Z值。具體算法如下：如果某個(gè)剖面數(shù)據(jù)點(diǎn)正好和某個(gè)網(wǎng)格節(jié)點(diǎn)重合，則該點(diǎn)的Z值取為網(wǎng)格節(jié)點(diǎn)的Z值；如果剖面數(shù)據(jù)點(diǎn)的位置不和任何網(wǎng)格節(jié)點(diǎn)重合，則需要對(duì)周圍的網(wǎng)格節(jié)點(diǎn)進(jìn)行線性插值或雙線性插值。

(4)特殊區(qū)域繪制

基于等值面對(duì)某些特殊區(qū)域(如特定等級(jí)區(qū)間、凸起區(qū)域或凹陷區(qū)域，如圖5)進(jìn)行繪圖并計(jì)算相關(guān)面積，是在等值線成果的基礎(chǔ)上進(jìn)一步進(jìn)行分析的高級(jí)應(yīng)用功能。特殊區(qū)域圖形繪制的實(shí)質(zhì)是基于網(wǎng)格數(shù)據(jù)的等值面重繪。特殊區(qū)域的繪圖和計(jì)算可為多領(lǐng)域中的實(shí)際工程提供圖形和數(shù)值依據(jù)，例如土地規(guī)劃和地基施工中對(duì)填挖方的計(jì)算、地下水資源評(píng)價(jià)中降落漏斗區(qū)的繪制等。

圖5　特殊區(qū)域繪制

(5)統(tǒng)計(jì)分析

統(tǒng)計(jì)分析以等值面圖為基礎(chǔ)，在疊加行政分區(qū)等分區(qū)因素的基礎(chǔ)上計(jì)算出等值面特征值及其分區(qū)面積或面積占比的量化數(shù)據(jù)。比如每個(gè)月水利部水情信息中心都要匯總各個(gè)省的地下水監(jiān)測(cè)數(shù)據(jù)，按照《地下水動(dòng)態(tài)月報(bào)》編制大綱要求，計(jì)算統(tǒng)計(jì)時(shí)段內(nèi)地下水水位變幅。在其業(yè)務(wù)分析與應(yīng)用系統(tǒng)以等值面技術(shù)二次開(kāi)發(fā)為基礎(chǔ)，結(jié)合平原區(qū)套省市分區(qū)實(shí)現(xiàn)此業(yè)務(wù)功能。

3結(jié)論

等值線在線修改系統(tǒng)的創(chuàng)新點(diǎn)為建立等值線數(shù)據(jù)流的連續(xù)性，并結(jié)合GIS技術(shù)實(shí)現(xiàn)了人機(jī)交互式的可視化的等值線自動(dòng)重繪，解決了等值線軟件和外部程序銜接和實(shí)用性差的問(wèn)題，利用網(wǎng)格計(jì)算等值線在水利信息化業(yè)務(wù)領(lǐng)域的應(yīng)用方面做出了新嘗試，例如剖面數(shù)據(jù)計(jì)算、分區(qū)面積計(jì)算、漏斗計(jì)算等功能。本系統(tǒng)的研發(fā)經(jīng)歷了上百個(gè)應(yīng)用系統(tǒng)的檢驗(yàn)，為今后的數(shù)據(jù)后處理工作提供了一種直觀便捷的等值線獲取方法。

參考文獻(xiàn)：

[1]趙敬和.基于矩形網(wǎng)格法的不規(guī)則區(qū)域的等值線生成與填充算法研究[D].北京：中國(guó)地質(zhì)大學(xué)，2013.

[2]宋麗娟，龔曉峰，鐘猛.基于網(wǎng)格法的等值線繪制方法[J].現(xiàn)代電子技術(shù)，2005(14)：65-67.

[3]成建梅，陳崇希，孫紅林.三角網(wǎng)格等值線自動(dòng)生成方法及程序?qū)崿F(xiàn)[J].水利學(xué)報(bào)，1998(10)：23-26.

[4]郭鵬，董蘭芳，夏澤舉.地質(zhì)數(shù)據(jù)的等值線繪制方法研究[J].計(jì)算機(jī)仿真，2009，26(9)：168-171.

[5]陳歡歡，李星，丁文秀.Surfer 8.0等值線繪制中的十二種插值方法[J].工程地球物理學(xué)報(bào)，2007，4(1)：52-57.

[6]蘇萬(wàn)益，劉衛(wèi)林，屈吉鴻.基于Matlab的地下水?dāng)?shù)值模擬計(jì)算的數(shù)據(jù)后處理[J].華北水利水電學(xué)院學(xué)報(bào)，2004，25(1)：50-53.

[7]孫海梅.繪圖軟件在數(shù)據(jù)處理中的應(yīng)用[J].運(yùn)城學(xué)院學(xué)報(bào)，2011，29(2)：50-52.

[8]許時(shí)光，鄭新奇.基于ArcGIS Engine的WebGIS空間分析功能拓展[J].測(cè)繪通報(bào)，2009(7)：25-27.

[9]李艷，龐振宇.利用MapGIS軟件繪制地質(zhì)等值線圖[J].現(xiàn)代科技：現(xiàn)代物業(yè)下旬刊，2010，9(5)：111-112.

[10] 劉勇奎，顏葉，石教英.一個(gè)有效的多邊形窗口的線裁剪算法[J].計(jì)算機(jī)學(xué)報(bào)，1999，22(11)：1209-1214.

The design and application of the online modification system of contour lines

XU Yingxue1，2，NIU Xinyan1，ZHAI Caixia1，2，XUE Wei1

(1.BeijingQingliuTechnologyCo.Ltd.，Beijing100073，China；2.BeijingQingliuGroundwaterResourcesInstitute，Beijing100073，China)

Abstract：Contour or isosurface is an important method of data processing, project presentation and forecasting-warning in water resources informatization. This article describes the self-development online modification system of contour lines (Conmas). The data stream technology is used to edit the original sample data that achieves the visualization of contour modified automatically, greatly improves the convenience of contour mapping. Combined the business needs of water resources informatization, it improves the function of cutting and expands the functions of the statistical analysis, cross-sectional drawing, special regional recognition etc.. The contour algorithm library after the second development has been applied to many projects of the Ministry of Water Resources Hydrologic Information Center successfully, such as groundwater business and analysis system, geological disaster monitoring and forecasting-warning platform.

Key words：interpolation; contour; online modification; data stream; Conmas

基金項(xiàng)目：北京市西城區(qū)科技計(jì)劃項(xiàng)目(2015IH12)

作者簡(jiǎn)介：徐映雪(1978-)，女，高級(jí)工程師，從事水文水資源及水利信息化研究工作。

中圖分類號(hào)：TV22;X830.3

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：2096-0506(2016)05-0012-05