數字地形分析及特征提取基本地形因子計算

坡度、坡向、曲率、粗糙度

面積、體積

地表日照度地形特征提取

溝脊值提取

水文分析

可視性分析主要內容地形因子的分類地形因子微觀因子宏觀因子坡度坡向坡度變率坡向變率平面變率剖面變率剖長剖型因子地形粗糙度地形起伏度高程變異系數地表切割深度地形因子姿態因子坡度坡向坡形因子宏觀上的坡形因子坡面變率因子坡面曲率因子坡長因子徑流坡長整體坡長坡位因子正地形坡位（溝澗地）負地形坡位（溝坡地、溝底地）坡面復雜度因子地形起伏度地表粗糙度地表切割深度高程變異系數地形因子的分類提取地形因子的基本方法先將地形的形態特征或各個坡面因子進行定量化描述，完成求導的數學模型，在此基礎上建立其以DEM為基本信息源進行提取的技術路線，并通過軟件形成一套易于計算機操作的方法因子提取方法應遵循地貌形態學與地貌成因學的基本理論DEM數據在表象上首先反映的是地面的海拔高程信息，更深層次的地形信息需要通過對高程的空間分布特征的提取，或更高層次地形要素的解算獲得。更復雜的還需通過統計分析、數模構建獲得微觀因子通常以DEM格網數據的空間矢量模型為基礎。通過空間向量的差分運算完成；宏觀因子一般通過移動分析窗口的方法完成。信息源的尺度特征會影響坡面因子的提取精度。坡度（Slope）：過地表一點的切平面與水平面的夾角。坡向（Aspect）：地表面上一點的切平面的法線在水平面的投影與該點的正北方向的夾角。聯系與區別：都是點位參數，只有理論意義，不具備地理意義。坡度坡向描述描述地表面在該點的傾斜程度。描述該點高程值改變量的最大變化方向。作用影響地表物質流動與能量轉換的規模與強度，制約生產力空間布局。決定地表面局部地面接收陽光和重新分配太陽輻射量的重要地形因子，直接造成局部地區氣候特征差異，影響各項農業生產指標。SlopeAspectxZY基本地形因子計算——坡度、坡向坡度指標地形表現<3°平坦平原、盆地中央部分、寬淺谷地底部、臺面3°~5°山前地帶、山前傾斜平原、沖積、洪積扇、淺丘、崗地、臺地、谷地等5°~15°山麓地帶、盆地周圍、丘陵15°~25°一般在200到1500米的山地中25°~30°大于1000米山地坡面上部（接近山頂部分）30°~45°大于1500米山體坡面上部>45°地理意義的垂直面坡度地理意義坡向地理坡向俗稱0°±22.5°北N陰坡45°±22.5°東北NE半陰坡315°±22.5°西北NW90°±22.5°東E270°±22.5°西W135°±22.5°東南SE半陽坡225°±22.5°西南SW180°±22.5°南S陽坡不存在不存在平地注：適用于南半球，北半球則正好相反N=0°NE=45°E=90°SE=135°S=180°SW=225°W=270°NW=315°坡向地理意義e5e2e6e1ee3e8e4e7坡度：Slope為坡度，Aspect為坡向，Slopewe為X方向上的坡度，Slopesn為Y方向上的坡度。Slopewe、Slopesn的計算如下：坡向：方法1：方法2：基于規則格網的坡度、坡向計算坡度：坡向：SlopeAspect基于不規則三角網的坡度、坡向計算坡度、坡向應用1．坡度選擇。坡度和坡向對果樹生長發育有一定影響。同一坡向不同坡度，對溫、熱、水分都有不同程度的差異。如南坡，10度坡太陽直接輻射量為平地的116％，20度坡為130％。表土的含水量，5度坡為52．38％，20度坡為34．78％，隨坡度增大而降低。土壤凍結深度，5度坡在20厘米以上，15度坡則為5厘米，但坡下低洼地冷空氣沉積，坡頂則寒冷，不易栽培。選擇葡萄栽培的適宜坡度5度～20度的斜坡地為好，15度坡最為合適。2．坡向選擇。坡向不同，光照、濕度、熱量、風量也不同。一般南坡、東南坡、西南坡，所獲得太陽光熱量大。北坡、東北坡，西北坡，則較冷涼。南坡與北坡近地面20厘米處氣溫平均相差0．4℃，80厘米深土層，南坡比北坡地溫高4℃～5℃。葡萄喜光、喜溫，以選擇南坡為宜。但南坡溫、濕度變化較大，水分蒸發量大，融雪、解凍比北坡早，因此必須加強水土保持工程。由于山地地勢非常復雜，南、北方氣候差異懸殊，但在中緯度的低山區，北坡水分蒸發量少，土壤墑情好，植被密生，土質較肥沃，土層較深厚，也能栽培葡萄樹。葡萄種植的地形條件1、緯度和海拔：世界上大部分葡萄園分布在北緯20-52度之間及南緯30-45度之間，絕大部分在北半球。海拔高度一般在400-600米。中國葡萄多在北緯30-43度之間海拔的變化較大，約200-1000米，河北懷來葡萄分布高度達1100米，山西清徐達1200米，西藏山南地區達1500米以上。緯度和海拔是在大范圍內影響溫度和熱量的重要因素。

