法律状态公告日
法律状态信息
法律状态
2019-04-26
授权
授权
2017-04-19
实质审查的生效 IPC(主分类):G06T5/00 申请日:20160921
实质审查的生效
2017-03-22
公开
公开
技术领域
本发明涉及计算机视觉技术领域,具体涉及一种散乱点云的保特征滤波方法。
背景技术
随着计算机视觉技术的快速发展,点云数据的研究得到了越来越广泛的关注。然而,尽管三维扫描设备的精度不断提高,但是由于设备自身的不足,人为误差及各种外部环境等因素造成扫描得到的数据中不可避免的含有噪声。因此,一种有效的点云去噪滤波算法对于后续的模型重建、优化等过程具有重要作用。
传统的采样模型去噪方法大致分为网格去噪算法和散乱点云去噪算法,然而由于扫描获取的数据不包含采样点间的连接结构,因此,散乱点云数据去噪算法应用更方便、广泛。传统的滤波方法只是基于采样点与邻域点间距离及其在法向上的投影对采样数据滤波,且没有对邻域点进行有效选取,不能在滤波的同时有效保持模型的细节。
发明内容
本发明目的在于提供一种滤波精度高、去燥效果好的散乱点云的保特征滤波方法。
为实现上述目的,采用了以下技术方案:本发明方法是获取采样点的k个邻域点中的有效邻域,根据采样点与其邻域点的曲率差异,设定阈值T,当曲率差异小于T时,该点为有效邻域点,否则为无效邻域点;该采样点所有有效邻域点构成有效邻域;采样点、邻域点及邻域点在采样点切平面的投影点构成的三角形的面积,与法向特征差异相结合作为滤波因子对散乱点云数据去噪。
所述方法的具体步骤如下:
步骤1,输入噪声点云集V={vi,i=1,2,…n},通过kdtree算法对邻域进行搜索并设定邻域点个数k,定义矩阵
其中,c为邻域中心;vij为采样点vi的第j个邻域点;对矩阵R进行特征分解,0≤λ1≤λ2≤λ3为矩阵R的三个非负特征值,e1,e2,e3为相应的特征向量,最小特征值λ1所对应的特征向量e1为该点的法向Ni,最大特征值λ3所对应的特征向量e3近似为该采样点的切向Ti,采样点的曲率
步骤2,根据采样点曲率ρi与邻域点曲率ρij的差异判断有效邻域,计算ρi-ρij,并设定阈值H,当ρi-ρij的值小于H时,该邻域点为有效邻域点,否则为无效邻域点,即定义二值函数η表示该点是否为有效邻域点;
步骤3,计算采样点vi与邻域点vij之间的位置差异||vij-vi||,qj为邻域点vij在采样点vi的切平面的投影点,计算邻域点vij与投影点qj间的距离<Ni,vij-vi>,采样点vi与投影点qj间的距离<Ti,vij-vi>,采样点vi、邻域点vij以及投影点qj构成的三角形的面积
步骤4,根据下式计算位置权重函数Wc和特征权重函数Ws
其中,σc为位置影响因子,σs为特征影响因子,Ni为采样点vi的法向,Nij为邻域点vij的法向;
步骤5,将Wc和Ws代入双边滤波因子h
其中,N(v)为采样点的邻域点集;
步骤6,将h代入下式中
vout=vi+Ni·h(6)
其中,vout为输出采样点;Ni为vi的法向;vi为原始待处理采样点;h为滤波因子;
步骤7,对于所有采样点重复以上步骤,得到去噪后的点云模型。
与现有技术相比,本发明具有如下优点:本发明对采样点的k个邻域点进行有效选取,避免了差异大的邻域点或不在同一面上的邻域点对采样点的影响。在有效邻域下,结合面积与法向特征差异作为滤波因子对点云数据去噪,使得,对于平面上的点,权重主要与距离呈反比,对于棱上的点,同一棱上的邻域点的权重更大,对于角点,各邻域点的权重都很小,因此,该方法能够兼顾平坦区域的光顺与高频区域细节特征的保持。
附图说明
图1是本发明方法的流程图。
图2为构造滤波因子过程示意图。
图3(a)为加噪声的fandisk模型图。
图3(b)为双边滤波算法去噪效果图。
图3(c)为本发明方法的去噪效果图。
具体实施方式
下面结合附图对本发明做进一步说明:
本发明方法是获取采样点的k个邻域点中的有效邻域,根据采样点与其邻域点的曲率差异,设定阈值H,当曲率差异小于H时,该点为有效邻域点,否则为无效邻域点;该采样点所有有效邻域点构成有效邻域;采样点、邻域点及邻域点在采样点切平面的投影点构成的三角形的面积,与法向特征差异相结合作为滤波因子对散乱点云数据去噪。
实施例一:
如图1所示,本发明方法的具体步骤如下:
步骤1,图3(a)为输入fandisk模型噪声点云集V={vi,i=1,2,…n},通过kdtree算法对邻域进行搜索并设定邻域点个数k,定义矩阵
其中,c为邻域中心;vij为采样点vi的第j个邻域点;对矩阵R进行特征分解,0≤λ1≤λ2≤λ3为矩阵R的三个非负特征值,e1,e2,e3为相应的特征向量,最小特征值λ1所对应的特征向量e1为该点的法向Ni,最大特征值λ3所对应的特征向量e3近似为该采样点的切向Ti,采样点的曲率
步骤2,在搜索到点vi的k个邻域中,根据采样点与邻域点间的曲率差异判断有效邻域,即计算ρi-ρij,其值越小,采样点与邻域点间的特征差异越小,设定阈值H,当ρi-ρji的值小于H时,该邻域点为有效邻域点,否则为无效邻域点,定义二值函数η表示该点是否为有效邻域点
由此,η=1的邻域点构成了点pi的有效邻域。
步骤3,计算采样点与邻域点之间的位置差异||vij-vi||,如图2,设qj为邻域点vij在采样点vi的切平面的投影点,计算邻域点与投影点间的距离<Ni,vij-vi>,也就是邻域点到采样点切平面的距离;采样点与投影点间的距离<Ti,vij-vi>,也就是采样点vi到过邻域点vij且以采样点的切向Ti为法向的平面的距离;采样点、邻域点以及投影点构成的三角形的面积
步骤4,根据下式计算位置权重函数Wc和特征权重函数Ws
其中,Ni为采样点vi的法向,Nij为邻域点vij的法向;σc为位置影响因子,它的值越大,平滑效果越好,但特征可能被过度平滑,σs为特征影响因子,其值越大,特征保持的越好,但可能使平坦区域不光滑。
步骤5,将Wc和Ws代入双边滤波因子h
其中,N(v)为采样点的邻域点集;
步骤6,将h代入下式中
vout=vi+Ni·h(6)
其中,vout为输出采样点;Ni为vi的法向;vi为原始待处理采样点;h为滤波因子。
步骤7,对于所有采样点重复以上步骤,得到去噪后的点云模型。
对fandisk加噪模型进行仿真实验,如图3,图3(a)为含噪声的fandisk模型图,图3(b)为双边滤波法去噪后的效果图,图3(c)为本发明去噪后的效果图,对比图3(b)和3(c)可以发现,在长方形内区域,图3(c)的棱线更为明显,实验表明,在同等条件下,本发明比双边滤波算法能够更好地保持模型的细节特征。
以上所述的实施例仅仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。
机译: 用于传输点云数据的设备,一种用于传输点云数据的方法,一种用于接收点云数据的设备和/或一种用于接收点云数据的方法
机译: 用于制备适合消费的饮料的保形垫,包括制垫机和这样的保形垫的饮料制备机的组装,一种用于制造保形垫的预成型垫的方法一种方法,一种制备饮料的方法,该饮料适合使用这样的保形垫或组合件以及包括多个保形垫的零件套件
机译: 用于制备适合消费的饮料的保形垫,包括制垫机和这样的保形垫的饮料制备机的组装,一种用于制造保形垫的预成型垫的方法一种方法,一种制备饮料的方法,该饮料适合使用这种保形垫或这样的组件以及包括多个这种保形垫的零件套件。