DEM分辨率對地形因子提取影響的研究

劉 娜，楊斌斌，房靖楠，郝 藝，楊海城，王維芳

(東北林業大學 林學院，哈爾濱 150000)

?

DEM分辨率對地形因子提取影響的研究

劉 娜，楊斌斌，房靖楠，郝 藝，楊海城，王維芳*

(東北林業大學 林學院，哈爾濱 150000)

以大興安嶺漠河地區1∶5萬比例尺的等高線數據為研究對象，生成不同分辨率的數字高程模型(DEM)數據，分析在不同分辯率條件下，由DEM提取的坡度、坡向及剖面曲率隨分辨率的變化情況，并利用信息熵公式進行不同分辨率DEM高程熵分析，得到高程熵隨分辨率的變化趨勢。結果表明，當DEM分辨率達到25 m左右時，各地形因子的變化情況趨于穩定，分辨率達到50 m以后變化明顯加劇。高程熵在分辨率增大到25 m時變化減慢，據此判斷，漠河地區1∶5萬比例尺等高線在生成DEM時最合適的分辨率應該是25 m。

數字高程模型(DEM)；地形因子；分辨率

0 引 言

數字高程模型(DEM)的概念于1958年由Miller首次提出，是用一組有序數值陣列形式表示地面高程的一種實體地面模型[1]。近年來用DEM提取地形因子(如高程、坡度、坡向、剖面曲率等)并對這些地形因子進行分析已經成為分析地形因子的主流方法，而通過數字高程模型提取的地形因子在水文模型建立、土壤侵蝕分析、水土流失監測、地貌形態模擬和生態環境研究等地學分析領域有著廣泛的應用[2]。DEM分辨率不同，所提取的地形因子的數值也會發生改變，從而對地形因子分析造成一定的影響。因此，對不同DEM分辨率下地形因子提取狀況的研究有著重要的意義。對DEM精度的研究自20世紀80年代開始取得突破，Brasington、Richards研究了各種模型及參數與DTM尺度的關系[3]，Brabb、Bolstad and Stower、Florinsky在探討DEM誤差對應用的影響時探討過DEM 誤差對坡度的影響[4-7]，Kienzle研究了DEM柵格分辨率對一階、二階及復合地形因子的影響[8]，Erskine 探討了DEM高程精度和柵格分辨率對地形屬性的影響[9]，湯國安提出了DEM 地形描述誤差的概念[10]，張磊等以應用目標、計算方法和尺度特征等為依據對地形因子進行了系統整合，為提高地形因子分析與應用的有效性提供了依據[11]，經過幾十年的發展，基于DEM的數字地形分析理論與技術方法正逐步走向成熟。本文以大興安嶺漠河地區作為研究區域。以該地區的高程、坡度、坡向和剖面曲率作為研究對象，研究不同分辨率DEM對提取地形因子的影響，確定最佳分辨率。

1 研究區域概況與研究方法

1.1 研究區域概況

漠河位于中國最北端，黑龍江省西北部，地理坐標為東經121°07′～124°20′，北緯52°10′～53°33′。東與塔河縣接壤，西與內蒙古額爾古納右旗交界，南與內蒙的額爾古納左旗為鄰，北與俄羅斯隔江相望。總面積18 233 km2，占大興安嶺地區總面積的21.6%，占全省總面積的3.9%。地勢南高北低，最大高差為1 350 m。

1.2 數據與方法

以大興安嶺漠河地區1∶5萬等高線為研究的基礎數據，通過ARCGIS10.2軟件將等高線生成不同分辨率的DEM，分辨率分別為5、10、15、20、25、50、75、100、150、200、250、300 m，分析由不同分辨率DEM生成的坡度，坡向以及剖面曲率圖的變化情況，并對各地形因子隨DEM分辨率變化的情況進行統計分析。

信息熵的概念最早由美國信息論創始人香農提出，其來源于信號通信理論，通過分析組成通訊信號的數字或符號的統計特征來定量表示信號通訊能力，即信息量的大小。本文采用信息熵對不同分辨率的DEM進行分析，以確定最佳分辨率。信息熵的表達式為：

