時(shí)慧潔

摘 要：簡單敘述了流域河網(wǎng)提取的兩種辦法，并以ArcGIS為操作平臺，介紹了基于坡面徑流漫流模型的數(shù)字河網(wǎng)提取原理和方法，其流程可簡單分為：DEM的預(yù)處理、水流流向的確定、匯流累積的計(jì)算以及河網(wǎng)生成。該文基于SRTM DEM數(shù)據(jù)以青龍河水域?yàn)檠芯繉ο筮M(jìn)行流域特征信息提取，提取結(jié)果與國家河網(wǎng)數(shù)據(jù)基本吻合，并就不同閾值劃分下的河網(wǎng)進(jìn)行比較。

關(guān)鍵詞：河網(wǎng) SRTM DEM ArcGIS 青龍河流域

中圖分類號：P33 文獻(xiàn)標(biāo)識碼：A 文章編號：1672-3791（2015）05（b）-0024-02

Abastract：The paper briefly describes two methods of extracting river basins.With the Hydrology Analysis tools of ArcGIS， the principle and method of ground water overflow modelare well introduced，the process of which includes the flow direction， flow accumulation， drainage network andthedrainagenetworkextraction. Based on the Qinglong River watershed SRTM DEM data， the fact that the river networksturn out with different threshold values is discussed. And the extraction results were verified and in accordance with the basic national river network data.

Key Words：drainagenetwork； SRTM DEM； ArcGIS； Qinglong River

流域河網(wǎng)水系是描述一個(gè)地區(qū)地理地貌和水文特征基本的重要條件，對流域整體規(guī)劃、水資源優(yōu)化配置、水利水電開發(fā)利用以及防洪抗旱、航運(yùn)、養(yǎng)殖、灌溉和旅游等都具有十分重大的意義[1]。目前，利用DEM數(shù)據(jù)提取流域數(shù)字河網(wǎng)主要有谷點(diǎn)提取模型和基于流向提取模型兩種[2]。谷點(diǎn)提取模型最初是由Greysukh，Pekuker， Douglas[3]等人提出，其原理是通過對比所研究網(wǎng)格的高程與相鄰網(wǎng)格的高程找出DEM中的谷點(diǎn)單元，再將谷點(diǎn)單元連接后形成河網(wǎng)[4]。該方法的問題是利用谷點(diǎn)單元形成的河網(wǎng)往往不連續(xù)，且在地形起伏較小的平原地區(qū)難以識別谷點(diǎn)單元[5]。而目前廣泛應(yīng)用的算法是基于OCallaghan和 Mark提出的坡面徑流模擬算法，其原理是先確定流域各柵格間的流向，再根據(jù)每個(gè)柵格點(diǎn)上游的計(jì)算得到匯流累積量，最后按不同閾值提取出河網(wǎng)。其中確定流向的方法多采用D8的算法，曹玲玲[6]、吳冰[7]、周國園[8]、沈中原[9]、熱紫燕[10]等人都在不同流域基于坡面徑流模擬算法進(jìn)行了河網(wǎng)提取。該文主要介紹利用DEM數(shù)據(jù)基于坡面徑流漫流模型的數(shù)字河網(wǎng)提取方法，并在ArcGIS環(huán)境下通過青龍河流域的實(shí)例進(jìn)行分析。

1 研究區(qū)背景概況

1.1 流域概況

青龍河是秦皇島市的重要水源地，青龍河流域是秦皇島市限制納污控制考核的重要水功能區(qū)。青龍河流經(jīng)灤河流域東側(cè)，地處東經(jīng)118°37′～119°37′，北緯39°51′ ～41°07′之間，流域面積6340 km2，占灤河流域面積的14%，是灤河水系年平均來水量最大的支流。該文擬以桃林口水壩址以上流域?yàn)槔?，進(jìn)行流域數(shù)字河網(wǎng)的提取研究。

1.2 數(shù)據(jù)來源

