基于DEM的古坊河流域信息提取及對比分析

金　晶，彭志榮，汪紅洲

(1.江西水利職業(yè)學院，江西 南昌 330013；2.江西省水利規(guī)劃設計研究院，江西 南昌 330029)

基于DEM的古坊河流域信息提取及對比分析

金晶1，彭志榮1，汪紅洲2

(1.江西水利職業(yè)學院，江西 南昌 330013；2.江西省水利規(guī)劃設計研究院，江西 南昌 330029)

摘要：以會昌縣古坊河流域為研究對象，基于數(shù)字高程模型DEM，通過對DEM的預處理、水流方向的提取、匯流累積量的提取、河網(wǎng)水系的提取等確定流域信息。將提取得到的結果與在CAD中人工判讀的結果進行比對分析。結果表明：流域面積與比降基本吻合，河長存在較大誤差，在今后的工作中可以直接利用GIS提取山洪溝流域面積和比降以減少實際的工作強度。

關鍵詞：DEM；流域信息；河網(wǎng)水系；GIS

古坊河系貢水二級支流、湘水一級支流，發(fā)源于江西省贛州市會昌縣安下村，自東北向西南流經(jīng)謝屋村、鄒屋村、禾塘村、小壩村，于狗脛村匯入湘水。古坊河流域范圍內(nèi)有兩座并聯(lián)水庫，分別是東鄉(xiāng)水庫和冬瓜坑水庫。東鄉(xiāng)水庫與冬瓜坑水庫之間平時補充灌溉用水，汛期兩水庫下泄洪量分別沿支流匯入古坊河。會昌縣古坊河由于受洪水的影響，流域內(nèi)溪溝山洪災害頻繁發(fā)生，直接威脅沿河兩岸的村莊、道路等安全，給集鎮(zhèn)人民生活帶來許多困難、工農(nóng)業(yè)生產(chǎn)帶來巨大損失，嚴重制約著區(qū)內(nèi)國民經(jīng)濟的發(fā)展和人民生活水平的進一步提高。因此，需要及時對其進行治理，建立起山洪災害防治體系。本文探討了對于山洪溝流域直接用DEM圖提取流域信息與人工判讀得出的結果的誤差，為山洪溝流域信息提取提供依據(jù)。

1流域地形處理與流域信息提取

1.1DEM的預處理

首先，根據(jù)古坊河的坐標和調(diào)查得出的古坊河流域內(nèi)的居民地以及古坊河水系的情況，經(jīng)過分析，在GIS中創(chuàng)建SHP圖形用于裁剪，然后在ERDAS中裁出能包括古坊河流域的邊界光滑的DEM圖。接著，對DEM圖進行填洼和平坦地區(qū)處理。之所以在DEM圖中出現(xiàn)洼地是因為真實的地形或DEM生成過程中的誤差導致的。從目前的研究情況來看，它形成的原因一般認為是因為誤差而導致的偽地形。水流流向的判讀會因為洼地的存在而帶來影響，可能會阻礙后續(xù)水流線的生成，使我們得到一些錯誤的信息。因此在ArcGIS中計算水流方向之前，需要對數(shù)字高程模型進行填洼處理。常用的填洼方法是0’callaghan和Mark提出的，直接在原始DEM中修正洼地內(nèi)的單元格的高程值[1]。填洼的原理是：如果遇到洼地，比較該柵格相鄰8個柵格的高程，用最小的值賦予這個柵格，對峰頂?shù)奶幚硪彩穷愃频姆椒╗2]。

運用ERDAS軟件中的Subset，裁剪出原始DEM圖見圖1。再利用ArcGIS水文分析中的填洼對包含古坊河流域的原始DEM圖進行填洼處理見圖2。

圖1　分辨率30 m的裁剪后的原始DEM圖

圖2　填洼處理后的DEM圖

1.2水流方向的提取

基于DEM水文分析的基礎是水流方向的提取。在目前的研究中，最廣泛采用的是單流向法——D8法[3]。它是基于3×3的網(wǎng)格窗口，假定每一個網(wǎng)格中的水只能流向與之相鄰的8個網(wǎng)格中的一個，從而對水流形態(tài)進行判別，如圖3所示。該方法是建立在最陡坡度的基礎上，認為網(wǎng)格以點源的形式產(chǎn)流，將千萬種可能的流向歸納成八種并用一維線來描述河道，從而確定網(wǎng)格流向[4]。這種方法簡單并且具有實用性從而被廣泛接受和使用。

圖3　D8法網(wǎng)格水流方向定義

將填洼處理后的DEM圖導入ArcGIS中，利用流向工具計算得出水流方向，結果見圖4。

圖4　流向工具處理后的水流方向圖

1.3匯流累計量的提取

