CN114663373A - 一种用于零件表面质量检测的点云配准方法及装置 - Google Patents

一种用于零件表面质量检测的点云配准方法及装置 Download PDF

Info

Publication number
CN114663373A
CN114663373A CN202210239157.4A CN202210239157A CN114663373A CN 114663373 A CN114663373 A CN 114663373A CN 202210239157 A CN202210239157 A CN 202210239157A CN 114663373 A CN114663373 A CN 114663373A
Authority
CN
China
Prior art keywords
point cloud
geometric
point
points
scanning
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
Application number
CN202210239157.4A
Other languages
English (en)
Other versions
CN114663373B (zh
Inventor
宋宝
李君�
刘永兴
周向东
唐小琦
王小柏
卢慧锋
张翔
高天赐
林志鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN202210239157.4A priority Critical patent/CN114663373B/zh
Publication of CN114663373A publication Critical patent/CN114663373A/zh
Application granted granted Critical
Publication of CN114663373B publication Critical patent/CN114663373B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/23Clustering techniques
    • G06F18/232Non-hierarchical techniques
    • G06F18/2321Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • G06T7/344Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods involving models
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10028Range image; Depth image; 3D point clouds

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Biology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Probability & Statistics with Applications (AREA)
  • Quality & Reliability (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Image Analysis (AREA)

Abstract

本发明涉及一种用于零件表面质量检测的点云配准方法及装置,本发明首先获取待测零件扫描点云集合及其对应的标准模型,并对点云集合进行预处理;将处理后的点云集合与标准模型进行初始配准;根据零件表面的几何特征将扫描点云集合划分成基础部分与几何特征部分,并构建匹配点对;利用匹配点对,计算协方差矩阵并进行奇异值分解,分别计算各部分对应的变换参数;为各部分分配不同的权重并计算总体变换参数,利用总体变换参数变换待测零件的扫描点云,重复迭代直至满足终止条件。本发明基于几何特征拟合的点云数据预处理及配准方法,对于存在大量背景及噪声的扫描点云能够提高配准稳定性,并提高配准精度及速度。

Description

一种用于零件表面质量检测的点云配准方法及装置
技术领域
本发明涉及图像处理技术领域,具体涉及一种用于零件表面质量检测的基于几何特征拟合的点云配准方法及装置。
背景技术
近几年来,零件表面质量三维检测技术已在工业制造领域广泛应用。其关键步骤为点云配准,即零件表面扫描得到的点云数据与通过三维软件创建的标准模型精确对齐。
基于三维检测技术进行零件表面质量检测一般是通过将零件扫描点云与标准模型配准后,两点云作差进行对比评估实现,因此,保证点云配准精度和稳定性极其重要。点云配准一般包括扫描点云预处理和配准两个阶段,在预处理阶段,通过体素栅格下采样减少点云数据量,实现保持点云结构的基础上降低数据处理难度;通过半径滤波、高斯滤波等滤波处理,实现点云降噪,去除离散点。在配准阶段,由于零件扫描点云与通过三维软件建模得到的标准模型处于两个不同的坐标系下,初始位置相隔较远,且点云方向未知,为保证配准精度,进行初始配准获得初始变换矩阵,为精配准提供良好的初始位置;然后在初始配准的基础上,采用ICP(Iterative Closest Point,迭代最近点)算法进行精配准,通过不断迭代优化点云中最近点距离,实现两点云的精确对齐。
由于零件的扫描环境一般较为复杂,通常存在辅助工具、地面等背景,扫描点云中除待测零件点云外可能会存在包含零件点云等多个高密度点云簇,常规滤波处理对离散噪点等干扰信息较为有效,但对于高密度点云簇的干扰,通常无法得到良好的预处理效果,从而也会影响点云配准精度与稳定性。
发明内容
本发明针对现有技术中存在的技术问题,利用零件具有柱面、球面等几何特征的特点,提出一种基于几何特征拟合的点云数据预处理及配准方法及装置。
本发明解决上述技术问题的技术方案如下:
第一方面,本发明提供一种用于零件表面质量检测的点云配准方法,包括:
S1,分别获取待测零件扫描点云集合P以及待测零件对应的标准模型Q,并对所述扫描点云集合P进行预处理得到点云集合Pt
S2,将点云集合Pt与标准模型Q进行初始配准;
S3,根据零件表面所具有的几何特征将点云集合Pt划分成基础部分P0与几何特征部分P1~Pn,并构建基础部分和几何特征部分中各点与标准模型Q中各点的匹配点对;其中基础部分为非几何特征部分,n为零件表面所具有的几何特征数量;
S4,利用基础部分和几何特征部分对应的匹配点对,计算协方差矩阵并进行奇异值分解,分别计算得到P0~Pn对应的变换参数
Figure BDA0003543593320000021
R为旋转矩阵,T为平移矩阵;
S5,根据零件表面具有的几何特征及其点云数据占比,为基础部分和几何特征部分分配不同的权重θk,并计算总体变换参数
Figure BDA0003543593320000022
S6,利用总体变换参数v变换待测零件的扫描点云,并重复步骤S3~S5,直至满足迭代终止条件,输出配准结果。
进一步的,将待测零件的扫描点云集合Pt与标准模型Q进行初始配准时,包括:
构建标准模型Q的Kd-tree结构,对于扫描点云集合Pt中的各点p,利用Kd-tree在标准模型Q中搜索最近点q,构成N个匹配点对。
进一步的,所述的几何特征的权重根据下式进行分配:
Figure BDA0003543593320000031
式中,N为匹配点对数量,N(Pk)为N个匹配点对中几何特征Pk的数量。
进一步的,步骤S5还包括:
1)判断变换后的扫描点云与标准模型的配准误差是否小于误差函数阈值;
2)判断迭代次数是否达到设定的最大迭代次数;
若1)或2)任一判断结果为是,则表示满足迭代终止条件。
进一步的,所述配准误差通过下式计算:
Figure BDA0003543593320000032
式中,pi表示第i个匹配点对中待测零件扫描点云中的点,qi表示第i个匹配点对中标准模型Q中的点。
进一步的,所述预处理,包括:
S601,对于未经预处理的扫描点云数据集P,根据点云密度将连续分布的点聚类,将扫描点云数据集划分成m个点云簇集合{C1,C2,...,Cm};
S602,根据待测零件的几何特征建立模型方程为:M=f(x,y,z,s,D),(x,y,z)为数据集P中点p的坐标值,s为求解模型所需的最少点的数量,该模型存在固有特征与模型表面上的点满足关系D;
S603,选取一个点云簇集合Cl,l=1~m,并从点云簇集合Cl中随机选择一个至少包含s个点的子集S,将子集S中各点的坐标根据模型方程M进行拟合得到模型M*
S604,对于点云簇集合Cl中其余点pt,计算其与模型M*的固有特征关系D*的偏离值评价误差erro(pt)=|D*-D|;若erro(pt)≤max erro,则保留该点,max erro为预设的误差阈值;当保留点的数量大于设定值时,则点云簇集合Cl即为待测零件对应的扫描点云Pt;删除除点云簇集合Cl外的其他点云簇集合。
进一步的,步骤S604还包括:
当保留点的数量大于设定值时,利用所保存的点根据最小二乘法拟合出新的几何模型并替换M*
第二方面,本发明提供一种用于零件表面质量检测的点云配准装置,包括:
数据获取及预处理模块,分别获取待测零件扫描点云集合P以及待测零件对应的标准模型Q,并对所述扫描点云集合P进行预处理得到点云集合Pt
初始配准模块,将点云集合Pt与标准模型Q进行初始配准;
点云划分及匹配点对构建模块,根据零件表面所具有的几何特征将点云集合Pt划分成基础部分P0与几何特征部分P1~Pn,并构建基础部分和几何特征部分中各点与标准模型Q中各点的匹配点对;其中基础部分为非几何特征部分,n为零件表面所具有的几何特征数量;
第一变换参数计算模块,利用基础部分和几何特征部分对应的匹配点对,计算协方差矩阵并进行奇异值分解,分别计算得到P0~Pn对应的变换参数
Figure BDA0003543593320000041
R为旋转矩阵,T为平移矩阵;
第二变换参数计算模块,根据零件表面具有的几何特征及其点云数据占比,为基础部分和几何特征部分分配不同的权重θk,并计算总体变换参数
Figure BDA0003543593320000042
利用总体变换参数v变换待测零件的扫描点云;
终止条件判断及输出模块,用于判断是否满足迭代终止条件,若满足则输出配准结果。
第三方面,本发明提供一种电子设备,包括:
存储器,用于存储计算机软件程序;
处理器,用于读取并执行所述计算机软件程序,进而实现本发明第一方面所述的一种用于零件表面质量检测的点云配准方法。
第四方面,本发明提供一种非暂态计算机可读存储介质,所述存储介质中存储有用于实现本发明第一方面所述的一种用于零件表面质量检测的点云配准方法的计算机软件程序。
本发明的有益效果是:本发明利用零件具有柱面、球面等几何特征的特点,提出一种基于几何特征拟合的点云数据预处理及配准方法,对于存在大量背景及噪声的扫描点云能够提高配准稳定性,并提高配准精度及速度。
附图说明
图1为本发明实施例提供的一种用于零件表面质量检测的点云预处理及配准方法流程示意图;
图2为本发明实施例提供的引入集合特征权重的ICP配准算法流程图;
图3为本发明实施例提供的一种用于零件表面质量检测的点云配准装置结构示意图;
图4为本发明实施例提供的一种电子设备实施例示意图;
图5为本发明实施例提供的一种非暂态计算机可读存储介质实施例示意图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
本发明实施例以锥齿轮零件为对象,对本方案进行说明。其中齿轮标准模型采用SolidWorks建模得到,并对模型表面均匀采样。
如图1所示,本发明提供一种用于零件表面质量检测的点云预处理及配准方法,主要包括预处理和配准两个部分。
1、预处理阶段:
由于零件的扫描环境一般较为复杂,通常存在辅助工具、地面等背景,扫描点云中除待测零件点云外可能会存在包含零件点云等多个高密度点云簇,常规滤波处理对离散噪点等干扰信息较为有效,但对于高密度点云簇的干扰,通常无法得到良好的预处理效果,从而也会影响点云配准精度与稳定性。
本发明实施例利用基于密度聚类以及几何特征拟合的预处理算法对待测零件的扫描点云数据集进行预处理。
为了降低噪声,去除离散群点,在获取待测零件扫描点云后利用点云数据间的距离差异,进行基于密度聚类的点云划分。
步骤1.1,对于点云数据集P,根据点云密度将连续分布的点聚类,使任一点的r*邻域内数量满足num>minnum,minnum为给定邻域最小数量,从而将点云数据集P划分成m个点云簇集合{C1,C2,...,Cm}。对于同一点云簇集合中的任意两点pA和pB满足distance(pA,pB)≤r*,即同一点云簇集合中任意两点间的距离小于等于r*。在本实施例中,设置近邻搜索半径r*=1.2mm,邻域内最少点数量minnum=10,为了滤掉噪声及背景部分,设定聚类需要的最少点数目为10k。
通过零件中所存在的平面、圆柱、球等几何特征拟合解决点云簇集合中目标点云的判断与提取问题。
根据待测零件的几何特征建立模型方程为:M=f(x,y,z,s,D),(x,y,z)为扫描点云数据集P中点p的坐标值,s为求解模型所需的最少点的数量,该模型存在固有特征与模型表面上的点满足关系D。
步骤1.2,由于锥齿轮具有圆柱几何特征,因此构建圆柱模型方程。以轴线作为圆柱模型的固有特征,则圆柱表面上任一点p(x,y,z)到圆柱轴线的距离满足关系D=r,r为圆柱底面半径。记圆柱轴线的单位向量(i,j,k),轴线上任一点p0(x0,y0,z0),则圆柱模型方程为:
(x-x0)2+(y-y0)2+(z-z0)2=r2+[i(x-x0)+j(y-y0)+k(z-z0)]2
求解圆柱模型至少需要7个点求出x0、y0、z0、i、j、k、r7个参数,即s=7。
步骤1.3,选取一个点云簇集合Cl,l=1~m,并从点云簇集合Cl中随机选择一个至少包含s个点的子集S,将子集S中各点的坐标根据模型方程M进行拟合得到模型M*
步骤1.4,对于点云簇集合Cl中其余点pt,计算其与模型M*的固有特征关系D*的偏离值评价误差erro(pt)=|D*-D|;若erro(pt)≤maxerro,则保留该点,maxerro为预设的误差阈值;当保留点的数量大于设定值时,则点云簇集合Cl即为待测零件对应的扫描点云Pt;删除除点云簇集合Cl外的其他点云簇集合。
由于通过erro(pt)≤maxerro保留的点,与模型M*之间存在一定误差,因此在保留的点的数量大于设定值时,利用最小二乘法重新拟合模型并更新M*。同时为了进一步保证拟合准确度,可以从点云簇集合中多次选择不同的子集S,根据满足误差的点的数量或调整误差阈值,确定最优几何模型。
2、配准阶段
本发明实施例在配准阶段在传统ICP算法中引入几何特征权重进行改进,即先采用传统ICP算法对经过预处理后的待测零件的扫描点云集合Pt与标准模型Q进行初始配准,再根据零件表面具有的几何特征及其点云数据占比,为每种几何特征分配不同的权重,然后求解变换矩阵R、T,使预处理后扫描点云Pt经过平移、旋转变换后与标准模型Q对齐。其算法流程如图2所示。
步骤2.1,针对经过预处理的待测零件对应的扫描点云Pt,根据零件表面所具有的几何特征将扫描点云集合Pt划分成基础部分P0与几何特征部分P1~Pn;其中基础部分为非几何特征部分,n为零件表面所具有的几何特征数量。
步骤2.2,为了提高点云搜索速度,构建标准模型Q的Kd-tree结构。
步骤2.3,对于扫描点云集合Pt中的各点,利用Kd-tree在Q中搜索最近点,构成N个匹配点对。
步骤2.4,利用N个匹配点对,计算协方差矩阵并进行奇异值分解,分别计算得到P0~Pn对应的变换参数
Figure BDA0003543593320000071
R为旋转矩阵,T为平移矩阵。
步骤2.5,根据零件表面具有的几何特征及其点云数据占比,根据下式为每种几何特征分配不同的权重θk
Figure BDA0003543593320000081
式中,N为匹配点对数量,N(Pk)为N个匹配点对中几何特征Pk的数量。
然后计算总体变换参数
Figure BDA0003543593320000082
利用总体变换参数v变换待测零件的扫描点云,并重复步骤2.3~S2.5,直至满足迭代终止条件,输出配准结果。
1)判断变换后的扫描点云与标准模型的配准误差是否小于误差函数阈值;
2)判断迭代次数是否达到设定的最大迭代次数;
若1)或2)任一判断结果为是,则表示满足迭代终止条件。
所述配准误差通过下式计算:
Figure BDA0003543593320000083
式中,pi表示第i个匹配点对中待测零件扫描点云中的点,qi表示第i个匹配点对中标准模型Q中的点。
本发明方法利用零件具有柱面、球面等几何特征的特点,提出一种基于几何特征拟合的点云数据预处理及配准方法,对于存在大量背景及噪声的扫描点云能够提高配准稳定性,并提高配准精度及速度。
基于上述方法,本发明实施例还提供一种用于零件表面质量检测的点云配准装置,如图3所示,包括:
数据获取及预处理模块,分别获取待测零件扫描点云集合P以及待测零件对应的标准模型Q,并对所述扫描点云集合P进行预处理得到点云集合Pt
初始配准模块,将点云集合Pt与标准模型Q进行初始配准;
点云划分及匹配点对构建模块,根据零件表面所具有的几何特征将点云集合Pt划分成基础部分P0与几何特征部分P1~Pn,并构建基础部分和几何特征部分中各点与标准模型Q中各点的匹配点对;其中基础部分为非几何特征部分,n为零件表面所具有的几何特征数量;
第一变换参数计算模块,利用基础部分和几何特征部分对应的匹配点对,计算协方差矩阵并进行奇异值分解,分别计算得到P0~Pn对应的变换参数
Figure BDA0003543593320000091
R为旋转矩阵,T为平移矩阵;
第二变换参数计算模块,根据零件表面具有的几何特征及其点云数据占比,为基础部分和几何特征部分分配不同的权重θk,并计算总体变换参数
Figure BDA0003543593320000092
利用总体变换参数v变换待测零件的扫描点云;
终止条件判断及输出模块,用于判断是否满足迭代终止条件,若满足则输出配准结果。
在此基础上,该系统还包括预处理模块,所述预处理模块包括:
聚类模块,对于未经预处理的扫描点云数据集P,根据点云密度将连续分布的点聚类,将扫描点云数据集划分成m个点云簇集合{C1,C2,...,Cm};
模型构建模块,根据待测零件的几何特征建立模型方程为:M=f(x,y,z,s,D),(x,y,z)为数据集P中点p的坐标值,s为求解模型所需的最少点的数量,该模型存在固有特征与模型表面上的点满足关系D;
拟合模块,选取一个点云簇集合Cl,l=1~m,并从点云簇集合Cl中随机选择一个至少包含s个点的子集S,将子集S中各点的坐标根据模型方程M进行拟合得到模型M*
评价筛选模块,对于点云簇集合Cl中其余点pt,计算其与模型M*的固有特征关系D*的偏离值评价误差erro(pt)=|D*-D|;若erro(pt)≤max erro,则保留该点,max erro为预设的误差阈值;当保留点的数量大于设定值时,则点云簇集合Cl即为待测零件对应的扫描点云Pt;保留点云簇集合Cl,并删除其他点云簇集合。
请参阅图4,图4为本发明实施例提供的电子设备的实施例示意图。如图4所示,本发明实施例提了一种电子设备500,包括存储器510、处理器520及存储在存储器520上并可在处理器520上运行的计算机程序511,处理器520执行计算机程序511时实现以下步骤:
S1,分别获取待测零件扫描点云集合P以及待测零件对应的标准模型Q,并对所述扫描点云集合P进行预处理得到点云集合Pt
S2,将点云集合Pt与标准模型Q进行初始配准;
S3,根据零件表面所具有的几何特征将点云集合Pt划分成基础部分P0与几何特征部分P1~Pn,并构建基础部分和几何特征部分中各点与标准模型Q中各点的匹配点对;其中基础部分为非几何特征部分,n为零件表面所具有的几何特征数量;
S4,利用基础部分和几何特征部分对应的匹配点对,计算协方差矩阵并进行奇异值分解,分别计算得到P0~Pn对应的变换参数
Figure BDA0003543593320000101
R为旋转矩阵,T为平移矩阵;
S5,根据零件表面具有的几何特征及其点云数据占比,为基础部分和几何特征部分分配不同的权重θk,并计算总体变换参数
Figure BDA0003543593320000102
S6,利用总体变换参数v变换待测零件的扫描点云,并重复步骤S3~S5,直至满足迭代终止条件,输出配准结果。
请参阅图5,图5为本发明实施例提供的一种计算机可读存储介质的实施例示意图。如图5所示,本实施例提供了一种计算机可读存储介质600,其上存储有计算机程序611,该计算机程序611被处理器执行时实现如下步骤:
S1,分别获取待测零件扫描点云集合P以及待测零件对应的标准模型Q,并对所述扫描点云集合P进行预处理得到点云集合Pt
S2,将点云集合Pt与标准模型Q进行初始配准;
S3,根据零件表面所具有的几何特征将点云集合Pt划分成基础部分P0与几何特征部分P1~Pn,并构建基础部分和几何特征部分中各点与标准模型Q中各点的匹配点对;其中基础部分为非几何特征部分,n为零件表面所具有的几何特征数量;
S4,利用基础部分和几何特征部分对应的匹配点对,计算协方差矩阵并进行奇异值分解,分别计算得到P0~Pn对应的变换参数
Figure BDA0003543593320000111
R为旋转矩阵,T为平移矩阵;
S5,根据零件表面具有的几何特征及其点云数据占比,为基础部分和几何特征部分分配不同的权重θk,并计算总体变换参数
Figure BDA0003543593320000112
S6,利用总体变换参数v变换待测零件的扫描点云,并重复步骤S3~S5,直至满足迭代终止条件,输出配准结果。
需要说明的是,在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详细描述的部分,可以参见其它实施例的相关描述。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式计算机或者其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包括这些改动和变型在内。

