李建飛

（遼寧省基礎(chǔ)測(cè)繪院，遼寧錦州 121003）

坡度圖是表示地面傾斜的地圖。主要用暈線或顏色在圖上直接表示出坡度的大小或陡緩，坡度傾角α的計(jì)算公式是：

式中，h表示高差，l表示水平距離。坡度數(shù)值的形式，通常用一個(gè)傾斜面對(duì)水平面間的夾角，即傾斜角度表示；也有用地面比降分?jǐn)?shù)式或百分率表示。

坡度圖的制作，通常利用DEM在地形圖上進(jìn)行坡度分級(jí)，其分級(jí)標(biāo)準(zhǔn)，多根據(jù)人們改造和利用自然實(shí)際需要的坡度的臨界極限，或各種地貌類型自然界限值進(jìn)行確定。坡度圖對(duì)農(nóng)業(yè)、林業(yè)、水利建設(shè)及軍事等方面均有重要使用價(jià)值。

一、坡度計(jì)算

1.坡度計(jì)算

坡度計(jì)算公式：

2.坡度計(jì)算模型

坡度計(jì)算模型主要有二階差分、三階不帶權(quán)差分、三階反距離平方權(quán)差分等。本文制作坡度圖應(yīng)用的是三階反距離平方權(quán)差分（圖1），即Horn算法：

圖1 DEM3×3局部移動(dòng)窗口

圖中G表示格網(wǎng)尺寸。e分別表示中心點(diǎn)e周圍格網(wǎng)點(diǎn)高程值。

3.坡度分級(jí)

按照耕地坡度分級(jí)要求，對(duì)耕地坡度分為五級(jí)，≤2°為1級(jí)，2°~6°為2級(jí)，6°~15°為3級(jí)，15°~25°為4級(jí)，＞25°為5級(jí)。

二、坡度圖制作方案

1.利用DEM制作坡度圖的技術(shù)方法與路線

DEM數(shù)據(jù)準(zhǔn)備→分幅DEM的鑲嵌→各縣區(qū)DEM數(shù)據(jù)裁切→柵格數(shù)據(jù)轉(zhuǎn)為坡度數(shù)據(jù)→坡度數(shù)據(jù)分級(jí)→分級(jí)后坡度數(shù)據(jù)轉(zhuǎn)為矢量數(shù)據(jù)→小面積圖斑融合→矢量數(shù)據(jù)后處理→質(zhì)量檢查與上交。

2.坡度圖制作

（1）DEM數(shù)據(jù)收集

遼寧省1:1萬DEM數(shù)據(jù)覆蓋遼寧省大部分地區(qū)，均已通過遼寧省測(cè)繪產(chǎn)品質(zhì)量監(jiān)督檢驗(yàn)站驗(yàn)收，作為全省縣級(jí)坡度圖生產(chǎn)制作的主要數(shù)據(jù)源。漏洞部分采用現(xiàn)有最新的1:5萬DEM數(shù)據(jù)作為補(bǔ)充數(shù)據(jù)源。

（2）DEM數(shù)據(jù)鑲嵌

本文DEM數(shù)據(jù)為遼寧省第三次全國(guó)國(guó)土調(diào)查生產(chǎn)的比例尺為1:1萬的數(shù)字高程模型數(shù)據(jù)，該數(shù)據(jù)分為40，41，42三個(gè)投影帶，中央經(jīng)線為120°，123°，126°。按照高斯-克呂格投影規(guī)則，遼寧省部分縣區(qū)涉及到跨帶，涉及跨帶的縣區(qū)，將中央經(jīng)線轉(zhuǎn)換成該縣區(qū)制作“三調(diào)”數(shù)據(jù)庫的投影帶。

（3）DEM數(shù)據(jù)轉(zhuǎn)為坡度數(shù)據(jù)

基于Arcgis10.3.1平臺(tái)下Spatial Analyst工具中的表面分析中的坡度工具，將合并后的DEM數(shù)據(jù)轉(zhuǎn)換成坡度柵格數(shù)據(jù)（圖3），其中像元大小定為5米。

圖3 坡度柵格數(shù)據(jù)

（4）坡度數(shù)據(jù)分級(jí)

將坡度柵格數(shù)據(jù)進(jìn)行分類，利用Arcgis10.3.1中3D Analyst工具中柵格重分類工具，本文選擇手動(dòng)分類的方法，類別為5類，分類標(biāo)準(zhǔn)按照《利用DEM確定耕地坡度分級(jí)技術(shù)規(guī)定》（圖4）。

圖4 坡度分級(jí)柵格數(shù)據(jù)

（5）柵格數(shù)據(jù)轉(zhuǎn)矢量數(shù)據(jù)

這一步利用Arcgis10.3.1中轉(zhuǎn)換工具下的由柵格轉(zhuǎn)出中柵格轉(zhuǎn)面工具實(shí)現(xiàn)，生成坡度分級(jí)矢量數(shù)據(jù)（圖5）。

圖5 坡度矢量數(shù)據(jù)

（6）消除小面積圖斑

本文所選用的消除小面積圖斑的方法為融合，將小于750m2的小面積圖斑進(jìn)按坡度級(jí)就低不就高原則并入鄰近圖斑，合并的次序是每次都要從此時(shí)最小面積圖斑開始進(jìn)行判斷，直到達(dá)到一定面積以下的小面積圖斑全部融合掉為止（圖6）。

圖6 融合后矢量數(shù)據(jù)

（7）矢量數(shù)據(jù)后處理

矢量數(shù)據(jù)后處理包括界線平滑、拓?fù)渲亟?、?shù)據(jù)裁切等處理，數(shù)據(jù)裁切按各縣區(qū)界線外擴(kuò)1km進(jìn)行裁切（圖7）。

圖7 坡度圖

3.存在問題及解決辦法

在坡度圖接邊時(shí)筆者遇到了嚴(yán)重的接邊問題，兩個(gè)相鄰縣區(qū)接邊出現(xiàn)錯(cuò)牙現(xiàn)象，導(dǎo)致接邊接不上（圖8），經(jīng)分析實(shí)驗(yàn)發(fā)現(xiàn)了制作過程中產(chǎn)生問題的原因，在坡度圖制作技術(shù)路線與方法中DEM數(shù)據(jù)的裁切是在DEM鑲嵌后就進(jìn)行的，然后用裁切過的柵格數(shù)據(jù)轉(zhuǎn)坡度數(shù)據(jù)、矢量數(shù)據(jù)，會(huì)導(dǎo)致矢量數(shù)據(jù)偏離正確的位置。也就是說不能用裁剪后的柵格數(shù)據(jù)制作坡度圖，需要按縣區(qū)界裁剪的，必須把這一步放到最后一步做，即不能裁剪柵格數(shù)據(jù)。

圖8 坡度圖接邊問題

三、結(jié)語

本文對(duì)土地坡度圖制作全過程進(jìn)行了富有成效的探索研究，在很多地方采用了非常好的技術(shù)處理和技術(shù)實(shí)現(xiàn)手段，在保證操作準(zhǔn)確處理基礎(chǔ)上，使運(yùn)行快捷、高效。坡度圖是遼寧省第三次全國(guó)國(guó)土調(diào)查的基礎(chǔ)數(shù)據(jù)和重要成果，是土地調(diào)查數(shù)據(jù)庫建設(shè)的重要組成部分，在未來該成果將在諸多相關(guān)領(lǐng)域取得更高更廣的應(yīng)用前景。