匯流累積分析是指通過提取水流流向柵格，從而生成匯流累積柵格的過程。匯流是一個不斷累積的過程，每一個匯流累積柵格的值，它表示上游匯流區(qū)域內(nèi)所有匯入該柵格的柵格單元總數(shù)[5]。簡而言之，就是假定最初每個柵格都代表一個單位水量，根據(jù)提取的水流流向，該點所通過的水量就是匯入該柵格的總個數(shù)。這是基于水流流向數(shù)據(jù)計算所得。其值越大，就表示該地區(qū)越容易產(chǎn)生地表徑流。

本研究利用水文分析工具集中的流量工具提取匯流累積量柵格，結果如圖5所示。

圖5　匯流累計量圖

1.4河網(wǎng)水系的提取

基于DEM開展水文分析，最主要的目的就是提取出流域的河流網(wǎng)絡。目前，在所有提取河網(wǎng)水系的方法中，最普遍的方法是地表徑流模型。先前我們已經(jīng)得到了流域的匯流累積量，也就是在水流方向上每一個柵格所流過的水量值。在匯流量超過一個臨界值時，地表徑流就會產(chǎn)生，該值被稱為積水面積閾值[6]。流域內(nèi)所有柵格與臨界值進行比較，超過臨界值的柵格所形成的路徑就是水流的路徑，即河網(wǎng)。然后，河網(wǎng)矢量化成河網(wǎng)形文件(即shape文件)輸出就可以將河網(wǎng)提取出來了。提取河網(wǎng)后，可以利用GIS的幾何計算的功能計算出河流的長度。

本研究利用ArcGIS的柵格計算功能，提取匯流量大于800的柵格單元，得到的柵格河網(wǎng)，再利用水文分析中的柵格河網(wǎng)矢量化功能得到河網(wǎng)水系圖，如圖6所示。

圖6　流域河網(wǎng)水系圖

1.5古坊河流域的確定

用于研究的DEM圖是根據(jù)分析古坊河河源河口的經(jīng)緯度坐標裁剪出來的。可以根據(jù)調(diào)查得出的古坊河流域內(nèi)的居民地以及古坊河水系的情況來確定古坊河的流域。

流域盆地是指被分水嶺分割而形成的匯水區(qū)域。它是基于水流方向數(shù)據(jù)，經(jīng)過分析確定出互相連接的同一流域盆地上的柵格。流域盆地的出水口在流域的邊界上，所以流域盆地可以通過分析流域邊界上的出水口位置來確定。確定之后，我們可以通過計算匯入該出水口的上游的柵格位置，最終確定流域盆地的集水區(qū)域。

以柵格為基礎，利用水文分析工具中的盆域分析，得出流域的柵格數(shù)據(jù)。再運用轉(zhuǎn)換工具柵格轉(zhuǎn)面，將流域柵格圖層轉(zhuǎn)換成矢量圖層。然后利用GIS的計算幾何功能計算出每個流域的流域面積，提取出來的古坊河流域。

1.6子流域的提取

操作ArcGIS時，應用河網(wǎng)連接工具對提取的柵格河網(wǎng)進行分析，將分析的結果作為小流域的出水口數(shù)據(jù)。這是因為通過河網(wǎng)連接分析得到的數(shù)據(jù)內(nèi)含有每一段河網(wǎng)弧段的連接信息。這些信息中，包含弧段的起點、終點等數(shù)據(jù)。一般情況下，小流域的出水口位置就在弧段的終點處。

子流域集水區(qū)域的生成原理與大流域的是一樣的。小流域的出水口位置被確定后，就等于確定了該集水區(qū)域的最低點，與水流方向結合在一起分析，該出水口上游所有流入該點的柵格位置就能被確定。這些柵格就組成了分割后的子流域區(qū)域。在ArcGIS中，應用水文分析工具中的河流連接，先利用計算得到的河網(wǎng)文件和流向柵格文件建立河流連接[7]，再結合現(xiàn)有的流向柵格，利用分水嶺工具提取出子流域。這之后，將子流域柵格文件轉(zhuǎn)換成矢量文件，并對其進行計算就可以得出子流域的面積，古坊河子流域圖見圖7。

圖7　古坊河子流域范圍圖

1.7河道平均比降的確定

河道平均比降J是指自分水嶺起根據(jù)沿流程的比降變化轉(zhuǎn)折點的高程。在量取河長的同時確定過河等高線，量取各等高線間距和相應高程，按下列計算公式計算平均比降[8]。

(1)

式中：hn為各轉(zhuǎn)折點至下游控制斷面的高程差，m；ln為各轉(zhuǎn)折點之間的距離，km；L為河長，km；J有效數(shù)字取三位，‰。

通過古坊河流域DEM圖，根據(jù)空間分析模塊中的等值線分析，提取出流域的等高線，見圖8。根據(jù)等高線，讀取若干比降變化特征點的高程和間距，并按照公式(1)計算出古坊河流域的平均比降。

