苗晶，王春梅，文佳昕，土祥，馬思煜

(西北大學城市與環(huán)境學院，陜西 西安，710127)

地面坡度是對地面傾斜程度的定量描述，也是一種基本的地貌形態(tài)指標，在流域水文分析與區(qū)域土壤侵蝕模擬等研究領域是重要的參數(shù)［1，2］。坡度數(shù)據(jù)一般基于中低分辨率遙感源 DEM 提?。?，4］。數(shù)字高程模型(Digital Elevation Model，DEM)是用一組有序數(shù)值陣列形式表示地面高程的一種實體地面模型，簡言之是一種地面形態(tài)的數(shù)字表現(xiàn)形式，基于數(shù)字高程模型(DEM)提取的各種地形因子如坡度、坡長、坡向、曲率等已經成為賴以進行地形分析的基本信息源和核心數(shù)據(jù)［5－7］。其中，覆蓋地球80%以上陸地表面的雷達影像數(shù)據(jù)SRTM(Shuttle Radar Topography Mission)［8］和覆蓋地球 99%陸地表面的 GDEM(Global Digital Elevation Model)［9］是目前應用較為廣泛的、可免費獲取的中低分辨率遙感源DEM數(shù)據(jù)。隨著分辨率的降低，基于中低分辨率DEM提取的坡度數(shù)據(jù)出現(xiàn)了明顯的衰減，即坡度的尺度效應問題［10，11］。這對中低分辨率遙感源DEM的高效使用造成了很大的局限，并限制了其在流域水文分析［12，13］與區(qū)域土壤侵蝕模擬等研究領域的應用。目前網上可免費獲取的30m分辨率1arc－GDEM和90m分辨率3arc－SRTM的高程及坡度數(shù)據(jù)迄今為止仍未對上述坡度尺度效應問題做處理。

已有研究表明，隨著DEM分辨率的降低，在其上求取的坡度不斷趨于平緩［14，15］。因而不能直接有效地表現(xiàn)地形的起伏狀況，也不能有效地提取與坡度有關的地貌和水文模型參數(shù)［16］，進而影響水文分析、土壤侵蝕模擬等的計算精度［17］。針對基于中低分辨率遙感源DEM提取坡度數(shù)據(jù)的尺度效應問題，研究者們做了大量的探索，認為對衰減了的坡度進行降尺度變換是可行的，但行之有效的科學方法正在探索中。湯國安提出通過對坡度圖譜的研究完成兩種比例尺之間各種坡度級別的統(tǒng)計學轉換，并開拓了一個十分活躍的研究領域——坡度圖譜研究［14，18，19］。楊昕針對地形濕度指數(shù)的提取也對坡度降尺度變換進行了比較系統(tǒng)的研究，基于對坡度遞減速度域和分形域的分析，在不同尺度范圍分別采用遞減速率方法(分辨率5～100 m)和分形法(分辨率100～1000 m)進行變換并取得較好效果［20］。楊勤科等提出了一種基于直方圖匹配［21－24］的方法，并實現(xiàn)了 100 m、50 m、25 m、10 m和5 m分辨率坡度之間的尺度轉換，取得了較好的效果。

本研究擬在直方圖匹配原理的基礎上，選取位于黃土高原丘陵溝壑區(qū)的延河流域為研究區(qū)，參照10m分辨率坡度數(shù)據(jù)，構建坡度降尺度變換模型，對衰減了的坡度進行降尺度變換，并在空間分布格局和統(tǒng)計意義上驗證變換精度，得到一批高質量、可共享的坡度數(shù)據(jù)和一套行之有效的中低分辨率遙感源DEM坡度數(shù)據(jù)降尺度變換方法。同時，本研究將為1arc－GDEM和3arc－SRTM數(shù)據(jù)在土壤侵蝕模擬、水文分析、地貌定量分析等領域的有效利用提供理論依據(jù)和技術方法。

1 研究方法

1.1 研究區(qū)

研究區(qū)位于延河流域，該流域地貌主要為黃土丘陵溝壑，梁峁丘陵相間分布，地形空間分異復雜。流域地勢西北高、東南低，地勢形態(tài)明顯表現(xiàn)為三種類型。河源至真武洞的上游，為峁梁丘陵溝壑區(qū)，梁多而峁小，河床比降大，植被稀少，侵蝕強烈;真武洞至甘谷驛的中游，為峁狀丘陵溝壑區(qū)，梁窄峁小，河谷寬闊，階地發(fā)育;甘谷驛至河口下游，為破碎塬區(qū)，塬面窄小，沖溝發(fā)育。據(jù)此將延河流域按上中下游劃分為三個分區(qū)(圖1)，在每個分區(qū)內選擇一個建模區(qū)和一個檢驗區(qū)(每個建模區(qū)與檢驗區(qū)均為面積約為40至70平方公里的完整小流域)。

