（毕业论文）基于不同算法的DEM地面曲率提取的比较分析.docVIP

基于不同算法的DEM地面曲率提取的比较分析 摘要：下文介绍了几种不同的曲率， 并简述了不同的DEM地面曲率提取的步骤，且通过南阳市DEM图用四次表面模型法和汤国安的SOS/SOA方法进行地面曲率的提取，得出了二种方法的不同，并进行了一定的比较分析，解决地面曲率 关键词:DEM; 曲率; 地面曲率 引言 ? 有效提取地形信息，特别是提取常规方法很难提取的地形要素，一直是许多测绘和地学工作者致力研究的目标。近年来，在 DEM 数据生产和分析方法方面取得了巨大进步， 但是从不同地形复杂度 DEM 提取地形信息，特别是研究几乎与坡度及 DEM 在各领域的广泛应用严重脱节。DEM 数据本身多尺度因素，加之地形、地貌特征具有宏观性与区域分异性的特点，直接的信息提取往往很难达到预期的目的。利用 DEM 制作图高效、省力， 但其精度有很大的不确定性，同时 DEM 制作过程中的误差传播、转移对信息的影响缺少系统的判断依据。选取实验样区， 在 DEM 生产的基础上提比较与统计分析 1.1 DEM数据获取 近年来,随着高新技术的运用,在DEM数据的建立DEM的方法有多种从数据源及采集方式讲有：(1)直接从地面测量,例如用GPS、全站仪、野外测量等;(2)根据航空或航天影像，通过摄影测量途径获取,如立体坐标仪观测及空三加密法、解析测图等等;(3)从现有地形图上采集，如格网读点法、数字化仪手扶跟踪及扫描仪半自动采集然后通过内插生成DEM等方法。DEM内插方法很多,主要有分块内插、部分内插和单点移面内插三种。目前常用的算法是通过等高线和高程点建立不规则的三角 网(Triangular Irregular Network, 简称TIN)。然后在TIN基础上通过线性和双线性内插建DEM。  地面曲率是对地形表面一点扭曲变化程度的定量化度量因子，地面曲率在垂直和水平两个方向上分量分别称为平面曲率和剖面曲率。: 首先是曲线的切线转过的角度( 简称切线转角) △α，△α 越大，则曲线的弯曲程度越大; 其次是曲线的长度，对于相同的切线转角△α，曲线越短，则曲线的弯曲程度越大． 因此，以单位弧长上曲线切线转角的值来衡量曲线的弯曲程度，即曲率 ．曲线的曲率是曲线切线倾角关于弧长变化率的绝对值． 曲面的曲率则是过该点的某一平面与曲面的交线的曲率. 2.1 剖面曲率(Kv) 剖面曲率是对地面坡度的沿最大坡降方向地面高程变化率量。 Kv)是地面坡度的变化率,是反映地形起伏变化 特征的重要指标之一,是确定坡形以及提取沟沿线、沟底线 等地形转折线等的重要定量指标 。 数学表达式 （1） 是X方向高程变化率； ，是Y方向高程变化率; ,是X方向高程变化的变化率; ，是x方向高程变化率在y方向的变化率; ,是y方向高程变化率的变化率. 2.2 正切曲率(kh) 正切曲率指在地形表面上，具体到任何一点 P，指切地形表面所得的曲线在该点的曲率值。曲率描述的是地表曲面沿水平方向的弯曲、 变化情况， 也就是该点所在的微小范围内坡向变化程度量。 地表曲面沿水平方向弯曲、变化情况，也就是该点所在的平面等高线的弯曲程度。从另一个角度讲，地形表面上一点的正切曲率也是对该点微小范围内坡向变化程度的度量。 正切曲率和剖面曲率都为负值时是凹形坡，都为正值时为凸形坡，坡型可以利用地表的曲率进行描述和量化，直线形和凸型斜坡在曲率上的体现是曲率≥0，凹型坡和阶梯型坡的曲率＜0数学表达式为： （2） 2.3 等高线曲率 在计算等高线曲率值时，通常采用样条函数进行等高线内插。待全区域所有等高线处理完毕，找出区域内的山顶点和谷底点并以这些点为起始点，根据一定的条件和搜索策略将已确定的地形特征线上点的备选点确定为各自所在的山脊线和山谷线。等高线曲率在确定地形特征线时，将地形特征线上点的判定与该点所属的地形特征线的确定分开来考虑，这将给后续地形特征线的跟踪带来很多麻烦。另外，在地形破碎地区或等高线不光滑的区域，地形特征的跟踪将十分困难。 等高线曲率是等高线方向的变化率，反映了坡向的变化。 其公式为： （3） 图1剖面曲率、正切曲率和等高线曲率 2.4 综合曲率 地面剖面曲率与正切曲率的和，显示地面曲率的综合指数。地面的综合曲率指数决定地面的复杂及地面切割的破碎程度。目前尚无统一的计算公式，可采用k=(正切曲率+剖面曲率)*tg(slope)计算地面综合区率。 2.5 三种地面曲率的内涵差异 Kv和Kh具有重要的理论和实践意义． Kv决定地表及土壤中物质移动的相对速度，Kv ＞ 0 表明移动加速，Kv ＜ 0 表明移动减速，反映了侵蚀或沉积的程度． Kh决定地表