李俊超，馬 倩，陶 鈞

（1.福建省測繪院，福建福州350000；2.內(nèi)蒙古自治區(qū)測繪院，內(nèi)蒙古呼和浩特010051；3.武進(jìn)規(guī)劃與測繪院，江蘇常州213000）

基于ArcGIS的水文流域分析及應(yīng)用

李俊超1，馬 倩2，陶 鈞3

（1.福建省測繪院，福建福州350000；2.內(nèi)蒙古自治區(qū)測繪院，內(nèi)蒙古呼和浩特010051；3.武進(jìn)規(guī)劃與測繪院，江蘇常州213000）

在理解基于DEM數(shù)據(jù)進(jìn)行水文分析的基本原理的基礎(chǔ)上，利用ArcGIS提供的水文分析工具對DEM數(shù)據(jù)水文流域進(jìn)行了分析和研究，并將其運(yùn)用于實(shí)際生產(chǎn)生活中，為險(xiǎn)情預(yù)報(bào)、土質(zhì)研究、降水分析等提供了幫助。

ArcGIS；水文分析；DEM

傳統(tǒng)的地形圖是將地面上的信息 (地貌、地物以及各種名稱)用圖形和注記的方式表示在圖紙上，其優(yōu)點(diǎn)在于直觀，便于攜帶與使用[1]。然而，隨著計(jì)算機(jī)技術(shù)和信息處理技術(shù)的飛速發(fā)展，紙質(zhì)地圖不能被計(jì)算機(jī)直接利用，無法滿足各種工程設(shè)計(jì)自動(dòng)化的研究要求。因此，地圖的數(shù)字化產(chǎn)品逐步得到開發(fā)應(yīng)用。數(shù)字高程模型（DEM）則是一種典型的數(shù)字化產(chǎn)品，是流域地形、地物識別的重要原始資料，具有廣泛的實(shí)際應(yīng)用價(jià)值：在測繪生產(chǎn)中可用于繪制等高線、坡度，制作立體景觀圖等；在工程項(xiàng)目中可用于體積、面積的計(jì)算，線路的設(shè)計(jì)等；在環(huán)境與規(guī)劃中可用于各種規(guī)劃、洪水險(xiǎn)情預(yù)報(bào)、水文研究等[2]。自20世紀(jì)60年代以來，在利用DEM提取流域水文特征，模擬地表水文過程等方面國內(nèi)外都開展了大量的研究[3]。本文基于ArcGIS軟件中水文分析工具(HydrologyModeling)，對水文流域分析進(jìn)行研究和比較，能起到對水文流域?yàn)?zāi)害的發(fā)生、發(fā)展 (時(shí)間、空間)并對其流向進(jìn)行區(qū)劃、預(yù)測預(yù)報(bào)作用，從而實(shí)現(xiàn)對水文流域的分析。ArcGIS提供的水文分析模塊主要用來建立地表水的運(yùn)動(dòng)模型，輔助分析地表水流從哪里產(chǎn)生以及要流向何處，再現(xiàn)水流的流動(dòng)過程[4]。同時(shí)，利用水文分析工具也可以有助于了解排水系統(tǒng)和地表水流過程的一些基本概念和關(guān)鍵過程，以及從DEM數(shù)據(jù)上獲取更多的水文信息。

1 基于ArcGIS的水文分析過程

水文分析的主要步驟如圖1所示。

由于被較高高程區(qū)域圍繞的洼地是進(jìn)行水文分析的一大障礙，因此在確定水流方向以前，必須先將洼地填充，獲得無洼地的DEM[5]。有些洼地是DEM生成過程中的數(shù)據(jù)錯(cuò)誤帶來的，但另外一些卻表示了真實(shí)的地形如采石場或巖洞等。因此在進(jìn)行洼地填充之前，必須計(jì)算洼地深度，然后在進(jìn)行洼地填充的過程中，設(shè)置合理的填充閾值。我們可以通過填充洼地(Fill Sinks)得到無洼地的DEM，具體步驟文中將作介紹。

圖1 利用DEM進(jìn)行水文分析的主要步驟圖

2 水文分析

2.1 水文流向分析原理

流向分析的原理如圖2所示，共有8個(gè)方向，分別是2的 n次方。水流的流向是通過計(jì)算中心柵格與鄰域柵格的最大距離權(quán)落差來確定的。距離權(quán)落差是指中心柵格與鄰域柵格的高程差除以兩柵格間的距離，柵格間的距離與方向有關(guān)，如果鄰域柵格對中心柵格的方向值為2、8、32、128，則柵格間的距離為SQRT (2)≈1.414，否則距離為1。如果高程差為正值，則為流出；為負(fù)值，則為流入。

