地形因子

第七章

1、本章主题编号

2、本章内容概述

（1）概述

● 坡面因子的分类及提取方法

● 确定坡面因子提取的算法基础

● 提取坡面因子的常用分析窗口

（2）坡度、坡向

● 坡度的提取

● 坡向的提取

（3）坡形

● 宏观坡形因子

● 地面曲率因子

● 地面变率因子

（4）坡长

（5）坡位

（6）坡面复杂度因子

3、本章内容

3.1 概述

（1）坡面因子的分类及提取方法

● 坡面因子的分类

按照坡面因子所描述的空间区域范围，可以将坡面因子划分为微观坡面因子与宏观坡面因子两种基本类型。常用的微观坡面因子主要有：坡度、坡向、坡长、坡度变率、坡向变率、平面曲率、剖面曲率等。常用的宏观坡面因子主要有：地形粗糙度、地形起伏度、高程变异系数、地表切割深度，以及宏观坡形因子（直线形斜坡、凸形斜坡、凹形斜坡、台阶形斜坡）等。

按照提取坡面因子差分计算的阶数，可以将坡面因子分为一阶坡面因子、二

阶坡面因子和高阶坡面因子。一阶坡面地形因子主要有坡度和坡向因子。二阶坡面因子主要有坡度变率、坡向变率、平面曲率、剖面曲率等因子。复合坡面因子有坡长、坡形因子、地形粗糙度、地形起伏度、高程变异系数和地表切割深度等。

按照坡面的形态特征，可将坡面因子进一步划分为：坡面姿态因子，坡形因子，坡位因子，坡长因子以及坡面复杂度因子五大类。

● 提取坡面因子的基本方法

首先将坡面的形态特征或各个坡面因子进行定量化描述，完成求导的数学模型，在此基础上，建立其以DEM为基本信息源进行提取的技术路线，并通过软件实现形成一套易于计算机操作的方法。

（2）确定坡面因子提取的算法基础

● DEM格网数据的空间矢量表达（如图7.1）

图7.1 DEM格网数据的空间矢量模型

● 基于空间矢量模型的差分计算

算法主要有数值分析方法、局部曲面拟合算法、空间矢量法、快速傅立叶变换等。其中数值分析方法包含有简单差分算法、二阶差分、三阶差分（带权或不带权）和Frame差分；局部曲面拟合又有线性回归平面、二次曲面和不完全四次曲面（据刘学军，2002）。

（3）提取坡面因子的常用分析窗口

● 窗口分析（领域分析）的基本原理是：对栅格数据系统中的一个、多个栅格点或全部数据，开辟一个有固定分析半径的分析窗口，并在该窗口内进行诸如极值、均值、标准差等一系列统计计算，或进行差分及与其它层面信息的复合分析等，实现栅格数据有效的水平方向扩展分析。

● 在坡面信息提取中，按照分析窗口的形状，可以将分析窗口划分为以下几类：

矩形窗口：以目标栅格为中心，分别向周围八个方向扩展一层或多层栅格。

圆形窗口：以目标栅格为中心，向周围作一等距离搜索区，构成一圆形分析窗口。

环形窗口：以目标栅格为中心，按指定的内外半径构成环形分析窗口。

扇形窗口：以目标栅格为中心，按指定的起始和终止角度构成扇形分析窗口。

矩形窗口最为常用，一般采用3×3位基本分析窗口，然而，按照分析的需要，分析窗口也可以扩大为5×5、7×7 或更大。

3.2 坡度、坡向

● 坡面姿态（坡度及坡向）是指局部地表坡面在空间的倾斜程度和朝向。

● 坡度表示了该局部地表坡面的倾斜程度，坡度大小直接影响着地表物质流动与能量转换的规模与强度，是制约生产力空间布局的重要因子。

● 坡向是决定地表面局部地面接收阳光和重新分配太阳辐射量的重要地形因子之一，直接造成局部地区气候特征的差异，同时，也直接影响到诸如土壤水分、地面无霜期以及作物生长适宜性程度等多项重要的农业生产指标。

（1）坡度的提取

