本发明涉及计算机三维数据处理，尤其涉及一种基于噪声分类与mls的点云数据泊松曲面重建方法。
背景技术：
近年来，通过三维扫描或激光扫描设备获取的三维点云数据具有物体表面坐标与属性信息，处理时无需考虑拓扑关系，简化数据结构与算法复杂度。隐式曲面重建算法恰好利用该优势，准确表达被测物体表面信息的同时降低曲面重建算法复杂度。与其他隐式曲面重建算法相比，泊松曲面重建算法通过结合全局和局部拟合方法的优点，可有效光顺散乱点云中的噪声并修复一定的数据缺失，减少点云孔洞。但利用泊松方程重建曲面需准确的点云法向量信息，这对点云质量提出了更高要求，使得点云质量较低的情况下，不能采用传统的泊松方程重建曲面。
技术实现要素：
发明目的：本发明针对现有技术存在的问题，提供一种基于噪声分类与mls的点云数据泊松曲面重建方法，该方法对点云质量要求不高，且重建后的视觉更逼近实际模型，曲面重建质量更高。
技术方案：本发明所述的基于噪声分类与mls的点云数据泊松曲面重建方法包括：
(1)对于待重建模型点云数据，采用限幅滤波去除偏离主体点云的第1类噪声点；
(2)采用mls方法，对去噪后的点云数据集进行孔洞修复，并估计修复后的点云的法向量；
(3)采用双边滤波对法向量估计中对与目标采样数据点混合的第2类噪声进行光顺，得到优化点云；
(4)对优化点云采用泊松方程重建曲面。
进一步的，步骤(1)具体包括：
(1.1)对于待重建模型点云数据集合p＝{pi|i＝1,2,...}，计算出每个点pi到其所有邻近点的平均距离di；
(1.2)将平均距离在预设范围之外的点，标记为第1类噪点，采用限幅滤波从点云数据集中p去除。
进一步的，步骤(2)具体包括：
(2.1)采用mls方法，对去噪后的点云数据集采用局部二元多项式逼近处理，得到逼近曲面，实现孔洞修复；
(2.2)通过带权协方差矩阵鲁棒估计逼近曲面的局部参考平面，并根据局部参考平面估算得到点云的法向量。
进一步的，步骤(2.1)具体包括：
(2.1.1)设局部二元多项式为：
g(u,v)＝au2+2buv+cv2+du+ev+f
式中，u,v为函数自变量，a,b,c,d,e,f为待求系数，[u2,uv,v2,u,v,1]构成逼近曲面基函数；
(2.1.2)采用最小二乘误差最小计算待求系数，即：
采用以上公式得到系数a,b,c,d,e,f，其中，fi为点云数据中点pi到局部参考平面h1的垂直高度，fi＝n·(pi-q)，(ui,vi)是pi投影到h1后的局部坐标，θ(d)为权函数，且
(2.1.3)对于去噪后的点云数据集，采用部双变量逼近多项式g:h1→r3来逼近，得到逼近曲面，根据逼近曲面得到新采样点集p'＝{p′i|i＝1,2,...}；其中，r3表示多项式。
进一步的，步骤(2.2)具体包括：
(2.2.1)搜索逼近曲面上新采样点p′i的k个最近邻点，并计算带权协方差矩阵c，其中：
式中，k为最近邻点个数，p'ij表示采样点p'i的k个最近邻点的集合，
(2.2.2)计算权协方差矩阵c的3个特征向量{e0,e1,e2}，并以e0,e1,e2为坐标轴，
(2.2.3)以
式中，，x表示局部参考平面h2上的点；
(2.2.4)搜索点p′i在h2的投影点的k个最近邻点，并计算带权协方差矩阵c'，其中：
式中，qi、qij、
(2.2.5)计算c'最小特征值的特征向量，将该特征向量估算为qi的法向量；
(2.2.6)返回执行(2.2.1)直至得到每个点的法向量。
进一步的，步骤(3)具体包括：
(3.1)根据点云数据的法向量计算双边滤波算子α，其中：
式中，n(·)是对应数据点的邻域点集合，p′i是逼近曲面上新采样点，||||与<>分别表示向量的模或绝对值与向量的内积，ni、nj分别是数据点p'i的法向量、邻近点p'j的法向量；wc、ws分别是数据点p'i到其邻域点的光顺滤波权因子、p'i到邻域点在p'i法向ni上的特征保持权因子；
(3.2)将法向量估计中对与目标采样数据点混合噪声称为第2类噪声，采用双边滤波算子对点云数据中点p'i进行去噪处理，得到优化点云的数据点p″i，其中：
式中，
其中，步骤(3.1)中wc、ws计算通式为：
wc(x)＝exp[-x2/(2σc2)]
ws(x)＝exp[-x2/(2σs2)]
其中:σc为p'i到其邻域点的距离对p'i的影响因子，具体为点p'i的邻域半径，σs为数据点p'j到邻近点的距离在其法向上的投影对数据点p'j的影响因子，具体为数据点p'j的邻域点标准差。
有益效果：本发明与现有技术相比，其显著优点是：
(1)本文算法具有良好适应性，能更好处理大范围噪声和混合噪声。对曲面光顺效果理想，同时有效保持了曲面复杂几何特征且对边缘轮廓有一定增强。
(2)可修补多个点云孔洞，且修补效果良好。通过空间基函数与权函数产生合理的深度图像与透视关系，实现了曲面融合。
(3)与传统点云边界提取技术相比，法向信息精确度更高，通过更准确法向信息的辅助，大大提高曲面重建的质量。提高曲面重建的质量。同时，改进后的点云边界提取技术性能稳定，效率较高。
附图说明
图1是本发明的一个实施例的流程示意图。
具体实施方式
本实施例提供了一种基于噪声分类与mls的点云数据泊松曲面重建方法，如图1所示，包括以下步骤：
(1)对于待重建模型点云数据，采用限幅滤波去除偏离主体点云的第1类噪声点。
三维扫描或激光扫描的测量误差通常会产生离群稀疏点，称为第1类噪声点。其具体包括距离主体大片点云中心较远、小而密集的点云和偏离主体点云且悬浮在主体点云上方的稀疏点。由于限幅方法去除明显离群点的效果良好且算法简单，采用限幅滤波去除第1类噪声。具体过程如下。(1.1)对于待重建模型点云数据集合p＝{pi|i＝1,2,...}，计算出每个点pi到其所有邻近点的平均距离di；(1.2)将平均距离在预设范围之外的点，标记为第1类噪点，采用限幅滤波从点云数据集中p去除。
(2)采用移动最小二乘法(movingleastsquares，mls)方法，对去噪后的点云数据集进行孔洞修复，并估计修复后的点云的法向量。
由于原点云数据集存在部分点云数据缺失，易致模型表面发生起伏并造成曲面孔洞。故采用mls方法，对去噪后的点云数据集采用局部二元多项式逼近处理，得到逼近曲面，实现孔洞修复；再通过带权协方差矩阵鲁棒估计逼近曲面的局部参考平面，并根据局部参考平面估算得到点云的法向量。法向量精度越高，曲面重建效果越好。
a、孔洞修复过程为：
(2.1.1)设局部二元多项式为：
g(u,v)＝au2+2buv+cv2+du+ev+f
式中，u,v为函数自变量，a,b,c,d,e,f为待求系数，[u2,uv,v2,u,v,1]构成逼近曲面基函数；
(2.1.2)采用最小二乘误差最小计算待求系数，即：
采用以上公式得到系数a,b,c,d,e,f，其中，fi为点云数据中点pi到局部参考平面h1的垂直高度，fi＝n·(pi-q)，(ui,vi)是pi投影到h1后的局部坐标，θ(d)为权函数，且
(2.1.3)对于去噪后的点云数据集，采用部双变量逼近多项式g:h1→r3来逼近，得到逼近曲面，根据逼近曲面得到新采样点集p'＝{p′i|i＝1,2,...}；其中，r3表示多项式。
5、根据权利要求3所述的基于噪声分类与mls的点云数据泊松曲面重建方法，其特征在于：步骤(2.2)具体包括：
(2.2.1)搜索逼近曲面上新采样点p′i的k个最近邻点，并计算带权协方差矩阵c，其中：
式中，k为最近邻点个数，p'ij表示采样点p'i的k个最近邻点的集合，
(2.2.2)计算权协方差矩阵c的3个特征向量{e0,e1,e2}，并以e0,e1,e2为坐标轴，
(2.2.3)以
式中，，x表示局部参考平面h2上的点；
(2.2.4)搜索点p′i在h2的投影点的k个最近邻点，并计算带权协方差矩阵c'，其中：
式中，qi、qij、
(2.2.5)计算c'最小特征值的特征向量，将该特征向量估算为qi的法向量；
(2.2.6)返回执行(2.2.1)直至得到每个点的法向量。
(3)采用双边滤波对法向量估计中对与目标采样数据点混合的第2类噪声进行光顺，得到优化点云。
其中，双边滤波在平滑滤波的同时能大量保留图像的边缘和细节特征，达到保边去噪的目的。具有简单、非迭代、局部的特点。故本实施例确定点云法向量后，采用双边滤波去除第2类点云噪声。第2类点云噪声包含与目标采样数据点混合的噪声。去噪具体步骤为：
(3.1)根据点云数据的法向量计算双边滤波算子α，其中：
式中，n(·)是对应数据点的邻域点集合，p′i是逼近曲面上新采样点，||||与<>分别表示向量的模或绝对值与向量的内积，ni、nj分别是数据点p'i的法向量、邻近点p'j的法向量；wc、ws分别是数据点p'i到其邻域点的光顺滤波权因子、p'i到邻域点在p'i法向ni上的特征保持权因子；wc、ws计算通式为：
wc(x)＝exp[-x2/(2σc2)]
ws(x)＝exp[-x2/(2σs2)]
其中:σc和σs为高斯滤波参数，反映计算任意采样点的双边滤波因子时切向和法向的影响范围。σc为p'i到其邻域点的距离对p'i的影响因子，它用来控制光顺程度，σc越大，说明选取的邻域点越多，则三维点云数据模型光滑性越好σs为数据点p'j到邻近点的距离在其法向上的投影对数据点p'j的影响因子，它用来控制特征保持程度，σs越大，表明改进双边滤波因子在平滑数据点p'j时在其法矢上移动的距离越长，则三维点云数据模型的特征保持性越好。通常情况下，σc可用数据点p'i的邻域半径来表示，σs可用数据点p'j的邻域点标准差来表示。当σc确定后，σs的不同取值对于三维点云数据模型平滑效果影响较明显，得到滤波后的采样数据，即为优化点云。
(3.2)采用双边滤波算子对点云数据中点p'i进行去噪处理，得到优化点云的数据点p″i，其中：
(4)对优化点云采用泊松方程重建曲面。
重建曲面方法为现有技术，不再赘述。
以上所揭露的仅为本发明一种较佳实施例而已，不能以此来限定本发明之权利范围，因此依本发明权利要求所作的等同变化，仍属本发明所涵盖的范围。