圖2 流向分析原理圖

2.2 原始資料數(shù)據(jù)的獲得

現(xiàn)以某地區(qū)1∶5萬DEM數(shù)據(jù)為例：在ArcMap中加載DEM數(shù)據(jù)，執(zhí)行工具條 [HydrologyModeling]中的菜單命令 [Hydrology]>>[Fill Sinks]，在出現(xiàn)的對話框中將 [InputSurface]參數(shù)指定為“DEM”，確定后得到無洼地的DEM數(shù)據(jù)：[Filled Sink1]（見圖3）。

圖3 無洼地的DEM數(shù)據(jù)

2.3 水文分析的具體步驟

1）流向分析：執(zhí)行工具條 [Hydrology Modeling]中的菜單命令 [Hydrology]>>[Flow Direction]，在出現(xiàn)的對話框中將 [Input Surface]參數(shù)指定為“Filled Sink1”；確定后得到流向柵格 [Flow Direction1]，流向柵格單元的數(shù)值表示的含義即是上述提到的流向分析原理中的值（見圖4）。

圖4 獲得流向柵格數(shù)據(jù)

2）計(jì)算流水累積量：執(zhí)行工具條[HydrologyModeling]中的菜單命令[Hydrology]>>[Flow Accumulation]，在出現(xiàn)的對話框中將 [Direction Raster]參數(shù)指定為“Flow Direction1”。

3）提取河流網(wǎng)絡(luò)：利用地圖代數(shù)來進(jìn)行河網(wǎng)提取[6]。目前常用的河網(wǎng)提取方法是采用地表徑流漫流模型計(jì)算：首先是在無洼地DEM上利用最大坡降的方法得到每一個(gè)柵格水流方向，然后利用水流方向計(jì)算出每一個(gè)柵格在水流方向上累積的柵格數(shù)，即匯流累積量，所得到匯流累積量則代表在一個(gè)柵格位置上有多少個(gè)柵格的水流方向流經(jīng)該柵格。

4）提取河流網(wǎng)絡(luò)柵格：打開 Arctoolbox，運(yùn)行工具 [SpatialAnalystTools]>>[Map Algebra]>>[單輸出地圖代數(shù)]。在 [地圖代數(shù)表達(dá)式]中輸入公式：con (Flow Accumulation1>800,1)[輸出柵格]指定為：Stream-Net。說明：通過此操作將流水累積量柵格 [Flow Accumulation1]中柵格單元值（流水累積量）大于800的柵格賦值為1，從而得到河流網(wǎng)絡(luò)柵格 [StreamNet]。

5）提取河流網(wǎng)絡(luò)矢量數(shù)據(jù)：執(zhí)行工具條[Hydrology Modeling]中的菜單命令[Hydrology]>>[Stream Network As Feature]，在出現(xiàn)的對話框中將 [Direction Raster]參數(shù)指定為"Flow Direction1"，[AccumulationRaster]參數(shù)指定為"Flow Accumulation1"，[M inimum NumberofCells for a Stream]參數(shù)指定為1 000，確定后得到河流網(wǎng)絡(luò)矢量數(shù)據(jù)。

6）平滑處理河流網(wǎng)絡(luò)：打開 [編輯器]工具欄，執(zhí)行工具欄中的命令 [編輯器]>>[開始編輯]，確保目標(biāo)圖層為河流網(wǎng)絡(luò)圖層[Shape1],通過打開[Shape1]屬性表，并選擇屬性表的所有行，選擇圖層 [Shape1]中的所有要素，也可以通過要素選擇按鈕選擇圖層中所有要素執(zhí)行 [編輯器]工具欄中的命令 [編輯器]>> [更多的編輯工具]>>[高級編輯]打開工具條：[高級編輯]，點(diǎn)擊其上的 [平滑]按鈕。在 [平滑]處理對話框中輸入?yún)?shù)[允許最大偏移]：3；得到平滑后的河流網(wǎng)絡(luò)矢量圖層，執(zhí)行命令：[編輯器]>>[停止編輯]，保存所做修改，如圖5所示。

圖5 平滑處理后的數(shù)據(jù)

