CN112396567A - 基于法向修正与位置滤波两步法的散乱点云去噪方法 - Google Patents
基于法向修正与位置滤波两步法的散乱点云去噪方法 Download PDFInfo
- Publication number
- CN112396567A CN112396567A CN202011350981.4A CN202011350981A CN112396567A CN 112396567 A CN112396567 A CN 112396567A CN 202011350981 A CN202011350981 A CN 202011350981A CN 112396567 A CN112396567 A CN 112396567A
- Authority
- CN
- China
- Prior art keywords
- normal
- point
- point cloud
- filtering
- correction
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 60
- 238000012937 correction Methods 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000005070 sampling Methods 0.000 claims abstract description 51
- 238000009499 grossing Methods 0.000 claims abstract description 11
- 230000014759 maintenance of location Effects 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000010606 normalization Methods 0.000 claims description 4
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 239000013598 vector Substances 0.000 claims description 3
- 238000000513 principal component analysis Methods 0.000 claims description 2
- 238000012847 principal component analysis method Methods 0.000 claims description 2
- 238000013178 mathematical model Methods 0.000 abstract 1
- 230000002146 bilateral effect Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000306 component Substances 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明为基于法向修正与位置滤波两步法的散乱点云去噪方法,包括以下步骤:初估点云法线;计算采样点与邻域点的空间欧式距离和法向偏差,进而求得空间欧氏距离核函数和法向偏差核函数;进行法向滤波修正,输出修正后的点云法向;将修正后的点云法向作为散乱点云的位置迭代方向,并结合滤波参数自适应的高斯核函数改进点云滤波数学模型,以采样点与邻域点的平均距离来控制平滑因子,采样点与邻域点的距离偏差在采样点上投影的标准差来控制特征保持因子,将修正后的法向、平滑因子和特征保持因子的值代入模型获得采样点沿修正后法向移动的距离,最后进行滤波迭代,输出去噪后的点云数据。实现了模型的保特征去噪,有效地避免了模型特征细节失真。
Description
技术领域
本发明涉及逆向工程中点云数据去噪预处理领域,具体涉及一种基于法向修正与位置滤波两步法的散乱点云去噪方法。
背景技术
航空发动机叶片作为航空发动机的核心部件之一,长期处于高温、高压等恶劣工作环境,极易造成损伤。研究叶片的修复技术能降低叶片的废弃率和维修成本,产生巨大的经济效益。通过逆向工程重建完整叶片过程中,获得高精度的点云数据则是极为关键的一步。然而,由于测量仪器、构件本身的缺陷以及测量环境等诸多因素的影响,使得测量所得海量点云数据中往往含有大量噪声,严重影响了利用点云模型完成曲面重建的效果。因此需要对获得的点云数据进行去噪预处理。
现有采样模型去噪算法分为两类,一类是网格模型去噪,另一类是散乱点模型去噪。散乱点云模型去噪光顺的基本原理是噪声点沿法向方向进行微小偏移,然而由于噪声的干扰使得采样点的法向估计存在较大误差,并且位置迭代滤波时滤波因子不能随采样点局部几何特征自适应变化,从而导致去噪后的模型特征失真。本发明主要针对散乱点云模型去噪光顺中存在的上述两种问题,提出一种基于法向修正与位置滤波两步法的散乱点云去噪方法,该方法在有效去除噪声的同时能够较好的保留模型的特征细节。
发明内容
针对传统散乱点云去噪模型易造成模型的过光顺和局部的失真问题,本发明提出了一种基于法向修正与位置滤波两步法的散乱点云去噪方法,本发明从法向修正和构建滤波参数自适应取值的高斯核函数两个方面对传统散乱点云滤波模型进行改进,实现了模型的保特征去噪,有效地避免了模型特征细节失真。
本发明解决所述技术问题的技术方案是:一种基于法向修正与位置滤波两步法的散乱点云去噪算法,包括以下步骤:
读取含噪声的点云数据P,设置邻域点数k、法向滤波修正次数T和滤波次数γ的值;
计算采样点pi与邻域点pij的空间欧式距离||pi-pij||和法向偏差获得距离偏差带宽σd和法向偏差带宽σn,进而求得空间欧氏距离核函数Wd和法向偏差核函数Wn;根据公式(3)进行法向滤波修正,当修正次数大于设置的法向滤波修正次数T时,停止迭代修正,输出修正后的点云法向
将修正后的点云法向作为散乱点云的位置迭代方向,并结合滤波参数自适应的高斯核函数改进点云滤波模型,
所述滤波参数自适应的高斯核函数为以采样点与邻域点的平均距离来控制平滑因子σc,采样点与邻域点的距离偏差在采样点上投影的标准差来控制特征保持因子σs,从而构建滤波参数自适应的高斯核函数中的空间距离权重核函数Wc和特征保持权重核函数Ws,构建新的滤波模型;
与现有技术相比,本发明的有益效果是:
本发明方法能够对邻域点法向进行双边加权,利用邻域法向加权之和更新采样点的法向,提高了模型尖锐特征处的法向精度,为散乱点云模型的位置滤波提供了精准的法向。构建了能随采样点局部几何特征自适应取值的滤波参数,使得位置滤波时模型的特征保持权重和光顺权重均能自适应改变。将修正后的法向和滤波参数自适应的高斯核函数相结合构建了新的点云滤波模型(公式(15)),实现了点云的保特征去噪,有效地避免了模型的过光顺和特征细节失真。
附图说明
图1为本发明基本流程示意图
图2为Fandisk模型的去噪结果。
图3为Head模型的去噪结果。
图4为Hand模型的去噪结果。
图5为修正前后法向精度对比图。
其中图2-图4中,(a)均为理想模型,(b)均为加噪声后的点云模型,(c)均为用双边滤波算法去噪结果,(d)均为本发明方法去噪结果。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,实施例是以本发明所述技术方案为前提进行的具体实施,给出了详细的实施方式和过程。但本申请的保护范围不限于下述的实施例描述。
本发明提供一种基于法向修正与位置滤波两步法的散乱点云去噪算法,该去噪算法实施包括如下步骤:
步骤1:读取含噪声的点云数据P,设置邻域点数k、法向滤波修正次数T和滤波次数γ的值。
步骤2:利用K-D树算法进行k邻域搜索,在局部k邻域内利用主成分分析法初估点云法线假定P=(pi,i=1,2,3,…,n)为原始点云数据,pi∈R3,pi1,pi2,pi3,…,pik为采样点pi的最近k邻域,主成分分析法得,pi的协方差矩阵C为
式中,为采样点pi的最近k邻域的质心,C为矩阵为3阶对称的半正定矩阵。它具有3个非负特征值λ1,λ2,λ3,其中,0≤λ1≤λ2≤λ3,e1,e2,e3分别为3个特征值所对应的特征向量。其中,最小特征值所对应的特征向量e1为采样点pi的初始估计法向,即公式(1)中T为转置。
步骤3:计算采样点与邻域点的空间欧式距离||pi-pij||和法向偏差根据采样点的局部几何特征自适应控制参数σd和σn的取值,利用采样点与邻域点的平均距离来控制距离偏差带宽σd,利用采样点法向与邻域点法向的标准偏差来控制法向偏差带宽σn,具体公式为公式(7)和(8),分别计算距离偏差带宽σd和法向偏差带宽σn。
式中,k为采样点的最近邻域点的个数,为采样点第t次修正后的法向,为采样点的第j个邻域点第t次修正后的法向,为采样点第t+1次修正后的法向,t为当前修正次数,Ki是归一化因子,Wd为空间欧氏距离核函数,Wn为法向偏差核函数,两者的具体表示如下:
本发明方法在进行散乱点云的位置滤波时,利用采样点与邻域点的平均距离来控制平滑因子σc,利用采样点与邻域点的距离偏差在采样点上投影的标准差来控制特征保持因子σs,通过对滤波参数自适应取值设置,使得特征保持权重和光顺权重能随采样点的几何特征自适应改变,有效避免了低频区去噪不完全和高频尖锐区特征失真,保留了原始点云模型的特征细节。将修正后的法向和滤波参数自适应的高斯核函数相结合构建了新的滤波模型(公式(16)),进而获得新的点云滤波模型(公式(15)),实现了点云的保特征去噪,有效地避免了模型的过光顺和特征细节失真。
图2-图4为对三种模型加入30dB的高斯白噪声,然后利用不同方法进行对比所得结果。本实施例中法向滤波修正次数T和滤波次数γ的值均为5。
三个图中的(a)均为原始模型,(b)为加入噪声后的模型,(c)为双边滤波算法处理后的结果(具体双边滤波算法参见文献Fleishman S,Drori I,Cohen-Or D.Bilateralmesh denoising[J].Acm Transactions on Graphics,2003,22(3):950-953.),(d)为本申请方法处理后的结果。
从部分放大图可以看出:本申请方法可以实现点云模型的低频平滑和高频保特征去噪。
表1是对去噪后的平均误差、均方根误差和最大误差的统计结果。从表中可以得出结论:本发明方法去噪后点云模型的三种误差均为最小,表明了本申请方法能更好地保留原始模型的特征细节,去噪结果更接近理想模型。
表1
图5是修正前后法向精度对比图,本发明方法提高了模型尖锐特征处的法向精度,为散乱点云模型的位置滤波提供了精准的法向。
利用本发明所述的技术方案,或本领域的技术人员在本发明技术方案的启发下,设计出类似的技术方案,而达到上述技术效果的,均是落入本专利的保护范围。
本发明未述及之处适用于现有技术。
Claims (5)
1.一种基于法向修正与位置滤波两步法的散乱点云去噪方法,包括以下步骤:
读取含噪声的点云数据P,设置邻域点数k、法向滤波修正次数T和滤波次数γ的值;
计算采样点pi与邻域点pij的空间欧式距离||pi-pij||和法向偏差获得距离偏差带宽σd和法向偏差带宽σn,进而求得空间欧氏距离核函数Wd和法向偏差核函数Wn;根据公式(3)进行法向滤波修正,当修正次数大于设置的法向滤波修正次数T时,停止迭代修正,输出修正后的点云法向
将修正后的点云法向作为散乱点云的位置迭代方向,并结合滤波参数自适应的高斯核函数改进点云滤波模型,
所述滤波参数自适应的高斯核函数为以采样点与邻域点的平均距离来控制平滑因子σc,采样点与邻域点的距离偏差在采样点上投影的标准差来控制特征保持因子σs,从而构建滤波参数自适应的高斯核函数中的空间距离权重核函数Wc和特征保持权重核函数Ws,构建新的滤波模型;
5.根据权利要求1所述的点云去噪方法,其特征在于,该方法用于逆向工程重建完整叶片过程中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011350981.4A CN112396567B (zh) | 2020-11-26 | 2020-11-26 | 基于法向修正与位置滤波两步法的散乱点云去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011350981.4A CN112396567B (zh) | 2020-11-26 | 2020-11-26 | 基于法向修正与位置滤波两步法的散乱点云去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112396567A true CN112396567A (zh) | 2021-02-23 |
CN112396567B CN112396567B (zh) | 2023-06-06 |
Family
ID=74604540
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011350981.4A Active CN112396567B (zh) | 2020-11-26 | 2020-11-26 | 基于法向修正与位置滤波两步法的散乱点云去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112396567B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114897026A (zh) * | 2022-05-24 | 2022-08-12 | 上海枢光科技有限公司 | 一种点云滤波的方法 |
CN116246121A (zh) * | 2023-05-12 | 2023-06-09 | 山东科技大学 | 一种基于改进最远点采样的点云数据处理方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102629367A (zh) * | 2012-01-17 | 2012-08-08 | 安徽建筑工业学院 | 基于KDTree的点云数据双边滤波去噪的方法 |
CN106709883A (zh) * | 2016-12-20 | 2017-05-24 | 华南理工大学 | 基于联合双边滤波和尖锐特征骨架提取的点云去噪方法 |
US20170193692A1 (en) * | 2015-12-30 | 2017-07-06 | Shenzhen Institutes Of Advanced Technology Chinese Academy Of Sciences | Three-dimensional point cloud model reconstruction method, computer readable storage medium and device |
CN110349094A (zh) * | 2019-06-12 | 2019-10-18 | 西安工程大学 | 基于统计离群和自适应双边混合滤波的3d点云去噪方法 |
-
2020
- 2020-11-26 CN CN202011350981.4A patent/CN112396567B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102629367A (zh) * | 2012-01-17 | 2012-08-08 | 安徽建筑工业学院 | 基于KDTree的点云数据双边滤波去噪的方法 |
US20170193692A1 (en) * | 2015-12-30 | 2017-07-06 | Shenzhen Institutes Of Advanced Technology Chinese Academy Of Sciences | Three-dimensional point cloud model reconstruction method, computer readable storage medium and device |
CN106709883A (zh) * | 2016-12-20 | 2017-05-24 | 华南理工大学 | 基于联合双边滤波和尖锐特征骨架提取的点云去噪方法 |
CN110349094A (zh) * | 2019-06-12 | 2019-10-18 | 西安工程大学 | 基于统计离群和自适应双边混合滤波的3d点云去噪方法 |
Non-Patent Citations (2)
Title |
---|
卢钰仁等: "《基于法向修正的双边滤波点云去噪处理》", 《仪表技术与传感器》 * |
袁小翠等: "《尖锐特征曲面散乱点云法向估计》", 《光学精密工程》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114897026A (zh) * | 2022-05-24 | 2022-08-12 | 上海枢光科技有限公司 | 一种点云滤波的方法 |
CN116246121A (zh) * | 2023-05-12 | 2023-06-09 | 山东科技大学 | 一种基于改进最远点采样的点云数据处理方法 |
CN116246121B (zh) * | 2023-05-12 | 2023-08-11 | 山东科技大学 | 一种基于改进最远点采样的点云数据处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112396567B (zh) | 2023-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109035163B (zh) | 一种基于深度学习的自适应图像去噪方法 | |
CN107705265B (zh) | 一种基于总曲率的sar图像变分去噪方法 | |
CN109242972B (zh) | 一种基于顶点特征的双法向网格模型光顺方法 | |
CN112396567B (zh) | 基于法向修正与位置滤波两步法的散乱点云去噪方法 | |
CN109671029A (zh) | 基于伽马范数最小化的图像去噪算法 | |
CN110930334A (zh) | 一种基于神经网络的网格去噪方法 | |
CN111028356A (zh) | 基于非凸非光滑二阶正则项和稀疏保真项的优化方法 | |
CN107341768B (zh) | 网格降噪 | |
CN110942077B (zh) | 基于权重局部变化度和l1中值优化的特征线提取方法 | |
Weijer et al. | Least squares and robust estimation of local image structure | |
Cuomo et al. | Numerical effects of the gaussian recursive filters in solving linear systems in the 3dvar case study | |
CN111105363A (zh) | 一种含噪高光谱图像快速解混方法 | |
CN108876711B (zh) | 一种基于图像特征点的草图生成方法、服务器及系统 | |
CN103761715A (zh) | 一种用于图像去噪的分数阶原始对偶方法 | |
JP4981313B2 (ja) | 3次元形状処理装置及び曲面作成プログラム並びに曲面作成方法 | |
CN115906688A (zh) | 用于跨声速风洞cfd计算的开孔壁壁面边界条件获取方法 | |
JP2011145999A (ja) | 連立一次方程式の計算方法及び装置 | |
CN114004952A (zh) | 基于统计流形曲率的带高密度噪声点云的数据处理方法 | |
CN107704724B (zh) | 基于Meridian分布的贝叶斯压缩感知的参数选取方法 | |
CN108563856B (zh) | 一种基于自由节点b样条建模的自适应采样方法 | |
CN110047046B (zh) | 图像分割和偏移场矫正方法、可读存储介质和终端 | |
Liu et al. | Toward designing convergent deep operator splitting methods for task-specific nonconvex optimization | |
CN113506212A (zh) | 一种改进的基于pocs的高光谱图像超分辨率重建方法 | |
CN109446480B (zh) | 一种基于b样条的阻尼振荡频谱数据预处理方法 | |
Volodine et al. | Smoothing of meshes and point clouds using weighted geometry-aware bases |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |