特征地理要素提取实用教案

1、8.1 8.1 地形(dxng)(dxng)特征点的提取地形特征点主要包括山顶(shn dn)点、凹陷点、山脊点、山谷点、鞍点、平地点等。第1页/共35页第一页，共35页。利用DEMDEM提取地形特征点，可通过3 3* *3 3或者更大的栅格窗口，通过中心格网点与8 8个邻域格网点的高程关系来进行判断和获取。即在一个局部区域内，用x x方向和y y方向上关于(guny)(guny)高程z z的二阶导数的正负组合关系来判断。第2页/共35页第二页，共35页。（1 1）加载DEMDEM第3页/共35页第三页，共35页。（2 2）生成(shn chn)(shn chn)间隔为1515米的等高线第4页

2、/共35页第四页，共35页。（3 3）生成(shn chn)(shn chn)间隔为7575米的等高线第5页/共35页第五页，共35页。（4 4）利用HillshadeHillshade模块提取光照晕渲图，作为等高线的背景。如果需要还要进行有效数据的提取，主要采用DEM=0DEM=0的办法进行低洼地的排除(pich)(pich)。然后按照1515米、7575米等高线、山脊阴影文件的顺序进行排序。第6页/共35页第六页，共35页。（5 5）利用Neighborhood StatisticsNeighborhood Statistics模块进行领域单元的统计，提取分析(fnx)(fnx)窗口最大值

3、。第7页/共35页第七页，共35页。（6 6）利用栅格计算模块，提取山顶点区域。主要(zhyo)(zhyo)采用maxpoint-DEM=0maxpoint-DEM=0。第8页/共35页第八页，共35页。8.28.2 山脊(shnj)(shnj)线、山谷线的提取8.2.18.2.1谷脊特征线提取方法概述p山脊线和山谷线构成了地形起伏变化的分界线，它对于地形研究有重要意义。p对于水文物理过程研究而言，山脊线和山谷线的提取实质上也是分水线与汇水线的提取。p已有的自动提取山脊线和山谷线的方法，按照使用的数据大致可以分为三种：p基于数字化等高线数据的方法；p基于规则格网DEMDEM数据的方法；p基于D

4、elaunayDelaunay三角网和VoronoiVoronoi数据的方法；p从算法设计(shj)(shj)原理上来分，大致可以分为以下四种：p基于图像处理技术的原理；p基于地形表面几何形态分析的原理；p基于地形表面流水物理模拟分析原理；p基于地形表面几何形态分析和流水物理模拟分析相结合的原理。第9页/共35页第九页，共35页。8.28.2 山脊线、山谷(shng)(shng)线的提取8.2.28.2.2谷脊特征线提取算法（1 1）基于图像处理技术的原理的算法采用各种滤波算子进行边缘提取。这种方法提取山脊山谷线的主要过程有两步：首先提取地形特征点（山脊点、山谷点、鞍点等）；将特征点连成地形特

5、征线（山脊线、山谷线）；该方法存在两个主要缺陷：提取特征点时必须排除DEMDEM中噪声的影响；将特征点连接成线时的算法设计较为(jio wi)(jio wi)困难。第10页/共35页第十页，共35页。（2 2）基于地形表面几何形态分析(fnx)(fnx)原理的算法典型算法就是断面极值法。其基本思想就是地形断面曲线上高程的极大值点就是分水点，而高程的极小值点就是汇水点。该方法的基本过程为：找出DEMDEM的纵向与横向的两个断面上的极大、极小值点，作为地形特征线上的备选点；根据一定的条件或准则将这些备选点划归各自所属的地形特征线。第11页/共35页第十一页，共35页。由于这种方法对地形特征线上点的

6、判定与其所属的地形特征线的判定是分开进行的，在确定地形特征线时，全区域采用一个相同的曲率阈值作为判定地形特征线上点的条件。因此它忽略了每条地形特征线必然存在的曲率变化现象。当阈值选择较大时，会丢失许多地形特征线上的点，导致后续跟踪的地形特征线间断且较短；如果选择过小，会产生地形特征线上点的误判，给后续地形特征线的跟踪带来困难。由于该方法只选择纵、横两个断面来去确定高程变化的极值点，因此它所确定的地形特征线具有一定的近似(jn s)(jn s)性，与实际的地形特征线有一定的差异，有时候还会出现遗漏。这种算法存在两个主要缺陷：第12页/共35页第十二页，共35页。（3 3）基于地形表面流水物理模拟

