基于散乱点高程采样数据的数字高程地形模型生成方法
【专利摘要】本发明公开了一种基于散乱高程采样数据的数字高程地形模型(DEM)生成方法。该方法从一组平面散乱采样点及其所对应的高程数据出发，首先确定DEM模型的采样区域Ω，通过为每个采样点定义邻近点集、支撑半径和支撑区域，进而为各采样点分别定义一个局部逼近函数和一个局部支撑函数；对于区域Ω内的任意一点，搜索支撑区域覆盖该点的各个采样点，并将它们的局部逼近函数用相应的局部支撑函数混合起来，定义覆盖区域Ω的全局高程函数；最后用全局高程函数计算平面上任意布局的DEM网格结点的高程值，生成所需的DEM模型。本发明方法能够过滤采样数据中的噪声，具有较好的内插能力，能够处理任意规模的散乱点采样数据。
【专利说明】基于散乱点高程采样数据的数字高程地形模型生成方法
【技术领域】
[0001] 本发明涉及一种基于散乱点高程采样数据的数字高程地形模型（DEM)生成方法， 属于数字化地形建模【技术领域】，尤其涉及一种由散乱点高程采样数据构造数字高程地形模 型（DEM)的方法。
【背景技术】
[0002] 建立DEM的方法有多种：从数据源及采集方式上讲有：（1)直接由地形测量得到， 例如用GPS、全站仪、野外测量等手段获取；（2)根据航空或航天影像，通过摄影测量途径获 取，如立体坐标仪观测及空三加密法、解析测图、数字摄影测量等等；（3)从现有地形图上 采集，如格网读点法、数字化仪手扶跟踪及扫描仪半自动采集，然后通过内插生成DEM等方 法。DEM内插方法很多，主要有整体内插、分块内插和逐点内插三种。整体内插（参见：周 兴华等，DEM内插方法与精度评定，测绘科学，2005, 30 (5) :86-88)在整个区域内用一个数 学函数式来表达地形曲面，该方法的拟合模型是由研究区内所有采样点的观测值（即采样 点的高程值）建立的。该方法得到的整体内插函数保凸性较差，不易得到稳定的数值解， 多项式系数的物理意义不明显，解算速度慢且对计算机容量要求较高，难以提供内插区域 的局部地形特征。分块内插（参见：陈国良等，一种由等高线模型生成规则格网模型的算 法，中国图象图形学报，2007, 12(6) : 1110-1113)将采样点所在空间分成若干大小相同的 块，对各分块使用不同的插值函数。其中，线性内插和双线性多项式内插是两种典型的分块 内插方法。线性内插法算法简单、易于实现、运算速度快。但由于此内插法的原理采用的是 选择邻近点内插，而不关心邻近点的特性，如果邻近点是极值点，那么内插出的结果与实际 情况的出入会较大。双线性内插法得到的模型表面比较平滑并且精度较高，缺点是地形的 特征点线面被网格平均化和平滑化，极值容易丢失。逐点内插（参见：李世平等，数字高程 模型的建立与分析应用，辽宁工程技术大学学报（自然科学版），2008, 5(27) :31-33)以待 插点为中心，定义一个局部函数去拟合周围的数据点，数据点的范围随待插位置的变化而 变化，因此又称为移动拟合法。逐点内插法对未知点高程值的获得只考虑距离因素，不考虑 距离之外的其它空间因素，也不关心未知点与周围其它点的拓扑关系和固有联系，该点的 高程值过多地受其邻近点高程值的影响，因此采样点坐标或采样点所在的坐标系的微小变 化都会使选点结果差别很大，结果可能造成数字高程模型表面的不连续。
[0003] 目前，DEM模型的拓扑结构有规则网络结构和不规则三角网（Triangular Irregular Network,简称TIN)结构两种。常用的DEM生成方法是：先在采样平面上构造网 格，然后再在所构造网格的基础上，通过线性和双线性内插建立DEM。用规则方格网高程数 据记录地表起伏信息的优点是：采样点的位置信息可隐含，无需全部作为原始数据存储，后 继的数据处理也比较方便，其缺点是数据采集量大，因为网格点不是特征点，一些微地形可 能没有被记录。TIN结构数据的优点是：能以不同层次的分辨率来描述地表形态。与规则 格网数据模型相比，TIN模型在某一特定分辨率下能用更少的空间和时间更精确地表示更 加复杂的地形表面。
【发明内容】
[0004] 本发明的目的是提供一种基于散乱点高程采样数据的数字高程地形模型（DEM) 生成方法，该方法能够过滤采样数据中的噪声，具有较好的内插能力，能够处理任意规模的 散乱点采样数据。
[0005] 为实现以上目的，本发明采用如下技术方案：
[0006] -种基于散乱点高程采样数据的数字高程地形模型（DEM)生成方法，其具体步骤 是：
[0007] a.输入平面采样点集合及相应的高程值集合：输入由η个平面采样点组成的采样 点集合，以及各个采样点处的高程值采样数据；
[0008] b.确定数字高程地形模型（DEM)的采样区域：遍历平面采样点集内的各点，确定 各采样点在平面上的分布范围；确定将要生成的数字地形模型（DEM)的采样区域；
[0009] c.确定任意采样点的邻点及支撑半径：对于每个采样点，计算与其距离最近的若 干个邻点；计算采样点到其各邻点的最大距离，并以此距离来定义其支撑半径；
[0010] d.构造各采样点的局部逼近函数：对于每个采样点，由其自身及其邻点、以及相 应的高程采样值，构造该采样点的局部高程逼近函数；
[0011] e.确定各个局部逼近函数的支撑函数：对于每个采样点，用该采样点的支撑半径 为其定义支撑区域；然后构造定义于支撑区域上的支撑函数；
[0012] f.构造全局高程计算函数：根据支撑函数的局部覆盖特性，将相应局部逼近函数 混合为一个全局高程计算函数，使得采样区域内的任一点均可由该函数计算出高程值；
[0013] g.构造 DEM数据：将采样区域分割为适当的网格拓扑结构，通过全局高程函数计 算网格各顶点处的高程值，生成数字化高程模型的DEM数据。
[0014] 本发明的方法考虑到了更多的采样数据之间的联系，顾及到了距离、位置、范围等 因素，同时可以有效地过滤采样高程数据中的噪声，既能生成规则网络结构的DEM模型，也 能生成不规则三角网（TIN)结构的DEM模型。此外，该方法能够具有较好的内插能力，能够 处理任意规模的散乱点采样数据。
【专利附图】
【附图说明】
[0015] 图1是本发明的方法流程图。
[0016] 图2是采样点的邻点与支撑半径示意图，其中，采样点Ptl共有8个邻点Pi (i = 1，2, Λ8)。
[0017] 图3是支撑函数示例，其中，（a)非均匀衰减的类高斯函数、(b)均匀衰减的线性函 数，随着自变量从〇变化到R，函数值均由1衰减到0。
[0018] 图4是支撑区域覆盖关系示意图，其中，Pl，p2, P3三个采样点的支撑区域覆盖点P。
[0019] 图5是DEM的网格拓扑结构示意图，其中，（a)是规则四边形网格、（b)均匀三角网 格、（c)具有细节层次的三角网格、（d)具有任意拓扑结构的三角网格。
【具体实施方式】
[0020] 下面结合附图和实施例来进一步说明本发明。
[0021] 如图1所示，基于散乱点高程采样数据的数字高程地形模型（DEM)生成方法，其具 体步骤是：
[0022] a.输入平面采样点集合及相应的高程值集合：输入由η个平面采样点组成的采样 点集合，以及各个采样点处的高程值采样数据。
[0023] b.确定数字高程地形模型（DEM)的采样区域：遍历平面采样点集内的各点，确定 各采样点在平面上分布的范围，进而在该范围基础上确定将要生成的数字地形模型（DEM) 的采样区域。
[0024] c.确定任意采样点的邻点及支撑半径：对于每个采样点，查找并记录与其距离最 近的若干个邻点（不同采样点的邻近点数可以不同）；计算其到各个邻点的最大距离d，并 以此距离进一步定义该采样点的支撑半径R，并使R > d。如图2所示，采样点Ptl的邻近点 共有8个，即Pi (i = 1，2, Λ 8)，其支撑半径取为采样点Ptl到采样点p6的距离。
[0025] d.构造各采样点的局部逼近函数：对于每个采样点，由其自身及其邻点、以及相 应的高程采样值，构造该采样点的局部高程逼近函数，使该函数在上述各采样点处的值逼 近相应的高程采样值。
[0026] e.确定各个局部逼近函数的支撑函数：对于每个采样点，用该采样点的支撑半径 定义一个平面区域，称其为该采样点的支撑区域；然后构造定义于支撑区域上的支撑函数， 支撑函数是一个径向单调非増函数，其值在该采样点处取为1，并沿径向外逐步衰减，到达 支撑区域边界时取值为0。图3给出了两个支撑函数的示例，其中采样点到支撑区域边界的 距离为R。
[0027] f.构造全局高程计算函数：根据支撑函数的局部覆盖特性，将相应局部逼近函数 混合为一个全局高程计算函数，使得采样区域内的任一点均可由其计算出高程值。如图4 所示，设P为DEM采样区域内的任一点，并且有且只有P 1、p2、p3三个采样点的支撑区域覆盖 着点P，那么P点处的高程值由 Pl、p2、P3的局部逼近函数与局部支撑函数共同确定。
[0028] g.构造 DEM数据：根据需要，将采样区域分割为适当的网格拓扑结构（如图5所 示，采样区域为一平面矩形区域；在区域上依据构造了规则四边形网格、均匀三角网格、具 有细节层次的三角网格、以及具有任意拓扑结构的三角网格），通过全局高程函数计算网格 各顶点处的高程值，生成数字化高程模型的DEM数据。
[0029] 实施例
[0030] 本发明方法的输入和输出分别是：
[0031] 输入：一组平面散乱采样点P ={ (X1, Y1)，（x2, y2)，Λ，（Xn, yn)}及其所对应的高程 值H = Oiuh2, Λ，hn}，其中：采样点（Xpyi)处的高程值为hi。
[0032] 输出：在采样点分布的确定区域范围内，由采样点高程值所确定的数字化高程地 形模型DEM。
[0033] 具体的实施步骤为：
[0034] 1.输入平面采样点集合及相应的高程值集合：输入由η个平面采样点组成的采样 点集合P= {Pi(Xi，yi) e R2, i = 1，Λ，η},以及这η个采样点处的高程值集合H= Oii e R, i =1, Λ , η}，并令采样点Pi (Xi, y)处的高程值为Iii (i = 1, Λ , η)。
[0035] 2.确定数字化地形DEM模型的采样区域：遍历平面点集P内的各个采样点，求出 各点在橫坐标和纵坐标上分布范围的最大值和最小值，即
【权利要求】
1.基于散乱点高程采样数据的数字高程地形模型生成方法，其特征在于，具体步骤 是： a. 输入平面采样点集合及相应的高程值集合；输入由n个平面采样点组成的采样点集 合，W及各个采样点处的高程值采样数据； b. 确定数字高程地形模型的采样区域；遍历平面采样点集内的各点，确定各采样点在 平面上的分布范围；确定将要生成的数字地形模型的采样区域； C.确定任意采样点的邻点及支撑半径：对于每个采样点，计算与其距离最近的若干个 邻点；计算采样点到其各邻点的最大距离，并W此距离来定义其支撑半径； d. 构造各采样点的局部逼近函数；对于每个采样点，由其自身及其邻点、W及相应的 高程采样值，构造该采样点的局部高程逼近函数； e. 确定各个局部逼近函数的支撑函数：对于每个采样点，用该采样点的支撑半径为其 定义支撑区域；然后构造定义于支撑区域上的支撑函数； f. 构造全局高程计算函数：根据支撑函数的局部覆盖特性，将相应局部逼近函数混合 为一个全局高程计算函数，使得采样区域内的任一点均可由该函数计算出高程值； g. 构造DEM数据；将采样区域分割为适当的网格拓扑结构，通过全局高程函数计算网 格各顶点处的高程值，生成数字化高程模型的DEM数据。
【文档编号】G06T11/00GK104463924SQ201410634465
【公开日】2015年3月25日 申请日期:2014年11月12日 优先权日:2014年11月12日
【发明者】庞明勇, 董晓芬 申请人:南京师范大学