一种基于点云特征分区的自适应点云精简方法技术

本发明专利技术涉及一种基于点云特征分区的自适应点云精简方法。首先对原始点云建立K

【技术实现步骤摘要】
一种基于点云特征分区的自适应点云精简方法


[0001]本专利技术属于点云处理与应用领域，特别是涉及一种基于点云特征分区的自适应点云精简方法。

技术介绍

[0002]随着虚拟匹配技术的发展与应用，由激光扫描仪扫描得到的点云模型质量的保证与提升越来越受到人们的关注，但点云数量的冗余极易造成点云处理效率低、重构模型质量低下等问题，所以点云精简处理就显得尤为重要。
[0003]常用的点云精简算法有随机采样法、包围盒采样法和曲率采样法等方法。随机采样法快速简单，其基本原理是通过随机数函数对点云数据进行选择并删除，该方法原理简单，但是精度低。包围盒采样法的精简算法是对点云进行最小立方体的构造，用立方体的中心点或重心点代替立方体中其他的点来进行下采样操作，精简效果与立方体尺寸的选取密切相关。基于曲率采样的点云精简算法原理是通过点云数据的曲率信息进行简化的曲率大的地方则数据点保留得多，曲率变化小的区域则数据点保留得少。
[0004]但采用任意一种精简算法对点云数据进行精简处理，无区别性地对点云进行下采样操作，极易造成点云模型的重要特征数据点误删除或过度简化等问题。所以Ji等人通过k近邻搜索设计了一种根据点权值简化算法，通过所设计的特征公式对点的重要性进行判断，最后采用八叉树对点云进行二次简化，方法思想即是通过特征对点云数据进行分类，而后再精简点云数据，这一思想在点云精简处理过程中越来越受到研究人员的关注。

技术实现思路

[0005]本专利技术针对现有技术的不足，提供了一种基于点云特征分区的自适应点云精简方法。即将点云先进行特征分类，依据点云特征对点云模型进行分区处理，有区别性地对点云数据进行简化操作，提高处理效率。其中点云轮廓线为强特征，包含了模型重要信息，一般保留不做精简，对于曲率变化较大的一般特征点云则采用曲率法进行简化处理，对于其他的非特征点云则采用包围盒法进行下采样操作。
[0006]为了达到以上目的，本专利技术提供点云精简技术方案包括以下步骤：
[0007]步骤1：基于K
‑
d树对目标点云建立索引；
[0008]步骤2：采用主成分分析法(Principal Component Analysis，PCA)计算点云法向量；
[0009]步骤3：采用弦法向量法计算点云主曲率k1、k2；
[0010]步骤4：通过主曲率计算点云平均曲率MC、点云整体平均曲率和标准偏差MC
std
；
[0011]步骤5：结合统计学原理设定阈值L与点云平均曲率对强特征点云进行识别并保留；
[0012]步骤6：计算点云K近邻法向量夹角平均值f
i
；
[0013]步骤7：依据法向量夹角阈值γ
i
进行一般特征点云区域与非特征点云区域进行区
分；
[0014]步骤8：采用曲率法对一般特征点云区域进行精简处理；
[0015]步骤9：采用包围盒法对非特征点云区域进行精简处理；
[0016]针对步骤1中对点云建立K
‑
D树，建立索引，提高点云后续处理效率。
[0017]针对步骤2中采用PCA算法计算点云法向量，计算过程如下：
[0018]对于点云中的每个扫描点p，搜索到与其近邻的k个相邻点，然后计算这些点的最小二乘意义上的局部平面P，此平面可以表示为：
[0019][0020]式中，为平面P的法向量，d为p到坐标原点的距离；p
i
为k邻域中任一搜索点到原点的方向向量，i为1～k中任一点编号。
[0021]可以认为由k个最近点拟合出的平面的法向量即当前扫描点的法向量。平面P的法向量可以由主成分分析(PCA)得到，由运算知p经过其K邻近域点的质心p0，且法向量满足先对式(2)中的协方差矩阵M进行特征值分解，求得M的各特征值，M最小特征值对应的特征向量即点p的法向量。
[0022][0023]通过PCA计算得到的点云法向量方向不确定，只能得到法向量所在的直线，无法确定法向量的最终指向，从而导致一个平面上点云法向指向不统一的问题，影响后续处理操作效果，本专利技术采用Hoppe提出的近似求解算法来进行法向量修正。
[0024]针对步骤3，点云表面上特定点的相邻点决定其局部形状，如果通过曲面拟合来估计曲率，可能会产生较大的误差，所以应该考虑法向量的贡献。所以本文采用弦法向量来估计一个点的曲率，其优点是只考虑一个相邻点，利用这一点来构造一条正截面曲线，进而构造法向截面圆，并根据两个点的位置和法向向量来估计特定点的法向曲率。具体计算过程如下：
[0025]对于点云中的每个点p，设p点的为法向量为N，使用点坐标和法向量来估计点p处的法向曲率：假设p点附近有m个近邻点，q
i
为点p第i个近邻点，q
i
的法向量为M
i
，设正交坐标{undefinedp，X，Y，Z}为点p的局部坐标系L，N表示p点的法向量。X和Y为正交的单位向量。在L中，p，q
i
，M
i
的坐标可以是(0，0，0)，q
i
为(x
i
，y
i
，z
i
)，M
i
为(n
x
，i，n
y
，i，n
z
，i)，则点p的法曲率可以通过其密切圆估计。
[0026]则p相对于q
i
的法曲率估计如下：
[0027][0028]式中，n为p点法向量，α是向量
‑
N和pq
i
之间的夹角，β是向量N和M
i
之间的夹角，等式(3)近似值也可表示为
[0029][0030]式中：
[0031]然后通过欧拉方程得出法向曲率与主曲率关系：
[0032][0033]式中：θ
i
+θ为点p过q
i
的法截线的切线与主方向的夹角，θ为向量e1和e2的夹角，θ
i
可以用点q
i
的局部坐标系计算，并且公式(5)可以写成一个优化问题：
[0034][0035]其中，m为p点附近点的个数，将公式(6)写成矩阵形式：
[0036][0037]其中M为矩阵：
[0038][0039]R为矩阵：
[0040][0041]对式(7)最小二乘拟合，可以获得μ的估计值，可以推断出主曲率k1、k2是矩阵：
[0042][0043]的特征值。把W的单位特征向量从局部坐标系L变换到全局坐标系，就得到主方向向量。
[0044]针对步骤4，点云平均曲率、整体平均曲率和标准偏差计算如公式(11)、(12)和(13)所示：
[0045][0046][0047][0048]其中，C
i
为k邻域中第i个点的曲率值。
[0049]针对步骤5，判断当前每个p
i...