Claims (10)

1.一种用于零件表面质量检测的点云配准方法,其特征在于,包括:
S1,分别获取待测零件扫描点云集合P以及待测零件对应的标准模型Q,并对所述扫描点云集合P进行预处理得到点云集合Pt
S2,将点云集合Pt与标准模型Q进行初始配准;
S3,根据零件表面所具有的几何特征将点云集合Pt划分成基础部分P0与几何特征部分P1~Pn,并构建基础部分和几何特征部分中各点与标准模型Q中各点的匹配点对;其中基础部分为非几何特征部分,n为零件表面所具有的几何特征数量;
S4,利用基础部分和几何特征部分对应的匹配点对,计算协方差矩阵并进行奇异值分解,分别计算得到P0~Pn对应的变换参数
Figure FDA0003543593310000011
R为旋转矩阵,T为平移矩阵;
S5,根据零件表面具有的几何特征及其点云数据占比,为基础部分和几何特征部分分配不同的权重θk,并计算总体变换参数
Figure FDA0003543593310000012
S6,利用总体变换参数v变换待测零件的扫描点云,并重复步骤S3~S5,直至满足迭代终止条件,输出配准结果。
2.根据权利要求1所述的方法,其特征在于,所述匹配点对的构建方法,包括:
构建标准模型Q的Kd-tree结构,对于扫描点云集合Pt中的各点p,利用Kd-tree在标准模型Q中搜索最近点q,构成N个匹配点对。
3.根据权利要求2所述的方法,其特征在于,所述的几何特征的权重根据下式进行分配:
Figure FDA0003543593310000013
式中,N为匹配点对数量,N(Pk)为N个匹配点对中几何特征Pk的数量。
4.根据权利要求2或3所述的方法,其特征在于,步骤S5还包括:
1)判断变换后的扫描点云与标准模型的配准误差是否小于误差函数阈值;
2)判断迭代次数是否达到设定的最大迭代次数;
若1)或2)任一判断结果为是,则表示满足迭代终止条件。
5.根据权利要求4所述的方法,其特征在于,所述配准误差通过下式计算:
Figure FDA0003543593310000021
式中,pi表示第i个匹配点对中待测零件扫描点云中的点,qi表示第i个匹配点对中标准模型Q中的点。
6.根据权利要求1所述的方法,其特征在于,所述预处理,包括:
S601,对于未经预处理的扫描点云数据集P,根据点云密度将连续分布的点聚类,将扫描点云数据集划分成m个点云簇集合{C1,C2,...,Cm};
S602,根据待测零件的几何特征建立模型方程为:M=f(x,y,z,s,D),(x,y,z)为数据集P中点p的坐标值,s为求解模型所需的最少点的数量,该模型存在固有特征与模型表面上的点满足关系D;
S603,选取一个点云簇集合Cl,l=1~m,并从点云簇集合Cl中随机选择一个至少包含s个点的子集S,将子集S中各点的坐标根据模型方程M进行拟合得到模型M*
S604,对于点云簇集合Cl中其余点pt,计算其与模型M*的固有特征关系D*的偏离值评价误差erro(pt)=|D*-D|;若erro(pt)≤maxerro,则保留该点,maxerro为预设的误差阈值;当保留点的数量大于设定值时,则点云簇集合Cl即为待测零件对应的扫描点云Pt;删除除点云簇集合Cl外的其他点云簇集合。
7.根据权利要求6所述的方法,其特征在于,步骤S604还包括:
当保留点的数量大于设定值时,利用所保存的点根据最小二乘法拟合出新的几何模型并替换M*
8.一种用于零件表面质量检测的点云配准装置,其特征在于,包括:
数据获取及预处理模块,分别获取待测零件扫描点云集合P以及待测零件对应的标准模型Q,并对所述扫描点云集合P进行预处理得到点云集合Pt
初始配准模块,将点云集合Pt与标准模型Q进行初始配准;
点云划分及匹配点对构建模块,根据零件表面所具有的几何特征将扫描点云集合Pt划分成基础部分P0与几何特征部分P1~Pn,并构建基础部分和几何特征部分中各点与标准模型Q中各点的匹配点对;其中基础部分为非几何特征部分,n为零件表面所具有的几何特征数量;
第一变换参数计算模块,利用基础部分和几何特征部分对应的匹配点对,计算协方差矩阵并进行奇异值分解,分别计算得到P0~Pn对应的变换参数
Figure FDA0003543593310000031
R为旋转矩阵,T为平移矩阵;
第二变换参数计算模块,根据零件表面具有的几何特征及其点云数据占比,为基础部分和几何特征部分分配不同的权重θk,并计算总体变换参数
Figure FDA0003543593310000032
利用总体变换参数v变换待测零件的扫描点云;
终止条件判断及输出模块,用于判断是否满足迭代终止条件,若满足则输出配准结果。
9.一种电子设备,其特征在于,包括:
存储器,用于存储计算机软件程序;
处理器,用于读取并执行所述计算机软件程序,进而实现权利要求1-8任一项所述的一种用于零件表面质量检测的点云配准方法。
10.一种非暂态计算机可读存储介质,其特征在于,所述存储介质中存储有用于实现权利要求1-7任一项所述的一种用于零件表面质量检测的点云配准方法的计算机软件程序。
CN202210239157.4A 2022-03-11 2022-03-11 一种用于零件表面质量检测的点云配准方法及装置 Active CN114663373B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210239157.4A CN114663373B (zh) 2022-03-11 2022-03-11 一种用于零件表面质量检测的点云配准方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210239157.4A CN114663373B (zh) 2022-03-11 2022-03-11 一种用于零件表面质量检测的点云配准方法及装置

