您的当前位置:首页正文

面向地形数据的点云简化算法

来源:化拓教育网
第35卷第3期 2015年6月 大地测量与地球动力学 Journal Of Geodesy and Geodynamies Vo1.35 No.3 Jun.,2015 文章编号:1671—5942(2015)03—0424 04 面向地形数据的点云简化算法 叶珉吕 花向红。 1佛山市城市规划勘测设计研究院,佛山市岭南大道北62号,528000 2武汉大学测绘学院,武汉市珞瑜路129号,430079 摘 要:针对地形点云数据量大、表面特征复杂多样等特点,提出面向地形数据的点云简化算法。基于K—D Tree搜索各点K邻域,构建点集空间拓扑关系;应用移动最小二乘法计算各点曲率,通过曲率的划分,在平缓 区域按距离进行简化,保证整个算法的效率;在突变区域根据曲率简化,确保曲率变化大的关键特征信息不丢 失,从而实现点云数据的简化。利用基于熵理论的定量评价方法,通过实例验证该方法的可行性和普适性。 关键词:三维激光扫描;点云;地形数据;简化算法;信息熵 中图分类号:P207 文献标识码:A 三维激光扫描技术获得的点云数据量越大、 点位越密集,对扫描物体的描述则越精确。然而, 效率高。 1.2点云数据曲率计算 庞大的点云数据会给后续处理及存储、传输与显 示等造成不利影响,因此需要对原始点云数据进 点云数据法向量、曲率等几何信息能够很好 地反映点云的表面特征,为了使简化结果能够较 行简化处理。目前,基于离散点云数据的简化方 法主要分为两类:均匀简化方法与基于曲率的简 化方法。均匀简化方法u 简单、高效,但是简化 过程没有顾及点云表面特征,容易造成点云特征 信息丢失,因此不适用于特征复杂的点云数据的 简化。基于曲率的简化方法能够较好地保留点云 特征,但是简化效率较低,且对于特征简单和曲率 好地保留点云特征,首先需要计算点云数据的几 何信息。点云几何信息是在点集拓扑结构的基础 上,通过曲面拟合的方法计算得到,点集的拓扑结 构一般通过K—D Tree、八叉树等方法搜索各点K 邻域而建立。点云数据中任意一点P的K邻域 为到该点的欧氏距离最近的K个点所组成的点 集,K—D Tree这种数据结构能够很好地提高空间 变化较小的区域容易过度简化,使得简化后的点 云数据出现分布不均匀的现象。对于数据量少、 搜索邻近点的效率,适用于海量的、空问分布不均 匀的离散点集。因此,本文基于K—D Tree进行 特征简单规则的模型、零件、建筑物等点云数据, 上述传统简化方法能够获得较好的简化效果,但 是对于海量的、表面特征复杂多样的地形点云数 据,目前并没有高效的、高精度的简化方法。对 此,本文提出了面向地形数据的点云简化算法,在 K邻域搜索。寻找任意一点P的邻域的基本流 程如下: 1)如果是根结点,从根结点开始搜索。 2)如果是叶子结点,采用递归法返回所有与 点P邻近的点,并计算最小距离值。 3)判断是否分割面另一边上有点更接近点 高保真、高精度的前提下快速对点云数据进行简 化,并利用定量的评价方法,通过实例验证算法的 可行性。 P,它是通过判断以点P为球心、最小距离值为半 径的球体与分割平面是否相交来判定的。如果相 l点云简化算法 1.1地形数据简化原则 交,分割面另一边上可能有邻近点,那么算法必须 向P点所在区域的上层区域回溯,从当前节点所 在K—D Tree分支移动到另一分支上去搜寻更近 点;如果不相交,那么算法继续沿着当前K—D Tree分支运行,当前节点的分割平面另一边整个 分支则被淘汰。 实用的点云简化算法应在保证失真较小的前 提下,最大限度地压缩点云数量,且点云的简化结 果能满足应用的精度要求,同时算法需简洁,执行 收稿日期:2014 06—04 项目来源:国家自然科学基金(4l174010);精密工程与工业测量国家测绘地理信息局重点实验室开放基金(PF201卜7)。 第一作者简介:叶珉吕,硕士,主要从事三维激光扫描点云数据处理研究,E mail:yeminlv@163.corn。 第35卷第3期 叶珉吕等:面向地形数据的点云简化算法 425 4)完成所有节点的搜寻,点P的邻域点集建 立。 所有点的邻域点集构建完成后,选取移动最 小二乘法计算点云数据的曲率信息 。移动最小 二乘法属于一种超平面估算方法,通过给定一个 局部高阶多项式,为采样点提供一个逼近或插值 参考平面,以这个超平面上该点处的曲率作为当 前采样点的曲率。具体步骤包括 : 1)向量场估算。法向量是构建点云移动最小 二乘曲面的前提。设Q为散乱点云数据,任意一 点 ∈Q,向量场 (z)可由下式给出: 一 器 ㈩ 其中,', 为点i法向量,Oix)一e—l—l z/hz为高斯 加权函数,h是与点间距离相关的一个常数。 2)移动最小二乘曲面投影。移动最小二乘法 通过拟合为采样点提供一个参考曲面,并将采样 点投影到曲面上。以曲面上该点处的曲率作为当 前采样点的曲率,MLS曲面隐式公式为_4]: gix)一∑ Q20(x)[ix— ) (z)](2) 3)曲率计算。MI S曲面确定后,在投影曲 面上利用隐式曲面的曲率计算公式计算采样点 的曲率: c一= ㈦ gix 一( ) ㈩ H(gix))一 ( gix))一 3gix) 3giz) 3g(x) 3x3x 3x3v azaz 3gix) 3g(z) 3g(z) aV3x 9y3y aV3z (5) 3gix) 9g(37) 3g( ) 3z3x 3V3z 3z3z 1.3简化算法实现步骤 在获取点云曲率信息后,首先通过设置一定 阈值(本文选取曲率平均值)将其划分为平缓区域 与突变区域。曲率小于阈值,则为平缓区域;反 之,则为突变区域。 对于特征简单的平缓区域,其所包含的特征 信息较少,可简单地根据距离原则进行简化,无需 顾及曲率信息,提高效率的同时避免了过度简化。 两点空间距离计算公式为: Dis—ff(37 一 )。+iy 一Y ) +( 一 ,) (6) 对于突变区域,特征信息较为丰富,具有各种 类型的地形特征点,其共同特点是与周围邻近点 的曲率差别较大,可以根据此特点进行简化,保留 点云的特征信息。两点曲率差值计算公式为: D。…一l C 一Cj『 (7) 其中,C 表示点i的曲率,C,表示点i邻近点 的曲率。 综上所述,本文简化算法的基本原理为:在平 缓区域按距离进行简化,保证整个算法的效率;在 突变区域根据曲率简化,确保曲率变化大的关键 特征信息不丢失。其具体实现步骤如下: 1)基于K—D Tree搜索各点K邻域,构建点 集空间拓扑关系; 2)在点集拓扑结构的基础上,通过移动最小 二乘法计算各点曲率; 3)根据阈值判断某点曲率的划分,如为平缓区 域,则执行步骤4);如为突变区域,则执行步骤5); 4)计算某点与其邻域内K个邻近点的距离, 删除与其距离最近的点; 5)计算某点与其邻域内K个邻近点的曲率 差值,删除与其曲率差值最小的点; 6)遍历每个点,输出简化后的点云数据,算法 结束。 2简化结果定量评价方法 目前,点云数据简化结果一般通过简化前后 的图像对比定性地进行评价,简化精度无法定量 分析,评价结果具有较大的主观性。文献[6]通过 熵理论描述点云数据的特征信息,指出某点熵越 大,该点所在局部区域的无序程度越高,该点提供 的信息量越大,越能精确反映扫描目标的细节。 因此,可以通过计算点云数据的熵值,定量评价简 化结果。曲率能够很好地反映点云的表面特征信 息,结合熵理论及点云数据曲率信息,计算某点熵 值的公式为 : H 一一P log2p 一 log2PJ (8) J=1 k P ===C /(Ci+∑Cj) (9) J=1 k P,一Cj/(C +∑Cj) (10) J=1 其中,C 表示点i的曲率,C 表示点i邻近点 的 曲率,P 与P 分别为点i及J的曲率概率分布。 则点云数据熵值等于各点所含的熵值之和: H一∑H (11) 2—1 点云数据熵值越大,其所含特征信息量越大, 第35卷第3期 叶珉吕等:面向地形数据的点云简化算法 427 型依然具有较高的精度,而且简化后的点云数据 大大减少了建模时间,提高了点云数据处理的效 率,从另一个角度证明了算法的正确性。 [43 Levin D. Mesh Independent Surface Interp01ation[c].Ge ometric Modeling for Scientific Visualization,Berlin,2004 [5] 杨荣华.地面三维激光扫描点云角度分辨率与数据处理模型 研究[I)_.武汉:武汉大学,2011(Yang Ronghua.Research on Point Cloud Angular Resolution and Processing Model of Ter— 4 结 语 针对地形点云数据量大、表面特征复杂多样 等特点,本文提出面向地形数据的点云简化算法, 并利用基于熵理论的定量评价方法,通过实例验 证该方法能在高保真、高精度的前提下高效地对 点云数据进行简化,具有较高的可行性与普适性。 restrial Laser Scanning[D].Wuhan:Wuhan University,201I) [6] 武剑洁.基于点的散乱点云处理技术的研究[D].武汉:华中 科技大学,2004(Wu Jianjie.Research of Point—Based Tech niques on Unorganized Point Cloud[D].Wuhan:Huazhong University of Science and Technology,2004) [7] 孟庆生.信息论[M].西安:西安交通大学出版社,1986 (Meng Qingsheng.Information Theory[M].Xi’an:Xi’an Jiaotong University Press.1 986) 该方法能够为后续的应用提供有效的数据信息, 节约后续工作的处理时间和硬件资源。 参考文献 [1]Weir D J,Milroy M J,Bradley C,et a1.Reverse Engineer ing Physical Models Employing Wrap—-around B——Spline Sur— -E8] 葛源坤.基于曲率特征信息的散乱点云数据预处理技术研 究[D].成都:西南交通大学,2012(Ge Yuankun.Research on Data Pre—Processing Technology of Scattered Point Cloud Based on Curvatures Feature[D].Chengdu:Xinan Jiaotong University,2012) faces and Quadrics[J].Journal of Engineering Manufac— ture,1996,210(2):147—157 [9] 尹婷.三维激光扫描数据处理技术的研究[D].武汉:武汉理 工大学,2010(Yin Ting.Research on 3D Scanning Data r2] Martin R R,Stroud I A,Marshal A D.Data Reduction for Processing Techniques[D].Wuhan:Wuhan University of Technology,2010) Reverse Engineering[R].Computer and Automation Insti— tute of Hungarian Academy of Science,1996 [1o] 贺美芳.基于散乱点云数据的曲面重建关键技术研究[D]. 南京:南京航空航天大学,2006(He Meifang.Research on Key Technologies of Surfaces Reconstruction Based on [3]郑德华.点云数据直接缩减方法及缩减效果研究[J].测绘 工程,2006,15(4):27—30(Zheng Dehua.The Data Reduc— tion of Point Cloud and Analysis of Reduction Effect[J]. Engineering of Surveying and Mapping,2006,15(4):27—30) Scattered Point Cloud Data[D].Nanjing:Nanjing Universi— ty 0f Aeronautics and Astronautics,2006) Point Cloud Simplification Algorithm for Terrain Data YE Minlii HUA Xianghong 1 Foshan Urban Planning Design and Surveying Research Institute,62 North-Lingnan Road,Foshan 528000,China 2 School of Geodesy and Geomatics,Wuhan University,129 Luoyu Road,Wuhan 430079,China Abstract:Large amounts of data will have an adversely impact on point cloud’s processing,storage,trans— mission and display.In order to improve the point cloud’s data processing efficiency,it is necessary to simply point cloud data and reduce data redundancy.Due to the large quantity and complex surface features of terrain point cloud data,the point cloud simplification algorithm for terrain data is proposed.We use K-D Tree to search the point’s K neighborhood and the moving least squares method to calculate the curvature of each point.Through the vision of curvature,point cloud data is simplified based on distance in the flat area to im— prove the efficiency of the algorithm and on the basis of curvature in the steep area to retain the feature infor— mation.According to the quantitative evaluation method,based on the entropy theory and verified by the ex— ample,it is indicated that the proposed new simplification algorithm for terrain data can simply the data in high-fidelity with high-precision.It is highly efficient feasible and universally applicable. Key words:terrestrial laser scanning;point cloud;terrain data;simplification algorithm;information entropy 

因篇幅问题不能全部显示,请点此查看更多更全内容