论文 DEM地形因子及其应用.docVIP

DEM地形因子及其应用 江 帆① 朱长青① 李小荣② ① 信息工程大学测绘学院 郑州 450052 ② 江西省景德镇陶瓷学院 景德镇 333000 摘要：地形因子能表示地形表面的基本特征，但很难用一种地形因子也很难准确具体地表达地形的特征。为此，可将各种地形因子综合考虑，在一定程度上来刻画地形表面的起伏变化，具有较好的效果。因此地形因子在研究地形变化时具有重要的应用价值。 关键字：微观地形因子 宏观地形因子 相关地形因子 1. 引言 地形表面是一个极不规则的曲面，DEM是地形的一个数学模型，从这个意义上讲，可将DEM看作一个或多个函数的和。实际上许多地形因子就是从这些函数中推导出来的。 如果对函数求一阶导数并进行组合，则可得到一系列的因子值如坡度、坡向、变差系数、变异系数等的函数；如果求二阶导数并进行组合则可得到坡度变化率、坡向变化率、曲率、凸凹系数等的函数。从理论上说，还可以继续求三阶、四阶等更高阶的导数直到无穷阶以派生更多的地形因子。但在实际应用中，对DEM进行高于二阶的求导意义已经很小，至少到目前为止还没有探讨过高于二阶的应用价值。 上述地形因子也称为地貌因子。用多种地表形态描述参数可以描述地表形态的一种或多种特征及其地形表面的复杂程度。本文中将地形因子可以划分为三类：微观地形因子、宏观地形因子、相关地形因子，对其进行归纳总结，并对其应用进行阐述。 2. 微观因子 空间信息的研究中，空间物体通常被抽象为点、线、面、体（曲面）等四大类，而除点以外的空间物体都具有形态特征。地形表面是一个极不规则的曲面，在地学研究中我们经常用基本地形因子的各种地表形态描述参数来描述地表形态的一种或多种特征以及地形表面的复杂程度。而微观因子包括表面积、体积、坡度、坡向、坡长、坡形、曲率等。 2.1. 表面积和体积 空间曲面表面积的计算与空间曲面拟合的方法以及实际使用的数据结构（规则格网或者三角形不规则格网）有关。对分块曲面拟合，曲面表面积由分块曲面片之和给出，因此问题的关键是要计算出曲面片的表面积。对于全局拟合的曲面，通常也是将计算区域剖分成若干规则单元，然后对每个单元计算出其面积，将其累积计算总面积。 考虑到实际地形表面的复杂性，可以利用数学曲面进行评价，数学曲面应最大限度的接近局部实际地行，简单曲面并不能很好的反映算法所具有的精度及相关参数确定。本文在给定的区域[0，100]×[0，100]内，用四个点对实际地形曲面逼近，并按不同的算法求解其曲面面积。对比几种算法效果,分别设定不同的DEM格网间距，选取函数。然后分别采用正方形格网划分为三角形格网、复化梯形公式、复化Simpson公式、复化Cotes公式计算曲面表面积，其计算结果见表1。为了比较计算效果，给出原始曲面精确值和基于双线性多项式的DEM表面精确值。 表1 公式 等 分数 DEM表面 精确值 正方形格网 划分为 三角形格网 复化 梯形公式 复化 Simpson 公式 复化 Cotes公式 10 2142944.6507 2143681.4923 2143132.6718 2142489.2096 1965798.8197 40 2142551.8155 2142599.8116 2142525.5103 2142483.6353 2098263.6390 80 2142490.8423 2142502.6176 2142494.0423 2142483.5279 2120368.8562 160 2142485.3455 2142488.2899 2142486.1461 2142483.5162 2131425.0050 200 2142484.6864 2142486.5707 2142485.1987 2142483.5157 2133636.5214 原始曲面精确值 2142483.5196 基于复化求积方法，可对DEM几种曲面函数表面运用不同的计算算法，得到了相应的函数关系公式，进而得到了其空间曲面表面积。实验表明，复化Simpson公式的结果较好，梯形公式和正方形格网划分为三角形格网的结果次之，而复化Cotes公式的结果较差。 体积是指空间曲面与一基准平面之间的空间的体积。与表面积计算一样，体积也与空间曲面拟合的方式、实际使用的数据结构（规则格网或者三角形不规则格网）有关。体积的计算通常也采用分块计算方法。 2.2. 坡度和坡向及其变化率 坡向与坡度是相互联系的两个参数，坡度反映斜坡的倾斜程度，坡向反映斜坡所面对的方向，两者是地形描述常用的参数。地面上某点的坡度是表示地表面在该点倾斜程度的一个矢量，其模等于地表曲面函数在该点的切平面与水平面夹角的正切，其方向等于在该切平面上沿最大倾斜方向的某