高一平

（太原理工大學 測繪科學與技術系，山西 太原 030024）

基于DEM的區(qū)域地形坡度分級圖制作

高一平

（太原理工大學 測繪科學與技術系，山西 太原 030024）

地形坡度分級是土地管理、地理信息服務的重要內容。基于DEM數據，按照設定的坡度分級標準，探討了地形坡度分級的數據內插、坡度分級生成分級多邊形合并方法。利用SRTM數據、GIS平臺，以太原市為實驗樣區(qū)，通過內插得到（30m×30m）數字高程模型?；跀底指叱棠Ｐ?，提取不同分辨率的太原市地形坡度信息，以90m分辨率的坡度分級圖為例，合并小圖斑到相鄰圖斑中，完成地形坡度分級圖。

SRTM；坡度分級；圖斑；合并

經過半個多世紀的發(fā)展，DEM在GIS中得到廣泛應用，成為地理信息系統(tǒng)中重要組成部分。DEM的應用主要有可視化地形圖制作、數字地形分析、以及DEM與其他地理要素的數據疊加分析。DEM是地球表面地形地貌的模擬數字表達，其中地形分析是以地形數字圖像資料、以及詳細野外地質調查等為研究基礎，進行地形地貌分析研究，主要包括高程剖面分析、坡度、坡向分析流域分析等，其實質是進行數值計算的過程。

1 試驗區(qū)概況及資料

1）數據來源：本次研究的柵格高程來自于SRTM3[1]數據，該高程數據覆蓋中國全境。其水平分辨率為90 m，高程基準是EGM-96的大地水準面，平面基準是WGS-84，標準絕對高程精度±16 m，絕對平面精度±20m。經太原市市區(qū)界線切割出太原市數字高程數據。

2）實驗區(qū)地理位置：太原市位于山西省境中央，全市整體地形北高南低、呈簸箕形。地理坐標為東經111°30′～113°09′，北緯37°27′～38°25′，市內最高海拔為2670m，最低點為760m。太原盆地的的北端，位于華北地區(qū)黃河流域中部，西、北、東三面環(huán)山；中、南部為河谷平原區(qū)域，輪廓呈蝙蝠形東區(qū)。

2 主要的技術流程

2.1 實驗區(qū)高程專題圖

SRTM的原始數據即高程數據，是DEM最基本的信息，可按實驗區(qū)實際，將高程分為5個級別。實驗數據以原始SRTM(90 m×90 m)和通過最鄰近內插方法得到(30 m×30 m)數據為例，圖1和圖2為高程分級結果圖，其高程分級按照自然間斷點分級法取整分級。

2.2 坡度計算

