基于法矢量修正和参数修正的点云双边滤波改进算法.pdf

第4 4卷 第3期2 0 2 3年9月制 导 与 引 信GU I D AN C E&F U Z E V o l.4 4 N o.3S e p.2 0 2 3文章编号:1 6 7 1-0 5 7 6(2 0 2 3)0 3-0 0 1 6-0 7基于法矢量修正和参数修正的点云双边滤波改进算法梁敏萱,李英博,强晶晶,谷琼琼,陈蕙心(上海无线电设备研究所,上海 2 0 1 1 0 9)摘 要:针对现有点云滤波算法在实现过程中存在的目标特征保留能力差和时间复杂度高的问题,提出了基于法矢量修正和参数修正的点云双边滤波改进算法。介绍了点云双边滤波算法原理,在该算法的基础之上分别修正了法矢量、距离权重因子和特征保持因子,提出了一种改进的点云双边滤波算法,并通过仿真分析进行了算法验证。与传统的点云半径滤波算法和点云双边滤波算法相比,所提算法能够在时间复杂度较低的情况下保留更多的目标细节特征,提高滤波算法实时性。关键词:点云双边滤波;法矢量;主成分分析;高斯核函数中图分类号:T N 9 5 7.5 1文献标志码:AD O I:1 0.3 9 6 9/j.i s s n.1 6 7 1-0 5 7 6.2 0 2 3.0 3.0 0 3收稿日期:2 0 2 3-0 3-2 3作者简介:梁敏萱(1 9 9 9),女,硕士研究生,主要从事信号处理技术研究。I m p r o v e d B i l a t e r a l F i l t e r i n g A l g o r i t h m f o r P o i n t C l o u d B a s e d o n N o r m a l V e c t o r C o r r e c t i o n a n d P a r a m e t e r C o r r e c t i o nL I ANG M i n x u a n,L I Y i n g b o,Q I ANG J i n g j i n g,G U Q i o n g q i o n g,CHEN H u i x i n(S h a n g h a i R a d i o E q u i p m e n t R e s e a r c h I n s t i t u t e,S h a n g h a i 2 0 1 1 0 9,C h i n a)A b s t r a c t:A i m i n g a t t h e p r o b l e m s o f p o o r f e a t u r e r e t e n t i o n a n d h i g h t i m e c o m p l e x i t y i n t h e i m p l e m e n t a t i o n o f e x i s t i n g f i l t e r i n g a l g o r i t h m s f o r p o i n t c l o u d,a n i m p r o v e d b i l a t e r a l f i l t e r i n g a l g o r i t h m f o r p o i n t c l o u d b a s e d o n n o r m a l v e c t o r c o r r e c t i o n a n d p a r a m e t e r c o r r e c t i o n w a s p r o p o s e d.T h e p r i n c i p l e o f b i l a t e r a l f i l t e r i n g a l g o r i t h m f o r p o i n t c l o u d w a s i n t r o d u c e d.O n t h e b a s i s o f t h e a l g o r i t h m,t h e n o r m a l v e c t o r,d i s t a n c e w e i g h t f a c t o r a n d f e a t u r e r e t e n t i o n f a c t o r w e r e m o d i f i e d r e s p e c t i v e l y,a n d a n i m p r o v e d f i l t e r i n g a l g o r i t h m f o r p o i n t c l o u d w a s p r o p o s e d.T h r o u g h s i m u l a t i o n a n a l y s i s,t h e p r o p o s e d a l g o r i t h m i s v e r i f i e d.C o m p a r e d w i t h t h e t r a d i t i o n a l r a d i u s f i l t e r i n g a l g o r i t h m a n d b i l a t e r a l f i l t e r i n g a l g o r i t h m f o r p o i n t c l o u d,t h e p r o p o s e d a l g o r i t h m c a n r e t a i n m o r e t a r g e t d e t a i l s i n t h e c a s e o f l o w t i m e c o m p l e x i t y,a n d i m p r o v e t h e r e a l-t i m e p e r f o r m a n c e o f f i l t e r i n g a l g o r i t h m.K e y w o r d s:p o i n t c l o u d b i l a t e r a l f i l t e r i n g;n o r m a l v e c t o r;p r i n c i p a l c o m p o n e n t a n a l y s i s;G a u s s i a n k e r n e l f u n c t i o n第3期梁敏萱,等:基于法矢量修正和参数修正的点云双边滤波改进算法0 引言 激光三维成像雷达可以获取三维立体图像,同时具有较强的目标识别能力及低功耗、低成本的优点,逐渐被应用于军事、民用和科学实验领域。但是在实际点云数据采集过程中,受天气、环境和设备本身测量误差的影响,会引入噪声。为了最大程度发挥激光三维成像雷达的优势,需要对其采集的点云数据进行滤波处理。点云滤波能够在消除噪声的同时很好地保留目标的特征信息。从2 0世纪9 0年代开始,国内外科研人员已经围绕点云滤波开展了大量的研究,近些年来更是精益求精,提出了很多先进的点云滤波技术。S CHA L L等1受均值漂移技术的启发,利用迭代策略将点移动到高概率的位置,采用基于核的聚类滤波方法实现了点云滤波,但是该方法会导致目标特征信息丢失严重。DUAN等2提 出 了 一 种 基 于 主 成 分 分 析(p r i n c i p a l c o m p o n e n t a n a l y s i s,P C A)的自适应聚类方法,该方法计算复杂度较低,在有效去除噪声的同时保留了环境细节特征,但是需要人为设置多个阈值。F L E I S HMAN等3和J ON E S等4将图像双边滤波应用到三维点云数据处理中,提出了点云双边滤波算法。该算法鲁棒性强,具有良好的降噪能力,但是对点云尖锐区域特征保留程度有限。焦亚男等5利用半径滤波处理离群孤立点,将改进的具 有 噪 声 基 于 密 度 的 聚 类(d e n s i t y-b a s e d s p a t i a l c l u s t e r i n g o f a p p l i c a t i o n s w i t h n o i s e,D B S C AN)算法和八叉树结构相结合,建立点云间的拓扑关系来滤除小尺度噪声,但是该算法计算复杂度非常高。黄五超等6提出基于r邻域的点云去噪算法,采用最大距离阈值区分信号点与噪声点。朱广堂等7基于移动最小二乘曲率特征进行点云去噪,通过判断曲率与设定阈值的关系来去除噪声点,但是该算法的鲁棒性不强。针对现有点云滤波算法对目标特征信息保留程度不高以及人为参与因素过多的问题,本文提出一种基于法矢量修正和参数修正的点云双边滤波改进算法。该算法利用加权的主成分分析法进行法矢量估计,并基于邻域点的影响对法矢量进行修正,最后对双边滤波因子中的高斯核函数参数进行修正。1 算法原理1.1 点云双边滤波 点云的双边滤波算法由图像的双边滤波算法发展而来。在二维图像处理中,双边滤波算法是用第i个数据点pi的k邻域点加权平均的灰度值来取代原始数据点的灰度值,从而达到去除噪声的目的。在点云中,双边滤波可以简单地理解为在 点 云 法 线 方 向 上 对 点 的 位 置 进 行 修 正。点pi修正后的位置矢量pi=pi+ni(1)式中:pi为点pi当前的位置矢量;为双边滤波因子;ni为点pi当前的法矢量。双边滤波因子的计算公式为=pi jk(pi)Wcpi-pi j Ws-1 pi jk(pi)Wcpi-pi j Ws-1 (2)其中Wcpi-pi j =e x p-pi-pi j2/(22c)(3)Ws-1 =e x p-1 2/(22s)(4)式中:pi j为点pi的第j个邻域点pi j的位置矢量;k(pi)为点pi邻域点的位置矢量集;Wc()为空间域权重函数,用以控制平滑程度;为矢量的二范数运算符;Ws()为特征域权重函数,用以表明邻域点间法矢量的变化;为矢量内积运算符;ni j为点pi的第j个邻域点的当前法矢量;c,s分别为距离权重因子与特征保持因子。1.2 加权主成分分析法估计法矢量 点云法矢量是点云表面形状的重要特征参量,是描述物质表面走向的矢量。因此准确估计点云法 矢 量 有 助 于 保 留 点 云 模 型 的 详 细 特71 制 导 与 引 信 第4 4卷征,并为 点 云 分 割 和 目 标 识 别 提 供 基 础 数 据。利用 主 成 分 分 析 法 估 计 点 云 法 矢 量 的 步 骤如下。步骤1,读取点云数据集合及邻域数据集合。设给定数据点pi的坐标为(xi,yi,zi),其中i=1,2,N,N为点云数据集合中数据点的总个数。则点云数据集合P可以表示为P=pi=(xi,yi,zi),i=1,2,N (5)以数据点pi为中心,建立其k邻域拓扑关系,数据点pi的k邻域数据集合Pi j=pi j=(xi j,yi j,zi j),j=1,2,k (6)式中:(xi j,yi j,zi j)为pi第j个邻域点pi j的三维坐标。步骤2,对于数据点pi,利用最小二乘法对其k邻域内的所有点拟合局部最优平面,即最小化邻域内所有的数据点到目标平面的空间距离。以此构建目标函数F(n,d)=a r g m i n(n,d)kj=1(npi j-d)2(7)式中:n为局部最优平面的法矢量;d为全局坐标原点到局部最优平面的空间距离;a r g m i n()为目标函数最小时自变量的取值函数。步骤3,用最小二乘法拟合平面求法矢量,可以转换为利用主成分分析法求协方差矩阵最小特征值对应的特征矢量,即将求解式(7)转化为求解协方差矩阵Ci。Ci的求解公式为Ci=1kkj=1(pi j-p-i)(pi j-p-i)TCiem=mem,m0,1,2 (8)式中:p-i表示数据点pi邻域内所有点位置的均值,即重心位置矢量;T为矩阵转置运算符;m和em分别表示协方差矩阵Ci的特征值和其对应的特征矢量,且012。最小特征值0对应的 特 征 矢 量e0就 近 似 为 数 据 点pi的 法矢量。传统的主成分分析法在利用拟合局部最优平面方法估计法矢量时,对k邻域内的所有点赋予同样的权重。该加权方式对平坦面法矢量的估计影响不大,但因忽略了复杂曲面等特殊条件下点云分布的差异,会导致复杂曲面的点云法矢量估计结果失真。R E N等8将邻域点权重引入主成分分析法中,使改进的主成分分析法能够准确估计点云的法矢量和特征值,从而保证后续滤波能够取得良好的效果。随机选取点云数据中的数据点pi,用dj表示点pi与其第j个邻域点pi j之间的距离,则pi j的权重wj可以表示为wj=d2jkj=1d2j -1(9)计算邻域点位置矢量的加权平均p-i=kj=1wjpi jkj=1wj(1 0)根据式(1 0),加权协方差矩阵Ci(wj)=1kkj=1wj(pi j-p-i)(pi j-p-i)T(1 1)空间距离越近的点,其几何特征越相似。这种加权方法更重视距离当前点pi更近的邻域点的影响,而减小了较远的邻域点对法矢量估计的干扰,提高了曲率变化较大的点云法矢量估计的鲁棒性。1.3 法矢量修正 基于加权主成分分析法估计点云法矢量具有二义性,该方法只能确定法矢量所在直线,而无法确定其方向。假设点云中点的数量足够多且采样平面处处光滑,那么相邻两点的法矢量会接近于平行9。令ni和nj分别为相邻两点pi和pj的法矢量,如果ninj0,说明两个法矢量的方向不一致,则将nj翻转,以此对求出的法矢量进行重定向。协方差矩阵各特征值对应的特征矢量为基矢量,而基矢量的各个维度下信息的丰富程度与特征值线性相关,所以可以用最大特征值对应的特征矢量来确定多维数据的主成分。主成分分析法就是通过删除一部分分辨率较低的数据来实现降维。虽然主成分分析法在计算法矢量时有一定的抗噪能力,但是如果点云中存在离群点,法矢量估计的准确性将会受到很大的影响,因此需要对计算得到的法矢量进行修正。将反映 物 体 特 征 的 点 称 为 特 征 点,反 之,则称为非 特征点。文献1 0 提 出基于双边 滤波的法矢量修正方法,修正的法矢量n i的表达式为81第3期梁敏萱,等:基于法矢量修正和参数修正的点云双边滤波改进算法n i=1kkj=1G1pi j-pi G2ni j-ni ni j(1 2)其中G1pi j-pi =e x p-pi j-pi2/(221)G2ni j-ni =e x p-ni j-ni2/(222)1=0.5,2(0.5,1.5)(1 3)式中:G1()和G2()为高斯函数;1和2为高斯核函数超参数。采用基于空间距离的加权主成分分析法可以降低离群点对法矢量估计的影响,但是该方法会影响特征点的法矢量估计。文献1 0 综合考虑了空间距离和法矢量差异来实现对点云法矢量的修正。如式(1 3)所示,G2ni j-ni 用于描述法矢量之间的差异,G2ni j-ni 的引入可以降低加权主成分分析法对特征点矢量估计平滑作用的影响。ni j-ni的值越小,意味着法线夹 角 越 小、两 个 点 的 几 何 特 征 越 相 近。但文献1 0 中的高斯核函数超参数是人为设置的。本文为了降低人为主观因素的影响,基于数据本身构造适应性的高斯核函数超参数,即21=kj=1pi j-pi2/(k-1)22=kj=1ni j-ni2/(k-1)(1 4)针对椅子部分点云数据,传统主成分分析法求得的未修正的法矢量如图1所示,加权主成分分析法求得的法矢量根据式(1 2)修正后的结果如图2所示。由图1和图2可以看出,“椅子”点云修正后的法矢量明显比未修正的法矢量有序,更有利于后续双边滤波正确改变点云位置。1.4 空间距离权重因子与特征保持因子的修正 空间距离权重因子c反映的是待处理的数据点pi与其k邻域中点的距离对滤波的影响。c越大,滤波后的模型越平滑,这也就意味着点云的细节特征将会被平滑。特征保持因子s反映的是数据点pi与其k邻域中点的距离在法向上的投影对滤波的影响。s越大,特征保持效果越好。图1“椅子”点云未修正的法矢量图2“椅子”点云修正后的法矢量目前传统点云双边滤波算法仍然需要根据点云数据,人为设定空间距离权重因子和特征保持因子。而本文利用局部特征权重和高斯核函数自适应设置双边滤波参数。文献1 1 中基于曲率信息计算数据点pi的k邻域局部特征权重wi的表达式为wi=1kkj=1|Ki j|-K-2+Ki-K-2(1 5)其中K-=1NNi=1Ki(1 6)式中:Ki j和Ki分别为邻域点pi j和当前点pi的平均曲率;K-为所有点云的平均曲率。如果点云数据中的当前点pi的k邻域局部特征权重wi大于给定的阈值w0,则认为该点为尖锐区域点;否则,则认为该点为平坦区域点。文献1 2 给出的权重阈值w0的表达式为91 制 导 与 引 信 第4 4卷w0=(wm i n+wm a x)/t(1 7)式中:wm i n和wm a x分别为所有点云数据的局部特征值中的最小值和最大值;t为通过实验得到的经验参数。对于平坦区域的点云,可以通过增大c来增加该区域的平滑效果,本文选取邻域点与当前点的最大距离作为c。此时s的作用不大,可以用当前点与邻域内各点法矢量的差异标量的标准差来表示。同理,对于尖锐区域的点云,可以通过减小c来提高对细节特征的保留程度,本文选取当前点与邻域内各点法矢量的差异标量的最小值作为s。此时c的作用不大,可以用当前点与邻域内各点距离的标准差来表示。综上所述,c和s可分别表示为c=m a x pi j-pi ,wiw0kj=1pi j-pi2/(k-1),wiw0 (1 8)s=kj=1-12/(k-1),wiw0m i n(-1),wiw0 (1 9)式中:m a x()为最大值函数;|为绝对值运算符;m i n()为最小值函数。综上所述,将经过修正的法矢量、空间距离权重因子和特征保持因子代入点云双边滤波算法中,即可得到改进的点云双边滤波算法。2 实验验证 实验采用斯坦福大学公开数据集中的兔子模型,利用M a t l a b 2 0 1 7进行程序的编写,邻域范围取k=4 0。2.1 空间距离权重因子与特征保持因子的选取 文献1 3 利用点云的平均曲率分布计算当前点的局部熵,用于描述当前点所代表的局部曲面变分的有序程度。局部熵越大,该点所在局部区域的无序程度越高,该点提供的信息量越大,表征的细节特征越多。定义当前点pi的k邻域局部熵Hi为Hi=-qil o g2qi-kj=1qi jl o g2qi j(2 0)其中qi=Ki/Ki+kj=1Ki j (2 1)qi j=Ki j/Ki+kj=1Ki j (2 2)式中:qi和qi j分别为当前点pi和其邻域点pi j的平均曲率的分布概率。若点云总数量为N,则所有点云的熵H=Ni=1Hi(2 3)参考文献1 2 和文献1 4,根据c和s的特性可知,c和s分别有三种可能的取值,如表1所示。表1 空间距离权重因子c和特征保持因子s的不同取值序号cskj=1pi j-pi2/(k-1)kj=1-12/(k-1)m a x pi j-pi ,wiw0kj=1pi j-pi2/(k-1),wiw0 kj=1-12/(k-1),wiw0m i n(-1),wiw0 m a x(pi j-pi),wiw01kkj=1pi j-pi,wiw0 kj=1ni j-ni2/(k-1),wiw01kkj=1-1,wiw0 文献1 4 还提出了一种利用法向离群因子改进双边滤波算法的方法。本文利用全局信息熵对各种组合进行判定,如表2所示。02第3期梁敏萱,等:基于法矢量修正和参数修正的点云双边滤波改进算法表2 各种组合下的全局信息熵csH有法向离群因子无法向离群因子3.3 1 4 9 7 1 7 0 7 0 4 0 5 0 81 091 43.3 1 5 2 4 6 8 4 1 6 5 0 5 3 81 093.3 1 4 9 5 9 3 9 8 6 3 2 4 9 11 093.3 1 5 2 6 3 9 5 5 9 5 4 4 1 71 093.3 1 4 9 5 9 2 2 3 9 3 8 2 9 81 091 23.3 1 5 2 6 3 7 2 2 6 5 7 6 4 51 093.3 1 4 9 7 1 8 4 1 8 7 5 2 2 81 093.3 1 5 2 4 6 8 3 7 2 0 1 7 4 31 093.3 1 4 9 5 9 5 3 3 5 3 4 8 7 61 093.3 1 5 2 6 4 8 8 2 3 8 7 6 8 11 09(本文)3.3 1 4 9 5 9 3 6 0 3 8 6 8 1 81 093.3 1 5 2 6 3 7 2 0 2 8 4 7 2 11 093.3 1 4 9 7 2 9 4 0 7 0 2 0 4 51 093.3 1 5 2 4 0 0 8 3 3 0 4 2 9 71 093.3 1 4 9 5 8 8 2 4 9 2 3 4 4 61 093.3 1 5 2 5 9 0 6 8 6 6 6 2 2 31 093.3 1 4 9 5 8 6 5 1 7 8 4 3 6 81 093.3 1 5 2 5 8 8 3 7 6 2 2 0 0 41 09 由表2可知,相比于文献1 2、文献1 4 以及其他排列组合,本文算法的熵更大,保留的目标特征信息更多。2.2 定量与定性评估 现有的滤波质量评价方法大多采用去噪前后点云数量的比值以及误差率作为评价参数。但是这种方法无法评估本身就带有噪声的源数据,因为其无法确定去除的点是否为噪声点,也无法判定目标细节特征的保留程度。文献1 4 提出了一种基于隐式移动最小二乘(i m p l i c i t m o v i n g l e a s t s q u a r e s,I ML S)的评估方法,利用当前点pi的邻域点构建隐式曲面,通过计算当前点到曲面的加权距离的平均值,即评估分数SI ML S,来评估滤波算法对细节特征的保留程度。细节特征保留程度越高,则评估分数SI ML S越低。SI ML S的计算公式为SI ML S=1NfNfi=1|dI ML S(i)|(2 4)其中dI ML S(i)=kj=1 Wcpi j-pi (pi j-pi)ni j)kj=1 Wcpi j-pi (2 5)式中:Nf为 经 过 滤 波 后 的 点 云 数 据 总 量;dI ML S(i)为利用当前点与邻域点之间的空间距离及其在邻域点法线方向上的投影共同得到的加权距离;为法向离群因子。采用传统双边滤波算法及本文算法对兔子模型进行滤波,滤波效果比较如图3所示。可以看出,无论是兔子的五官、背上的花纹,还是脚趾的形状,采用本文算法滤波后兔子模型能被清晰地刻画出来,滤波效果明显优于传统双边滤波算法。图3 传统双边滤波算法与本文算法滤波效果比较 采用I ML S评估方法对不同滤波算法的性能进行评价,结果如表3所示。可以看出,文献1 5提到的半径滤波算法对目标细节特征的保留程度最差,而且耗时最长。除此之外,使用半径滤波算12 制 导 与 引 信 第4 4卷法需要不断调整参数以达到最优效果,人为因素干扰过多。与双边滤波算法相比,本文算法时间复杂度略高,但是保留了更多的目标细节特征。上述分析证明本文算法能够应用到实际点云处理中,实现有效滤波。表3 不同滤波算法性能评估表滤波算法SI ML S 耗时/s半径滤波6.2 5 3 1 6 9 2 7 4 8 6 1 9 9 71 0-67 3.1 3 4 双边滤波2.7 5 3 0 7 5 8 0 1 2 5 0 4 5 41 0-63 0.1 9 5 本文算法9.7 4 1 5 5 9 7 3 9 6 7 3 1 8 81 0-74 9.1 8 6 3 结论 针对现有点云滤波算法目标特征保留程度不高以及人为参与因素过多的问题,本文提出一种基于法矢量修正和参数修正的点云双边滤波改进算法。利用熵对所有空间距离权重因子与特征保持因子的组合进行评估,选出性能最佳的组合作为高斯核函数的参数,降低主观因素影响。利用滤波效果图和I ML S评估方法分别定性和定量地比较双边滤波改进前后的滤波效果。实验结果证明,本文算法能够详细刻画目标特征。算法的有效性得到了验证。参考文献1 S CHA L L O,B E L YA E V A,S E I D E L H P.R o b u s t f i l t e r i n g o f n o i s y s c a t t e r e d p o i n t d a t aC/P r o c e e d i n g s E u r o g r a p h i c s/I E E E VG T C S y m p o s i u m P o i n t-b a s e d G r a p h i c s.P i s c a t a w a y,N J:I E E E P r e s s,2 0 0 5:7 1-7 7.2 DUANY,YAN G C C,C H E N C,e t a l.L o w-c o m p l e x i t y p o i n t c l o u d f i l t e r i n g f o r L i D A R b y P C A-b a s e d d i m e n s i o n r e d u c t i o nJ.O p t i c s C o mm u n i-c a t i o n s,2 0 2 1,4 2 8:1 2 6 5 6 7.3 F L E I S HMAN S,D R O R I I,C OHE N-OR D.B i l a t e r a l m e s h d e n o i s i n gJ.A CM T r a n s a c t i o n s o n G r a p h i c s,2 0 0 3,2 2(3):9 5 0-9 5 3.4 J ON E S T R,DUR AN D F,D E S B R UN M.N o n-i t e r a t i v e,f e a t u r e-p r e s e r v i n g m e s h s m o o t h i n gJ.A CM T r a n s a c t i o n s o n G r a p h i c s,2 0 0 3,2 2(3):9 4 3-9 4 9.5 焦亚男,马杰,钟斌斌.一种基于尺度变化的点云并行去噪方法J.武汉大学学报(工学版),2 0 2 1,5 4(3):2 7 7-2 8 2.6 黄五超,韩玲,黄勃学,等.一种运用P C A的深度学习激光点云分类方法J.计算机与现代化,2 0 2 1(5):3 1-3 7.7 朱广堂,叶珉吕.基于曲率特征的点云去噪及定量评价方法研究J.测绘通报,2 0 1 9(6):1 0 5-1 0 8.8 R E N Y J,L I T Z,XU J K,e t a l.O v e r a l l f i l t e r i n g a l g o r i t h m f o r m u l t i s c a l e n o i s e r e m o v a l f r o m p o i n t c l o u d d a t aJ.I E E E A c c e s s,2 0 2 1,9:1 1 0 7 2 3-1 1 0 7 3 4.9 邢正全,邓喀中,薛继群.基于K-近邻搜索的点云初始配准J.测绘科学,2 0 1 3,3 8(2):9 3-9 5.1 0 Z T I R E L I A C,GU E NN E B AU D G,G R O S S M.F e a t u r e p r e s e r v i n g p o i n t s e t s u r f a c e s b a s e d o n n o n-l i n e a r k e r n e l r e g r e s s i o nJ.C o m p u t e r G r a p h i c s F o r u m,2 0 0 9,2 8(2):4 9 3-5 0 1.1 1 程效军,贾东峰,程小龙.海量点云数据处理理论与技术M.上海:同济大学出版社,2 0 1 4:1 8-8 8.1 2 刘静静.三维点云重建中的去噪算法研究D.北京:北京交通大学,2 0 1 9.1 3 武剑洁.基于点的散乱点云处理技术的研究D.武汉:华中科技大学,2 0 0 4.1 4 陈亚超,樊彦国,禹定峰,等.考虑法向离群的自适应双边滤波点云平滑及I ML S评价方法J.图学学报,2 0 2 3,4 4(1):1 3 1-1 3 8.1 5 鲁冬冬,邹进贵.三维激光点云的降噪算法对比研究J.测绘通报,2 0 1 9(增刊2):1 0 2-1 0 5.22