高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于主成分分析与曲面拟合的激光点云滤波去噪

柳赟 孙淑艳

引用本文:
Citation:

基于主成分分析与曲面拟合的激光点云滤波去噪

    作者简介: 柳赟(1986-), 女, 硕士, 工程师, 研究方向为电子信息工程。E-mail:kevin1320622@126.com.
  • 基金项目:

    中央高校基本科研业务费专项资金资助项目 2019MS005

  • 中图分类号: TN249

Laser point cloud denoising based on principal component analysis and surface fitting

  • CLC number: TN249

  • 摘要: 为了消除激光点云采集时点云中的噪声点,避免噪声尤其是一些孤立离群点对点云数据质量的影响,将散乱的、含有噪声点云变成规则的、高精度的点云,采用了基于主成分分析与曲面拟合进行点云去噪的方法,首先提出了点云区域的主成分分析计算方法,在主成分分析的法向量进行粗去噪,而后去噪后的点云进行曲面拟合,最后根据点到曲面的距离进行了点云的滤波,得到滤波后的点云。结果表明,该方法去噪效果精度高,尤其针对散乱点云,去噪效果明显,最佳滤波性能误差仅为0.018mm。该研究为散乱激光点云的去噪滤波提供了参考。
  • Figure 1.  Point cloud chair

    a—before filter b—after filter

    Figure 2.  Point cloud gargo

    a—before filter b—after filter

    Table 1.  Performance comparison of several filter method

    filter method average error max error
    chair gargo chair gargo
    method in reference [24] 0.026 0.038 0.033 0.042
    method in reference [18] 0.023 0.037 0.031 0.040
    method in reference [28] 0.024 0.035 0.031 0.041
    proposed method 0.018 0.028 0.023 0.032
    下载: 导出CSV

    Table 2.  Effect of surface fitting window T to filter

    point cloud average error max error
    T=16 T=22 T=28 T=31 T=37 T=35 T=16 T=22 T=28 T=31 T=37 T=35
    chair 0.032 0.029 0.019 0.026 0.031 0.031 0.031 0.027 0.025 0.028 0.032 0.035
    gargo 0.041 0.032 0.027 0.031 0.037 0.042 0.039 0.039 0.037 0.034 0.041 0.046
    下载: 导出CSV

    Table 3.  Effect of surface window M to filter

    point cloud average error max error
    M=10 M=15 M=20 M=25 M=30 M=35 M=10 M=15 M=20 M=25 M=30 M=35
    chair 0.028 0.026 0.018 0.023 0.027 0.031 0.029 0.025 0.024 0.025 0.028 0.031
    gargo 0.033 0.031 0.028 0.030 0.034 0.037 0.039 0.035 0.032 0.036 0.039 0.042
    下载: 导出CSV
  • [1]

    TAUBIN G. A signal processing approach to fair surface design[C]//ACM Proceedings of the 22nd Annual Conference on Computer Graphics and Interactive Techniques. Palo Alto, USA: Stanford University Press, 1995: 351-358. 
    [2]

    TAUBIN G. Linear anisotropic mesh filters[C]//IBM Research Report RC22213(W0110-051). New York, USA: Computer Science, 2001: 110-115.
    [3]

    JONES T R, DURAND F, DESBRUN M. Non-iterative, feature-preserving mesh smoothing[C]//ACM Transactions on Graphics. New York, USA: ACM, 2003, 22(3): 943-949. 
    [4]

    HILDEBRANDT K, POLTHIER K. Anisotropic filtering of non-linear surface features [J]. Computer Graphics Forum, 2004, 23(3):391-400. doi: 10.1111/j.1467-8659.2004.00770.x
    [5]

    ZHANG W, DENG B, ZHANG J, et al. Guided mesh normal filtering[J]. Computer Graphics Forum, 2015, 34(7): 23-34. doi: 10.1111/cgf.12742
    [6]

    YAGOU H, OHTAKE Y, BELYAEV A. Mesh smoothing via mean and median filtering applied to face normals[C]//IEEE Geometric Modeling and Processing. New York, USA: IEEE, 2002: 124-131. 
    [7]

    FAN H, YU Y, PENG Q. Robust feature-preserving mesh denoising based on consistent subneighborhoods[J]. IEEE Transactions on Visualization and Computer Graphics, 2010, 16(2): 312-324. doi: 10.1109/TVCG.2009.70
    [8]

    FLEISHMAN S, DRORI I, COHEN-OR D. Bilateral mesh denoising[J]. ACM Transactions on Graphics, 2003, 22(3): 950-953. doi: 10.1145/882262.882368
    [9]

    CHOUDHURY P, TUMBLIN J. The trilateral filter for high contrast images and meshes[C]//Eurographics Workshop on Rendering Techniques. New York, USA: ACM, 2003: 186-196.
    [10]

    DUGUET F, DURAND F, DRETTAKIS G.Robust higher-order filtering of points[J]. Rapport de Recherche Inria, 2004, 17(4):12-18. 
    [11]

    TASDIZEN T, WHITAKER R, BURCHARD P, et al. Geometric surface smoothing via anisotropic diffusion of normals[C]// IEEE Proceedings of the Conference on Visualization. New York, USA: Computer Society, 2002: 125-132. 
    [12]

    MEDEROS B, VELHO L, de FIGUEIREDO L H. Robust smoothing of noisy point clouds[C]//Slam Conference on Geometric Design & Computing. New York, USA: IEEE, 2003: 405-416. 
    [13]

    CLARENZ U, RUMPF M, TELEA A. Fairing of point based surfaces[C]//IEEE Computer Graphics International. New York, USA: IEEE, 2004: 600-603. 
    [14]

    PAULY M, KOBBELT L, GROSS M. Multiresolution modeling of point-sampled geometry[M].Zurich, Swiss: Swiss Federal Institute of Technology Zurich Press, 2002:42-49.
    [15]

    WANG R, CHEN W, ZHANG S, et al. Similarity-based denoising of point-sampled surfaces[J]. Journal of Zhejiang University, 2008, 9(6): 807-815. doi: 10.1631/jzus.A071465
    [16]

    FAN H Q, MIAO L F, ZHANG X, et al. Robust high-order bilateral denoising filter[J]. Journal of Computer-aided Design & Computer Graphics, 2009, 21(6):813-818 (in Chinese). 
    [17]

    FU Y, ZHAI J L. Research on scattered points cloud denoising algorithm [C]// IEEE Intemationai Conference on Signal Communications and Computing. New York, USA: IEEE, 2015: 1-5. 
    [18]

    SONG D H, LI Zh K, WANG Zh, et al. A feature-preserving algorithm of point cloud smoothing[J]. China Mechanical Engineering, 2014, 25(11):1498-1501 (in Chinese). 
    [19]

    NARVÁEZ E A L, NARVÁEZ N E L. Point cloud denoising using robust principal component analysis[C]//Proceedings of the First International Conference on Computer Graphics Theory and Applications. New York, USA: IEEE, 2006: 51-58. 
    [20]

    LANGE C, POLTHIER K. Anisotropic smoothing of point sets[J]. Computer Aided Geometric Design, 2005, 22(7): 680-692. doi: 10.1016/j.cagd.2005.06.010
    [21]

    HU G, PENG Q, FORREST A R. Mean shift denoising of point-sampled surfaces[J]. The Visual Computer, 2006, 22(3):147-157. doi: 10.1007/s00371-006-0372-0
    [22]

    LEVIN D. Mesh-independent surface interpolation[M]. Berlin, Germany: Springer Press, 2004: 37-49.
    [23]

    GUENNEBAUD G, GROSS M. Algebraic point set surfaces[J]. ACM Transactions on Graphics, 2007, 26(3): 231-236. 
    [24]

    WANG J, YU Z. Feature-preserving mesh denoising via anisotropic surface fitting[J]. Journal of Computer Science and Technology, 2012, 27(1): 163-173. doi: 10.1007/s11390-012-1214-3
    [25]

    ZHENG Y, FU H, AU O K C, et al. Bilateral normal filtering for mesh denoising[J]. IEEE Transactions on Visualization and Computer Graphics, 2011, 17(10): 1521-1530. doi: 10.1109/TVCG.2010.264
    [26]

    SHI L, GUO T, PENG Ch, et al. Segmentation of laser point cloud and safety detection of power lines [J].Laser Technology, 2019, 43(3): 341-3346 (in Chinese). 
  • [1] 王云云唐菲菲王章朋肖敏唐天俊王铜川 . 植被茂密地区点云双重滤波方法研究. 激光技术, 2022, 46(2): 233-238. doi: 10.7510/jgjs.issn.1001-3806.2022.02.014
    [2] 李靖钱建国王伟玺李晓明李游 . 基于机载激光点云的电力线自动提取方法. 激光技术, 2021, 45(3): 362-366. doi: 10.7510/jgjs.issn.1001-3806.2021.03.017
    [3] 唐菲菲阮志敏刘星 . 城区机载LIDAR数据滤波方法研究. 激光技术, 2011, 35(4): 527-530. doi: 10.3969/j.issn.1001-3806.2011.04.021
    [4] 闫德立高尚李韶华霍萌 . 基于激光雷达的道路不平度及可行驶区域检测. 激光技术, 2022, 46(5): 624-629. doi: 10.7510/jgjs.issn.1001-3806.2022.05.007
    [5] 吴兆勇杜正春 . 基于误差椭球的激光测量系统的不确定度分析. 激光技术, 2017, 41(1): 29-33. doi: 10.7510/jgjs.issn.1001-3806.2017.01.007
    [6] 屈文俊谢卉胡强高孙莉萍 . 高斯光束经衍射光栅的滤波特性研究. 激光技术, 2012, 36(4): 540-543. doi: 10.3969/j.issn.1001-806.2012.04.027
    [7] 苍桂华岳建平 . 基于加权总体最小二乘法的点云平面拟合. 激光技术, 2014, 38(3): 307-310. doi: 10.7510/jgjs.issn.1001-3806.2014.03.005
    [8] 汪杰君杨杰张文涛王新强叶松 . 雾天偏振成像影响分析及复原方法研究. 激光技术, 2016, 40(4): 521-525. doi: 10.7510/jgjs.issn.1001-3806.2016.04.014
    [9] 董军彭翰生魏晓峰胡东霞周维赵军普程文雍刘兰琴 . 线性啁啾脉冲频谱干涉特性的模拟研究. 激光技术, 2009, 33(3): 232-235.
    [10] 刘兆亮杜龙龙王荣扬赵文超 . 基于激光点云技术的机床夹具设计与改进. 激光技术, 2024, 48(4): 597-602. doi: 10.7510/jgjs.issn.1001-3806.2024.04.020
    [11] 史洪云虢韬王迪王时春赵健刘欣龙新 . 基于激光点云的电力线悬挂点定位方法. 激光技术, 2020, 44(3): 364-370. doi: 10.7510/jgjs.issn.1001-3806.2020.03.017
    [12] 王德旺王改利 . 自适应中值滤波在云雷达数据预处理的应用. 激光技术, 2012, 36(2): 217-220,224. doi: 10.3969/j.issn.1001-3806.2012.02.019
    [13] 徐梁刚时磊陈凤翔王时春龙新王迪 . 基于激光点云的输电线路杆塔倾斜检测算法. 激光技术, 2022, 46(3): 390-396. doi: 10.7510/jgjs.issn.1001-3806.2022.03.015
    [14] 林承达谢良毅韩晶胡方正 . 基于激光点云的农田玉米种植株数数目识别. 激光技术, 2022, 46(2): 220-225. doi: 10.7510/jgjs.issn.1001-3806.2022.02.012
    [15] 曾旭陈伯建潘磊李诚龙江波 . 基于机载激光点云的电网绝缘子识别方法. 激光技术, 2023, 47(1): 80-86. doi: 10.7510/jgjs.issn.1001-3806.2023.01.012
    [16] 吴华刘海燕丁高峰曹飞 . 复杂环境中电力线激光点云的自动提取. 激光技术, 2020, 44(4): 509-514. doi: 10.7510/jgjs.issn.1001-3806.2020.04.019
    [17] 时磊虢韬彭赤石书山杨立任曦胡伟 . 电力线激光点云的分割及安全检测研究. 激光技术, 2019, 43(3): 341-346. doi: 10.7510/jgjs.issn.1001-3806.2019.03.010
    [18] 柳长安孙书明赵丽娟 . 基于激光点云实现杆塔提取的轻量级网络. 激光技术, 2021, 45(3): 367-372. doi: 10.7510/jgjs.issn.1001-3806.2021.03.018
    [19] 吴俊河林松施向丰 . 层次化点云边界快速精确提取方法研究. 激光技术, 2021, 45(5): 571-575. doi: 10.7510/jgjs.issn.1001-3806.2021.05.006
    [20] 王果王成王宏涛张成龙杨福芹 . 利用SLAM点云的玉米株数自动识别. 激光技术, 2024, 48(1): 140-144. doi: 10.7510/jgjs.issn.1001-3806.2024.01.022
  • 加载中