地面某點坡度(slope)是過該點的切平面與水平地面的夾角，是高度變化的最大值比率，是對地面傾斜程度的定量描述，它是重要的地形因子之一。常見的坡度計算方法有：空間矢量分析法、擬合平面法、曲面擬合法、直接解法等，其中為最佳的方法為曲面擬合法[2,3](窗口微分分析法)見圖3，其原理是通過在3×3的DEM柵格窗口，中心像元及其相鄰的八個像元的值確定水平增量(dz/dx)和垂直增量(dz/dy)。這些相鄰的像元使用字母a至i標識，其中e表示當前正在計算坡向的像元。計算坡度的基本算法是:tan(p)=√（dz/dx）2+（dz/dy）2.像元e在x方向上的變化率使用以下算法計算：[dz/dx]=((c+2f+i)-（a+2d+g）/（8 *x_cellsize）.像元e在y方向上的變化率使用以下算法計算:

圖1 太原市高程圖(90 m×90 m)

圖2 太原市高程圖(30 m×30 m)

其中，dz/dx、dz/dy分別表示x，y方向的偏導數，P為坡度。坡度值越小，地形越平坦；坡度值越大，地形越陡。

圖3 DEM3×3移動窗口

2.3 坡度分級方法的選擇

根據不同的目的，完成坡度信息提取的基礎上，進行合理的坡度分級，是完成最終的坡度專題圖的又一關鍵。湯國安[4]對各種分級方法比較后歸納為三種：主管分級法、臨界分級法、模式分級法；其中主管分析法有較強的主觀性和隨機性；臨界分級法需要大量的實驗數據基礎；模式分級法中常見的有手動分類、等間距、分位數、等面積、標準差、自然斷裂點等分級方法。

本次實驗是以坡耕地坡度分級標準[5]為依據，將坡度范圍分為五級，即小于等于2。為1級，大于2。小于等于6°為2級，大于6。小于等于15°的為3級，大于15°小于等于25°為4級，大于25°的為5級。見圖4和圖5。

由圖4和圖5看出，通過最鄰近內插方法插值后生成的坡度圖，其圖斑的破碎程度大大增加，坡度變化的相對連續(xù)性也下降。

這是因為曲面擬合法(窗口微分分析法)是利用中心像元與周圍八個像元的高程來計算坡度。分辨率高的DEM可以反映地表坡度的細微變化，分辨率低的主要是反映整體的坡度變化。

圖4 坡度分級圖(90m×90m)

圖5 坡度分級圖(30m×30m)

表1 不同分辨率的DEM生成的平均坡度

對坡度進行分級后，不可避免會出現一些破碎圖斑，需對其進行處理。根據不同應用要求，設置不同閾值，只有面積達到一定閾值時才被單獨提取。此次試驗以90m分辨率的坡度圖為例，按實地面積不小于75000m2的要求進行綜合取舍，將圖斑單元小于9個的圖斑并入臨近圖斑。具體實現方法為：基于以上坡度分級圖借助ArcGIS中柵格計算功能，對連通性分析后面積仍小于75 000 m2（約9個柵格像元數）融合到周圍大圖斑中。見圖6。

圖6 圖斑合并后的坡度圖

2.4 地形坡度分級圖

ArcGIS中，完成柵格轉矢量的操作后，還有可能出現碎小圖斑，計算圖斑面積，將小于閾值的圖斑提取，利用ArcGIS中Eliminate命令，將小于75 000 m2小圖斑合并在相鄰面積最大的圖斑中。

自動生成矢量圖的過程中，圖斑邊界一般均存在非圓滑的折線，對制圖而言，為了使區(qū)域表現更加細膩，需對圖斑邊界進行平滑處理。ArcGIS中基于paek算法、200 m容差限值平滑之后，平滑前后的對比圖見圖7。

圖7 平滑前后的對比圖

3 結束語

（1）利用SRTM-DEM（90 m）對太原區(qū)域地形坡度進行提取，坡度相對集中在25。以下，通過設定閾值合并碎小圖斑，有效減少圖斑的數量；但此次試驗的基礎數據為SRTM-DEM數據，其原始數據分辨率較大，存在一定的數據精度問題；若能采用數據精度較高的高程數據進行地形坡度分級，會得到更客觀的成果。（2）地形坡度信息是影響地表物質流動與能量轉換的規(guī)模和強度重要因素，也是制約生產力空間布局的重要因素。地形坡度信息與其他專題圖進行疊加，對不同領域進行應用研究及定量化分析均有重要意義。

[1] 汪凌.美國航天飛機雷達地形測繪使命簡介[J].測繪通報,2000,12(3).38-40.

[2]陳楠,王欽敏,湯國安.基于DEM的坡向提取算法對比分析[J].遙感應用,2007(1):70-75.

[3]陳楠，王欽敏，湯國安.自DEM由不同算法提取坡度的對比分析[J].測繪工程，2006，15(1)：6-9.

[4]湯國安，宋佳.基于DEM坡度圖制圖中坡度分級方法的比較研究[J].水土保持學報，2006，20(2)：157-192.

[5]TDT1014--2007第二次全國土地調查技術規(guī)程[S].

Regional Terrain Slope Grading Map based on the DEMData

GAO Yi-ping
(Department of Surveying and Mapping,Taiyuan University of Technology,Taiyuan Shanxi 030024)

Terrain slope grading is essential for land management and geographic information service. Based on the DEMdata and given slope grading criteria,the paper discusses interpolation,slope grading generation,and polygons merging method.With SRTMdata and GIS platform,taking Taiyuan area as a sample,Digital Elevation Model(30m*30m)is achieved through the interpolation.On the Digital Elevation Model,terrain slope information of Taiyuan with different resolutions is extracted.Taking the slope grading map (resolution is90m)as an example,small polygons are merged with adjacent polygons so that the grading map is accomplished.

SRTM;slope grading;polygons;merge

P285

A

1672-5050（2012）09-0044-03

2012-05-20

高一平（1985—），女，內蒙古巴彥淖爾人，碩士研究生，從事空間數據采集方法和數據處理的研究工作。

劉新光