圖1 研究區(qū)及分區(qū)示意圖

1.2 數(shù)據(jù)基礎

本研究基礎數(shù)據(jù)包括:延河流域3arc－SRTM數(shù)據(jù)(分辨率約為90 m)、延河流域1arc－GDEM數(shù)據(jù)(分辨率約為30 m)(作為坡度降尺度變換研究的基礎數(shù)據(jù))以及延河流域建模區(qū)與檢驗區(qū)六個小流域的1∶50，000地形圖(基于該數(shù)據(jù)構建10 m分辨率DEM并用于坡度降尺度變換模型的構建與驗證)，收集整理的基礎數(shù)據(jù)見表1。

表1 基礎數(shù)據(jù)

1.3 研究方法

1.3.1 數(shù)據(jù)預處理

對下載的延河流域1arc－GDEM與3arc－SRTM的高程數(shù)據(jù)進行空值填充、濾波除噪，并提取坡度數(shù)據(jù)。對1∶50，000地形圖掃描、數(shù)字化，在ANUDEM軟件環(huán)境下建立10 m分辨率DEM并提取坡度，得到檢驗區(qū)與建模區(qū)高分辨率坡度數(shù)據(jù)。

1.3.2 坡度降尺度模型構建

在各建模區(qū)基于經預處理所得1arc－GDEM與3arc－SRTM坡度數(shù)據(jù)，以10 m分辨率DEM坡度數(shù)據(jù)為參照，根據(jù)直方圖匹配原理［25］，建立坡度降尺度變換模型。最后，在各檢驗區(qū)進行坡度降尺度變換及精度驗證，運用構建的模型對中低分辨率遙感源DEM坡度數(shù)據(jù)做降尺度變化，與高分辨率坡度數(shù)據(jù)在統(tǒng)計意義上作比較，同時分析變換后坡度的空間分布格局。技術路線見圖2。

圖2 技術路線

1.3.3 坡度降尺度變換精度驗證

本研究從統(tǒng)計特征和空間分布格局兩個方面驗證精度:

(1)基于統(tǒng)計特征的驗證

在各模型檢驗區(qū)，將較低分辨率1arc－GDEM與3arc－SRTM坡度數(shù)據(jù)進行降尺度變換，并繪制直方圖與10 m高分辨率坡度直方圖對比判斷精度，以相交直方圖相似度指數(shù)［25］(即直方圖的交Histogram Intersection，HI)作為評價指標。

式中X、Y分別表示不同分辨率坡度數(shù)據(jù)的頻率曲線，xi和yi表示在第i個分級間隔上X和Y頻率曲線對應的頻率值。HI的取值范圍為0≤HI(X，Y)≤1，其值越大，表明變換后的較低分辨率頻率曲線與較高分辨率頻率曲線重合度越大，精度越高。

(2)基于空間分布格局的驗證

通過比較降尺度變換后的1arc－GDEM、3arc－SRTM坡度表面與10m高分辨率坡度表面空間分布格局的差異性，討論模型的適用性。

2 結果與分析

2.1 坡度尺度效應分析

隨著分辨率降低，基于中低遙感源DEM提取的坡度數(shù)據(jù)較高分辨率坡度發(fā)生明顯衰減，不能準確表現(xiàn)地表起伏度，圖3為分區(qū)Ⅰ的建模區(qū)Ⅰ_m在不同分辨率下的坡度表面，可以明顯看出隨著分辨率的下降，坡度所表現(xiàn)的地形越來越粗糙，坡度的層次減少，地貌結構簡化，很多細碎的地形特征難以表達。從表2的統(tǒng)計特征值可以看出隨著分辨率的降低，坡度數(shù)據(jù)的最大值、均值和標準差都發(fā)上了不同程度的減小。最大值由67.82減小到63.74和40.51:平均值在10 m分辨率時為26.28，在30 m分辨率時減小到16.86，到90 m分辨率時坡度為13.15:標準差也從10.90減少到8.78和6.55。

圖3 不同分辨率坡度表面差異

表2 Ⅰ_m區(qū)基于不同分辨率DEM的坡度統(tǒng)計表

2.2 坡度降尺度變換

2.2.1 坡度降尺度變換模型構建

本研究將研究區(qū)分為建模區(qū)與檢驗區(qū)兩個部分，表3是不同研究區(qū)的坡度降尺度變換模型。

2.2.2 坡度降尺度變換結果分析

