吳強(qiáng)建中+徐珍+蔣婷+賀春明+黎武

摘 要：該文應(yīng)用水文分析方法提取南充境內(nèi)河流網(wǎng)絡(luò)。采用ARCGIS10.0作為軟件平臺(tái)，基于南充地區(qū)影像圖得到無洼地DEM。從無洼地DEM提取出水流方向、匯流累積量，最終得到矢量河網(wǎng)。實(shí)驗(yàn)發(fā)現(xiàn)閾值的設(shè)置對(duì)提取到的河網(wǎng)起決定性作用。

關(guān)鍵詞：DEM；南充市；河網(wǎng)提取

中圖分類號(hào) P333 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 1007-7731（2017）11-0033-04

River Network Extraction in Nangchong Based on DEM

Wu Qiangjianzhong et al.

（School of Land and Resources，China West Normal University，Nanchong 637009，China）

Abstract：This paper applicated the method of hydrology for extraction of river network in Nanchong County of Sichuan based on DEM.Using ARCGIS10 as software platform，based on satellite imagery received non_depression DEM.Extract direction of flow，flow accumulation from non_depression DEM and finally received river network.The setting of threshold plays a decisive role for extracting river network.

Key words：DEM；Nanchong City；River network extraction

DEM是描述地球表面形態(tài)的多種信息空間分布的有序數(shù)字陣列，是地球表面地形地貌的一種離散表達(dá)，現(xiàn)已經(jīng)成為國(guó)家空間數(shù)據(jù)基礎(chǔ)設(shè)施的“4D”產(chǎn)品之一，是水文分析的主要數(shù)據(jù)[1]。利用DEM可以快捷的得到數(shù)字流域信息[2]。數(shù)字流域信息能夠反映河流分布的規(guī)律并能模擬河流分布的形態(tài)，廣泛應(yīng)用于水文預(yù)測(cè)、城市規(guī)劃等領(lǐng)域，是地理數(shù)據(jù)的核心內(nèi)容之一，也是分布式水文模型的基礎(chǔ)數(shù)據(jù)[3]。數(shù)字流域信息對(duì)地區(qū)規(guī)劃設(shè)計(jì)，洪澇災(zāi)害監(jiān)測(cè)都有非常重要的意義。因此，基于DEM的流域分析已經(jīng)非常常見[4-8]，例如，郭超穎研究不同水平分辨率對(duì)流域水文特征的影響發(fā)現(xiàn)，相同閾值下分辨率越低，地形及水流就越簡(jiǎn)單，相反則越詳細(xì)[9]；李霜林等借助ArcGIS平臺(tái)，在DEM上進(jìn)行流域分析，提取了長(zhǎng)江三峽地區(qū)地形的河網(wǎng)、集水區(qū)域等重要水文信息，方便了三峽地區(qū)水資源利用，地質(zhì)災(zāi)害的監(jiān)測(cè)等[10]。本文基于DEM地表水文分析的方法，運(yùn)用ArcGIS10.0作為軟件平臺(tái)，提取了南充段嘉陵江流域的水流方向、匯流累積量等河流信息，最終得到矢量河網(wǎng)。

1 研究區(qū)域概況

南充市地處四川盆地東北部，介于30°35′N～31°51′N、105°27′E～106°58′E；南北跨度165km（圖1），東西跨度143km。東鄰達(dá)州市，南連廣安市，西與遂寧市、綿陽市接壤，北與廣元市、巴中市毗鄰，幅員面積12479.96km2。嘉陵江流經(jīng)南充市的閬中、南部、儀隴、蓬安4縣市和高坪、順慶、嘉陵3區(qū)，全長(zhǎng)298km，沿途接納了東河、構(gòu)溪河、西河、西充河等支流，流域面積9753km2，占南充市總面積的78%。北與廣元、巴中市毗鄰，南與廣安市相接，東臨達(dá)州市，西接遂寧、綿陽市。南充水系比較發(fā)達(dá)，除嘉陵江外南充還有20余條支小河流和24座大中小型通航水庫(kù)，較大的支流有東河、構(gòu)溪河、西河、西充河等，與嘉陵江交匯形成樹枝狀水系。嘉陵江水系河道蜿蜒曲折，多曲流，河道彎曲系數(shù)一般在2以上，在南充青居牛肚壩形成了堪稱“世界地貌奇觀”的大曲流，河道回旋接近360°；嘉陵江流量大，季節(jié)分配不均，年際變化大，夏秋多，冬春少；嘉陵江含沙量大，以細(xì)沙粒物質(zhì)為主。

