黃土塬區(qū)DEM水文分析中消除地面?zhèn)螠瞎鹊姆椒?/h1>
2011-07-26 10:58:16趙明偉湯國(guó)安李發(fā)源
中國(guó)水土保持 2011年8期 關(guān)鍵詞：實(shí)驗(yàn)
趙明偉，湯國(guó)安，李發(fā)源
(南京師范大學(xué)地理科學(xué)學(xué)院，江蘇南京210046)
1 問(wèn)題的提出
溝壑密度是反映黃土高原地區(qū)地面受侵蝕程度重要的區(qū)域水土流失因子，在該地區(qū)的土壤侵蝕與侵蝕產(chǎn)沙研究中得到了廣泛的應(yīng)用［1－3］。目前，不同學(xué)者提出了不同的提取溝壑網(wǎng)絡(luò)的方法［4－7］，其中基于DEM的坡面流模擬方法因?yàn)閿?shù)據(jù)來(lái)源廣泛、方法簡(jiǎn)單、算法穩(wěn)定性高而得到普遍應(yīng)用。這一方法根據(jù)水文學(xué)坡面流模型來(lái)判別水流路徑，首先對(duì)地面洼地進(jìn)行填充，構(gòu)建無(wú)洼地的連續(xù)地面［8］，并根據(jù)每一個(gè)柵格單元與相鄰單元之間的最陡坡度，識(shí)別求取沿水流路徑的匯流累積量值，再選擇適當(dāng)?shù)膮R流累積量閾值來(lái)確定河網(wǎng)。該方法可以直接生成連續(xù)的河網(wǎng)。但是，在計(jì)算格網(wǎng)單元流量累積量值時(shí)，本不會(huì)產(chǎn)生地面徑流的平坦地區(qū)(例如黃土塬面等)也將產(chǎn)生很大的匯流累積量值，從而造成提取溝谷網(wǎng)絡(luò)時(shí)會(huì)有大量的偽溝谷出現(xiàn)［9－10］。以上問(wèn)題的出現(xiàn)，給黃土高原地區(qū)溝壑密度的正確提取以及水文分析的有效進(jìn)行，都帶來(lái)了很大的困難。
圖1為實(shí)驗(yàn)樣區(qū)水文分析中存在的偽溝谷示意圖，(a)為實(shí)驗(yàn)樣區(qū)的地形等高線圖，(b)為該地區(qū)采用坡面流方法提取的溝谷網(wǎng)絡(luò)，對(duì)比可以發(fā)現(xiàn)，原始提取的溝谷網(wǎng)絡(luò)中存在大量的偽溝谷。
圖1 偽溝谷示意
影響提取溝谷網(wǎng)絡(luò)結(jié)果的關(guān)鍵因素是匯流累積量閾值的選取［11］，一般說(shuō)來(lái)，當(dāng)該閾值變大時(shí)，所提取的溝谷網(wǎng)絡(luò)結(jié)果數(shù)量會(huì)變少。李俊、陳濤等以陜北黃土高原為研究樣區(qū)，研究了溝壑網(wǎng)絡(luò)提取中匯流累積量閾值的選取依據(jù)［12－13］。但是，對(duì)于黃土塬、黃土殘塬等地區(qū)，由于大量平地的存在，試圖僅僅通過(guò)調(diào)整匯流累積量閾值來(lái)消除偽溝谷是不可行的。圖2為DEM數(shù)據(jù)的山體陰影暈渲圖，可以看出A點(diǎn)處于溝谷中，而B(niǎo)點(diǎn)處于塬面平坦區(qū);對(duì)應(yīng)于水流路徑匯流累積量圖3，B點(diǎn)的匯流累積量值卻大大高于A點(diǎn)。由此可見(jiàn)，在黃土塬、黃土殘塬等存在大面積平坦區(qū)域的地區(qū)，不可能存在合適的閾值，使得提取的溝谷網(wǎng)絡(luò)既包含較小的溝谷，又不存在偽溝谷。
如上所述，在平原臺(tái)地地區(qū)(如黃土塬面)，由于大量平坦區(qū)域的存在，使得當(dāng)采用坡面流模擬方法進(jìn)行溝谷提取時(shí)會(huì)產(chǎn)生大量偽溝谷，且無(wú)法通過(guò)調(diào)整匯流累積量閾值來(lái)消除。因此，急需尋求一種簡(jiǎn)單而有效的解決方法。
圖2 山體陰影暈渲圖
圖3 匯流累積量圖
2 實(shí)驗(yàn)方法
2.1 實(shí)驗(yàn)樣區(qū)與數(shù)據(jù)
實(shí)驗(yàn)樣區(qū)位于陜西省淳化縣城西北部涇河中游地區(qū)，黃土塬及殘塬為主要地貌類(lèi)型，塬面地形平緩，但被諸多大溝谷深切割裂，平均海拔1045.04 m，相對(duì)高差390 m，塬面平均坡度為8.1°，面積100 km2(10 km×10 km)。此外，實(shí)驗(yàn)中為了對(duì)比文中提出的實(shí)驗(yàn)方法消除偽溝谷的效果，還選取了地貌類(lèi)型屬于黃土臺(tái)塬、黃土殘塬的樣區(qū)進(jìn)行附加實(shí)驗(yàn)。附加實(shí)驗(yàn)樣區(qū)的基本情況如表1所示。
表1 附加實(shí)驗(yàn)樣區(qū)概況
實(shí)驗(yàn)數(shù)據(jù)為國(guó)家測(cè)繪部門(mén)生產(chǎn)的5 m分辨率的DEM數(shù)據(jù)，數(shù)據(jù)精度滿(mǎn)足實(shí)驗(yàn)要求。
2.2 實(shí)驗(yàn)思路與方法
如前所述，基于DEM提取溝谷網(wǎng)絡(luò)是在得到水流路徑匯流累積量矩陣的基礎(chǔ)上，根據(jù)一定的閾值提取溝谷網(wǎng)絡(luò)。水流路徑匯流累積量矩陣的產(chǎn)生示意圖見(jiàn)圖4，可以看出，某點(diǎn)匯流累積量值的確定既受水流方向的控制，同時(shí)也與分析區(qū)域面積緊密相關(guān)。因此，在大面積的平原臺(tái)地地區(qū)，雖然地勢(shì)變化緩慢，但是仍然會(huì)產(chǎn)生很大的匯流累積量值，因而采用常規(guī)方法提取水系時(shí)會(huì)產(chǎn)生很多偽溝谷。
圖4 匯流累積量矩陣計(jì)算示意
根據(jù)以上分析，可以在計(jì)算匯流累積量矩陣時(shí)對(duì)每個(gè)柵格點(diǎn)增加權(quán)重值，使得平原臺(tái)地地區(qū)匯流累積量減小，從而減小偽溝谷出現(xiàn)的可能性。實(shí)驗(yàn)發(fā)現(xiàn)，在黃土高原地區(qū)，將地形表面劃分為正地形與負(fù)地形之后，溝谷網(wǎng)絡(luò)屬于負(fù)地形區(qū)域。因此，本實(shí)驗(yàn)根據(jù)地表面的正負(fù)地形屬性來(lái)設(shè)置權(quán)重值，正地形區(qū)域的匯流累積量權(quán)重值為0，負(fù)地形區(qū)域的匯流累積量權(quán)重值為1;將權(quán)重矩陣與原始匯流累積量矩陣做點(diǎn)乘運(yùn)算得到新的匯流累積量矩陣，再設(shè)置閾值提取溝谷網(wǎng)絡(luò)。綜合以上分析，在平原臺(tái)地地區(qū)基于DEM數(shù)據(jù)提取正確的溝谷網(wǎng)絡(luò)的步驟為:
(1)基于DEM數(shù)據(jù)得到匯流累積量矩陣［14］。首先，對(duì)原始DEM數(shù)據(jù)進(jìn)行洼地填充處理，然后，利用GIS水文分析，得到提取區(qū)域的水流方向矩陣及匯流累積量矩陣。
(2)基于DEM數(shù)據(jù)獲取實(shí)驗(yàn)樣區(qū)地表的正負(fù)地形屬性［15］，見(jiàn)圖5。將正地形區(qū)域柵格屬性設(shè)為0，負(fù)地形區(qū)域柵格屬性設(shè)為1，作為權(quán)重矩陣。
(3)將步驟(1)得到的匯流累積量矩陣與權(quán)重矩陣做點(diǎn)乘運(yùn)算，得到新的匯流累積量矩陣。
(4)基于新的匯流累積量矩陣設(shè)置閾值，得到正確的溝谷網(wǎng)絡(luò)，其中閾值的選取要保證達(dá)到所要提取的最小溝谷的匯流累積量。
2.3 實(shí)驗(yàn)結(jié)果與分析
實(shí)驗(yàn)結(jié)果見(jiàn)圖6，圖7是基于該地區(qū)的DOM數(shù)據(jù)手動(dòng)勾繪的溝谷網(wǎng)絡(luò)。實(shí)驗(yàn)結(jié)果與圖1(b)的原始結(jié)果比較可以看出，本文提出的消除偽溝谷的方法可以有效去除平坦塬面上的偽溝谷，并且細(xì)小的溝谷以及平坦的大溝谷都得到了很好的保留;根據(jù)處理結(jié)果計(jì)算的溝壑密度為1.06 km/km2，而根據(jù)原始數(shù)據(jù)計(jì)算的溝壑密度為2.05 km/km2，二者的差異是相當(dāng)明顯的。由手動(dòng)勾繪得到的溝谷網(wǎng)絡(luò)，相應(yīng)的溝壑密度為0.92 km/km2，其值略小于偽溝谷消除后的結(jié)果。從圖中可以看出，手動(dòng)勾繪的溝谷大部分是直線，而消除偽溝谷后的結(jié)果中溝谷有一定的彎曲，相對(duì)于手動(dòng)勾繪的溝谷網(wǎng)絡(luò)，偽溝谷消除以后的溝谷網(wǎng)絡(luò)更符合實(shí)際情況。
為了定量描述實(shí)驗(yàn)方法對(duì)偽溝谷的消除效果，如果將消除偽溝谷后的溝谷網(wǎng)絡(luò)近似看作實(shí)際區(qū)域的溝谷網(wǎng)絡(luò)，則定義偽溝谷消除率為偽溝谷總長(zhǎng)度與實(shí)際溝谷總長(zhǎng)度的比值，即
式中:R是溝谷消除率;∑L真是真實(shí)溝谷的總長(zhǎng)度，這里把消除偽溝谷后的溝谷網(wǎng)絡(luò)作為真實(shí)溝谷網(wǎng)絡(luò)的近似值;∑L偽是偽溝谷的總長(zhǎng)度，在計(jì)算時(shí)可以用原始提取的溝谷總長(zhǎng)度減去真實(shí)溝谷的總長(zhǎng)度。
表2為實(shí)驗(yàn)樣區(qū)消除偽溝谷之前的溝壑密度(DS1)、消除偽溝谷后的溝壑密度(DS2)以及偽溝谷消除率(R)的計(jì)算值。
表2 實(shí)驗(yàn)前后溝壑密度以及偽溝谷消除率計(jì)算值
從表中可以看出，在各種地貌類(lèi)型的實(shí)驗(yàn)樣區(qū)中，黃土塬的偽溝谷消除率最高，黃土臺(tái)塬次之，黃土殘塬的偽溝谷消除率最低。這是因?yàn)辄S土塬區(qū)存在大面積的平坦區(qū)域，因此在采用坡面流模擬方法進(jìn)行溝谷網(wǎng)絡(luò)提取時(shí)，這些平坦區(qū)域的匯流累積量值很大，造成了大量的偽溝谷;黃土臺(tái)塬的情況類(lèi)似，只是塬面平坦區(qū)域面積有所減小，因此偽溝谷消除率較之黃土塬區(qū)也有所下降;而黃土殘塬地區(qū)由于破碎程度較高，不存在大面積的平坦區(qū)域，所以這些地區(qū)在提取溝谷網(wǎng)絡(luò)時(shí)產(chǎn)生的偽溝谷數(shù)量較少，偽溝谷消除率在三類(lèi)實(shí)驗(yàn)樣區(qū)中也是最小的。
3 結(jié)論與討論
本文根據(jù)黃土塬、黃土殘塬以及黃土臺(tái)塬等地區(qū)的地形特點(diǎn)，以該區(qū)域的正負(fù)地形屬性作為權(quán)重矩陣，與原始匯流累積量矩陣做點(diǎn)乘運(yùn)算產(chǎn)生新的匯流累積量矩陣，以此匯流累積量矩陣設(shè)置閾值所提取的溝谷水系有效地去除了原始方法結(jié)果中的偽溝谷。該方法的偽溝谷消除效果對(duì)于大面積平坦地區(qū)最為有效，隨著樣區(qū)地形破碎程度的增加，偽溝谷消除率呈遞減的趨勢(shì)。
［1］白占國(guó).從地貌空間結(jié)構(gòu)特征預(yù)測(cè)土壤侵蝕的研究——以神府、東勝煤田區(qū)為例［J］.中國(guó)水土保持，1993(12):23－24.
［2］盧金發(fā).黃河中游流域地貌形態(tài)對(duì)流域產(chǎn)沙量的影響［J］.地理研究，2002，21(2):171 －178.
［3］廖義善，蔡強(qiáng)國(guó)，秦奮，等.基于DEM的黃土丘陵溝壑區(qū)流域地貌特征及侵蝕產(chǎn)沙尺度研究［J］.水土保持學(xué)報(bào)，2008，22(1):1 －6.
［4］O’Callaghan J F，Mark D M.The extraction of drainage networks from digital elevation data［J］.Computer Vision，Graphics，and Image Processing，1984，28(3):323 －344.
［5］韋中亞，周貴云，羅萬(wàn)勤.一種基于數(shù)學(xué)形態(tài)學(xué)的溝壑密度提取算法［J］.地理學(xué)與國(guó)土研究，2001，17(2):24 －27.
［6］陳永良，劉大有，虞強(qiáng)源.從DEM中自動(dòng)提取自然水系［J］.中國(guó)圖象圖形學(xué)報(bào)，2002，7(1):91 －96.
［7］李昌峰，馮學(xué)智，趙銳.流域水系自動(dòng)提取的方法和應(yīng)用［J］.湖泊科學(xué)，2003，15(3):205 －212.
［8］謝順平，都金康，羅維佳，等.基于DEM的復(fù)雜地形流域特征提取［J］.地理研究，2006，25(1):96－102.
［9］祝士杰，湯國(guó)安，賈旖旎.基于DLG水系的DEM修正方法研究［J］.現(xiàn)代測(cè)繪，2009，32(6):3 －4.
［10］徐濤，胡光道.基于數(shù)字高程模型自動(dòng)提取水系的若干問(wèn)題［J］.地理與地理信息科學(xué)，2004，20(5):11 －14.
［11］吳良超.基于DEM的黃土高原溝壑特征及其空間分異規(guī)律研究［D］.西安:西北大學(xué)，2005.
［12］李俊，湯國(guó)安，張婷，等.利用DEM提取陜北黃土高原溝谷網(wǎng)絡(luò)的匯流閾值研究［J］.水土保持通報(bào)，2007，27(2):75－78.
［13］陳濤，朱金兆，馬浩，等.黃土高原溝壑區(qū)溝道提取時(shí)閾值范圍的研究［J］.長(zhǎng)江科學(xué)院院報(bào)，2008，25(3):28 －31.
［14］周貴云，劉瑜，鄔倫.基于數(shù)字高程模型的水系提取算法［J］.地理學(xué)與國(guó)土研究，2000，16(4):77 －81.
［15］周毅.基于DEM的黃土正負(fù)地形特征研究［D］.南京:南京師范大學(xué)，2008.
猜你喜歡
我做了一項(xiàng)小實(shí)驗(yàn)作文·小學(xué)低年級(jí)(2025年2期)2025-02-13 00:00:00 記住“三個(gè)字”，寫(xiě)好小實(shí)驗(yàn)小雪花·小學(xué)生快樂(lè)作文(2024年11期)2024-12-31 00:00:00 我做了一項(xiàng)小實(shí)驗(yàn)作文·小學(xué)低年級(jí)(2024年2期)2024-04-29 00:00:00 我做了一項(xiàng)小實(shí)驗(yàn)作文·小學(xué)低年級(jí)(2023年3期)2023-04-29 00:00:00 記一次有趣的實(shí)驗(yàn)小獼猴智力畫(huà)刊(2022年9期)2022-11-04 02:31:42 有趣的實(shí)驗(yàn)小主人報(bào)(2022年4期)2022-08-09 08:52:06 微型實(shí)驗(yàn)里看“燃燒”中學(xué)生數(shù)理化·中考版(2022年11期)2022-02-16 07:01:20 做個(gè)怪怪長(zhǎng)實(shí)驗(yàn)小哥白尼(趣味科學(xué))(2019年6期)2019-10-10 01:01:50 NO與NO2相互轉(zhuǎn)化實(shí)驗(yàn)的改進(jìn)發(fā)明與創(chuàng)新(2016年38期)2016-08-22 03:02:52 實(shí)踐十號(hào)上的19項(xiàng)實(shí)驗(yàn)太空探索(2016年5期)2016-07-12 15:17:55
趙明偉，湯國(guó)安，李發(fā)源
(南京師范大學(xué)地理科學(xué)學(xué)院，江蘇南京210046)
1 問(wèn)題的提出
溝壑密度是反映黃土高原地區(qū)地面受侵蝕程度重要的區(qū)域水土流失因子，在該地區(qū)的土壤侵蝕與侵蝕產(chǎn)沙研究中得到了廣泛的應(yīng)用［1－3］。目前，不同學(xué)者提出了不同的提取溝壑網(wǎng)絡(luò)的方法［4－7］，其中基于DEM的坡面流模擬方法因?yàn)閿?shù)據(jù)來(lái)源廣泛、方法簡(jiǎn)單、算法穩(wěn)定性高而得到普遍應(yīng)用。這一方法根據(jù)水文學(xué)坡面流模型來(lái)判別水流路徑，首先對(duì)地面洼地進(jìn)行填充，構(gòu)建無(wú)洼地的連續(xù)地面［8］，并根據(jù)每一個(gè)柵格單元與相鄰單元之間的最陡坡度，識(shí)別求取沿水流路徑的匯流累積量值，再選擇適當(dāng)?shù)膮R流累積量閾值來(lái)確定河網(wǎng)。該方法可以直接生成連續(xù)的河網(wǎng)。但是，在計(jì)算格網(wǎng)單元流量累積量值時(shí)，本不會(huì)產(chǎn)生地面徑流的平坦地區(qū)(例如黃土塬面等)也將產(chǎn)生很大的匯流累積量值，從而造成提取溝谷網(wǎng)絡(luò)時(shí)會(huì)有大量的偽溝谷出現(xiàn)［9－10］。以上問(wèn)題的出現(xiàn)，給黃土高原地區(qū)溝壑密度的正確提取以及水文分析的有效進(jìn)行，都帶來(lái)了很大的困難。
圖1為實(shí)驗(yàn)樣區(qū)水文分析中存在的偽溝谷示意圖，(a)為實(shí)驗(yàn)樣區(qū)的地形等高線圖，(b)為該地區(qū)采用坡面流方法提取的溝谷網(wǎng)絡(luò)，對(duì)比可以發(fā)現(xiàn)，原始提取的溝谷網(wǎng)絡(luò)中存在大量的偽溝谷。
圖1 偽溝谷示意
影響提取溝谷網(wǎng)絡(luò)結(jié)果的關(guān)鍵因素是匯流累積量閾值的選取［11］，一般說(shuō)來(lái)，當(dāng)該閾值變大時(shí)，所提取的溝谷網(wǎng)絡(luò)結(jié)果數(shù)量會(huì)變少。李俊、陳濤等以陜北黃土高原為研究樣區(qū)，研究了溝壑網(wǎng)絡(luò)提取中匯流累積量閾值的選取依據(jù)［12－13］。但是，對(duì)于黃土塬、黃土殘塬等地區(qū)，由于大量平地的存在，試圖僅僅通過(guò)調(diào)整匯流累積量閾值來(lái)消除偽溝谷是不可行的。圖2為DEM數(shù)據(jù)的山體陰影暈渲圖，可以看出A點(diǎn)處于溝谷中，而B(niǎo)點(diǎn)處于塬面平坦區(qū);對(duì)應(yīng)于水流路徑匯流累積量圖3，B點(diǎn)的匯流累積量值卻大大高于A點(diǎn)。由此可見(jiàn)，在黃土塬、黃土殘塬等存在大面積平坦區(qū)域的地區(qū)，不可能存在合適的閾值，使得提取的溝谷網(wǎng)絡(luò)既包含較小的溝谷，又不存在偽溝谷。
如上所述，在平原臺(tái)地地區(qū)(如黃土塬面)，由于大量平坦區(qū)域的存在，使得當(dāng)采用坡面流模擬方法進(jìn)行溝谷提取時(shí)會(huì)產(chǎn)生大量偽溝谷，且無(wú)法通過(guò)調(diào)整匯流累積量閾值來(lái)消除。因此，急需尋求一種簡(jiǎn)單而有效的解決方法。
圖2 山體陰影暈渲圖
圖3 匯流累積量圖
2 實(shí)驗(yàn)方法
2.1 實(shí)驗(yàn)樣區(qū)與數(shù)據(jù)
實(shí)驗(yàn)樣區(qū)位于陜西省淳化縣城西北部涇河中游地區(qū)，黃土塬及殘塬為主要地貌類(lèi)型，塬面地形平緩，但被諸多大溝谷深切割裂，平均海拔1045.04 m，相對(duì)高差390 m，塬面平均坡度為8.1°，面積100 km2(10 km×10 km)。此外，實(shí)驗(yàn)中為了對(duì)比文中提出的實(shí)驗(yàn)方法消除偽溝谷的效果，還選取了地貌類(lèi)型屬于黃土臺(tái)塬、黃土殘塬的樣區(qū)進(jìn)行附加實(shí)驗(yàn)。附加實(shí)驗(yàn)樣區(qū)的基本情況如表1所示。
表1 附加實(shí)驗(yàn)樣區(qū)概況
實(shí)驗(yàn)數(shù)據(jù)為國(guó)家測(cè)繪部門(mén)生產(chǎn)的5 m分辨率的DEM數(shù)據(jù)，數(shù)據(jù)精度滿(mǎn)足實(shí)驗(yàn)要求。
2.2 實(shí)驗(yàn)思路與方法
如前所述，基于DEM提取溝谷網(wǎng)絡(luò)是在得到水流路徑匯流累積量矩陣的基礎(chǔ)上，根據(jù)一定的閾值提取溝谷網(wǎng)絡(luò)。水流路徑匯流累積量矩陣的產(chǎn)生示意圖見(jiàn)圖4，可以看出，某點(diǎn)匯流累積量值的確定既受水流方向的控制，同時(shí)也與分析區(qū)域面積緊密相關(guān)。因此，在大面積的平原臺(tái)地地區(qū)，雖然地勢(shì)變化緩慢，但是仍然會(huì)產(chǎn)生很大的匯流累積量值，因而采用常規(guī)方法提取水系時(shí)會(huì)產(chǎn)生很多偽溝谷。
圖4 匯流累積量矩陣計(jì)算示意
根據(jù)以上分析，可以在計(jì)算匯流累積量矩陣時(shí)對(duì)每個(gè)柵格點(diǎn)增加權(quán)重值，使得平原臺(tái)地地區(qū)匯流累積量減小，從而減小偽溝谷出現(xiàn)的可能性。實(shí)驗(yàn)發(fā)現(xiàn)，在黃土高原地區(qū)，將地形表面劃分為正地形與負(fù)地形之后，溝谷網(wǎng)絡(luò)屬于負(fù)地形區(qū)域。因此，本實(shí)驗(yàn)根據(jù)地表面的正負(fù)地形屬性來(lái)設(shè)置權(quán)重值，正地形區(qū)域的匯流累積量權(quán)重值為0，負(fù)地形區(qū)域的匯流累積量權(quán)重值為1;將權(quán)重矩陣與原始匯流累積量矩陣做點(diǎn)乘運(yùn)算得到新的匯流累積量矩陣，再設(shè)置閾值提取溝谷網(wǎng)絡(luò)。綜合以上分析，在平原臺(tái)地地區(qū)基于DEM數(shù)據(jù)提取正確的溝谷網(wǎng)絡(luò)的步驟為:
(1)基于DEM數(shù)據(jù)得到匯流累積量矩陣［14］。首先，對(duì)原始DEM數(shù)據(jù)進(jìn)行洼地填充處理，然后，利用GIS水文分析，得到提取區(qū)域的水流方向矩陣及匯流累積量矩陣。
(2)基于DEM數(shù)據(jù)獲取實(shí)驗(yàn)樣區(qū)地表的正負(fù)地形屬性［15］，見(jiàn)圖5。將正地形區(qū)域柵格屬性設(shè)為0，負(fù)地形區(qū)域柵格屬性設(shè)為1，作為權(quán)重矩陣。
(3)將步驟(1)得到的匯流累積量矩陣與權(quán)重矩陣做點(diǎn)乘運(yùn)算，得到新的匯流累積量矩陣。
(4)基于新的匯流累積量矩陣設(shè)置閾值，得到正確的溝谷網(wǎng)絡(luò)，其中閾值的選取要保證達(dá)到所要提取的最小溝谷的匯流累積量。
2.3 實(shí)驗(yàn)結(jié)果與分析
實(shí)驗(yàn)結(jié)果見(jiàn)圖6，圖7是基于該地區(qū)的DOM數(shù)據(jù)手動(dòng)勾繪的溝谷網(wǎng)絡(luò)。實(shí)驗(yàn)結(jié)果與圖1(b)的原始結(jié)果比較可以看出，本文提出的消除偽溝谷的方法可以有效去除平坦塬面上的偽溝谷，并且細(xì)小的溝谷以及平坦的大溝谷都得到了很好的保留;根據(jù)處理結(jié)果計(jì)算的溝壑密度為1.06 km/km2，而根據(jù)原始數(shù)據(jù)計(jì)算的溝壑密度為2.05 km/km2，二者的差異是相當(dāng)明顯的。由手動(dòng)勾繪得到的溝谷網(wǎng)絡(luò)，相應(yīng)的溝壑密度為0.92 km/km2，其值略小于偽溝谷消除后的結(jié)果。從圖中可以看出，手動(dòng)勾繪的溝谷大部分是直線，而消除偽溝谷后的結(jié)果中溝谷有一定的彎曲，相對(duì)于手動(dòng)勾繪的溝谷網(wǎng)絡(luò)，偽溝谷消除以后的溝谷網(wǎng)絡(luò)更符合實(shí)際情況。
為了定量描述實(shí)驗(yàn)方法對(duì)偽溝谷的消除效果，如果將消除偽溝谷后的溝谷網(wǎng)絡(luò)近似看作實(shí)際區(qū)域的溝谷網(wǎng)絡(luò)，則定義偽溝谷消除率為偽溝谷總長(zhǎng)度與實(shí)際溝谷總長(zhǎng)度的比值，即
式中:R是溝谷消除率;∑L真是真實(shí)溝谷的總長(zhǎng)度，這里把消除偽溝谷后的溝谷網(wǎng)絡(luò)作為真實(shí)溝谷網(wǎng)絡(luò)的近似值;∑L偽是偽溝谷的總長(zhǎng)度，在計(jì)算時(shí)可以用原始提取的溝谷總長(zhǎng)度減去真實(shí)溝谷的總長(zhǎng)度。
表2為實(shí)驗(yàn)樣區(qū)消除偽溝谷之前的溝壑密度(DS1)、消除偽溝谷后的溝壑密度(DS2)以及偽溝谷消除率(R)的計(jì)算值。
表2 實(shí)驗(yàn)前后溝壑密度以及偽溝谷消除率計(jì)算值
從表中可以看出，在各種地貌類(lèi)型的實(shí)驗(yàn)樣區(qū)中，黃土塬的偽溝谷消除率最高，黃土臺(tái)塬次之，黃土殘塬的偽溝谷消除率最低。這是因?yàn)辄S土塬區(qū)存在大面積的平坦區(qū)域，因此在采用坡面流模擬方法進(jìn)行溝谷網(wǎng)絡(luò)提取時(shí)，這些平坦區(qū)域的匯流累積量值很大，造成了大量的偽溝谷;黃土臺(tái)塬的情況類(lèi)似，只是塬面平坦區(qū)域面積有所減小，因此偽溝谷消除率較之黃土塬區(qū)也有所下降;而黃土殘塬地區(qū)由于破碎程度較高，不存在大面積的平坦區(qū)域，所以這些地區(qū)在提取溝谷網(wǎng)絡(luò)時(shí)產(chǎn)生的偽溝谷數(shù)量較少，偽溝谷消除率在三類(lèi)實(shí)驗(yàn)樣區(qū)中也是最小的。
3 結(jié)論與討論
本文根據(jù)黃土塬、黃土殘塬以及黃土臺(tái)塬等地區(qū)的地形特點(diǎn)，以該區(qū)域的正負(fù)地形屬性作為權(quán)重矩陣，與原始匯流累積量矩陣做點(diǎn)乘運(yùn)算產(chǎn)生新的匯流累積量矩陣，以此匯流累積量矩陣設(shè)置閾值所提取的溝谷水系有效地去除了原始方法結(jié)果中的偽溝谷。該方法的偽溝谷消除效果對(duì)于大面積平坦地區(qū)最為有效，隨著樣區(qū)地形破碎程度的增加，偽溝谷消除率呈遞減的趨勢(shì)。
［1］白占國(guó).從地貌空間結(jié)構(gòu)特征預(yù)測(cè)土壤侵蝕的研究——以神府、東勝煤田區(qū)為例［J］.中國(guó)水土保持，1993(12):23－24.
［2］盧金發(fā).黃河中游流域地貌形態(tài)對(duì)流域產(chǎn)沙量的影響［J］.地理研究，2002，21(2):171 －178.
［3］廖義善，蔡強(qiáng)國(guó)，秦奮，等.基于DEM的黃土丘陵溝壑區(qū)流域地貌特征及侵蝕產(chǎn)沙尺度研究［J］.水土保持學(xué)報(bào)，2008，22(1):1 －6.
［4］O’Callaghan J F，Mark D M.The extraction of drainage networks from digital elevation data［J］.Computer Vision，Graphics，and Image Processing，1984，28(3):323 －344.
［5］韋中亞，周貴云，羅萬(wàn)勤.一種基于數(shù)學(xué)形態(tài)學(xué)的溝壑密度提取算法［J］.地理學(xué)與國(guó)土研究，2001，17(2):24 －27.
［6］陳永良，劉大有，虞強(qiáng)源.從DEM中自動(dòng)提取自然水系［J］.中國(guó)圖象圖形學(xué)報(bào)，2002，7(1):91 －96.
［7］李昌峰，馮學(xué)智，趙銳.流域水系自動(dòng)提取的方法和應(yīng)用［J］.湖泊科學(xué)，2003，15(3):205 －212.
［8］謝順平，都金康，羅維佳，等.基于DEM的復(fù)雜地形流域特征提取［J］.地理研究，2006，25(1):96－102.
［9］祝士杰，湯國(guó)安，賈旖旎.基于DLG水系的DEM修正方法研究［J］.現(xiàn)代測(cè)繪，2009，32(6):3 －4.
［10］徐濤，胡光道.基于數(shù)字高程模型自動(dòng)提取水系的若干問(wèn)題［J］.地理與地理信息科學(xué)，2004，20(5):11 －14.
［11］吳良超.基于DEM的黃土高原溝壑特征及其空間分異規(guī)律研究［D］.西安:西北大學(xué)，2005.
［12］李俊，湯國(guó)安，張婷，等.利用DEM提取陜北黃土高原溝谷網(wǎng)絡(luò)的匯流閾值研究［J］.水土保持通報(bào)，2007，27(2):75－78.
［13］陳濤，朱金兆，馬浩，等.黃土高原溝壑區(qū)溝道提取時(shí)閾值范圍的研究［J］.長(zhǎng)江科學(xué)院院報(bào)，2008，25(3):28 －31.
［14］周貴云，劉瑜，鄔倫.基于數(shù)字高程模型的水系提取算法［J］.地理學(xué)與國(guó)土研究，2000，16(4):77 －81.
［15］周毅.基于DEM的黃土正負(fù)地形特征研究［D］.南京:南京師范大學(xué)，2008.
猜你喜歡
作文·小學(xué)低年級(jí)(2025年2期)2025-02-13 00:00:00
小雪花·小學(xué)生快樂(lè)作文(2024年11期)2024-12-31 00:00:00
作文·小學(xué)低年級(jí)(2024年2期)2024-04-29 00:00:00
作文·小學(xué)低年級(jí)(2023年3期)2023-04-29 00:00:00
小獼猴智力畫(huà)刊(2022年9期)2022-11-04 02:31:42
小主人報(bào)(2022年4期)2022-08-09 08:52:06
中學(xué)生數(shù)理化·中考版(2022年11期)2022-02-16 07:01:20
小哥白尼(趣味科學(xué))(2019年6期)2019-10-10 01:01:50
發(fā)明與創(chuàng)新(2016年38期)2016-08-22 03:02:52
太空探索(2016年5期)2016-07-12 15:17:55