DEM的數(shù)據(jù)來源是SRTM DEM（Shuttle Radar Topography Mission）數(shù)據(jù)，是2000年2月11日美國發(fā)射的“奮進(jìn)”號航天飛機(jī)上搭載的SRTM系統(tǒng)采集的地表面數(shù)據(jù)。SRTM DEM數(shù)據(jù)采用WGS1984投影參數(shù)，通過中國科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國際科學(xué)數(shù)據(jù)鏡像網(wǎng)站可下載空間分辨率為 3″（約90m）的中國區(qū)域數(shù)據(jù)[11]。由于下載的圖像數(shù)據(jù)不一定剛好在一個(gè)圖幅范圍內(nèi)，需對下載的數(shù)據(jù)進(jìn)行剪裁及修正，從而得到青龍河流域部分。

2 河網(wǎng)提取方法

2.1 DEM預(yù)處理

在利用SRTM數(shù)據(jù)提取流域的數(shù)字河網(wǎng)之前，首先需要對SRTM數(shù)據(jù)進(jìn)行預(yù)處理，以消除SRTM中存在的無效數(shù)據(jù)區(qū)域。此外，由于原始的DEM圖像中數(shù)據(jù)缺陷或是一些真實(shí)地形的存在，使得DEM表面存在著一些凹陷的區(qū)域，使水流無法按照其本身流向流出。所以為避免水流出現(xiàn)逆流或中斷的現(xiàn)象，首先要對DEM進(jìn)行無洼地的填充。洼地填充的基本過程是先利用水流方向數(shù)據(jù)計(jì)算出DEM數(shù)據(jù)中的洼地區(qū)域，并計(jì)算其洼地深度，再依據(jù)這些洼地深度設(shè)定填充閾值進(jìn)行洼地填充。

2.2 流向確定

柵格單元的流向是指水流出該柵格單元的方向。每個(gè)柵格中的流向是通過計(jì)算該柵格周圍最大的落差坡度來確定的，這種計(jì)算方法稱為D8法。D8法將柵格中的水流方向分為8個(gè)方向，即與該柵格相鄰的8個(gè)柵格分別代表一個(gè)方向。在此基礎(chǔ)上計(jì)算中心網(wǎng)格與各相鄰網(wǎng)格間的距離權(quán)落差，取距離權(quán)落差最大的網(wǎng)格與中心網(wǎng)格的連線作為水流流出方向。水流方向的確定過程就是在柵格周圍尋找坡度變化最快的方向。

2.3 匯流累積量計(jì)算

區(qū)域匯流累積量的計(jì)算需要建立在水流方向的基礎(chǔ)上，實(shí)際是曲面上水流過程的模擬。水流方向確定后便可計(jì)算每個(gè)柵格單元的上游集水區(qū)范圍，并生成與之對應(yīng)的水流方向矩陣，計(jì)算該矩陣中某個(gè)柵格的累積流量，正是根據(jù)該柵格周圍的流向數(shù)值編碼計(jì)算而來的。

2.4 河網(wǎng)生成

通過匯流累積量的計(jì)算，完成柵格之間流量傳遞的過程，匯流累積量代表著柵格在流量傳遞過程中累計(jì)的柵格數(shù)。每一個(gè)柵格都攜帶著一份流量，當(dāng)柵格攜帶的流量在同一位置累計(jì)到一定值時(shí)，便會(huì)產(chǎn)生水流，這些水流形成的路徑就是河網(wǎng)，累計(jì)到形成水流的值稱為河網(wǎng)閾值。其原理是把所有大于或等于給定閾值的匯流柵格單元標(biāo)記為1，其余柵格則沒有數(shù)據(jù)，遂將賦值為1的單元格提取出來。

3 實(shí)驗(yàn)與結(jié)論