7、分析原理的算法这种算法的基本思想是：按照流水从高至低的自然规律，顺序计算每一栅格点上的汇水(hushu)(hushu)量，然后按汇水(hushu)(hushu)量单调增加的顺序，由高到低找出区域中的每一条汇水(hushu)(hushu)线。根据得到的汇水(hushu)(hushu)线，通过计算找出各自汇水(hushu)(hushu)区域的边界线，就得到了分水线。这一算法采用了DEMDEM的整体追踪分析的思路与方法，分析结果具有系统性好，还便于进行相应的径流成因分析。第13页/共35页第十三页，共35页。第14页/共35页第十四页，共35页。p该方法也存在以下两个明显的缺陷：p由于该算法所计算的

8、汇水量与汇水量小，高程值小的地形特征线上的点的汇水量大。因此，可能导致低处非地形特征线上的点的汇水量也较大而被误认为地形特征线上的点；而位于高处的地形特征线上的点会因为汇水量小而被排除(pich)(pich)；这就造成用该算法所确定的地形特征线（汇水线）的两端效果很差。p由于该算法降格汇水区域的公共边界视为分水线，因此它所确定的分水线均为闭合曲线，这与实际的地形特征线（山脊线）不符。第15页/共35页第十五页，共35页。（4 4）基于地形表面几何形态和流水物理模拟分析的算法基本思路是：首先采取较稀疏的DEMDEM格网数据，按流水物理模拟算法去提取区域内概略的地形特征(tzhng)(tzhng)

9、线；然后用其引导，在其周围邻近区域对地形进行几何分析，来精确的确定区域的地形特征(tzhng)(tzhng)线。关键步骤：求出已提取的概略的地形特征(tzhng)(tzhng)线与DEMDEM格网线的交点，在该交点附近的一个小区域内，对DEMDEM数据进行几何分析，即找出该区域内与概略的地形特征(tzhng)(tzhng)线正交方向地形断面上高程变化的极值点，该点即为地形特征(tzhng)(tzhng)线的精确位置。基本过程：1 1）概略DEMDEM的建立；2 2）地形流水物理模拟；3 3）概略地形特征(tzhng)(tzhng)线提取；4 4）地形几何分析；5 5）地形特征(tzhng)(t

10、zhng)线精确确定。第16页/共35页第十六页，共35页。（5 5）平面曲率与坡位组合法即首先利用DEMDEM数据提取地面(dmin)(dmin)的平面曲率及地面(dmin)(dmin)的正负地形，取正地形上平面曲率的大值即为山脊，负地形上平面曲率的大值为山谷。该种方法提取的山脊、山谷的宽度可由选取平面曲率的大小来调节，方法简便效果好。第17页/共35页第十七页，共35页。8.38.3 沟沿线的提取沟沿线是黄土高原地貌最重要的一条地形特征线。利用(lyng)DEM(lyng)DEM提取沟沿线方法的核心是利用(lyng)(lyng)沟间地和沟坡地在形态与成因上的差异性进行提取。坡度：沟间地和沟

11、坡地在坡度上存在明显差异。剖面曲率：沟沿线与大剖面曲率值的黑色条带呈很好的吻合。第18页/共35页第十八页，共35页。p沟壑：利用DEM首先提取沟坡地上的切沟与冲沟，并以此(y c)构建缓冲区，再作等值连接，即形成沟沿线。第19页/共35页第十九页，共35页。p高程变异：基于沟沿线的数学特征，引入边缘检测算子，检测反映沟沿线两侧高程变化的候选点，借助形态学方法把这些(zhxi)候选点连接成线段，再通过滤除细碎线段，最终生成沟沿线。第20页/共35页第二十页，共35页。8.48.4 水系的提取水系提取的有关概念：集水流域、集水出口(ch ku)(ch ku)、子流域、分水岭、水流网络第21页/共

12、35页第二十一页，共35页。p基于地表径流(jngli)(jngli)漫流模型的水系提取方法p四个基本问题：p洼地的处理p平地的处理p水流方向及水流累积量的确定p水道起始位置的确定第22页/共35页第二十二页，共35页。第23页/共35页第二十三页，共35页。p水流(shuli)(shuli)方向及水流(shuli)(shuli)累积量的确定第24页/共35页第二十四页，共35页。第25页/共35页第二十五页，共35页。第26页/共35页第二十六页，共35页。8.58.5 流域(liy)(liy)的提取第27页/共35页第二十七页，共35页。第28页/共35页第二十八页，共35页。第29页/共35页第二十九页，共35页。8.68.6 可视性分析(fnx)(fnx)两点之间的可视性第30页/共35页第三十页，共35页。第31页/共35页第三十一页，共35页。第32页/共35页第三十二页，共35页。p可视性分析最基本的用途p可视查询p地形可视结构(jigu)(jigu)计算（即可视域的计算）p水平可视计算第33页/共35页第三