2 數(shù)據(jù)與方法

2.1 數(shù)據(jù) 本研究的數(shù)據(jù)主要為南充及周邊地區(qū)的影像地圖。為了契合小流域范圍內(nèi)的水域分析，選擇下載有符號(hào)整型像素深度16位的單波段TIFF格式影像數(shù)據(jù)，投影類型D_WGS_1984_UTM_Zone_48N，線性單位為：meter，分辨率為30m×30m。實(shí)驗(yàn)中下載的DEM數(shù)據(jù)共有4幅。將四幅影像圖在ArcGIS中合并，設(shè)置合并后的數(shù)據(jù)基本信息為有符號(hào)整型、像素深度16位、波段數(shù)為1。運(yùn)用按掩膜提取工具去除南充轄區(qū)以外的區(qū)域，得到只包含南充地區(qū)的DEM數(shù)據(jù)。圖幅范圍（上：35421199.13288，左：499973.201333，右：692964.201333，下：3318729.13288）；灰度值范圍：43～183；可以看出圖中灰度值較高的區(qū)域?yàn)樗w，以河流為主。

2.2 研究方法 地表河流網(wǎng)提取是DEM水文分析的主要內(nèi)容之一，目前河網(wǎng)提取方法主要采用地表徑流漫流模型：首先在無洼地DEM上擬用最大坡降法計(jì)算每一個(gè)柵格的水流方向，然后依據(jù)自然水流由高處流往低處的自然規(guī)律，計(jì)算出每一個(gè)柵格在水流方向上累積的柵格數(shù)，既匯流累積量。假設(shè)每個(gè)柵格攜帶一份水流，那么柵格的匯流累積量就代表該柵格的水流量?；谶@種思維，當(dāng)匯流量達(dá)到一定值就會(huì)產(chǎn)生地表水流，所有匯流量大于臨界值的柵格就存在潛在的水流路徑，這些水流路徑構(gòu)成的河流網(wǎng)絡(luò)就稱為河網(wǎng)。本實(shí)驗(yàn)從原始DEM數(shù)據(jù)處理，最后得到矢量河網(wǎng)，主要可以分為提取無洼地DEM的水流方向、計(jì)算匯流累積量、設(shè)置臨界積水面積閾值生成河網(wǎng)、河網(wǎng)矢量化等幾個(gè)步驟[11]。具體流程如圖2所示。

DEM是一種地形表面模型，其中包含有真實(shí)洼地和假的洼地既由數(shù)據(jù)精度不夠高所造成的水流積聚地。洼地填平的主要功能是排除由于DEM的精度不夠高所產(chǎn)生的假的洼地。這些區(qū)域的存在會(huì)讓分析得不到我們想要的效果甚至出現(xiàn)錯(cuò)誤的水流方向，因此做水文分析之前需要對(duì)原始DEM數(shù)據(jù)進(jìn)行洼地填充生成無洼地DEM。具體步驟如下：

2.2.1 提取原始DEM的水流方向 對(duì)于DEM中的某一個(gè)格網(wǎng)，水流方向表征的是該格網(wǎng)表面的降水受重力作用在地形表面約束下的離開該格網(wǎng)的方向。在ArcGIS中通過對(duì)中心柵格的8個(gè)鄰域柵格進(jìn)行編碼，中心柵格的某一方向便可由其中的某一值來確定。例如中心柵格的水流方向向右，這水流方向被賦值為1，方向值以2的冪值指定是因?yàn)榇嬖诟窬W(wǎng)水流方向不確定的情況，此時(shí)需要將數(shù)個(gè)方向值相加，這樣在后續(xù)處理中從相加結(jié)果便可以確定相加時(shí)中心格網(wǎng)的領(lǐng)域格網(wǎng)的情況如圖3所示。ArcGIS的水文分析以D8算法為基礎(chǔ)的。D8算法的基本思想是假設(shè)單個(gè)柵格中的水流只能流入與之相鄰的8個(gè)柵格中。采用最陡坡度法來確定水流的方向，在3×3的DEM柵格上，計(jì)算中心柵格單元與各相鄰的8個(gè)柵格單元間的距離權(quán)落差（即柵格中心點(diǎn)落差除以柵格中心點(diǎn)之間的距離），取距離權(quán)落差最大的柵格為中心柵格的流出柵格既水流方向，并且規(guī)定一個(gè)柵格的水流方向用一個(gè)特征碼表示。被處理柵格單元同相鄰8個(gè)柵格單元之間坡降的算法為：