利用ArcGIS的水文分析模塊，運(yùn)用坡面徑流模擬方法，提取桃林口書庫壩址上游流域的河網(wǎng)。對區(qū)域設(shè)置不同的閾值，所繪制出的河網(wǎng)疏密程度不同（圖2）。當(dāng)閾值較大時(shí)，形成水流所需的累計(jì)柵格數(shù)較多，不易形成河網(wǎng)；當(dāng)設(shè)置閾值較小時(shí)，形成水流所需要的柵格數(shù)較少，容易形成河網(wǎng)。設(shè)定匯流累積閾值時(shí)，要充分考慮水系的實(shí)際情況，并通過觀察和與地形圖相比較等方法適當(dāng)?shù)倪x定閾值，保證徑流模擬準(zhǔn)確性和實(shí)際性。如當(dāng)面積閾值從25.8164 km2（柵格數(shù)為3000）調(diào)整為129.0822 km2（柵格數(shù)為15000）時(shí)，河流分支數(shù)量相應(yīng)從197個(gè)河道減少到37個(gè)河道。在確定上游集水面積閾值為15000個(gè)柵格單元時(shí)的提取結(jié)果與1：25萬的國家五級河網(wǎng)疊加分析可知（圖2（a）），圖中河流水系與實(shí)際河網(wǎng)基本吻合能夠滿足水文模擬與生產(chǎn)需要。

該文介紹了基于ArcGIS平臺支持下河網(wǎng)特征的提取原理，以青龍河流域?yàn)榘咐龑ζ溥M(jìn)行DEM的預(yù)處理、水流流向的確定、匯流累積的計(jì)算，并基于流域不同集水面積閾值提取不同疏密程度的水系特征。選用空間分辨率為90m的SRTM DEM數(shù)據(jù)作為數(shù)據(jù)來源避免了傳統(tǒng)地面測繪容易出現(xiàn)的數(shù)據(jù)插值誤差從而導(dǎo)致平坦地區(qū)河網(wǎng)提取效果不佳，基于SRTM DEM數(shù)據(jù)進(jìn)行河網(wǎng)提取可以最大程度地保證提取的河網(wǎng)的準(zhǔn)確性。該文成功提取青龍河水域數(shù)字河網(wǎng)，這對青龍河流域的水資源規(guī)劃及配置提供技術(shù)條件并對實(shí)現(xiàn)水資源現(xiàn)代化管理有著重要意義。

參考文獻(xiàn)

[1] 梁科.基于SRTM數(shù)據(jù)的流域水系提取與三維可視化[D].長沙：中南大學(xué)，2008.

[2] 鄒鳴.基于DEM的數(shù)字河網(wǎng)提取及其應(yīng)用[J].水資源研究，2013（2）：14-15.

[3] Puecker TK， Doughlas DH.Detection of surface-specific points by local parallel processing of discrete terrain elevation data[J].Computer Graphics and Image Processing，1975（4）：375-387.

[4] Lawrence WM， Jurgen Garbrecht.Automated recognition of valley lines and drainage networks from grid digital elevation models：a review and a newmethod[J].Journal of hydrology，1995（167）：393-396.

[5] 葉愛中，夏軍，王綱勝，等.基于數(shù)字高程模型的河網(wǎng)提取及子流域生成[J].水利學(xué)報(bào)，2005（5）：531-537.

[6] 曹玲玲，張秋文.基于SRTM的數(shù)字河網(wǎng)提取及其應(yīng)用[J].人民長江，2007（8）：150-152.

[7] 吳冰，李昌文.基于SRTM DEM的盤龍河流域河網(wǎng)提取研究[J].水資源與水工程學(xué)報(bào)，2013（6）：157-162.

[8] 周國園，侯志華.基于DEM的汾河流域水系提取試驗(yàn)研究[J].科技創(chuàng)新與生產(chǎn)力，2014（11）：60-63.

[9] 沈中原，李占斌，李鵬，等.基于DEM的流域數(shù)字河網(wǎng)提取算法研究[J].水資源與水工程學(xué)報(bào)，2009（1）：20-23，28.

[10] 熱紫燕，阿不都沙拉木，穆艾塔爾·賽地，等.烏魯木齊河流域數(shù)字河網(wǎng)提取分析[J].新疆水利，2012（6）：11-15.

[11] Koch A，Heipke C.Quality assessment of digital surfacz models derived from the Shuttle Radar Topography Mission（SRTM）[C].Sydney：Proc.of the IEEE 2001 International Geoscience and Remote Sensing Symposium，2001.