基于局部曲面变化因子的散乱点云数据精简处理方法
【技术领域】
[0001] 本发明属于逆向工程中计算机=维建模技术领域，为了提高=维建模效率，主要 针对由产品数字化测量设备得到的=维散乱点云数据，进行的散乱点云数据精简处理。
【背景技术】
[0002] 光学非接触式测量设备有数据获取快，数据丰富的优点，不仅能精确的包含被测 物的全部形貌特征数据，还包含法线方向、反身强度等信息。在逆向工程的应用邻域中得到 越来越广泛的应用。
[0003] 随着光学=维测量技术的发展，测量设备的分辨率越来越高，一次测量采集可得 到多达上百万个点数据，但是并不是所有的点数据都可W被用于后续的曲面重建等工作， 如在特征平缓的地方，很多点数据都是冗余的；同时，过多的点数据会使得在计算机的储存 和计算等操作中消耗过多的时间，在曲面重建等工作中的效率会大大降低，因此在不影响 曲面重建等工作的基础上，需要首先对点云数据进行精简处理。针对点云精简处理问题，国 内外研究学者做了大量的研究，并产生了很多优秀的研究成果。目前，其处理方法主要分为 两大类：均匀精简的方法和特征精简的方法。
[0004] 1.均匀精简的方法
[0005] 均匀精简方法不考虑细节曲率等特征，基于距离等方式实现均匀的精简点云数 据，比较有代表性的方法有均匀网格法和包围盒法。均匀网格法是用投影的方式，将点云投 影到一个小栅格里，求解靠近网格中屯、位置的点所对应的=维点，W此代替网格中其他的 投影点多对应的=维点，从而实现点云精简。精简的比例取决网格的大小，网格越大，精简 的比例也越大。包围盒法是分别寻找散乱点云中=个坐标的最大值最小值，得到能包围所 有点云的最小包围盒的边长，并W此按照一定的空间划分策略，将此包围盒划分成大小相 等的小立方体栅格，并将点云分配到每一个小立方栅格中，保留每个小立方体栅格中靠近 中屯、的=维点数据，或取小立方体栅格中点云的重屯、点，来代替小立方体栅格中的其他点 云数据，从而达到精简点云数据的目标。
[0006] 均匀精简的效果对于均匀的点云数据效果比较好，但是对于复杂的自由曲面模型 的点云数据，其曲率变化较大，如果将栅格划分过小，则保留的数据点较多，不能达到精简 目的；如果将栅格划分过大，则保留的数据点较少，会破坏棱角特征，甚至会出现孔洞。
[0007] 2.特征精简的方法
[000引特征精简的方法是充分考虑点云数据模型的细节特征，在实现点云数据精简的情 况下，保护细节特征部分不被删除，主要的方法有非均匀网格法、基于特征点的点云精简法 和基于曲率特征的点云数据精简方法。其中基于曲率特征精简点云最为成熟的一种方法； 首先需要散乱点云拓扑关系的建立，即建立散乱点云的k邻域信息，然后通过一定的方法 来估算k邻域的曲率大小，根据估算得到的曲率的大小变化，对点云数据进行精简，在曲率 小的地方，保留较少的点，曲率大的地方，保留较多的点，从而达到基于特征精简点云数据。
[0009] 该种基于曲率特征的精简方法可W在最大量保存实物模型的细节特征的前提下， 尽可能的增加精简效率，但是存有w下不足：
[0010] 一方面是在建立散乱点云k邻域信息方面，主要可W分为=类方法：一是利用采 样点集Voronoi图实现k个最近点捜索，但主要应用于二维散乱点集，且其计算量仍然很 大，在=维应用领域中，并不能适合所有类型的点云数据；二是基于树的层级结构捜索k邻 域，如KD-Tree、八叉树等方法，但是该些方法在捜索不同层的相邻叶节点的时候就会变得 特别复杂，若邻域点和采样点处于不同的层级结构时，该时会导致捜索范围扩大，致使捜索 效率大大降低；S是基于立方体小栅格捜索k邻域，利用不同的划分空间策略，将点云空间 分成许多大小相同的立方体栅格，按照一定的捜寻扩展方式，快速捜索k邻域，且该些方法 都是按照固定的k值捜寻，容易使其中的捜索循环陷入死循环，捜索效率较慢，且在一些边 界位置，用固定的k邻域处理点云数据会导致更大的错误，因而会影响整体的点云数据精 简的速度和精度。
[0011] 另一方面是在估算曲率方面，主要有两种方法；一种方法，是建立其他模型替代原 先的S维点云数据来估算曲率，如局部坐标系中的二次参数曲面模型、S角网格模型等，该 种方法过程复杂，且不能保证转换后的模型的准确性，另一种直接采用原始点云模型估算 曲率，大部分局限在规则点云的条件之上，对散乱点云模型直接估算曲率计算复杂，设及很 多大参数和矩阵计算，因而也会严重影响点云数据精简的整体效率。
【发明内容】
[0012] 针对现有技术存在的问题，本发明克服了现有的点云数据精简方法中的不足之 处，提出一种基于局部曲面变化因子的散乱点云数据精简处理方法，通过建立动态的点云k 邻域，快速建立散乱点云的拓扑关系，W局部曲面变化因子为规则精简点云数据，实现基于 特征保护的点云数据精简方法，从而提高点云数据处理、曲面重建的效率和精度。
[0013] 为了解决上述技术问题，本发明提出的一种基于局部曲面变化因子的散乱点云数 据精简处理方法，包括W下步骤：
[0014] 步骤一、读取测量得到的点云数据；
[00巧]步骤二求取点云的中屯、点，包括：
[0016] (2-1)根据点云的密度P，W及点云数据总数N，中屯、点最大的邻域点数K、点云数 据的数据范围建立立方体栅格，具体内容如下：
[0017] 根据散乱点云中坐标中的最大值和最小值，得到最小包围盒的边长及其体积V，由 式（1)求得初次划分的立方体栅格边长L。：
【主权项】
1. 一种基于局部曲面变化因子的散乱点云数据精简处理方法，其特征在于，包括以下 步骤： 步骤一、读取测量得到的点云数据； 步骤二、求取点云的中心点，包括： (2-1)根据点云的密度P，以及点云数据总数N，中心点最大的邻域点数K、点云数据的 数据范围建立立方体栅格，具体内容如下： 根据散乱点云中坐标中的最大值和最小值，得到最小包围盒的边长及其体积V，由式 (1)求得初次划分的立方体栅格边长L。：
式⑵和式⑶中，含有数据点的栅格总数，a为调控因子；由式⑵和式（3) 得到二次划分的立方体栅格边长L，并以该二次划分的立方体栅格边长L划分立方体栅格， 控制立方体栅格中心点的个数以此保证中心点邻域覆盖全部点云数据； (2-2)求得立方体栅格中心的点，并作为立方体栅格中心点； 步骤三、基于立方体栅格搜寻中心点的动态k个邻域点，从而建立散乱点云的拓扑关 系，具体内容如下： (3-1)搜寻并扩展搜寻：首先，初始搜寻二次划分的立方体栅格，并计算二次划分的立 方体栅格的中心点距六个面的距离d[i]，iG[〇, 5]，并从小到大排序d[0]?d[5]; 如果在该二次划分的立方体栅格中没有搜寻到设定的k个邻域点，则以当前搜寻的中 心点为圆心，依次以递增的radius为半径作一个球，以与该球相交的栅格作为新的搜索区 域进行扩展搜寻，直到搜寻到设定的k个邻域点为止； 其中，radius的递增原则是：radius=d[i]+LXn(4) 式（4)中，0彡i彡5,0彡n彡5;当n= 0时，1彡i彡5;当1彡n彡5时，0彡i彡5，i,n均属于自然数， (3-2)建立动态k邻域点：根据点云数据特性、精度要求和运算时间，确定k的阈值为 [M,K]； 步骤四、采用协方差分析法计算中心点k邻域的曲面变化因子，具体内容如下： 设Pi为点云数据中心点，其k邻域为Knn(p) = {q^，qi;1. . .qi，k}，通过分析k邻域Knn(Pi)的统计特性，求得点云数据中心点Pi的法向量、曲面变化量的属性， 点云数据中心点a.和k邻域Knn(n,_)的质心为：
k邻域拟合的平面过该质心L从而得到关于点云数据中心点？1的3X3协方差矩阵 Cov：
从矩阵Cov得到点云数据中心点Pi处的统计特性，点云数据中心点pk邻域Knn(pJ 关于所述质心^的距离平方差的偏离程度； 所述3X3协方差矩阵Cov的特征向量有： Cov?a』=入』?a』，jG{〇, 1，2} (7) 由式（7)得到三个非负特征值A。，A2，其中，彡AA2;与所述三个非负特 征值--对应的特征向量为a&aa2; 其中，所述三个非负特征值反映了k邻域Knn(Pi)到所述质心&的偏离程度，表述为：
其中，三个非负特征值特征值的比例反映点云数据中心点Pi处的曲面变化情况，将所 述三个非负特征值特征值的比例定义为曲面变化因子I，计算公式为：
曲面变化因子I的值越大，点云数据中心点Pi处的曲面变化程度越大； 步骤五、确定每个立方体栅格中心点k邻域的精简率Rotio，并在k邻域范围内均匀精 简：
式（10)中，r为精简率调控因子用以控制点云整体的精简率。
【专利摘要】本发明公开了一种基于局部曲面变化因子的散乱点云数据精简处理方法，包括：1)读取测量得到的点云数据；2)求取点云的中心点；3)基于立方体栅格搜寻中心点的动态k个邻域点,从而建立散乱点云的拓扑关系；4)采用协方差分析法计算中心点k邻域的曲面变化因子；5)确定每个立方体栅格中心点k邻域的精简率，并在k邻域范围内均匀精简。本发明通过建立散乱点云动态k邻域点信息，建立散乱点云的拓扑关系。用曲面变化因子替代复杂的曲率计算。通过曲面变化因子ξ调整精简率Rotio，实现在k邻域范围内均匀精简，不仅能保护曲率大的细节特征，在精简程度较大时，也保护了曲率较小的平面特征。使点云数据处理、曲面重建的效率和精度得以提高。
【IPC分类】G06T17-00
【公开号】CN104616349
【申请号】CN201510050544
【发明人】林滨, 盛金月, 亓振良
【申请人】天津大学
【公开日】2015年5月13日
【申请日】2015年1月30日