2、坡向和坡度：在大地形條件相似情況下，不同坡向的小氣候有明顯差異。通常以南向（包括正南向、西北向和東南向）的坡地受光熱較多，平日氣溫較高。坡地的增溫效應與其坡度密切相關。一般坡地向南每傾斜1度，相當于推進1緯度。受熱最多的坡地角度約為20-35度（在北緯40-50度范圍）。葡萄因較耐干旱和土壤瘠薄，可以在相對不大范圍內發育根系，所以比其它果樹更適宜在坡地上栽培，然而坡度越大水土流失越嚴重，因此，在種植葡萄時應優先考慮坡度在20-25度以下的土地。

3、水面的影響：海洋、湖泊、江河、水庫等大的水域，由于吸收的太陽輻射能量多，熱容量較大，白天和夏季的溫度比陸地低，而夜間和冬季的溫度比內陸高。因此，臨近水域沿岸的氣候比較溫和，無霜期較長。臨近大水面的葡萄園由于深水反射出大量的藍紫光和紫外線，漿果著色和品質好，所以選擇葡萄園時盡量靠近大的湖泊、河流與海洋的地方。

坡度、坡向應用坡度是評價耕地質量的主要指標，也是衡量土地利用是否合理的一個關鍵因子。借助地理信息系統、3D技術，準確分析耕地坡度信息，對土地調查具有重要意義。第二次全國土地調查耕地坡度分類圖制作流程讀入源數據按縣級行政區拼接是否剪切？分塊DEMDEM生成坡度坡度圖重新分類參數：坡度分級表坡度分級柵格圖轉換為矢量格式坡度分級矢量圖矢量數據編輯,去除細碎多邊形,邊界光滑二調：精度參數坡度分級圖土地利用圖斑數據耕地坡度分級圖圖幅整飾中間成果成果

定義：地形曲面在各個界面方向上的形狀，凹凸變化的反映，是平面點位的函數。反映了地形結構和形態。影響著土壤有機物含量的分布，在地表過程模擬、水文、土壤等領域有著重要的應用價值和意義。基本地形因子計算——曲率定義：反映地表起伏變化和侵蝕程度的指標。一般定義為地表單元的曲面面積與其在水平面上的投影面積之比。粗糙度是能夠反映地形的起伏變化和侵蝕程度的宏觀地形因子，是衡量地表侵蝕程度的重要量化指標。基本地形因子計算——粗糙度S為曲率，AreaS為格點面元趨勢面面積，AreaH為格點面元面積。粗糙度計算定義：實現步驟：實際應用時，當分析窗口為3*3時，可采用以下公式近似計算。1、根據DEM提取坡度因子S；2、粗糙度：R=1/cos(S)

基本規則格網的粗糙度計算由于基于不規則格網的地表面積計算較為簡單，可以將格網DEM的每個格網分解為三角形，從而簡化規則格網地表面積計算。計算三角形的表面積使用海倫公式：Di表示第i（1≤i≤3）對三角形兩頂點之間的表面距離，S表示三角形的表面積，P表示三角形周長的一半。整個DEM的表面積則是每個三角形表面積的累加。基本地形因子計算——地表面積投影面積指的是任意多邊形在水平面上的面積，主要有兩類計算方法：1、直接采用海倫公式進行計算，只要將表面積計算公式中的距離改為平面上兩點的距離即可；2、根據梯形法則，如果一個多邊形由順序排列的N個點（Xi，Yi，i=1，…，N）組成并且第N點與第1點相同，則水平投影面積計算公式為：如果多邊形頂點按順時針方向排列，則計算的面積值為負；反之，為正。基本地形因子計算——投影面積基于數字地面模型的體積計算實質上是計算地形表面和給定參考面之差，在實際應用中，這種差值計算往往被歸類于土方計算，應用在工程施工、水利工程規劃等方面。對DEM進行挖或填后，土方量可由原始DEM體積V0減去新的DEM體積V1，乘以相應的物質密度求得。

V=V0–V1當V>0時，表示挖方；當V<0時，表示填方；當V=0時，表示既不挖方也不填方。基本地形因子計算——體積基于規則格網的體積計算可以將格網單元視為平面，計算多個立方體的體積總和。

基本規則格網的體積計算由于每一個空間三角形與其在指定計算高程平面上的投影都形成一個三棱柱和一個三棱錐，所以基于不規則三角網的體積計算最終歸結為若干三棱柱和三棱錐的體積運算：對每個三棱柱計算其填、挖方量并加以累積，就可以得到最終結果。基本不規則三角網的體積計算土方計算結果

由于地面的起伏，地面各點所接受的太陽輻照度是不相同的，其計算方法為：el為太陽高度角，αz為太陽方位角，α為當前點的坡角，β為當前點的坡向。基本規則格網的地表日照度計算太陽高度角=90度太陽高度角=45度基本地形因子計算——地表日照度

