《基于不同算法的剖面曲率提取研究-毕业论文》由会员分享，可在线阅读，更多相关《基于不同算法的剖面曲率提取研究-毕业论文（18页珍藏版）》请在金锄头文库上搜索。
1、 南阳师范学院20XX届毕业生毕业论文（设计） 题 目： 基于不同算法的剖面曲率提取研究 完 成 人： 班 级： 学 制： 专 业： 地理信息系统 指导教师： 完成日期： 目 录摘要(1)1引言(1)2剖面曲率的含义(1)3剖面曲率的计算方法(2)3.1汤国安提出的SOS算法(3)3.2四次表面模型(3)3.3基于空间矢量模型的差分计算(3)4由不同算法提取的剖面曲率(4)4.1研究区概况(4)4.2由汤国安提出的SOS算法提取的剖面曲率(6)4.3由四次表面模型算法提取的剖面曲率(7)4.4由基于空间矢量模型的差分计算算法提取的剖面曲率(8)5试验结果及分析(9)5.1不同算法提取剖面曲率的
2、统计指标分析(9)5.2基于三种算法提取的剖面曲率的优缺点、误差及制图效果分析(11)6结论(14)参考文献(15)Abstract(15)基于不同算法的剖面曲率提取研究 摘要：本文在不考虑DEM误差的情况下，通过对研究区（南阳西峡县局部地区）的DEM用ArcGIS软件做图分析，基于目前主要的三种剖面曲率的计算算法，即汤国安的SOS算法、四次表面模型、基于空间矢量模型的差分计算，分别提取剖面曲率，并对它们进行归纳分析，以对比它们的优缺点、误差、及制图效果，为以后求出最适合的地形剖面曲率提供基础。研究结果表明，在计算地形起伏不是很大的地区的剖面曲率时，汤国安的SOS算法最为精确，四次表面模型算法
3、次之。关键词：剖面曲率；算法；SOS；1 引言地形因子是为定量表达地貌形态特征而设定的具有一定数学意义的数学参数或指标1。剖面曲率（Profile Curvature）KV 是地形因子之一，基于提取算法的坡面因子分，剖面曲率属于二阶坡面因子。在地理信息系统（GIS）中，剖面曲率的提取算法有三种：汤国安提出的SOS/SOA算法、四次表面模型算法和基于空间矢量模型的差分计算算法。数字高程模型（DEM）作为地理信息系统( GIS)空间数据库的核心数据库之一，是进行二维地形空间模拟和地形图生产的基础，DEM 能派生出各种地形因子，如坡度、坡向、平面曲率、剖面曲率、汇水面积等1。其中，剖面曲率是最重要的
4、地形因子之一，地面曲率影响着，也是制约生产力空间布局及自然地理研究的重要因子 。尽管利用GIS 平台在DEM 上提取坡度的算法已经成熟，但是由于算法会影响坡度的测算，所以需要比较不同算法对提取坡度的影响。与坡度、坡向等基本地形参数一样，剖面曲率在数学定义上是明确的和唯一的，具有确定的表达式。尽管其理论定义是明确的，然而DEM是地形曲面的微分模拟，算法设计必然存在各种各样的假设，不同假设和前提导致不同的剖面曲率计算模型和结果，这虽然对地形特征的可视化和地形分类的影响不大，但对以数值计算为主的地学分析模型的影响却是非常显著的2。2 剖面曲率的含义剖面曲率是地面曲率在垂直方向上的分量，是对地面坡度的
5、沿最大坡降方向的地面高程变化率的度量，是研究地形表面的重要因子之一。剖面曲率因子提取算法的基本原理为：将离散的高程数值拟为一个连续的曲面，基于微分几何的思想，模拟曲面上每一个点所处的平行于水平面的曲线，利用曲线曲率的求算方法推导得出各个曲率因子的公式，进而求算出每一点的曲率值1。剖面曲率KV 具有重要的理论和实践意义。KV 决定地表及土壤中物质移动的相对速度，KV 0表明移动加速，KV 0表明移动减速，反映了侵蚀或沉积的程度。KV 通常可以揭示出谷脊线呈突起状或阶梯状，还可以用在很多领域，如在大尺度区域中模拟土壤湿度、土层厚度、土壤侵蚀、山崩石头的散布、植被等；也可用于不同尺度的地貌研究，研究
6、区域尺度和次大陆尺度地质构造、地貌和土壤的空间相关规律；还用于识别显示断层、地质断陷线等3。笛卡尔坐标系中，x、y 表示空间位置，z = f( x，y) 表示地形曲面函数，假定函数具有单值性、连续性和二阶可导性，则剖面曲率的计算公式为：Kv=-p2r+2pqs+q2tp2+q2(1+p2+q2)3/2 （1）其中，p=zx ，是x方向高程变化率；q=zy ，是y方向高程变化率；r=2zx2 ，是x方向高程变化率的变化率；s=2zxy ，是x方向高程变化率在y方向的变化率；t=2zy2 ，是y方向高程变化率的变化率剖面曲率提取算法一般都是根据P.A.Borrough提出的窗口微分法遍历整个图幅进
7、行计算的。即是在33的DEM栅格窗口中进行的，窗口在DEM数据矩阵中连续移动后而能完成整幅图的计算工作4。Z1Z2Z3Z4Z5Z6Z7Z8Z9图1 33地形窗口(Z1Z9分别代表各个栅格的高程值)3 剖面曲率的计算方法GIS中剖面曲率的计算大多是基于DEM数据的，由离散的数据计算剖面曲率有不同算法，在考虑到算法的通用性和使用性的情况下，这里选择了ArcGIS采用的汤国安提出的SOS/SOA算法、四次表面模型算法和基于空间矢量模型的差分计算算法1三种常用的剖面曲率算法进行分析。3.1 汤国安提出的SOS算法剖面曲率为地面高程变化的二次导数，即坡度在垂直方向的最大变化率；Tang已经证明5，在Ar
8、c/Info等GIS软件中，地面剖面曲率的数字矩阵可以直接通过对DEM数据的求取地面坡度的坡度（Slope of Slope,SOS）而获得。SOS就是坡度的变率，是在所提取的坡度值的基础上对地面每一点再求算一次坡度，坡度是地面高程变化率的求解，因此，坡度变率表示了地球表面高程相对于水平面变化的二阶导数3。这里的“坡度”并非是坡度，而是按照坡度的算法实现了对坡度数据阵列中其变化率的量化提取。所获得的剖面曲率值在0,90内，并非是真正的曲率值，但却真实的地反映了坡度和等高线的变化率，可提取山谷、山脊等地形因子，而且因其算法简单易于实现，所以也称剖面曲率，有些教材中称之为地面变率因子3。3.2 四
9、次表面模型Zeverbergen 和 Thorne(1987)提出了一种计算算法，构造一个含有9个参数的四次方程（x、y的最高次数为2）6，见公式（2）.函数的所有参数由33窗口（如图1）内的单元值唯一确定，函数经过每个数据点，包括中心单元。系数A、B、C、I在计算曲率时未用到，可以忽略，系数DH由公式（3）（7）计算，则中心格网的剖面曲率的表达式变为表面函数系数组成的式子，见公式（8）3。Z = Ax2y2 + Bx2y + Cxy2 + Dx2 + Ey2 + Fxy + Gx + Hy + I. (2)D =（Z4 + Z6）/2-Z5/L2 . (3)E =（Z2 + Z8）/2-Z5
10、/L2 . (4)F =(- Z1 + Z3 +Z7 Z9)/4L2. (5)G =(- Z4 + Z6)/2L. (6)H =(Z2 Z8)/2L . (7)剖面曲率：KV = 2(DG2 + EH2 + FGH)/(G2 + H2). (8)3.3 基于空间矢量模型的差分计算 在DEM数据的基础上，从微分几何的思想出发，求出剖面曲率公式(1)中的各个导数pt,即可求出各曲率值，差分计算常用的是三阶反距离平均权差分。对每个栅格点都确定一个33分析窗口，如图2（1），ai分别代表各个栅格的高程值，L表示栅格大小，p、q可由以下公式求出3。p=zx=c+2f+i-(a+2d+g)8L （9）p=
11、zx=a+2b+c-(g+2h+i)8L (10)abcaxbxcxaybxcxdefdxexfxdxexfxghigxhxixgxhxix （1） （2） （3）图2 提取地面曲率的分析窗口 求出所有 p、q 值，则组成 p 值栅格图（ 见图1（ 2 ） ） 和 q 值栅格图（ 见图1（ 3 ） ） 。 以此数据为基础，进行差分运算，可求出r、s、t8,见公式( 10 ) 、（ 11 ） 和（ 12 ） 。r=2zx2=cx+2fx+ix-(ax+2dx+gx)8L （11）s=2zxy=ax+2bx+cx-(gx+2hx+ix)8L （12）t=2zy2=ay+2by+cy-(gy+2hy
12、+iy)8L （13）4 由不同算法提取的剖面曲率由任何一种算法得出的剖面曲率，我们都要评估其精确性和制图效果，并以此来作为我们求算地形的剖面曲率的依据。因此，一个科学、合理、准确的剖面曲率分析是目前研究剖面曲率的首要任务9。在此，以西峡县的某一地区为研究区，我将依据三种剖面曲率算法所作出的图像和数据进行比较分析，得到柱状图和各统计量，为找到不同高程域的适合剖面曲率提取算法提供依据。4.1 研究区概况研究区地处伏牛山腹地，豫、鄂、陕交界地区，东临南阳盆地，西嵌秦岭，南襟荆襄平原，北负八百里伏牛山。境内地形复杂，北部是海拔高、坡度大的中低山地，南部是鹳河谷地，两侧是起伏大的低山丘岭。海拔181米，自然坡降为33%。研究区内的高程分布如表1所示。由表1知，分布在120405之间的高程数值最多，为3197380个，约占总数值个数的27.9%，分布在405690之间的高程数值约占总数值个数的22.5%，分布在690991之间的高程数值约占总数值个数的21.5%，