H=-∑pi(x)ln[pi(x)]。

由地面高程的分布近似于正態分布N(μ，σ2)可得高程信息熵公式[12]：

H=1/2+lnσ+lnσ(2π)/2-lnΔx。

其中σ為高程的標準差，Δx為DEM分辨率。采用隨機抽樣的方法，通過ArcGIS 創建隨機點的方法實現隨機點的選取，共選取799個隨機點，再通過值提取至點來獲取每個隨機點的高程值。

2 結果與分析

2.1 坡度分析

坡度一般定義為地表水平面和實際地形表面之間的夾角的正切值[13]。坡度是描述地表形態及獲取其他地形因子的重要基礎數據，坡度的變化會引起一系列的相關變化。因此，研究不同DEM分辨率下坡度的變化情況具有一定的實際意義。

2.1.1 不同坡度占總面積比的分析

將坡度按0°～3°，3°～5°，5°～10°，10°～15°，15°～20°，20°～25°，25°以上進行分級，分別統計不同分辨率DEM生成的坡度分類圖中各坡度分級的面積和百分比，并繪制折線圖，如圖1所示。由圖1可以看出，當分辨率為5～25 m時，各坡度級別所占的比例基本上保持不變，當分辨率達到50 m時，3°～5°這一級所占的比例開始顯著的增加，當分辨率達到100 m時，0°～3°，10°～15°和15°～20°這幾級所占的比例都開始明顯降低，這說明，當DEM分辨率為25 m時，對提取的坡度值是穩定的。研究區大部分面積的坡度都在10°以下，地形比較平緩。

圖1 不同坡度占總面積的比隨 DEM分辨率變化的折線圖Fig.1 Change ratio of different slopes to total area with the DEM resolution

2.1.2 坡度最大值、最小值、平均值分析

不同分辨率DEM生成的坡度數據中，坡度最小值均為零，最大值和平均值隨分辨率的降低而減小，說明隨分辨率的粗略化，地形起伏更加概括，地勢整體趨于平坦。這種變化在分變率在5～50 m之間時并不明顯，但分辨率在50～300 m時，坡度的最大值和平均值變化情況都變得明顯。如圖2和圖3所示。

圖2 坡度最大值隨DEM分辨率變化的折線圖Fig.2 Change of the maximum value of the slope with the DEM resolution

圖3 坡度平均值隨DEM分辨率的變化Fig.3 Change of the average value of the slope with the DEM resolution

2.2 坡向分析

坡向是指斜坡的朝向。坡向不同，坡面所接受的太陽輻射不同會影響到濕度、溫度條件、坡面植被生長狀況、地面組成物質的風化速率、風化類型等[14]。

由不同分辨率的DEM生成坡向分布圖，并按北(0°～22.5°，337.5°～360°)，東北(22.5°～67.5°)，東(67.5°～112.5°)，東南(112.5°～157.5°)，南(157.5°～202.5°)，西南(202.5°～247.5°)，西(247.5°～292.5°)，西北(292.5°～337.5°)，平地(-1)進行重分類，計算不同坡向的面積，及占總面積的百分比，如圖4所示。從圖4可以看出，當分辨率為5～25 m時，各坡向所占比例沒有較大變化，當分辨率達到50 m時，平面所占比例急劇下降，其它各坡向也開始比較顯著地增加，說明當DEM分辨率為25 m時，對提取的坡向值是穩定的。主要原因可能是隨著分辨的降低，地貌細部特征被簡化，使總的平坦區面積比例減小，而其他各類坡向面積比例變化較慢。這可能是由于坡向是切平面上沿最大傾斜方向的某個矢量在水平面上的投影方向，而DEM 分辨率的粗略化使得對地形描述簡化，但對坡向影響不如坡度大。

圖4 不同坡向占總面積的比隨DEM分辨率的變化Fig.4 Change of the ratio of different slope aspects to total area with the DEM resolution