將構建的三個樣區(qū)坡度降尺度變換模型應用在延河流域對應三個分區(qū)，得到延河流域經降尺度變換后的1arc－GDEM與3arc－SRTM坡度數(shù)據(jù)，圖4是三個分區(qū)中低分辨率坡度表面、降尺度變換后的坡度表面，圖5為三個分區(qū)拼接后整個延河流域的中低分辨率坡度及降尺度變換坡度表面。從表4坡度降尺度變換前后的統(tǒng)計特征表達上看，經過降尺度變換，坡度衰減現(xiàn)象得到較好的改善，坡度值整體升高，空間格局更符合研究區(qū)實際情況。

表3 降尺度變換模型

圖4 三個分區(qū)中低分辨率坡度、降尺度變換坡度

圖5 延河流域中低分辨率坡度、降尺度變換坡度

表4 降尺度變換前后各分區(qū)及延河流域坡度基本統(tǒng)計特征

2.3 坡度降尺度變換結果精度驗證

表5是檢驗區(qū)坡度變換前后基本統(tǒng)計特征數(shù)據(jù)，統(tǒng)計結果表明，變換后的坡度較變換前的低分辨率坡度數(shù)據(jù)的最大值、平均值及標準差都增大，各統(tǒng)計值已經接近較高分辨率的DEM坡度的統(tǒng)計值。從圖6不同分辨率的頻率曲線可以看出，變換后坡度的頻率曲線相比變換前中低分辨率坡度的頻率曲線峰值向高坡度范圍移動，與高分辨率坡度頻率曲線分布更加接近，變換后頻率曲線相交部分的面積即直方圖的交HI的值都在0.85以上(表6)，表明坡度降尺度變換模型的精度相對較高。

表5 坡度降尺度變換前后基本統(tǒng)計特征

圖6 頻率曲線

表6 坡度降尺度變換前后HI統(tǒng)計表

圖7 三個檢驗區(qū)中低分辨率坡度及降尺度變換后的坡度表面

圖7是三個檢驗區(qū)中低分辨率坡度及降尺度變換后的坡度表面。圖中可以明顯看出降尺度變換后比變換前的坡度表面更接近于高分辨率坡度表面，較高坡度所占比例增加。

3 結論與討論

3.1 結論

1arc－GDEM＆3arc－SRTM數(shù)據(jù)是目前可免費獲取、覆蓋范圍廣泛的中低分變率遙感源DEM數(shù)據(jù)，然而從中提取的坡度數(shù)據(jù)出現(xiàn)坡度衰減問題限制了其應用，本研究通過構建數(shù)學模型對其進行降尺度變換，生成一個新的坡度表面，有效地解決了該問題，并通過驗證將模型應用于較大區(qū)域。主要研究結論如下:

1)相比10m分辨率坡度，基于1arc－GDEM和3arc－SRTM的坡度數(shù)據(jù)發(fā)生明顯衰減，平均值與標準差均降低，這種衰減影響了中低分辨率遙感源坡度數(shù)據(jù)的應用;

2)本文通過劃分降尺度變換單元，在延河流域各變換單元基于直方圖匹配原理，通過在各變換單元選擇典型小流域構建模型，構建的坡度降尺度變換模型經過在各檢驗樣區(qū)的驗證，其模型精度較高;

3)本研究構建的分區(qū)坡度降尺度變換模型可應用于延河流域的坡度降尺度變換，可將延河流域基于1arc－GDEM與3arc－SRTM的坡度數(shù)據(jù)進行降尺度變換，使其在統(tǒng)計信息上接近10m分辨率坡度，從而可提供一套延河流域質量較好、可公開使用的遙感源坡度數(shù)據(jù)，為土壤侵蝕、水文等領域提供數(shù)據(jù)支持。

3.2 討論

本研究對于流域水文分析與區(qū)域土壤侵蝕模擬中地形參數(shù)求取和中低分辨率DEM高效利用等，均具有重要意義。然而本研究還存在以下不足，以期在今后研究中進行進一步探討:

1)本研究中應用的1arc－GDEM和3arc－SRTM數(shù)據(jù)，對于地表植被等非地面高程信息未進行處理，今后將探討該方面較為有效的方法進行改進。

2)經過降尺度變換坡度的衰減問題在統(tǒng)計特征和整體空間格局上得到了改善，然而在對地形頻率層次的表達上有無可能達到10 m分辨率的水平依然有待進一步研究。

［1］王中根，劉昌明，吳險峰.基于DEM的分布式水文模型研究綜述［J］.自然資源學報，2003，18(2):168－173.