圖8　古坊河流域等高線示意圖

2結果分析

基于DEM圖利用GIS的水文分析功能提取出的流域信息與根據(jù)現(xiàn)有的地形圖在CAD中勾繪得出的流域信息進行對比的結果見表1。

對于ArcGIS與CAD計算結果的相對偏差的絕對值認為在2%以內(nèi)是合理的。根據(jù)表1可知，基于DEM利用GIS計算得出的流域面積與根據(jù)地形圖在CAD中人工勾繪出的相對偏差為0.82%，計算結果很接近。因此，對于山洪溝所在的中小流域，運用GIS水文分析功能提取出來的流域面積是合理的。在GIS中對應CAD中的主河道的比降相對偏差的絕對值為1.75%，我們也認為是可信的。但是，GIS提取出的水系對應CAD的主河道的長度偏差較大，不適合使用。造成這種結果是由于以下幾點原因造成的：第一，它們的數(shù)據(jù)源不同，這里使用的CAD中的地形圖是根據(jù)現(xiàn)場測量得出的，而使用的DEM圖是航拍得到的。第二，數(shù)據(jù)的時效性，這里使用的DEM圖是比較早期生成的，古坊河流域內(nèi)后來新建兩座水庫等一系列的工程措施，這引起了流域的變化。第三，圖像的比例尺的不同，CAD中的地形圖是1∶1萬，而DEM圖是1∶25萬的，這引起了平面位置精度誤差和垂直高度精度誤差。

表1　流域信息對比表

3結語

綜上所述，本文利用EARDS根據(jù)測量的經(jīng)緯度裁剪DEM后，用GIS對DEM圖進行水文模擬分析，將模擬得出的流域信息與CAD中根據(jù)地形圖人工判讀得出的流域信息進行對比分析，結果得出流域面積和比降基本吻合，但河長存在差異。故對于山洪溝流域可以基于DEM圖根據(jù)GIS進行水文模擬分析得出流域面積以減少工作量。由于數(shù)據(jù)源不同和數(shù)據(jù)的時效性，河長存在較大的差異。雖然比降基本吻合，但是運用兩種軟件計算的方法是一樣的，不會明顯減少工作量。因此，對于山洪溝流域的研究中，對于無地形資料的山洪溝流域水文計算應用GIS提取出流域信息是不合適的。對于有地形資料的較復雜的山洪溝流域，可以應用GIS進行水文模擬分析計算得出流域面積和比降，再通過CAD人工判讀出河長。這樣對于今后的山洪溝研究工作，可以在一定程度上減少工作量，提高工作效率。

參考文獻：

[1]Callaghan F,Mark D M.The extraction of drainage networks from digital elevation data[J].Computer Vision，Graphics and Images Processing，1984，28(3):323-344.

[2]林金輝.基于GIS的數(shù)字化渭河流域地理特性研究[D].西安：西安建筑科技大學，2007.

[3]Garbrecht J,Martz L W.The assignment of drainage direction over flat surfaces in raster digital elevation models[J].Journal of Hydrology，1997,5(193):204-213.

[4]劉杰.基于GIS的水庫產(chǎn)匯流模擬研究[D].南京：河海大學，2013.

[5]陳健飛,張筱林.地理信息系統(tǒng)導論[M].北京:科學出版社，2010.

[6]王玉富,王翰釗.ArcGIS 環(huán)境下基于 DEM 的流域特征提取[J].湖北民族學院學報，2010，28(4):439-442.

[7]劉剛,趙榮.基于DEM的瀾滄江流域水文信息提取方法的研究[J].地理信息世界，2007(1):56-59.

[8]戴榮,王正發(fā).資料匱乏地區(qū)中小流域設計洪水計算方法的研究[J].西北水電，2012(2):1-5.

Gufang River Basin information extraction and analysis based on DEM

JIN Jing1, PENG Zhirong1, WANG Hongzhou2

(1.Jiangxiwaterconservancyvocationalcollege,Nanchang330013,China;2.Jiangxiprovincialwaterconservancyplanningdesignandresearchinstitute,Nanchang330029,China)

Abstract：Take the Gufang River of Huichang county as a research subject, through the pretreatment of DEM, information extraction of water direction and river network to confirm basin information. The results extracted were compared and analyzed against the results of artificial read in CAD. The results showed that: the drainage area basically agrees well with hydraulic slope, but?the length of river exist notable error. We can extract catchment area and slope of mountain gullies directly by using GIS in the future work to reduce the actual work intensity.

Key words：DEM;watershed information; river network; GIS

作者簡介：金晶(1991-)，女，助理講師，主要從事水文水資源方面的研究。

中圖分類號：P333.9;TP311.52

文獻標志碼：A

文章編號：2096-0506(2016)05-0040-04