Publications (2)

Publication Number Publication Date
CN114663373A true CN114663373A (zh) 2022-06-24
CN114663373B CN114663373B (zh) 2024-08-20

Family

ID=82028527

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210239157.4A Active CN114663373B (zh) 2022-03-11 2022-03-11 一种用于零件表面质量检测的点云配准方法及装置

Country Status (1)

Country Link
CN (1) CN114663373B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114882087A (zh) * 2022-07-12 2022-08-09 武汉瀚迈科技有限公司 针对具有不完整基本图元的三维扫描点云实时配准方法
CN116109685A (zh) * 2023-04-11 2023-05-12 成都飞机工业(集团)有限责任公司 一种零件点云配准方法、装置、设备及介质
CN116817771A (zh) * 2023-08-28 2023-09-29 南京航空航天大学 一种基于圆柱体素特征的航天零件涂层厚度测量方法
CN118196104A (zh) * 2024-05-17 2024-06-14 深圳大学 钢构件质量智能检测方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20190073744A1 (en) * 2017-09-01 2019-03-07 Tata Consultancy Services Limited Systems and methods for shape constrained 3d point cloud registration
CN110009667A (zh) * 2018-12-19 2019-07-12 南京理工大学 基于罗德里格斯变换的多视点云全局配准方法
CN110220476A (zh) * 2019-06-03 2019-09-10 西安飞机工业(集团)有限责任公司 一种基于三维激光扫描的航空零件自动快速测量方法
CN113267122A (zh) * 2021-05-12 2021-08-17 温州大学瓯江学院 一种基于3d视觉传感器的工业零件尺寸测量方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20190073744A1 (en) * 2017-09-01 2019-03-07 Tata Consultancy Services Limited Systems and methods for shape constrained 3d point cloud registration
CN110009667A (zh) * 2018-12-19 2019-07-12 南京理工大学 基于罗德里格斯变换的多视点云全局配准方法
CN110220476A (zh) * 2019-06-03 2019-09-10 西安飞机工业(集团)有限责任公司 一种基于三维激光扫描的航空零件自动快速测量方法
CN113267122A (zh) * 2021-05-12 2021-08-17 温州大学瓯江学院 一种基于3d视觉传感器的工业零件尺寸测量方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114882087A (zh) * 2022-07-12 2022-08-09 武汉瀚迈科技有限公司 针对具有不完整基本图元的三维扫描点云实时配准方法
CN114882087B (zh) * 2022-07-12 2022-10-04 武汉瀚迈科技有限公司 针对具有不完整基本图元的三维扫描点云实时配准方法
CN116109685A (zh) * 2023-04-11 2023-05-12 成都飞机工业(集团)有限责任公司 一种零件点云配准方法、装置、设备及介质
CN116109685B (zh) * 2023-04-11 2023-08-04 成都飞机工业(集团)有限责任公司 一种零件点云配准方法、装置、设备及介质
CN116817771A (zh) * 2023-08-28 2023-09-29 南京航空航天大学 一种基于圆柱体素特征的航天零件涂层厚度测量方法
CN116817771B (zh) * 2023-08-28 2023-11-17 南京航空航天大学 一种基于圆柱体素特征的航天零件涂层厚度测量方法
CN118196104A (zh) * 2024-05-17 2024-06-14 深圳大学 钢构件质量智能检测方法及系统