［2］姚志宏.基于GIS的區(qū)域水土流失過程模擬研究［D］.中國科學院研究生院:教育部水土保持與生態(tài)環(huán)境研究中心，2010.

［3］楊勤科，李銳，徐濤，姚志紅，崔琰.區(qū)域水土流失過程及其定量描述的初步研究［J］.亞熱帶水土保持，2006，18(2):20－23.

［4］姚志宏，楊勤科，吳喆，崔琰.區(qū)域尺度侵蝕產沙估算方法研究［J］.中國水土保持科學，2007，5(4):13－17.

［5］Hutchinson M F.A new procedure for gridding elevation and stream line data with automatic removal of spuriouspits［J］.Journal of Hydrology，1989，106(3):211－232.

［6］楊勤科，師維娟，Tim R.McVicar，Tom G.Van Niel，Lingtao Li.水文地貌關系正確DEM的建立方法［J］.中國水土保持科學，2007，5(4):1－6.

［7］Yang Q K，McVicar T R，Van Niel T G，Hutchinson M F，Li L T，Zhang X P.Improving a digital elevation model by reducing source data errorsand optimising interpolation algorithm parameters:an example in the Loess Plateau，China.［J］.International Journal ofApplied Earth Observation and Geoinformation，2007，9(3):235－246.

［8］Yang L P，Meng X M，Zhang X Q.SRTM DEM and its application advances［J］.International Journal of Remote Sensing，2011，32(14):3875－3896.

［9］Tachikawa T，Kaku M，Iwasaki A.ASTER GDEM Version 2 Validation Repotr［R］.Report to the ASTER GDEM Version 2 Validation Team，2011.

［10］劉新華，楊勤科，湯國安.中國地形起伏度的提取及在水土流失定量評價中的應用［J］.水土保持通報，2001，21(1):57－59，62.

［11］楊勤科，趙牡丹，劉詠梅，郭偉玲，王雷，李銳.DEM與區(qū)域土壤侵蝕地形因子研究［J］.地理信息世界，2009，7(1):25－31+45.［12］Lacroix M P，Martz L W，Kite G W，Garbrecht J.Using digital terrain analysis modeling techniques for the parameterization of a hydrologic model［J］.Environmental Modelling＆software，2002，17(2):125－134.

［13］閆寶偉，郭生練，郭靖，張娜.多變量水文分析計算方法的比較［J］.武漢大學學報(工學版)，2009，42(1):10－15.

［14］湯國安，趙牡丹，李天文，劉詠梅，謝元禮.DEM提取黃土高原地面坡度的不確定性［J］.地理學報，2003，58(6):824－830.

［15］Wang C M，Yang Q K，Guo W L，Liu H G，Jupp D，Li R.Zhang H M.Influence of resolution on slope in areas with different topographic characteristics［J］.Computers＆Geosciences，2012，41:156－168.

［16］胡鵬，楊傳勇，吳艷蘭等.新數(shù)字高程模型:理論、方法、標準和應用［M］.北京:測繪出版社.2007.

［17］任希巖，張雪松，郝芳華，程紅光，楊志峰.DEM分辨率對產流產沙模擬影響研究［J］.水土保持研究.2004，11(1):1－4，26.

［18］湯國安，劉學軍，閭國年.數(shù)字高程模型及地學分析的原理與方法［M］.北京:科學出版社.2005:338－354.

［19］李發(fā)源.黃土高原地面坡譜及空間分異研究［D］.中國科學院成都山地災害與環(huán)境研究所:中國科學院.水利部成都山地災害與環(huán)境研究所，2007.

［20］楊昕，湯國安，劉學軍，李發(fā)源，祝士杰.數(shù)字地形分析的理論、方法與應用［J］.地理學報.2009，64(9):1058－1070.

［21］師維娟.基于DEM和GIS的坡度變換方法研究［D］.西北農林科技大學.2007.

［22］Yang Q K，Jupp D，Li R，et al.Re－ scaling lower resolution slope by histogram matching［C］/Zhou Q M，Lees B G，Tang G A.Advances in Digital Terrain Analysis2.Springer－Verlag Berlin and Heidelberg GmbH＆ Co.K，2008:193－210.

［23］郭蘭勤，楊勤科，胡潔.基于分形的中低分辨率坡度降尺度變換方法研究［J］.西北農林科技大學學報(自然科學版).2011，39(12):173－180.

［24］王春梅.坡度尺度效應與轉換及其對土壤侵蝕評價影響研究［D］.中國科學院研究生院:教育部水土保持與生態(tài)環(huán)境研究中心.2012.

［25］牛亮，楊勤科.DEM尺度變換中直方圖相似度計算與應用［J］.水土保持研究.2010，17(3):120－125.