2.3 剖面曲率分析

剖面曲面是地面任意一點的地表坡度的變化率，或稱為高程變化的二次導數。剖面曲率作為重要的地形因子，能夠反映局部地形結構，在地表過程模擬、水土保持等領域的研究中有重要的價值[15]。剖面曲率的精度受到DEM分辨率的制約，如何平衡剖面曲率的精度與DEM的分辨率之間的矛盾關系，成為一個需要研究的問題。

利用不同分辨率的DEM數據分別生成剖面曲率分布圖。計算剖面曲率最大值、最小值、平均值以及標準差，來研究剖面曲率與DEM分辨率的關系，如圖5和圖6所示。從圖中可以看出，隨著DEM分辨率的降低，剖面曲率的最小值與平均值逐漸增大，最大值與標準差逐漸減小。當分辨率在5～25 m時，變化極其明顯，而當分辨率的值大于50 m時，曲線趨于平緩，這說明當分辨率的值大于50 m時，地形信息遺失很嚴重，這是因為隨著DEM分辨率的粗略化，對地形復雜程度的描述愈來愈簡單，地形變化趨于平緩。而當分辨率增大到一定值時，地形信息的遺失維持在一定的水平。這也說明，25 m是一個比較合適的分辨率值。

圖5 剖面曲率最大值最小值隨DEM 分辨率變化的折線圖Fig.5 Change of the minimum and maximum values of profile curvature with the DEM resolution

圖6 剖面曲率平均值標準差隨DEM 分辨率變化的折線圖Fig.6 Change of the standard deviations of average value of profile curvature with the DEM resolution

2.4 高程信息熵分析

信息熵表示信源輸出后，每個柵格所提供的地形信息量的總和[16]。隨著分辨率的改變，柵格的分布相應發生變化，其所包含的地形信息量也相應發生改變，信息熵是度量不同分辨率的DEM所含信息容量的重要指標之一[12]。利用信息熵公式，計算不同分辨率DEM的信息熵，繪制圖7所示的高程信息熵隨DEM分辨率變化的曲線。由圖7可見，隨著DEM分辨率的增大，高程信息熵逐漸增大，尤其是在分辨率由300 m變化到25 m這一階段信息熵增加較快，從25 m到5 m繼續增加，但速度有所放緩。因此可以認為，隨著DEM分辨率的增大，高程信息熵也逐漸增在，但25 m是個拐點。

圖7 高程熵隨DEM分辨率變化的折線圖Fig.7 Change of elevation entropy with the DEM resolution

3 結 論

DEM分辨率在地形因子的提取中起著非常重要的作用，一般情況下，分辨率越高，包含的信息容量也越多，相應地要占用更多的存儲空間，需要較長的運行時間[17]。本文利用大興安嶺漠河地區1∶5萬等高線數據，研究不同分辨率DEM在提取相關地形因子時的變化情況。從研究的結果來看，隨著DEM分辨率的逐漸增加，高程所包含的信息熵值也逐漸增大，說明分辨率越高，DEM中的信息容量也越大。通過對不同分辨率DEM提取的坡度、坡向和剖面曲率分析，當DEM分辨率在25 m左右時，各種指標逐漸趨于平穩，當DEM分辨率大于50 m時地形信息損失嚴重。因此，在對該地區進行地形因子的分析研究時，選擇DEM分辨率25 m比較合適。可以認為漠河地區由1∶5萬等高線生成DEM時，最合理的分辨率大約為25 m。

在確定DEM分辨率等級時，直接由25 m變動到50 m，跨度較大，最佳分辨值應該位于25～50 m，因此還需要更細致的研究進一步來確定具體的最佳分辨率值，以縮小范圍。近年來，許多研究者采用不同的方法，取得了關于最佳分辨率的研究成果，但是不同研究區域的地貌差異較大，所以研究成果是否具有通用性，還有等于進一步的研究。

[1]張欣方，黃敏兒.數字高程模型簡介[J].科技致富導向，2011(20)：95.