Also Published As

Publication number Publication date
CN114663373B (zh) 2024-08-20

Similar Documents

Publication Publication Date Title
CN114170279B (zh) 一种基于激光扫描的点云配准方法
CN114663373A (zh) 一种用于零件表面质量检测的点云配准方法及装置
CN111299815B (zh) 一种用于低灰度橡胶垫视觉检测与激光切割轨迹规划方法
CN110807781B (zh) 一种保留细节与边界特征的点云精简方法
CN109887015A (zh) 一种基于局部曲面特征直方图的点云自动配准方法
CN110599506B (zh) 一种复杂异形曲面机器人三维测量的点云分割方法
CN113628263A (zh) 一种基于局部曲率及其近邻特征的点云配准方法
CN111696210A (zh) 一种基于三维点云数据特征轻量化的点云重构方法及系统
CN108376408A (zh) 一种基于曲率特征的三维点云数据快速加权配准方法
CN108388902B (zh) 结合全局框架点与局部shot特征的复合3d描述子构建方法
CN114972459B (zh) 一种基于低维点云局部特征描述符的点云配准方法
CN108550166B (zh) 一种空间目标图像匹配方法
CN116402866A (zh) 基于点云的零件数字孪生几何建模与误差评定方法及系统
CN116452583B (zh) 点云缺陷检测方法、装置、系统及存储介质
CN115601574A (zh) 一种改进akaze特征的无人机影像匹配方法
CN115861397A (zh) 一种基于改进fpfh-icp的点云配准方法
CN115797418B (zh) 一种基于改进icp的复杂机械零件测量点云配准方法及系统
CN116452644A (zh) 基于特征描述子的三维点云配准方法、装置及存储介质
CN111783722B (zh) 一种激光点云的车道线提取方法和电子设备
CN117132630A (zh) 一种基于二阶空间兼容性度量的点云配准方法
CN117495891B (zh) 点云边缘检测方法、装置和电子设备
CN114898118A (zh) 基于多源点云的输电线路房屋拆迁量自动统计方法及系统
CN117274339A (zh) 一种基于改进的iss-3dsc特征结合icp的点云配准方法
CN109345571B (zh) 一种基于扩展高斯图像的点云配准方法
CN117593204B (zh) 一种基于超体素邻近聚类的点云实例分割方法

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