【技术保护点】

【技术特征摘要】
1.一种基于点云特征分区的自适应点云精简方法，其特征在于，包括如下步骤：步骤1：基于K
‑
d树对目标点云建立索引；步骤2：采用主成分分析法计算点云法向量；步骤3：采用弦法向量法计算点云主曲率k1、k2；步骤4：通过主曲率计算点云平均曲率MC、点云整体平均曲率和标准偏差MC
std
；步骤5：结合统计学原理设定阈值L与点云平均曲率对强特征点云进行识别并保留；步骤6：计算点云K近邻法向量夹角平均值f
i
；步骤7：依据法向量夹角阈值γ
i
对一般特征点云区域与非特征点云区域进行识别与区分；步骤8：采用曲率法对一般特征点云区域进行精简处理；步骤9：采用包围盒法对非特征区域进行精简处理。2.如权利要求1所述的一种基于点云特征分区的自适应点云精简方法，其特征在于：所述步骤1对原始点云模型建立K
‑
d树索引，为后续点云模型相关参数计算做准备，提高点云处理效率。3.如权利要求1所述的一种基于点云特征分区的自适应点云精简方法，其特征在于：所述步骤2采用主成分分析法PCA计算点云法向量，为后续计算点云曲率与法向量夹角做准备，具体计算过程如下：对于点云中的每个扫描点p，搜索到与其近邻的k个相邻点，然后计算这些点的最小二乘意义上的局部平面P，此平面表示为：式中，为平面P的法向量，d为p到坐标原点的距离；p
i
为k邻域中任一搜索点到原点的方向向量，i为1～k中任一点编号；由k个最近点拟合出的平面的法向量即当前扫描点的法向量；平面P的法向量由主成分分析得到，由运算知p经过其K邻近域点的质心p0，且法向量满足先对式(2)中的协方差矩阵M进行特征值分解，求得M的各特征值，M最小特征值对应的特征向量即点p的法向量：通过PCA计算得到的点云法向量方向不确定，采用Hoppe提出的近似求解算法来进行法向量修正，得到最终的法向量。4.如权利要求1所述的一种基于点云特征分区的自适应点云精简方法，其特征在于：所述步骤3采用弦法向量计算点云主曲率k1、k2，计算过程如下：对于点云中的每个点p，设p点的为法向量为N，使用点坐标和法向量来估计点p处的法向曲率：假设p点附近有m个近邻点，q
i
为点p第i个近邻点，q
i
的法向量为M
i
，设正交坐标{undefinedp,X,Y,Z}为点p的局部坐标系L，N表示p点的法向量，X和Y为正交的单位向量；在L中，p，q
i
，M
i
的坐标是(0,0,0)，q
i
为(x
i,
y
i
,z
i
)，M
i
为(n
x
,i，n
y
...

【专利技术属性】
技术研发人员：田应仲，侯玉琴，李龙，
申请(专利权)人：上海大学，
类型：发明
国别省市：