[2]王 妍，劉洪斌，寧茂歧，等.柵格DEM的水平分辨率對地形信息的影響分析[J].測繪與空間地理信息，2005，28(6)：41-47.

[3]Brasington J，Richards K.Interactions between Model Predictions，Parameters and DTM Scales for TOP-MODEL[J].Computer and Geosciences，1998，24(4)：299-314.

[4]Bolstad P V，Stowe T.An evaluation of DEM accuracy：Elevation，slope，and aspect [J].Photogrammetric Engineering and Remote Sensing.1994，60(11)：1327-1332.

[5]Brabb E E.Innovation approaches to landslide hazard and risk mapping[C]，Proc.of 4th ISL，Toronto，1984：307-323.

[6]Folrinsky I V.Combined analysis of digital terrain models and remotely sensed data in landscape investigations[J].Progress in Physical Geography，1998，22(l)：33-60.

[7]Folrinsky I V，Kuryakova G A.Determination of grid size for digital terrain modeling in landscape investigations exemplified by soil moisture distribution at a microscale[J].International Journal of Geographical Information Science，2000，14(8)：815-832.

[8]Kienzle S.The effect of DEM raster resolution on first order，second order and compound terrain derivatives[J].Transactions in GIS，2004，8(1)：83-111.

[9]Erskine R H，Green T R，Ramirez J A.Digital elevation accuracy and grid cell size：effects on estimated terrain attributes [J].Soil Science Society of America，2007，71(4)：1371-1380.

[10]王 峰，王春梅.地形因子與DEM分辨率關系的初步研究[J].水土保持研究，2009，16(4)：225-229.

[11]張 磊.基于核心地形因子的黃土地貌形態空間格局研究[D].南京：南京師范大學，2013.

[12]朱 偉，王東華，周曉光.基于信息熵的DEM最佳分辨率確定方法研究[J].遙感信息，2008(8)：79-82.

[13]湯國安，楊勤科，張 勇，等.不同比例尺DEM提取地面坡度的精度研究[J].水土保持通報，2001，21(1)：53-56.

[14]申成磊，馬勁松.柵格尺度對DEM分析中地形因子的影響[J].現代測繪，2007，30(6)：3-6.

[15]陳 楠.剖面曲率精度變化與DEM分辨率關系[J].中國礦業大學學報，2013，42(1)：147-151.

[16]張朝忙，劉慶生，劉高煥，等.中國地區SRTM3 DEM高程精度質量評價[J].測繪工程，2014，23(4)：14-19.

[17]安辛克，李琳琳，張運剛.不同分辨率SRTM-DEM的數字地形分析比較[J].測繪與空間地理信息，33(6)：73-75.

Study on the Effect of DEM Resolutions Extracting on the Landform Factors

Liu Na，Yang Binbin，Fang Jingnan，Hao Yi，Yang Haicheng，Wang Weifang

(College of Forestry，Northeast Forestry University，Harbin 150000)

The 1：50000 contours data of Mohe area Daxingan Mountains was used as the objective to generate different resolution Digital Elevation Model (DEM) and analyze the changes of slope，aspect and profile curvature under condition of the different resolution DEM.The change trend of DEM height-entropy with different resolution was also analyzed.The results showed that the changes of each topography factor became stable when the DEM resolution came to around 25 meters and the change became apparently intensive when the resolution came to over 50 meters.Height-entropy changed slowly when the resolution enlarged to 25 meters.It was concluded from above results that the most appropriate resolution should be 25 meters in the Mohe area when 1：50，000 contours was used to generate DEM.

Digital Elevation Model(DEM)；landform factors；resolutions

2016-05-06

國家級大學生創新訓練計劃項目(20150225067)

劉 娜，本科生。研究方向：地理信息系統。

*通信作者：王維芳，博士，副教授。研究方向：森林經理。E-mail：weifangwang@126.com

劉 娜，楊斌斌，房靖楠，等.DEM分辨率對地形因子提取影響的研究[J].森林工程，2016，32(6)：57-60.

P217

A

1001-005X(2016)06-0057-04