图(2) / 表(3)
计量
  • 文章访问数:  10493
  • HTML全文浏览量:  7036
  • PDF下载量:  41
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-10-18
  • 录用日期:  2019-12-12
  • 刊出日期:  2020-07-25

基于主成分分析与曲面拟合的激光点云滤波去噪

    作者简介: 柳赟(1986-), 女, 硕士, 工程师, 研究方向为电子信息工程。E-mail:kevin1320622@126.com
  • 华北电力大学 电气与电子工程学院,北京 102206
基金项目:  中央高校基本科研业务费专项资金资助项目 2019MS005

摘要: 为了消除激光点云采集时点云中的噪声点,避免噪声尤其是一些孤立离群点对点云数据质量的影响,将散乱的、含有噪声点云变成规则的、高精度的点云,采用了基于主成分分析与曲面拟合进行点云去噪的方法,首先提出了点云区域的主成分分析计算方法,在主成分分析的法向量进行粗去噪,而后去噪后的点云进行曲面拟合,最后根据点到曲面的距离进行了点云的滤波,得到滤波后的点云。结果表明,该方法去噪效果精度高,尤其针对散乱点云,去噪效果明显,最佳滤波性能误差仅为0.018mm。该研究为散乱激光点云的去噪滤波提供了参考。

English Abstract

    • 在进行激光点云采集时,受仪器本身影响及现场环境干扰,点云中会混有噪声点。点云去噪是激光雷达数据处理的基础。激光点云数据分为两种,有序点云及无序点云。有序点云是指激光点距离整体点云偏差小。无序点云是指比较孤立的、离群较远的噪声点。针对不同类型的点云数据,研究人员采用不同的去噪滤波方法。

      针对有序点云数据,通常采用的去噪算法包括中值滤波、高斯滤波、拉普拉斯滤波、法线方向平滑、双边滤波、曲面拟合去噪等。最早对有序点云进行去噪处理的是中值滤波和高斯滤波,中值滤波的中心思想是将需要处理的采样点按照位置及邻域关系进行重新排序,将周围领域点计算得到的均值点作为当前点的位置,该算法简单快捷,但其缺点是会产生新的噪声点,并且整体细节特征会比较模糊,这会使得目标失真,影响到点云的后期处理。高斯滤波是对邻域点进行加权处理、求取平均值,加权是按照噪声点与邻域点的距离,根据距离分配权重,距离近,权值大,距离远,权值相对较小,这种方法能够在一定程度上保留目标的细节特征。TAUBIN将拉普拉斯算法用于点云去噪[1],能够达到平滑的效果,但是,由于拉普拉斯算子在各个方向的扩散速度相同,该算法在细节特征保留方面较弱。TAUBIN根据拉普拉斯算法对点云去噪,同时采用最小二乘算法去噪,可以减小漂移以及外形收缩问题[2]。在此之后,JONES提出了基于特征保持来进行滤波的平滑算法[3],通过邻域点进行计算,采用加权处理来进行采样点中噪声的滤波,但由于法线计算误差,滤波出现了顶点漂移。在通过法向量优化滤波中,HILDEBRANDT将离散的各向异性平均曲率流用于网格曲面模型的去噪中,精度随噪声大小而变化,当噪声较大时,该算法的稳定性降低[4]。ZHANG通过1阶邻域法向量确定导向,然后进行滤波[5]。YAGOU等人首先对法向采用均值滤波,然后对点云沿着法线方向进行滤波,可以取得较高精度,但当噪声较大时,点云易发生变形[6]。针对这种现象,FAN等人对法向进行分类,并对法线方向进行平滑,虽然改善了漂移问题,但导致了复杂性过高[7]。与此同时,FLEISHMAN等人[8]先对法线方向进行平滑,然后采用双边滤波方法,对点云进行处理,经过迭代得到最终的目标点云。在此基础上,CHOUDHURY等人提出了三边滤波方法,在双边滤波理论基础上,根据法线向量确定滤波核函数的倾斜度,提高了细节特征的精度[9],该方法的复杂度也较高。为了进一步提高去噪精度,DUGUET提出局部二次曲面拟合[10],二次曲面拟合能够改善点云体积收缩问题,但其带来的顶点漂移易导致点云发生扭曲变形。TASDIZEN等人采用了水平曲面模型,提高了点云去噪的精度[11]

      针对无序点云数据,很多学者都提出了处理方法,也有很多研究成果,主要包括传统滤波算子、双边滤波及曲面滤波等。MEDEROS等人采用平滑算子进行点云滤波,点云的部分细节特征能够保留[12]。CLARENZ将基于偏微分方程用于散乱点云去噪中,但在计算过程中,由于需要多次迭代,复杂度高[13]。PAULY采用先加权后拉普拉斯变化的方法,反复迭代实现点云去噪[14]。WANG等人分析了点云几何密度与特征的关系,将二者结合对点云进行滤波去噪[15]。在点云去噪滤波中,双边滤波也具有较好的性能,双边滤波通过采用邻域信息计算当前点的滤波值,其特点是既具有较强的平滑性能,也能保留较高的细节特征。但是,针对较大噪声的点云,双边滤波处理的能力降低。针对具有复杂细节特征的目标,双边滤波也很难有较好的处理效果。除此之外,双边滤波的复杂度较高。FAN等人将双边滤波算法和鲁棒理论相结合,对点云数据进行滤波去噪,保证了点云外形的特征[16]。基于双边滤波,这些改进的去噪滤波方法能够不同程度上对点云的噪声进行处理,有的能够防止高频细节特征被过度的平滑,有的能够防止点云外形的收缩扭曲变形[17]。SONG等人综合考虑距离、法向及曲率等特征的相对强弱关系,计算出滤波参量,进而进行去噪滤波,效果较好[18]。NARVÁES等人通过分配权重因子,采用线性局部拟合来减小点云的收缩[19]。LANGE等人采用了各向异性平均曲率流进行去噪滤波,算法的精度受到曲率变化的影响,特别针对复杂的场景,去噪效果不明显[20]。HU等人[21]首先采用均值漂移算法对采样点分类,而后采用三边滤波算法进行点云滤波,精度较高,但复杂度也稍高。曲面拟合滤波也是一种思路,LEVIN将点云拟合成光滑的曲面,而后通过采样点到曲面的投影来去噪[22]。GUENNEBAUD等人给予曲面投影技术,采用代数球面拟合进行去噪[23]。与双边滤波相比,采用各向异性的网格曲面模型进行滤波去噪对处理较低噪声场景,效果较好,但噪声变大时,去噪效果降低,容易产生失真、变形、伪特征[6]。针对这种情况,WANG等人提出双边滤波和曲面拟合相结合的方法,首先根据每个点的法向加权距离进行邻域的调整,而后采用双边滤波方法对法向进行去噪,最后,去除噪声后的法向信息与曲面信息进行拟合来进行去噪[24]。除此之外,ZHENG等人提出了另外一种思路,首先对计算出来的法向量去噪处理,然后提出了法向场的去噪方法,该方法的提出使得在强噪声的点云中,仍然能够有效地完成去噪,对于复杂的点云曲面,该方法仍具有较好地去噪性能[25]。SHI等人提出了一种基于密度的空间聚类方法的电力线激光点云单条分割提取算法[26],用于得到单条电力线激光点云数据。

      从上述分析可以看到,不同的去噪滤波方法能够处理不同类型的点云噪声问题,但是尚未有一种滤波去噪方法,能够兼而有之的对有序点云数据和无序点云数据进行处理。作者提出了基于主成分分析与曲面拟合进行点云去噪,该方法适用于有序及无序点云噪声,能够较高精度的恢复目标原始细节。

    • 点云数据的目的是将噪声点云进行处理,去除噪声干扰,使其尽最大可能回归到应有的位置,将散乱的、含有噪声的点云变成规则的、高精度的点云,从而使得目标中的平面区域更加平顺,细节部分更加精细,在此同时,尽最大可能保留目标应该有的细节信息,避免在数据去噪过程中产生细节消失、扭曲、伸缩等现象,得到目标的真实外形点云, 为后续的点云处理做好准备。为了实现这个目的,采取了基于主成分分析的去噪、曲面拟合、滤波3个步骤。其中主成分分析的目的是找到曲面的前两个主成分,从而方便在第3个主成分进行去噪滤波;曲面拟合将在前两个主成分的基础上进行拟合,求出每个曲面;最后将基于拟合出来的曲面和第3个主成分方向进行去噪处理。具体描述如下。

    • 对点云,以一定半径窗口进行滑动,每个窗口内,求取点云的xyz坐标均值(ox, oy, oz)。求取协方差矩阵:

      $ \left[ \begin{array}{l} \ \ \ \ \ \ \ \sum {{{({x_i} - {o_x})}^2}\ \ \ \ \ \ \ \ \ \ \sum {({x_i} - {o_x})({y_i} - {o_y})} \ \ \ \ \ \sum {({x_i} - {o_x})({z_i} - {o_z})} } \\ \sum {({y_i} - {o_y})({x_i} - {o_x})}\ \ \ \ \ \ \ \ \ \ \ \sum {{{({y_i} - {o_y})}^2}}\ \ \ \ \ \ \ \ \ \ \ \ \ \sum {({y_i} - {o_y})({z_i} - {o_z})} \\ \sum {({z_i} - {o_z})({x_i} - {o_x})}\ \ \ \ \sum {({z_i} - {o_z})({y_i} - {o_y})} \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \sum {{{({z_i} - {o_z})}^2}} \end{array} \right] $

      (1)

      式中, (xi, yi, zi)为窗口内每个点的坐标。求取矩阵的两个主方向($\overrightarrow x ,\overrightarrow y $),既为窗口点云的主方向,两个主方向差乘,得到的方向为点云噪声方向$\overrightarrow z $。在两个主方向上对进行局部平面拟合,然后求出每个点(xi, yi, zi)与平面沿噪声方向$\overrightarrow z $的距离${d_{\overrightarrow z }}$。

      每个点在噪声方向进行位移,位移距离d为:

      $ d = {d_{\overrightarrow z }} \times \alpha $

      (2)

      式中, α为系数。

      该步骤可以去除较大的噪声,尤其离群点噪声。

    • 去除较大噪声后,对点云进一步进行曲面拟合的方程为:

      $ \begin{array}{c} f(x,y) = \sum\limits_{i = 0}^n {\sum\limits_{j = 0}^n {{a_{ij}}{x^i}{y^i}} } \\ (i,j = 0,1,2, \cdots ,n) \end{array} $

      (3)

      式中, aij为系数,为了求得最优aij,定义误差函数为:

      $ E({a_{ij}}) = \sum\limits_{i = 1}^n {{{[{z_i} - f({x_i},{y_i})]}^2}} $

      (4)

      将所有点代入,转为多项式求解,通过整体误差最小,求得系数:

      $ \frac{{\partial E({a_{ij}})}}{{\partial {a_{ij}}}} = 0 $

      (5)

      整幅点云范围较大,为了方便处理,在计算过程中,采取窗口扫面方式对点云进行拟合处理。针对每个窗口点云,采用(4)式和(5)式进行处理。当窗口遍历完所有点云后,得到整体点云曲面。

      曲面大小的计算方法为距离待处理像素(x0, y0, z0)的欧式距离在一定范围T内的点云,即:

      $ \begin{array}{c} \{ {x_i},{y_i},{z_i}\} \in \\ \sqrt {{{({x_i} - {x_0})}^2} + {{({y_i} - {y_0})}^2} + {{({z_i} - {z_0})}^2}} \end{array} $

      (6)

      式中, {xi, yi, zi}为在(x0, y0, z0)周围进行曲面拟合的所有点云。

    • 滤波算法为:

      $ f(x,y,z) = \frac{{\sum {{\sigma _{x,y,z}}} \times (x,y,z) \times \widetilde d}}{{\sum {{\sigma _{x,y,z}}} }} $

      (7)

      式中, $\widetilde d$为该点到曲面的距离,x, y, z为坐标值。

      在(7)式中,σx, y, z为窗口内标准差,计算方法如下:

      $ \begin{array}{c} {\sigma _{x,y,z}} = \frac{1}{{(2M + 1)(2M + 1)}} \times \\ \sum {[(} x,y,z) - ({o_x},{o_y},{o_z}){]^2} \end{array} $

      (8)

      式中, 窗口大小为(2M+1)(2M+1)。随不同区域起伏不同,M值也随着调整。M的计算过程为:

      $ M = k \times f $

      (9)

      式中,k为给定的系数,f为范围T内的点云的起伏情况,f的计算方法为:

      $ f = {\rm mean}({\widetilde d_i}) $

      (10)

      即每个点到曲面距离绝对值的均值,其中${\widetilde d_i}$为每个点到曲面距离绝对值。

    • 实验结果分别如图 1图 2所示。可以看出,经本文中提出的方法滤波后,消除了较大的噪声,有效地保留模型中的尖锐、边缘特征及细节特征。从图 1中可以看到,椅子边缘的噪声被滤波去掉,保留的椅子边缘细腻。从图 2中可以看到,经过滤波后,gargo的细节能够比较清晰地展现出来。

      Figure 1.  Point cloud chair

      Figure 2.  Point cloud gargo

      表 1中对4种去噪算法进行了滤波效果比较,从中可以发现,与3种滤波算法相比,提出的算法在进行滤波处理后,点云平均误差更小,通过统计发现,整体上与原始模型更为接近;同时,提出的滤波算法,相对其它3种算法比较,最大误差较更小,这表明所提出的算法能够更好地将目标的高频特征或细节特征保留。

      Table 1.  Performance comparison of several filter method

      filter method average error max error
      chair gargo chair gargo
      method in reference [24] 0.026 0.038 0.033 0.042
      method in reference [18] 0.023 0.037 0.031 0.040
      method in reference [28] 0.024 0.035 0.031 0.041
      proposed method 0.018 0.028 0.023 0.032

      表 2中给出了曲面拟合窗口的大小T对整体滤波精度的影响。从表中可以看到,T对去噪效果有一定影响。T过大及过低,对滤波效果都有坏处,合适的滤波窗口,将带来最佳的滤波效果。在T=28时,能够得到最低的平均误差和最大误差,可见T=28为最佳的曲面拟合窗口。随着滤波目标的形状不同,T会发生适当的变化。

      Table 2.  Effect of surface fitting window T to filter

      point cloud average error max error
      T=16 T=22 T=28 T=31 T=37 T=35 T=16 T=22 T=28 T=31 T=37 T=35
      chair 0.032 0.029 0.019 0.026 0.031 0.031 0.031 0.027 0.025 0.028 0.032 0.035
      gargo 0.041 0.032 0.027 0.031 0.037 0.042 0.039 0.039 0.037 0.034 0.041 0.046
    • 在激光对目标进行3维扫描过程中,由于环境中的粉尘干扰、周围杂散光干扰、物体表面反射材料、雨衰、测量中人操作产生的误差等因素,必然会包含一定的噪声点。噪声点会导致展现的目标形状和原实物相差很多,尤其是点云中存在一些孤立噪声点、离群点直接影响点云数据的质量和后期处理等问题。这会严重影响后续点云的处理,尤其在一些实际应用中,比如自动驾驶中的目标提取计算、3维点云立体重构中的重构精度等。因此,在对点云处理工作时,必须要进行点云去噪滤波工作。

      从前人的研究论文中可以发现,经历了多年的探索研究,有众多的去噪滤波算法,有些技术分支已经日趋成熟。众多的思路基本都是需要根据点云之间的相对位置关系,进行综合处理。但是在实际应用过程中发现,由于环境干扰,实际扫描的点云所含噪声都较大,或者噪声较多时候,点云之间的相对关系易被破坏,点云滤波的精度难以保证,滤波的处理复杂度不断增加,此时,前人的滤波方法对处理效果没有明显的增强,滤波过程中还容易产生局部细节损失现象。根据这种情况,进行了研究,从实验结果可以看到,提出的方法能解决上述问题, 得到较高精度的去噪滤波结果。

      进一步,表 3中给出了滤波窗口的大小M对整体滤波精度的影响。从表中可以看到,去噪效果随着滤波窗口M的变化而变化。M增大时,点云去噪最大误差呈减小趋势,随着M进一步增大,点云去噪最大误差逐步平缓,而后M越大,点云去噪最大误差增大。这是由于M控制滤波窗口点云采样个数,采样个数的适当增加,对滤波效果有好处,但是当增加到一定程度的时候,效果反而变差,因此,滤波窗口系数的合理选择,对滤波去噪效果有很大的影响。

      Table 3.  Effect of surface window M to filter

      point cloud average error max error
      M=10 M=15 M=20 M=25 M=30 M=35 M=10 M=15 M=20 M=25 M=30 M=35
      chair 0.028 0.026 0.018 0.023 0.027 0.031 0.029 0.025 0.024 0.025 0.028 0.031
      gargo 0.033 0.031 0.028 0.030 0.034 0.037 0.039 0.035 0.032 0.036 0.039 0.042

      同时从表 1中曲面拟合窗口大小和表 3中滤波窗口大小对处理精度比较来看,窗口大小对精度的影响略有不同,在曲面拟合过程中,窗口大一些好,在滤波过程中,窗口稍小一些好。这是因为在曲面拟合过程中,窗口适当大一些,拟合的精度高,不会变形,而在滤波过程中,小窗口对精度有提升。

      从上述实验结果及分析可以看出,该方法通过拟合运算,滤波精度高。有效处理了噪声点、离群点。对于陡峭、平坦的区域,该方法具有较好的滤波效果,适用于各种目标点云。该算法结构简单,同时能够有效保留细节信息。

      由于传统的移动曲面的最小二乘滤波只在较小区域时候有效,而且受所取曲面周围区域的大小制约,将主成分分析基础上进行曲面拟合,而后去噪,取得了较为理想的效果。与采用双边滤波算法、改进双边滤波、三边滤波或Fleishman方法等算法进行比较,提出的方法能够在较大范围内、较高噪声情况下,进行大规模散乱点云处理,去噪效果明显,细节特征保留,同时光顺区域误差也比较小。

    • 近些年来,随着3维视觉技术的飞速发展,人们对3维信息愈发感兴趣,研究也日渐深入。由于3维数据可以描述实物场景,且其结构简单,能广泛应用于军事、石化、水利、航天和考古等诸多邻域。在采样数据去噪的过程中,关键之一在于点云在去噪的时候,细节特征能够完整保留,平滑区域更加光顺。对此,国内外学者们进行了广泛深入的研究,并且取得了很多研究成果。针对点云的噪声点,需要对点云进行去噪滤波工作。将离群点去除,将散乱的、含有噪声的点云变成规则的、高精度的点云。作者提出了基于主成分分析与曲面拟合进行点云去噪。首先针对区域点云进行主成分分析计算方法,在主成分分析的法向量进行粗去噪,而后针对粗去噪声后的点云进行曲面拟合,最后根据点到曲面的距离进行了点云的高精度滤波。实验表明去噪效果精度高。下一步工作中,如何结合点云的噪声特性来自适应地优化将是下一步工作的重点。另外,基于深度学习的点云平滑去噪,也是将是下一步工作的重点之一。

参考文献 (26)

目录

    /

    返回文章
    返回