● 严格地讲，地表面任一点的坡度是指过该点的切平面与水平地面的夹角。坡度表示了地表面在该点的倾斜程度，在数值上等于过该点的地表微分单元的法

矢量与z轴的夹角（如图7.2所示），即：

Slope

= （7.1）

图7.2 地表单元坡度示意图

● 基于DEM的坡度提取通常在3×3的DEM栅格分析窗口中，采用几何平面来拟合或差分计算的方法进行。分析窗口在DEM数据矩阵中连续移动完成整个区域的计算工作。

（2）坡向的提取

● 坡向定义为：地表面上一点的切平面的法线矢量在水平面的投影与

过该点的正北方向的夹角（如表7.1中的坡向示意图所示，x轴为正北方向）。其数学表达公式

为：（7.2）

● 对于地面任何一点来说，坡向表征了该点高程值改变量的最大变化方向。坡向值有如下规定：正北方向为0度，按顺时针方向计算，取值范围为0°～360°。

● 坡向可在DEM数据中用式7.2直接提取。但应注意，由于式7.2求出坡向有与x轴正向和x轴负向夹角之分，此时就要根据fx和fy的符号来进一步确定坡向值（如表7.1所示）。

表7.1 坡向值的判断

α =

注：上述情况假定所建立的DEM数据从南向北获取的，且x轴与正北方向重合，否则上述公式求得的坡向值，还应加上x轴偏离正北方向的夹角值。

3.3 坡形

● 坡形是指局部地表坡面的曲折状态。宏观上讲，一般可分为直线形斜坡、凸形斜坡、凹形斜坡和台阶形斜坡四种基本类型。从微观角度上，一般可采用地面曲率因子和地面变率因子度量地面表面一点的弯曲变化程度。

（1）宏观坡形因子

● 直线形斜坡：从分水岭到斜坡底部地面坡度基本上不变。

● 凸形斜坡：地面坡度随着距分水岭距离增加而增加。邻近分水岭附近的地面平缓，以后随坡长的增加，坡度亦增加。

● 凹形斜坡：斜坡上半部坡度较陡，下半部坡度较缓。此种坡形常以沉积为主，较多分别在山区与阶地平原接壤处或河谷的两岸。

● 台阶形斜坡：台阶形斜坡是斜坡与阶地相间的复式，可以看作是凸形坡与凹形坡的组合。

（2）地面曲率因子

● 地面曲率是对地形表面一点扭曲变化程度的定量化度量因子，地面曲率在垂直和水平两个方向上分量分别称为平面曲率和剖面曲率。

● 剖面曲率是对地面坡度的沿最大坡降方向地面高程变化率的度量。数学表达式为：

（7.3）

● 平面曲率指在地形表面上，具体到任何一点P，指用过该点的水平面沿水平方向切地形表面所得的曲线在该点的曲率值。平面曲率描述的是地表曲面沿水平方向的弯曲、变化情况，也就是该点所在的微小范围内坡向变化程度的度量。数学表达式为：

（7.4）

● 曲率因子的提取算法的基本原理为：在DEM数据的基础上，根据其离散的高程数值，把地表模拟成一个连续的曲面，从微分几何的思想出发，模拟曲面上每一点所处的垂直于和平行于水平面的曲线，利用曲线曲率的求算方法的推导得出各个曲率因子的计算公式。

（3）地面变率因子

● 地面变率描述的是地表局部范围内坡度、坡向两个基本的地形指标的变化情况，它包括坡度变率、坡向变率两个基本因子。

● 地面坡度变率，是地面坡度在微分空间的变化率，是依据坡度的求算原理，在所提取的坡度值的基础上对地面每一点再求算一次坡度。即坡度之坡度（Slope of Slope, 简称SOS）。

● 地面坡向变率，是指在地表的坡向提取基础之上，进行对坡向变化率值的二次提取，亦即坡向之坡度（Slope of Aspect, SOA）。地面坡向变率在所提取的地表坡向矩阵的基础上沿袭坡度的求算原理，提取地表局部微小范围内坡向的最大变化情况。

3.4 坡长

● 坡长通常是指在地面上一点沿水流方向到其流向起点间的最大地面距离在水平面上的投影长度。其数学表达为：