式（1）中：FlowDiretion為2個(gè)柵格之間的距離權(quán)落差；Dz為2個(gè)柵格單元之間的高程差；Di為2個(gè)柵格單元中心之間的距離。格網(wǎng)間距離與有方向有關(guān)，若領(lǐng)域柵格對(duì)中心柵格的方向值為2、8、16、128，則柵格間的距離為[2]倍的柵格大小否則距離為1。該步驟主要用流向提取工具，柵格數(shù)據(jù)為原始DEM，得到水流方向柵格flowdir如圖4所示。

2.2.2 洼地計(jì)算 雙擊Hydrology工具集中的Sink工具，選擇計(jì)算的水流方向數(shù)據(jù)Flowdir作為輸入數(shù)據(jù)，指定輸出數(shù)據(jù)為Sink，得到洼地?cái)?shù)據(jù)sink。

2.2.3 洼地填充 計(jì)算得知原始DEM上有洼地，需要進(jìn)行洼地填充。雙擊Hydrology工具集中的Fill工具，進(jìn)行原始DEM挖地填充。在Input surface raster中選擇原始DEM數(shù)據(jù)，指定輸出文件名為filldem，得到無洼地DEM如圖5所示。

2.2.4 計(jì)算無洼地DEM的匯流累積量 使用Flow Accumulation工具，柵格數(shù)據(jù)就是上面獲得的數(shù)據(jù)flowdirfill，得到無洼地的DEM數(shù)據(jù)flowacc。無洼地DEM水流方向見圖6，匯流累積量見圖7。

2.2.5 河網(wǎng)的生成 河網(wǎng)生成的思路是根據(jù)需要設(shè)定恰當(dāng)?shù)拈撝?，依?jù)設(shè)定的閾值逐個(gè)對(duì)柵格單元進(jìn)行判斷，若單元格匯流量大于閾值其屬性為1，小于或等于閾值的柵格設(shè)置為無數(shù)據(jù)（NoDate）。具體操作使用spatial analyst中的柵格計(jì)算器，將匯流累積量計(jì)算結(jié)果作為輸入數(shù)據(jù)，提取公式為“Con（"flowacc">100，1）”。設(shè)置的閾值是100、500、2000、5000，得柵格河網(wǎng)str100、str500、str2000、str5000圖。在圖中像元值只有0、1，其中1表示是河網(wǎng)，0是非河網(wǎng)區(qū)域。

2.2.6 柵格河網(wǎng)矢量化 使用Stream To Feature工具，將是柵格河網(wǎng)作為輸入數(shù)據(jù)，將flowdirfill作為水流方向，得到矢量化后的河網(wǎng)。

3 結(jié)果與分析

本文運(yùn)用ArcGIS10.0軟件平臺(tái)對(duì)南充流域的DEM數(shù)據(jù)進(jìn)行了一系列水文分析，設(shè)置不同的閾值提取得到了不同的矢量河網(wǎng)?，F(xiàn)對(duì)比分析閾值分別為5000、2000、500、100……時(shí)提取的河流局部圖。圖7是相鄰的4個(gè)區(qū)域不同閾值提取得到的河網(wǎng)對(duì)比圖；圖8是同一區(qū)域不同閾值提取的河網(wǎng)疊加在一起形成的對(duì)比圖。對(duì)比發(fā)現(xiàn)設(shè)置閾值越小河流越詳細(xì)，設(shè)置閾值越大河流信息越簡(jiǎn)單，但在某些河流段主要是一些較大的主干河流，閾值為500與閾值為100時(shí)提取的河網(wǎng)完全相同；既閾值較大的河網(wǎng)與閾值較小的河網(wǎng)的主干部分完全重合。在河流的主干部分閾值越大得到的河網(wǎng)越理想，但會(huì)讓河流的部分分支尤其是在源頭處的細(xì)小河流被忽略；在河流的分支尤其是河流的源頭閾值越小得到的河流越理想，但閾值較小會(huì)在河流的主干較寬的部分產(chǎn)生很多平行的偽河道[12]。

