三角网数字地面模型的构建PPT课件

.,1,第七章数字地面模型,.,2,1概述一、数字地面模型的定义和特点,1.数字地面模型的定义数字地面模型是描述地球表面形态多种信息空间分布的有序数据阵列。用公式表示为：式中：Kp为第p号地面点上第k类地面特征信息的取值（如高程值、土壤类型等）；up、vp为第p号地面点的二维坐标（地图投影平面坐标、经纬度、矩阵行列号）m为地面特征信息类型的数目,n为地面点的个数。,.,3,2.数字地面模型的内容地貌信息：高程、坡度、坡向等基本地物信息：居民地、水系、道路境界线等主要的自然资源和环境信息：土壤、植被、地质、气候等主要的社会经济信息：人口分布、工农业产值、收入等,.,4,3.数字高程模型（）数字高程模型是对地球表面起伏的一种离散的数字表达，是数字地面模型中最基本的部分。可以表示成区域上三维向量序列：当该序列中各向量的平面点位是规则格网排列时，则平面坐标可省略，简化为一维向量序列：,.,5,4.DEM的优点：数字形式，可直接输入计算机，便于各种计算机辅助系统利用；容易以多种形式显示地形信息；精度不会损失；便于存储、更新、传播。,.,6,5.DEM的应用作为国家地理信息的基础数据数字线划图DLG、数字高程模型DEM、数字正射影像DOM、地理名称、数字栅格图，前三为我国国家空间数据基础设施（NSDI）的框架数据。测绘领域工程军事规划、环境遥感,.,7,6.DEM的生产流程数据采集、数据处理、DEM应用,.,8,二、DEM的形式,离散的点规则格网（优缺点）不规则三角网TIN（triangulatedirregularnetwork）混合形式,.,9,2DEM的数据获取,一、数据来源1.影像航测一直是地形图测绘和更新的最有效和最主要的手段。航摄影像是高精度、大范围DEM生产最有效的数据源。遥感影像(SPOT)也是快速生产DEM的有效手段，但精度较低。目前高分辨率卫星影像（Quickbird）、雷达图像也已用于生产DEM。,由SRTMC波段获取DEM再与TM图像叠加的结果,.,11,2.地形图地形图是DEM的另一个主要数据来源。存在的问题：地图符号的数字化、现势性差。3.地面本身GPS、全站仪、经纬仪计算机，用于小范围、大比例尺（1/2000）数字测图、土方计算，但工作量大。,.,12,二、数据采集方法1.现有地形图常用、低廉、快速，精度较低、现势性差。,.,13,地形图数字化生产DEM,.,14,2.摄影测量方法常用：等高线法、规则格网法、选择采样法、渐进采样法、混合采样法、自动采集法（全数字摄影测量）沿等高线采样法适合于地形复杂地区、不适合平坦地区，有等距离间隔和等时间间隔两种方式。,.,15,规则格网采样采样点的平面坐标由计算机控制，是一种简单、效率较高的方法，缺点是不能反映地形特征点、线。渐进采样法是一种使采样点分布合理、减少数据量的方法。,.,16,渐进采样法,此法先以较稀的间隔采样，得到p1、p3、p5点，再判断p1、p3间是否需要加密。条件是，两点间中点的二次插值与线性插值之差是否超限。如超限，则加密格网点。h2=1/8(6h3+3h1-h5)h2=1/2(h3+h1)h2=1/8(2h3-h1-h5)如h2T则内插p2、p4两点。,.,17,选择性采样为了准确反映地形，根据地形特征选择采样，沿山脊线、山谷线、断裂线和离散特征点（山顶）采点。优点是，少量的点即能反映地面的主要特征；缺点是，效率低，速度慢。混合采样混合采样是把规则格网采样（渐进采样）与选择性采样结合起来，以提高精度。,.,18,全数字摄影测量系统采集数据点由数字摄影测量系统，根据影像匹配原理，经自动相对定向、绝对定向后，建立数字立体模型，采集规则格网的数字地面模型。,.,19,.,20,3DEM数据预处理,DEM数据预处理是数据处理的工作之一，是DEM数据内插的前的准备工作。预处理包括：数据格式的转换；坐标系统的变换；数据的编辑；栅格数据的矢量化；数据的分块。,.,21,一、格式转换由于数据采集的软件、硬件不同，所采集的数据各有自己的格式，要被某一专业软件处理，建立DEM，必须首先进行数据格式的转换。如GeoTin软件只能处理God格式的数据，其它格式（txt/dxf/E00）的原始数据，都要转成God形式的。,.,22,二、坐标变换非地面坐标系变换为国家同一的地面坐标系或局部坐标系。三、数据编辑将采集的数据绘图输出或显示，以检查数据是否有错、遗漏和重复。四、栅格数据的矢量化可用自动矢量化软件实现。,.,23,五、数据分块目的是提高数据处理（内插）的速度，避免检索不必要的数据点。,.,24,4DEM的内插,目的DEM的内插就是根据已知数据点的高程，计算其它待定点的高程。原因由于DEM产品是一种密集的规则格网形式表示地面高程，而采集的数据点，通常是稀疏格网或沿等高线采样或沿地形线采样或任意采样。方法由已知点的高程，建立一个数学模型（求待定系数），求出待定点的高程。常用方法有整体函数内插、局部函数内插和逐点内插三种。,.,25,一、基于矩形格网的DEM多项式内插用已知规则格网点坐标内插出更密的规则格网点（或其它点）的高程。方法是由4点或更多的周边的相邻点拟合一个曲面（平面），然后计算待定点的高程。曲面必须连续甚至光滑。,.,26,1.双线性多项式内插根据最近的4个数据点，确定一个双线性多项式：或用矩阵形式表示为：利用4个已知点求出4个系数，再由求出的系数和待定点的平面坐标（X,Y）内插出该点的高程。,.,27,当4个已知点是规则格网点时，可由双线性多项式的定义、直积运算和阵列代数法推导出内插公式。直接由定义推导同理,.,28,把、代入得,.,29,二、双三次多项式内插三次曲面方程为,.,30,用阵列代数表示为,.,31,对规则格网的数据点，将坐标原点移到待定点所在格网的左下角。为了简单，令格网的长度为1。由于双三次多项式有16个系数，因而除了所在格网的4个角点的过程外，还需要已知这些点处的一阶偏导数和二阶混合偏导数。,.,32,四个格网点的一阶偏导数和二阶混合导数为：,.,33,因为,且,所以,.,34,因此,16个系数为,双三次多项式能保证相邻曲面的连续和光滑。,.,35,三、移动曲面拟合法,移动曲面拟合法是逐点内插法，它以每一待定点为中心，定义一个局部函数去拟合周围的数据点。此法灵活简单、精度较高、占内存少，但一次只算一点，速度较慢，是一种广泛采用的方法。,.,36,计算过程：1.对某待定点，找出周围的几个分块格网中的数据点，并将坐标原点移到该DEM格网点（xP,yP）,对任一数据点i有:,.,37,2.以待定点为中心，为半径作圆，落在圆内的数据点即被选用，来拟合一个二次曲面，数据点个数6，的大小可变，以便点数适中。,.,38,3.列出误差方程式，通常用二次曲面，则对第数据点可列出误差方程式,.,39,4.计算每一数据点的权（相关程度），权的大小与该点距待定点的距离有关，愈大，权愈小。计算公式有如下几种：,.,40,5.法化求解由于所以系数就是待定点的高程Zp。,.,41,实际应用中，也可用加权平均法，它是移动拟合法的简化。,.,42,5三角网数字地面模型,不规则三角网（TINTriangulatedIrregularNetwork）数字地面模型能很好地反映地貌的特征点、线。,.,43,一、三角网数字地面模型的构建,三角网数字地面模型的构建应基于最佳三角形条件，尽可能保证每个三角形是锐角三角形，避免出现过大的钝角或过小的锐角。1.角度判别法建立TIN该方法是当已知三角形的两个顶点后，利用余弦定理计算备选第三点为角顶点的三角形内角的大小，选择最大者对应的点为该三角形的第三顶点。步骤如下：,.,44,1.角度判别法建立TIN,（1）将原始数据分块，以便搜索所处理三角形邻近的点，而不必检索全部数据。（2）确定第一个三角形从离散数据点中任取一点（数据文件中第一点或左下角第一点），在其附近选取距离最近的一点作为三角形的第二点。然后对附近的点计算角,A,B,Ci,C1,c,bi,ai,.,45,1.角度判别法建立TIN,若则为该三角形的顶点（3）三角形的扩展由第一个三角形向外扩展，将全部点构成，并要保证三角网中没有重复和交叉的三角形。方法是依次对每个已生成的三角形新增加的两边，按角度最佳原则向外扩展，并进行是否重复的检查。,.,46,1.角度判别法建立TIN,向外扩展处理如已建P1、P2、P3三角形，如边扩展，应取异侧的点。P1P2方程为若备选点的坐标满足则在的异侧，作为备选扩展顶点。,P1,P2,P3,P,.,47,重复与交叉的检测由于任意一边最多只能是两个三角形的公共边，因此只需给每一边记下扩展的次数，如超过2，则该扩展无效。当所有生成的三角形新边都经扩展后，全部离散的数据点即连成了一个不规则三角网DEM。,1.角度判别法建立TIN,.,48,2.泰森多边形与狄络尼三角网,若区域上有n个离散点Pi(xi,yi)(i=1,2.n),若将用一组直线段分成n个相互邻接的多边形，满足下列条件（1）每个多边形内含且仅含一个离散点；（2）中任一点P(xi,yi)若位于所在多边形内，则满足：（3）若P(xi,yi)在Pi、Pj所在的两多边形公共边上，则,.,49,泰森多边形,.,50,2.泰森多边形与狄络尼三角网,则这些多边形称为泰森多边形。用直线连接每两个相邻多边形内的离散点而成的三角网称为狄络尼三角网。头以上定义知，泰森多边形是唯一的，每个多边形是凸多边形，任意两个多边形无公共区域。建立的三角网在均匀分布点的情况下，可避免不良三角形。,.,51,二、三角网数字地面模型的存储,1.直接表网点邻接关系的结构,特点是，存储量少，编辑方便，但计算量大，不便检索和显示,.,52,(2)直接表示三角形和邻接关系的结构,特点：检索网点拓扑关系效率高，但存储量大，不便编辑,.,53,（3）混合表示网点和三角形邻接关系的结构,.,54,5等高线的自动绘制,由规则格网自动绘制等高线有两个步骤1.等高线的跟踪由规则格网点高程内插等高线点，并将这些点按顺序排列。2.等高线的光滑利用这些顺序排列的等高线点的平面坐标进行插补，进一步加密等高线点并绘制成光滑的曲线。,.,55,一、等高线的跟踪,有两种方法方法一：对每条等高线边内插边排序；方法二：同一高程的等高线先内插所有等高线点，再逐一排列每条等高线的点。,.,56,一、等高线的跟踪,1.按每条等高线的走向顺序插点此法逐条按等高线的走向边搜索边插点，内插和排列同时进行。（1）确定等高线的高程先根据DEM中最低点的高程Zmin和最高点的高程Zmax计算最低等高线的高程zmin和最高等高线的高程zmax:,.,57,1.按每条等高线的走向顺序插点,上式中，Z为等高距，当时其它等高线高程为,.,58,（2）计算状态矩阵为了记录第条等高线通过DEM格网的情况，设置两个状态矩阵（）和（），分别表示等高线穿过DEM格网水平边和竖直边的状态。,1.按每条等高线的走向顺序插点,.,59,1.按每条等高线的走向顺序插点,.,60,1.按每条等高线的走向顺序插点,高程为zk的等高线通过格网（i,j）水平边的条件是：等高线高程介于此水平边两端点高程之间，即同理，高程为zk的等高线通过格网（i,j）垂直边的条件是：,也就是,.,61,则状态矩阵H(k)与V(k)的元素为为了避免上述判别式为零的情况，可将所有等于等高线高程的格网点上的高程加上（减）一个微小值0,1.按每条等高线的走向顺序插点,.,62,1.按每条等高线的走向顺序插点,（3）搜索等高线的起点与边界相交的等高线为开曲线，不与边界相交的等高线为闭曲线。通常先跟踪开曲线，即沿DEM的四边搜索，所有的元素都对应着一条开曲线的起点（终点）。,.,63,1.按每条等高线的走向顺序插点,在搜索到一个开曲线的起点后，要将其相应的状态矩阵元素置零，以免重复搜索。处理完开曲线后，再处理闭曲线。此时可按先列后行（先行后列）的次序搜索内部格网的竖直边和水平边，搜索到的第一个等高线通过的边即闭曲线的起点边，也是终点边，对应的矩阵元素仍保留原值1。,.,64,（4）内插等高线点等高线点的坐标一般采用线性内插。格网（i,j）水平边上等高线点坐标（xp,yp）为,1.按每条等高线的走向顺序插点,i,j,xi,p,yj,x,zk,xp,q,i+1,j,.,65,格网（i,j）竖直边上等高线点坐标（xp,yp）为（5）搜索下一个等高线点当等高线起点找到后，继续搜索该等高线与格网边的下一个交点。为此将每一格网边编号为1，2，3，4，则等高线的进入边号IN有4种可能。,1.按每条等高线的走向顺序插点,.,66,1.按每条等高线的走向顺序插点,设IN1，按一定方向（如逆时针）搜索等高线的离去边号OUT。当vi+1,j1时，OUT2，并令vi+1,j0；否则，搜索下一边。找到后，进入下一格网，继续以上算法，直到状态矩阵H(k)与V(k)的元素全为零，此时高程为zk的等高线就全部被搜了索出来。,1,2,3,4,i,j,i+1,j,.,67,1.按每条等高线的走向顺序插点,当格网中有两条高程为zk的等高线通过时,.,68,1.按每条等高线的走向顺序插点,当等高线过特征线时，必须用特征线与格网交点与相应的格网点内插等高线点。,.,69,1.按每条等高线的走向顺序插点（方法一）,（6）搜索等高线的终点对开曲线，当一个点是边界上的点时，该点即为终点；对于闭曲线，当一个点也是该等高线的第一点是，该点即为终点。,.,70,2.整体内插出整个DEM所有等高线点，然后逐条排列、存储,将离散的等高线点顺序地排列起来，可按两个条件进行。（1）方向条件：要求从已经排好