7）流域分析：執(zhí)行工具條 [Hydrology Modeling]中的菜單命 [Hydrology]>>[Watershed]，在出現(xiàn)的對話框中將[Direction Raster]參數(shù)指定為"Flow Direction1"，[AccumulationRaster]參數(shù)指定為"Flow Accumulation1"，[M inimum Number of Cells for a Stream]參數(shù)指定為2 000，確定后得到流域柵格 [Watershed1]。打開 [空間分析]工具欄，執(zhí)行命令:[空間分析]>>[轉(zhuǎn)換]>> [柵格到要素]，將流域柵格轉(zhuǎn)換成為矢量圖層，得到矢量數(shù)據(jù)：[WaterShed.shp]。設(shè)置圖層 [Watershed]屬性[符號]后得到類似圖6所示的效果。

圖6 流域柵格轉(zhuǎn)換為矢量圖層

3 水文分析的應(yīng)用——降水分析

執(zhí)行命令[Hydrology]>>[InteractiveProperties]選中“RainDrop”按鈕，在DEM上任意位置點(diǎn)擊（表示這里有降雨），則顯示一條路徑表示此處的降水的最終流向。圖面非常直觀、易理解，這給實(shí)際降水降雨提供了一個(gè)很好的平臺。降水流向分析見圖7。

圖7 降水流向分析

4 結(jié)語

基于ArcGIS的水文流域分析可以對水文流域?yàn)?zāi)害的發(fā)生、發(fā)展 (時(shí)間、空間)以及其流向進(jìn)行區(qū)劃、預(yù)測預(yù)報(bào)，對模型庫管理具有模型存取、查詢、維護(hù)、操作、新建模型等功能[7]。利用ArcGIS進(jìn)行水文流域分析提高了水文三維可視化的質(zhì)量和速度，能夠帶來巨大的社會效益[8]。但是，目前的 ArcGIS也不盡完美，當(dāng)用戶分析內(nèi)容過大時(shí)，將會帶來數(shù)據(jù)量的增加，系統(tǒng)的運(yùn)行將會變得遲緩，過大的模型導(dǎo)入到程序中進(jìn)一步操作也會出現(xiàn)類似的問題。為了更合理地體現(xiàn)出流域徑流，在DEM數(shù)據(jù)中提取水文信息的傳統(tǒng)方法基礎(chǔ)上，加入子流域界數(shù)據(jù)和河流位置數(shù)據(jù)，利用這 3種數(shù)據(jù)共同分析，將會達(dá)到更好的效果。這樣能確保在與實(shí)際情況相符的前提下，更加精確地反映出河道附近格網(wǎng)的水流方向。這也將是作者今后研究的方向。

[1] 程雄，王紅.GIS軟件應(yīng)用[M].武漢：武漢大學(xué)出版社，2004

[2] 劉志平，張素華，杜啟勝，等.基于ArcGIS的DEM生成方法及應(yīng)用[J].地理空間信息，2009,7（5）：69-71

[3] 張瑞軍.數(shù)字高程模型(DEM)的構(gòu)建及其應(yīng)用[J].工程勘察，2005，7(5)：61-64

[4] 蘇姝，李霖，劉慶華.基于DEM對都市流域的水文分析[J].武漢理工大學(xué)學(xué)報(bào)，2005，27（11）：59-62

[5] 郭立華，郭生練.基于DEM的數(shù)字河網(wǎng)生成方法的探討[J].長江科學(xué)院院報(bào),2003(4):14-17

[6] 程朋根，劉少華，王偉，等.三維地質(zhì)模型構(gòu)建方法的研究及應(yīng)用[J].吉林大學(xué)學(xué)報(bào)：地球科學(xué)版，2004，4(2)：309-313

[7] 何政偉，黃潤秋，許強(qiáng)，等.基于ArcGIS的地質(zhì)災(zāi)害防治信息與決策支持系統(tǒng)的研制[J].吉林大學(xué)學(xué)報(bào):地球科學(xué)版，2004，34（4）：601-606

[8] 侯海濱，侯海江.基于MapGIS的縣級土地利用數(shù)據(jù)庫的建設(shè)[J].測繪與空間地理信息，2010(1)：136-139

Analysis and Application of Hydrological Basin Based on ArcGIS

by LI Junchao

Understanding the basic principle ofhydrologic analysis based on the DEM data,we utilized hydrologic analysis tools provided by ArcGIS to analysis the data for practical life production.These tools can make contribution todanger forecasting,soil research,and precipitation analysis.

ArcGIS，hydrologic analysis，DEM

2012-02-02

P208

B

1672-4623(2012)06-0121-03

李俊超,研究方向?yàn)楣こ虦y量、城市空間信息。