地形特征提取地形特征點：山峰、谷底、鞍部點等地形特征線：山脊線、山谷線等水系特征提取分水線：類似山脊線合水線：類似山谷線基于格網DEM的特征提取基于圖像處理基于地形曲面幾何分析基于地形曲面流水物理模擬分析

2）和3）分析方法的綜合地形特征提取溝脊值提取方法：1、提取特征點；2、將提取的特征點按照某種方式連成線。基本規則格網的地表溝脊提取

水文分析與計算是對所研究的水文變量或過程，作出盡可能正確的概率描述，對防止水旱災害和開發、利用、保護水資源的工程或非工程措施的規劃、設計、施工以及管理運用有著重要的意義，也是DEM數據應用的一個重要方面，主要用于研究與地表水流有關的各種自然現象比如：洪水水位及泛濫情況或者可以劃定受污染源影響的地區、以及預測當改變某一地區的地貌時對整個地區造成的后果等。在城市和區域規劃、農業及森林等許多領域對地球表面形狀的理解很重要。這些領域需要知道水流怎樣經過某一地區，以及這個地區的地貌的改變會以什么樣的方式影響水流的流動。地形特征提取——水文分析水系Watershed流域(Basin,Catchment,

Contributingarea)WatershedBoundaries

(DrainageDivides)流域邊界匯流點PourPoints(Outlets)

DEM水平和垂向的分辨率、DEM生成過程的內插、以及內插精度和網格單元內高程信息取平均等原因造成DEM中存在凹陷洼地和平地。

洼地：區域地形的集水區域,洼地底點(谷底點)的高程通常小于其相鄰近點(至少八鄰域點)的高程。原始DEM數據處理后數據水文分析——數據預處理洼地填平算法——Moran和Vezina算法

首先用一個極大高程水面數據將原始地面的DEM數據表面淹沒，然后移除DEM上多余的水，最后得到的高程數據就是填洼處理的高程數據。Ａ、對DEM數據初始化，除了邊界外，用一極大臨界水面高程數據將原始地面上DEM數據表面淹沒。Ｂ、通過迭代方法將W值逐步減小，最后W會收斂于Wf,Wf代表的高程點數據既是填洼后的數據。平地墊高算法——Martz和Garbrecht算法

用高程增量疊加算法處理平地。對平地范圍內的單元格增加一微小增量，每個單元格的增量大小是不一樣的，就可以消除平地。數據預處理方法在經過填充洼地后的DEM（FilledDem），流水可以暢通無阻地流至區域地形的邊緣。在經過填充洼地后的DEM是流向分析的基礎預處理結果D8算法：A、如果最大落差值小于0，則賦以負值以表明此格網方向未定；B、如果最大落差值大于或等于0，且最大值只有一個，則將對應此最大值的方向值作為中心格網處的方向值；C、如果最大落差等于0，且有一個以上的0值，則以這些0值所對應的方向值相加。在極端情況下，如果8個鄰域高程值都與中心格網高程值相同，則中心格網方向值賦以255；D、如果最大落差值大于0，且有一個以上的最大值，則任選一個方向作為水流方向。被處理柵格單元同相鄰8個柵格單元之間坡降的算法為：

Slope=Dz/Di式中:Slope為兩個柵格之間的坡降；Dz為兩個柵格單元之間的高程差；Di為兩個柵格單元中心之間的距離。326412816K1842水文分析——水流方向78726971584974675649465069534437384864585522312468614721161974533412111222244822244811248412812812482214441111416326412816K1842原始DEM矩陣水流方向矩陣水流方向提取結果水流方向提取結果

區域地形每點的流水累積量，它可以用區域地形表面的流水模擬方法獲得。規則格網表示的數字地面高程模型每點處有一個單位的水量，按照自然水流從高處流到低處的自然規律，根據區域地形的水流方向數字矩陣計算每點處所流過的水量數值，便可以得到該區域匯流累積量。水文分析——匯流累積量22244822244811248412812812482214441111416326412816K1842水流方向矩陣匯流累積量矩陣匯流累積量提取結果000000011220037540000200100012400247352326412816K1842水流方向矩陣匯流累積量矩陣匯流累積量提取結果匯流累積量提取結果河網的提取采用的是地表徑流漫流模型，通過模擬地表徑流的流動來產生水系。A、在無洼地的DEM上計算出每一個柵格的水流方向矩陣；B、根據自然水流由高處流到低處的自然規律利用水流方向矩陣計算出匯流累積量；C、當匯流累積量達到一定值得時候，就會產生地表水流，所有匯流累積量大于臨界值的柵格就是潛在的水流路徑，由這些水流路徑構成的網絡就是河網。

在實現的過程中要確定匯流能力閾值，不低于該閾值的單元格標記為水系組成部分，該方法簡單，直接產生連續的流線段。水文分析——提取柵格河網單閾值>20單閾值>100提取柵格河網結果河流的分級一般依據河流的流量