地形因子计算详解
DEM坡面地形因子提取与分析
DEM坡面地形因子提取与分析
一、实验目的
了解基于DEM坡面地形因子提取的原理；掌握坡度、坡向、坡面曲率因子的提取方法及坡度分级图的制作；能够利用坡面地形因子与其它空间分析方法相结合以解决实际应用问题。
二、实验数据
一幅5m分辨率的黄土地貌DEM数据
三、实验步骤
3.1坡度
1、提取坡度
图1、坡度数据
2、3度等间距分级
3、土地利用模式分级
3.2坡向
图4、坡向图
3.3提取地面变率因子
1、坡度变率（SOS）
2、坡向变率（SOA）
3.4提取DEM层的反地形；
图7、DEM反地形提取
3.5DEM的平均曲率
图9、平均曲率计算
图10、平均曲率图与光照晕渲图叠加
四、思考
对于平均地形因子的提取，我认为可以通过邻域分析，计算窗口的最大最小值，求其平
均值，再不断扩大，知道窗口扩大到整个ＤＥＭ。
实验五DEM坡面地形因子提取
实验五DEM坡面地形因子提取实验目的：通过数字高程模型（DEM）数据提取坡度和坡向地形因子，以分析地形特征对水文过程和土地利用分布的影响。
实验步骤：1.数据准备a) 获取高分辨率的地形DEM数据，可以选择使用Lidar数据或者采用其他方式获取DEM数据。
b)进行数据预处理，包拟合DEM数据，去除噪声和突出值等。
2.坡度计算a)在DEM上采样，计算每个像元上的坡度。
b)坡度计算可以通过以下公式进行计算：Slope(i,j) = arctan(sqrt((dz/dx)^2 + (dz/dy)^2))其中，Slope(i,j)代表坡度， dz/dx代表DEM在x方向的梯度，dz/dy代表DEM在y方向的梯度。
3.坡向计算a)在DEM上采样，计算每个像元上的坡向。
b)坡向计算可以通过以下公式进行计算：Aspect(i,j) = arctan(dz/dx / dz/dy)其中，Aspect(i,j)代表坡向， dz/dx代表DEM在x方向的梯度，dz/dy代表DEM在y方向的梯度。
4.地形指数计算a)根据坡度和坡向的计算结果，可以进一步计算其他地形指数，例如地形湿度、地形开阔度等。
b)地形湿度可以通过计算每个像元周围的流通路径长度来估算。
c)地形开阔度可以通过计算每个像元周围的可见面积来估算。
5.结果分析a)可视化坡度和坡向地形因子，以了解地形特征。
b)利用地形指数，可以分析地形特征对水文过程和土地利用分布的影响。
实验结果分析：通过提取DEM的坡度和坡向地形因子，可以分析出地形特征，进而对水文过程和土地利用分布进行预测和分析。
例如，通过分析坡度可以了解一个地区的地势起伏程度，从而对洪水灾害的发生概率进行预测。
通过分析坡向可以了解水流在地表的流向，从而对土壤侵蚀和水资源分布进行预测。
此外，通过计算其他地形指数，还可以分析地形湿度和地形开阔度对生态环境的影响，为环境管理和规划提供数据支持。
总结：本实验通过DEM数据的处理和分析，提取了坡度和坡向地形因子，并通过计算其他地形指数，以分析地形特征对水文过程和土地利用分布的影响。
6 地形分析坡面因子提取
3.2 山脊线、山谷线的提取
求出已提取的概略地形特征线与DEM 格网 线的交点，在该交点附近的一个小区域， 对DEM 数据进行几何分析，即找出该区域 内与概略的地形特征线正交方向地形断面 上高程变化的极值点，该点即为该条地形 特征线的精确位置。
3.2 山脊线、山谷线的提取
平面曲率与坡位组合法 : 利用DEM数据提取地面的平面曲率及地面 的正负地形，取正地形上平面曲率的大值 即为山脊，负地形上平面曲率的大值为山 谷。 提取的山谷和山脊可以通过曲率的值 来进行调节。
1，表示谷点 1，表示脊点 VR i，j 2，表示鞍点 0，表示其他点
(i+1,j-1)
(i+1,j)
(i+1,j+1)
差分算法示意图
山顶点
鞍பைடு நூலகம்点
图例
等高线
山顶点
鞍部
利用ArcView GIS 软件及DEM数据提取的山顶、鞍部
3.2 山脊线、山谷线的提取
山谷线 山脊线与山谷线
2.3 坡面复杂度因子
坡面复杂因子是宏观的地形信息因子，包括地 形起伏度、地形粗糙度、地表切割深度和沟壑 密度等 ；
地形起伏度：是在所指定的分析区域内所有栅 格中最大高程与最小高程的差，
RFi H max H min
.
中国1：100万DEM 提取的中国陆地区域地形起伏度图
2.3 坡面复杂度因子
3.2 山脊线、山谷线的提取
基本思想 首先用较稀的DEM格网 数据用地形流水物理分 析方法提取区域内概略 的地形特征线，然后用 其引导，在其周围邻近 区域对地形进行几何分 析来精确确定区域的地 形特征线。
概略DEM建立 地形流水物理模拟
地形因子计算详解
第七章1、本章主题编号2、本章内容概述（1）概述●坡面因子的分类及提取方法●确定坡面因子提取的算法基础●提取坡面因子的常用分析窗口（2）坡度、坡向●坡度的提取●坡向的提取（3）坡形●宏观坡形因子●地面曲率因子●地面变率因子（4）坡长（5）坡位（6）坡面复杂度因子3、本章内容3.1概述（1）坡面因子的分类及提取方法●坡面因子的分类按照坡面因子所描述的空间区域范围，可以将坡面因子划分为微观坡面因子与宏观坡面因子两种基本类型。
常用的微观坡面因子主要有：坡度、坡向、坡长、坡度变率、坡向变率、平面曲率、剖面曲率等。
常用的宏观坡面因子主要有：地形粗糙度、地形起伏度、高程变异系数、地表切割深度，以及宏观坡形因子（直线形斜坡、凸形斜坡、凹形斜坡、台阶形斜坡）等。
按照提取坡面因子差分计算的阶数，可以将坡面因子分为一阶坡面因子、二阶坡面因子和高阶坡面因子。
一阶坡面地形因子主要有坡度和坡向因子。
二阶坡面因子主要有坡度变率、坡向变率、平面曲率、剖面曲率等因子。
复合坡面因子有坡长、坡形因子、地形粗糙度、地形起伏度、高程变异系数和地表切割深度等。
按照坡面的形态特征，可将坡面因子进一步划分为：坡面姿态因子，坡形因子，坡位因子，坡长因子以及坡面复杂度因子五大类。
●提取坡面因子的基本方法首先将坡面的形态特征或各个坡面因子进行定量化描述，完成求导的数学模型，在此基础上，建立其以DEM为基本信息源进行提取的技术路线，并通过软件实现形成一套易于计算机操作的方法。
（2）确定坡面因子提取的算法基础● DEM格网数据的空间矢量表达（如图7.1）图7.1 DEM格网数据的空间矢量模型●基于空间矢量模型的差分计算算法主要有数值分析方法、局部曲面拟合算法、空间矢量法、快速傅立叶变换等。
其中数值分析方法包含有简单差分算法、二阶差分、三阶差分（带权或不带权）和Frame差分；局部曲面拟合又有线性回归平面、二次曲面和不完全四次曲面（据刘学军，2002）。
ArcGIS 第7章 坡面地形因子提取
区域水土流失地形因子定量指标的制定
本章结束！
• L指坡长，m指地表面沿流向的水流长度， • θ指水流地区的地面坡度值。
5. 坡 位
• ������ 坡位是指坡面所处的地貌部位。 • ������ 正地形、负地形 • ������ 沟间地(沟壑区<25⁰，塬墚区<15⁰) • 沟谷地 • 沟坡地
6. 坡面复杂度因子
地表粗糙度
• ������ R = S曲面/ S水平 R = 1/cos(S)
坡面因子分类体系
DEM误差分类体系
坡面因子分类体系
提取坡面因子的基本方法
• • • • ������ ������ ������ ������ 明确各个坡面因子的数字特征 建立解译模型 研究基于DEM的提取方法 软件实现
坡面因子提取的算法基础
• ������ DEM格网数据的空间矢量表达
坡面因子提取的算法基础
地面曲率因子
• 地面曲率是对地形表面一点扭曲变化程度的定量化度量因 子，地面曲率在垂直和水平两个方向上的分量分别称为 • 平面曲率和剖面曲率
地面变率因子
• ������ • ������ 坡度变率（SOS） 坡向变率（SOA）
4. 坡 长
• 坡长通常是指在地面上一点沿水流方向到其流向 起点间的最大地面距离在水平面上的投影长度。
简化的差分公式
• fx是X方向高程变化率，fy是Y方向高程变化率
坡 向
坡向的计算
3. 坡 形
指局部地表坡面的曲折程度，宏观上讲，可分为 直线形斜坡、凸形斜坡、凹形斜坡和台阶形斜坡。 从微观角度上看， 一般可采用地面曲 率因子和地面变率 因子度量地面一点 的弯曲变化程度。
• 宏观坡形因子
DEM坡面地形因子提取技术文档
DEM坡面地形因子提取技术文档一、引言坡面地形因子（Terrain Factors）是描述地形地貌特征的一种指数，它在地质、地形、水文及环境科学研究中扮演着重要的角色。
地形因子通常由数字高程模型（DEM）数据中提取而得，其中包括坡度、坡向、高程等。
在本文中，我们将介绍一种提取DEM坡面地形因子的技术。
二、技术原理1.DEM数据预处理首先，需要对DEM数据进行预处理。
预处理包括裁剪、填充、平滑等操作，以去除无效数据和噪声干扰。
这样可以得到一份清洁、准确的DEM数据供后续分析使用。
2.坡度计算坡度是地形表面在一个给定点处的曲率。
坡度可以通过计算DEM中两个相邻像元之间的高度差来获得。
大致可以使用以下公式计算坡度：坡度= arctan(√(（∂z/∂x）^2 + (∂z/∂y)^2))其中，z是DEM中其中一像元的高程，x和y是该像元与其相邻像元的水平位置。
通过计算所有像元的坡度，即可获得整个地形表面的坡度分布。
3.坡向计算坡向是地表倾斜的指向，即地面水流流向的方向。
坡向可以通过计算DEM中每个像元的局部水平面斜率及其方向来获得。
常用的计算方法有以下两种：-最大坡向：将DEM视为一个等高线，计算累积坡度最大的方向作为坡向。
-朗巴特坡向：根据DEM的高程变化来计算坡向。
该方法利用光学效应的原理，将DEM分成若干小块，分别计算每个块中的坡向，再通过插值方法将坡向合并为整体。
4.高程计算高程是地表在垂直方向上的绝对高度。
在DEM数据中，高程信息已经包含在每个像元的值中。
因此，只需简单地读取DEM数据中的高程值即可获得地形表面的高程分布。
三、技术流程1.获得并预处理DEM数据，去除无效数据和噪声干扰。
2.计算坡度：计算DEM中每个像元的坡度值。
3.计算坡向：根据所选择的坡向计算方法，计算DEM中每个像元的坡向值。
4.计算高程：读取DEM数据中每个像元的高程值。
四、技术应用坡面地形因子对地质、地形、水文及环境科学研究具有广泛的应用。
江西省地形因子(LS)分析
4）.得到洼池填平图
• 3.水流方向计算
• 操作步骤： • 1）打开Arctoolbox，在Arctoolbox列表框中 点击Spatial Analysis Tool\Hydrology，双击 Flow Direction，打开Flow Direction操作的对 话框。 • 2）在Input surface raster下拉列表框中选择 经过填充的DEM图层，在Output flow direction raster文本框中确定输出文件的名 称及路径，点击选中Force all edge cells to flow outward前的复选框
• 1）点击Spatial Analyst 工具条上的Spatial Analyst/ Raster Calculator,打开Raster Calculator对话框。
• 2）利用公式 LS Am (sin ) n 来计算地形 因子LS。其中，A是汇流累积量（FlowAcc）， β是坡度（slope，单位为角度），m=1.6， n=1.3。
• 主要因子包括：坡度、坡向、平面、水流 长度。大多数的地形因子由地形表面的方 向导数计算而成。它们可通过一个二阶差 分格式或者满足DEM的一个二元插值函数 z=f(x, y)计算获得，之后计算函数的导数。 或许，我们想先计算出DEM的低洼地，就 必须指出一个或多个规则来决定排水方向 和每个像元间的联系，以便计算水流长度 和上坡汇水区域。目标是能够利用计算出 的因子来描述形态测量学、汇水位置以及 山坡和组成流域通道的地表因子。 1. 《地理信息系统在流域生态水文过程模 拟研究中的应用》，夏佰成、胡金明、宋 新山，水土保持研究，2004，11 • 2. 《不受DEM 空间分辨率影响的地形指数 计算》，徐静,任立良,程媛华,袁飞，清华大 学学报(自然科学版)，2008年第48卷第6期
pva计算公式
pva计算公式PVA(Porosity Coefficient) 计算公式是土壤力学中一个重要的参数，用于描述土壤孔隙结构的状况，其大小直接影响土壤的承载压力和水土保持性能。
下面是常见的 PVA 计算公式及其解释:1. 常规 PVA 计算公式常规 PVA 计算公式为:PVA = (1 - A%)/100其中，PVA 为土壤孔隙率 (Porosity Coefficient),A%为土壤固体颗粒含量 (Adhesion Coefficient)。
该公式基于土壤孔隙率和固体颗粒含量之间的关系，通过将土壤孔隙率降低到 100% 以下，可以将土壤孔隙率表示为土壤固体颗粒含量的百分比。
2. 地形因子 PVA 计算公式地形因子 PVA 计算公式为:PVA = (1 - B%)/100其中，PVA 为地形因子 (Topographic Coefficient),B%为土壤粘粒含量 (Bogosity Coefficient)。
该公式基于土壤孔隙率和粘粒含量之间的关系，考虑了地形对土壤孔隙率的影响。
在地形因子的计算中，通常会考虑到地形高度、坡度和地形复杂度等因素。
3. 水文 PVA 计算公式水文 PVA 计算公式为:PVA = (1 - C%)/100其中，PVA 为水文 PVA(Hydrologic Porosity Coefficient),C%为土壤孔隙中水分含量 (Water Content Coefficient)。
该公式基于土壤孔隙率和水分含量之间的关系，考虑了土壤孔隙对水分的容纳能力。
在水文 PVA 的计算中，通常会考虑到土壤的饱和度、湿度、孔隙大小等因素。
以上是常见的 PVA 计算公式，这些公式可以帮助工程师和研究人员更好地理解和预测土壤的力学性质和水文特性。
同时，PVA 计算公式也是土壤力学和水土保持领域中重要的研究内容之一。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第七章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坡向值的判断fy fxα =Aspect 坡向示意=0 >0 ∕90 =0 ∕0 <0 ∕270>0 >0 0～90 α=0 0 0<0 -90～0 360+α<0 >0 -90～0 180+α=0 0 180 <0 0～90 180+α注：上述情况假定所建立的DEM数据从南向北获取的，且x轴与正北方向重合，否则上述公式求得的坡向值，还应加上x轴偏离正北方向的夹角值。
注：上述情况假定所建立的DEM数据从南向北获取的，且x轴与正北方向重合，否则上述公式求得的坡向值，还应加上x轴偏离正北方向的夹角值。
3.3坡形●坡形是指局部地表坡面的曲折状态。
宏观上讲，一般可分为直线形斜坡、凸形斜坡、凹形斜坡和台阶形斜坡四种基本类型。
从微观角度上，一般可采用地面曲率因子和地面变率因子度量地面表面一点的弯曲变化程度。
（1）宏观坡形因子●直线形斜坡：从分水岭到斜坡底部地面坡度基本上不变。
●凸形斜坡：地面坡度随着距分水岭距离增加而增加。
邻近分水岭附近的地面平缓，以后随坡长的增加，坡度亦增加。
●凹形斜坡：斜坡上半部坡度较陡，下半部坡度较缓。
此种坡形常以沉积为主，较多分别在山区与阶地平原接壤处或河谷的两岸。
●台阶形斜坡：台阶形斜坡是斜坡与阶地相间的复式，可以看作是凸形坡与凹形坡的组合。
（2）地面曲率因子●地面曲率是对地形表面一点扭曲变化程度的定量化度量因子，地面曲率在垂直和水平两个方向上分量分别称为平面曲率和剖面曲率。
●剖面曲率是对地面坡度的沿最大坡降方向地面高程变化率的度量。
数学表达式为：（7.3）●平面曲率指在地形表面上，具体到任何一点P，指用过该点的水平面沿水平方向切地形表面所得的曲线在该点的曲率值。
平面曲率描述的是地表曲面沿水平方向的弯曲、变化情况，也就是该点所在的微小范围内坡向变化程度的度量。
数学表达式为：（7.4）●曲率因子的提取算法的基本原理为：在DEM数据的基础上，根据其离散的高程数值，把地表模拟成一个连续的曲面，从微分几何的思想出发，模拟曲面上每一点所处的垂直于和平行于水平面的曲线，利用曲线曲率的求算方法的推导得出各个曲率因子的计算公式。
（3）地面变率因子●地面变率描述的是地表局部范围内坡度、坡向两个基本的地形指标的变化情况，它包括坡度变率、坡向变率两个基本因子。
●地面坡度变率，是地面坡度在微分空间的变化率，是依据坡度的求算原理，在所提取的坡度值的基础上对地面每一点再求算一次坡度。
即坡度之坡度（Slope of Slope, 简称SOS）。
●地面坡向变率，是指在地表的坡向提取基础之上，进行对坡向变化率值的二次提取，亦即坡向之坡度（Slope of Aspect, SOA）。
地面坡向变率在所提取的地表坡向矩阵的基础上沿袭坡度的求算原理，提取地表局部微小范围内坡向的最大变化情况。
3.4坡长●坡长通常是指在地面上一点沿水流方向到其流向起点间的最大地面距离在水平面上的投影长度。
其数学表达为：(7.5)式中L指坡长，m指地表面沿流向的水流长度，θ指水流地区的地面坡度值。
●自然条件下，水流向低处流动，遇到洼地，首先将其填满，然后再从该洼地的某一最低出口流出。
但在一个连续的栅格中，地形洼地的存在，导致依据水流方向矩阵所提取的排水网络不连续，使自然水流不能畅通无阻地流至区域地形的边缘。
因此，对已有的DEM数据，首先要进行洼地填充，生成无洼地DEM。
DEM中的洼地可分为凹陷型洼地和阻挡型洼地。
一般情况下，对于阻挡型洼地，可降低阻挡物存在处的高程，使水流穿过障碍物；对于凹陷型洼地，采用常规的将洼地内所有栅格单元垫高至洼地周围最低栅格单元高程的方法。
坡长提取流程如图7.3。
● DEM数据中的平地，包括原始DEM中的平地和洼地填平产生的平地。
平地区域的存在使得该区于水流方向的确定出现不确定性，因此需要对平地区域进行处理。
基本的处理方法思想是对平地范围内的单元格增加一微小增量，每个单元格的增量大小是不一样的，这样每个单元格就有一个明确的水流方向，以便能够产生合理的汇流水系。
●水流方向是水流离开此格网时的指向。
确定水流方向的算法有很多，每一种算法的假设前提都不尽不同，因此得到的结果也有差异。
基本的算法根据其基本思想可以大致可以分为：单流向算法（SFD）和多流向算法（MFD）以及其他算法。
单流向算法是将某单元格上产生的径流都流向一个最低的相邻单元格，多流向算法将径流按一定的比例流向若干相对较低的相邻单元格。
根据流向对每个格网进行追踪，并记录其追踪路径的距离，即坡长结果。
3.5坡位●坡位是指坡面计算单元所整个大坡面的地貌部位。
例如，位于正地形还是负地形等；处于沟间地还是沟坡地。
从一定意义上讲，DEM单元栅格坡位的提取核心的工作是地貌部位、地貌类型的划分。
以黄土高原沟间地、沟谷地的划分为例来说明。
●黄土高原地区地形复杂，具有独特的堆积地貌特征，从侵蚀地貌形态的坡度组合特征、地貌成因的侵蚀特征和土地利用特征等三个主要方面，可以将该地区的区域地貌在垂直方向上划分为沟间地、沟坡地和沟底地三种基本类型。
沟间地是指黄土高原地区位于地形结构线（沟沿线）以上相对比较平缓的地形部分。
●基于DEM提取沟间地的基本方法是先从沟间地的形态特征和成因原理出发，根据坡度特征提取缓坡图层，再利用沟谷缓冲分析，获得沟谷图层，两层相减，得到沟间地。
3.6坡面复杂度因子（1）地形起伏度●地形起伏度是指，在所指定的分析区域内所有栅格中最大高程与最小高程的差。
●可表示为如下公式： (7.6) 式中，RFi指分析区域内的地面起伏度，Hmax指分析窗口内的最大高程值，Hmin 指分析窗口内的最小高程值。
●地形的起伏是反映地形起伏的宏观地形因子，在区域性研究中，利用DEM 数据提取地形起伏度能够直观的反映地形起伏特征。
在水土流失研究中，地形起伏度指标能够反映水土流失类型区的土壤侵蚀特征，比较适合区域水土流失评价的地形指标。
（2）地表粗糙度●地表粗糙度是反映地表的起伏变化和侵蚀程度的指标，一般定义为地表单元的曲面面积S曲面与其在水平面上的投影面积S水平之比。
●用数学公式表达为： R = S曲面 / S水平(7.7)●地表粗糙度能够反映地形的起伏变化和侵蚀程度的宏观地形因子。