在河網(wǎng)的屬性表中記錄了有關(guān)河流的信息，現(xiàn)按照閾值從大到小排序，如表1所示。從表1可以看出，隨著閾值的變小，河流總條數(shù)急劇增加，其原因是小河溝也加入了進(jìn)來，導(dǎo)致總數(shù)激增。隨著閾值的變小，最長(zhǎng)河流在變短，其原因是閾值變小的過程中，小河流的加入將最長(zhǎng)的那條河流分割成兩條甚至更多，所以最長(zhǎng)河流在發(fā)生變化。隨著閾值的變小，河流的總長(zhǎng)度也在激增，同樣也是由于在閾值變小的過程中，加入了小河流的緣故。

4 結(jié)論與討論

本文主要基于DEM地表水文分析的方法，對(duì)南充段嘉陵江流域進(jìn)行了研究。采用ARCGIS10.0作為軟件平臺(tái)，提取了研究區(qū)域矢量河網(wǎng)，實(shí)驗(yàn)發(fā)現(xiàn)預(yù)知的設(shè)置對(duì)河網(wǎng)起至關(guān)重要的作用。閾值設(shè)置越小河流越詳細(xì)，閾值設(shè)置越大河流越稀疏，且閾值較大時(shí)得到的河網(wǎng)與閾值較小時(shí)得到的河網(wǎng)的主干部分完全重合。隨著閾值的變小，許多小河流被納入河網(wǎng)中，使得河流總條數(shù)急劇增加，河流的總長(zhǎng)度也在激增。

在整個(gè)研究過程中，由于閾值的設(shè)定不是很科學(xué)，如果采用科學(xué)梯度法可能得到更為理想的河網(wǎng)柵格數(shù)據(jù)。閾值的設(shè)置對(duì)提取的河網(wǎng)起決定性作用，閾值過大會(huì)造成關(guān)鍵河流的丟失，閾值過小會(huì)產(chǎn)生很多多余的數(shù)據(jù)，都得不到我們想要的數(shù)據(jù)。怎樣確定合理的閾值？閾值為多少時(shí)得到的河網(wǎng)與現(xiàn)實(shí)中河流最為接近？這些都是以后的研究需要提高和完善的，該方面的研究還需要更高分辨率的數(shù)據(jù)，更加科學(xué)和專業(yè)的研究方法。

參考文獻(xiàn)

[1]李志林.數(shù)字地形分析[J].地理學(xué)報(bào)，2006，61（12）：1326-1326.

[2]王云，梁明，汪桂生.基于ArcGIS的流域水文特征分析[J].西安科技大學(xué)學(xué)報(bào)，2012（5）：581-585.

[3]王敏嬙.基于DEM的水文要素提取及水文拓?fù)浣Y(jié)構(gòu)建立的研究[D].西安：西北大學(xué)，2012（2）：64-71.

[4]宋向陽，吳發(fā)啟，趙龍山，等.基于DEM的延河流域水文特征提取與分析[J].干旱地區(qū)農(nóng)業(yè)研究，2012（4）：200-206.

[5]原曉平.基于DEM的格爾木河流域數(shù)字地形分析[D].北京：中國(guó)地質(zhì)大學(xué)（北京），2013（2）：56-75.

[6]朱海玲，楊曉暉，張學(xué)培，等.基于DEM的密云水庫(kù)上游流域特征提取與分析[J].中國(guó)水土保持科學(xué)，2013（3）：66-72.

[7]常直楊，王建，白世彪，等.基于DEM的岷江上游流域構(gòu)造活動(dòng)強(qiáng)度分析[J].地球信息科學(xué)學(xué)報(bào)，2014（4）：568-574.

[8]常直楊，王建，白世彪，等.基于DEM的白龍江流域構(gòu)造活動(dòng)定量分析[J].第四紀(jì)研究，2014（2）：292-301.

[9]郭超穎.基于DEM的馬蓮河流域數(shù)字地形分析[D].北京：北京林業(yè)大學(xué)，2008.

[10]李霜林，楊武年，梁英杰.基于DEM的長(zhǎng)江三峽地區(qū)流域分析[J].測(cè)繪，2012（4）：147-150.

[11]湯國(guó)安，楊昕.地理信息系統(tǒng)空間分析實(shí)驗(yàn)教程[M].北京：科學(xué)出版社，2006.

[12]黃付友，賴格英，許為一.基于ArcHydro Tools對(duì)流域特征自動(dòng)提取的若干問題[J].許昌學(xué)院學(xué)報(bào)，2008（2）：43-46.

（責(zé)編：張宏民）