CN115526922B - 基于分块标签描述子的航空发动机叶片点云数据配准方法 - Google Patents
基于分块标签描述子的航空发动机叶片点云数据配准方法 Download PDFInfo
- Publication number
- CN115526922B CN115526922B CN202211498139.4A CN202211498139A CN115526922B CN 115526922 B CN115526922 B CN 115526922B CN 202211498139 A CN202211498139 A CN 202211498139A CN 115526922 B CN115526922 B CN 115526922B
- Authority
- CN
- China
- Prior art keywords
- cloud data
- point cloud
- point
- blade
- block
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 102
- 230000000903 blocking effect Effects 0.000 title claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims abstract description 53
- 230000009466 transformation Effects 0.000 claims abstract description 53
- 238000005070 sampling Methods 0.000 claims abstract description 13
- 239000013598 vector Substances 0.000 claims description 47
- 238000000638 solvent extraction Methods 0.000 claims description 18
- 238000005259 measurement Methods 0.000 claims description 13
- 230000003287 optical effect Effects 0.000 claims description 9
- 230000001154 acute effect Effects 0.000 claims description 5
- 238000000354 decomposition reaction Methods 0.000 claims description 5
- 230000001131 transforming effect Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 13
- 238000000513 principal component analysis Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000005192 partition Methods 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 238000003754 machining Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000011324 bead Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
- G06T7/344—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods involving models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- 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
-
- 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/20—Special algorithmic details
- G06T2207/20021—Dividing image into blocks, subimages or windows
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开基于分块标签描述子的航空发动机叶片点云数据配准方法,通过采样得到叶片场景点云数据和叶片模型点云数据,并分别计算点云数据的坐标轴,对点云数据进行分块并附标签,得到点云数据分块标签,对点云数据分块上的点进行关键点检测,采用分块标签和检测得到的关键点对叶片场景点云数据和叶片模型点云数据进行分块对应查找和关键点组合配对,得到多个关键点对,从每个点云数据分块中选出一个关键点对,组成点云数据的四点基,并求解变换矩阵,遍历点云数据分块上的多个关键点对,得到多个变换矩阵,结合距离均方差值最小化法得到最优变换矩阵,并对叶片场景点云数据进行坐标变换。该方法能够为精配准提供良好的初始值。
Description
技术领域
本发明属于叶片加工误差检测领域,具体涉及一种基于分块标签描述子的航空发动机叶片场景与叶片模型点云数据配准方法。
背景技术
航空发动机叶片被誉为“工业皇冠上的明珠”,其高压、高温等恶劣的工作环境对发动机叶片加工精度提出了严格的要求,目前针对于发动机叶片的检测技术主要是基于三坐标机的测量,而该测量方式效率低且需要人工干预,随着视觉技术的发展,光学测量技术被引入用于发动机叶片的加工检测。
结合光学测量技术,点云配准技术被发展用于处理获取的点云数据。为了进行加工误差检测,首先要将多视角场景点云数据进行拼接,然后将拼接完成的场景点云数据与模型点云数据进行配准,使两点云处于相同位姿与坐标系下,最后根据两点云的距离关系进行误差估计。在光学测量过程中需要对被测物体粘贴标记点或用夹具平台对叶片底部夹持进行测量,这样的测量方式导致获取的点云数据存在部分缺失,即标记点位置缺失或被夹持部分缺失,这对点云配准方法提出了更高的要求。现有的点云配准方法中常用的有以下几种:一种是ICP配准,该配准方法要求两点云的距离值较小,常被用于精配准过程;一种是特征点结合描述子配准,该配准方法通过提取特征点并结合特征点附近点集信息形成该点描述子,但该方法需要根据不同形状的物体预先进行参数调试,以及针对于自由曲面的叶片配准效果不佳;还有一种针对完整点云配准的方法,是采用PCA(Principal ComponentAnalysis,主成份分析)方法进行粗配准,该方法通过计算两个点云的坐标轴进行坐标轴对齐求解变换矩阵,但用于光学测量设备获取的部分缺失的点云数据,该方法会失效,为了得到好的结果,针对于夹具测量方式,可以通过在共有平面手动选择对应点的方式避免这个问题,但该方法不满足自动检测要求。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供基于分块标签描述子的航空发动机叶片场景点云与叶片模型点云数据配准方法。
基于分块标签描述子的航空发动机叶片点云数据配准方法,包括以下:
S1、通过采样得到点云数据,并计算得出点云数据的坐标轴,点云数据包括叶片场景点云数据和叶片模型点云数据;
S2、根据点云数据的坐标轴对点云数据进行分块,得到点云数据分块,对点云数据分块附标签,得到对应的点云数据分块标签;
S3、对点云数据分块上的点进行关键点检测,采用点云数据分块标签对检测得到的关键点进行描述,得到分块标签描述子,通过分块标签描述子对叶片场景点云数据和叶片模型点云数据进行分块对应查找和关键点组合配对,得到点云数据分块上的多个关键点对;
S4、从每个点云数据分块中选出一个关键点对,组成点云数据的四点基,通过四点基求解变换矩阵,遍历点云数据分块上的多个关键点对,得到多个变换矩阵,结合距离均方差值最小化法得到最优变换矩阵,根据最优变换矩阵对叶片场景点云数据进行坐标变换,得到变换后的叶片模型点云数据坐标系下的叶片场景点云数据。
优选地,S1具体为:
S11、通过光学测量设备采集叶片多视角场景点云数据,对其进行拼接和采样,得到叶片场景点云数据;
S12、对叶片CAD模型数据进行采样,得到叶片模型点云数据;
S13、采用PCA方法分别对叶片场景点云数据和叶片模型点云数据进行坐标轴提取,得到叶片场景点云数据和叶片模型点云数据的坐标轴。
优选地,S13具体为:
S131、定义点云数据;
S132、采用PCA方法对点云数据求解协方差矩阵;
S133、提取协方差矩阵的特征值和特征向量,并依据特征值的大小对相应的特征向量进行排列;
优选地,S2具体为:
S23、对点云数据分块附标签,并根据夹角的角度对点云数据分块标签赋值,得到点云数据分块标签。
优选地,S3中对点云数据分块上的点进行关键点检测,具体为:
S33、选取点云数据分块上距离最大的点作为检测到的点云数据分块上的关键点。
优选地,步骤S32中计算点云数据分块上的点到坐标平面和的距离,包括计算叶片场景点云数据分块上的点到叶片场景点云数据的坐标平面和的距离,和计算叶片模型点云数据分块上的点到叶片模型点云数据的坐标平面和的距离,具体公式为:
优选地,S3中通过分块标签描述子对叶片场景点云数据和叶片模型点云数据进行分块对应查找和关键点组合配对,得到点云数据分块上的多个关键点对,具体为:
S34、通过分块标签描述子对叶片场景点云数据和叶片模型点云数据进行分块对应查找,得到对应的点云数据分块;
S35、将叶片场景点云数据分块上的关键点和叶片模型点云数据分块上的关键点按对应分块进行组合配对,得到点云数据分块上的多个关键点对。
优选地,S4中从每个点云数据分块中选出一个关键点对,组成点云数据的四点基,通过四点基求解变换矩阵,具体为:
S41、从每个点云数据分块中选出一个关键点对,组成点云数据的四点基;
S42、对叶片场景点云数据的四点基和叶片模型点云数据的四点基求解协方差矩阵;
S43、对协方差矩阵采用奇异值分解法求解,得出点云数据的变换矩阵。
优选地,S4中结合距离均方差值最小化法得到最优变换矩阵,具体为:
S44、通过变换矩阵对叶片场景点云数据进行变换,得到叶片模型点云数据坐标系下的叶片场景点云数据;
S45、查找叶片模型点云数据坐标系下的叶片场景点云数据中的点与叶片模型点云数据中对应的最近点;
S46、计算得到的所有叶片场景点云中的点与叶片模型点云数据中对应的最近点的距离并求解均方差值,均方差值最小时对应的变换矩阵为最优变换矩阵。
上述基于分块标签描述子的航空发动机叶片点云数据配准方法,解决了现有配准方法中需要预先调整参数、针对自由曲面现有配准方法的配准效果不佳,以及需要手动选择对应点等问题,能够为精配准提供良好的初始值。
附图说明
图1是本发明一实施例中基于分块标签描述子的航空发动机叶片点云数据配准方法流程图;
图2是本发明一实施例中基于分块标签描述子的航空发动机叶片点云数据配准方法的整体框图;
图3是本发明一实施例中关键点检测过程中的距离求解示意图,其中,(a)为计算叶片场景点云P上的点到坐标平面xoz的距离的示意图,(b)为计算叶片场景点云P上的点到坐标平面yoz的距离的示意图,(c)为计算叶片模型点云Q上的点到坐标平面xoz的距离的示意图,(d)为计算叶片模型点云Q上的点到坐标平面yoz的距离的示意图;
图4为本发明一实施例中采用PCA方法计算得到的坐标轴中y轴和z轴的存在情形,其中(a)为相对于点云本身,坐标轴z轴向上,y轴向右,(b)为相对于点云本身,坐标轴z轴向上,y轴向左,(c)为相对于点云本身,坐标轴z轴向下,y轴向右,(d)为相对于点云本身,坐标轴z轴向下,y轴向左。
具体实施方式
为了使本技术领域的人员更好地理解本发明的技术方案,下面结合附图对本发明作进一步的详细说明。
在一个实施例中,基于分块标签描述子的航空发动机叶片点云数据配准方法,具体包括以下:
S1、通过采样得到点云数据,并计算得出点云数据的坐标轴,点云数据包括叶片场景点云数据和叶片模型点云数据;
S2、根据点云数据的坐标轴对点云数据进行分块,得到点云数据分块,对点云数据分块附标签,得到对应的点云数据分块标签;
S3、对点云数据分块上的点进行关键点检测,采用点云数据分块标签对检测得到的关键点进行描述,得到分块标签描述子,通过分块标签描述子对叶片场景点云数据和叶片模型点云数据进行分块对应查找和关键点组合配对,得到点云数据分块上的多个关键点对;
S4、从每个点云数据分块中选出一个关键点对,组成点云数据的四点基,通过四点基求解变换矩阵,遍历点云数据分块上的多个关键点对,得到多个变换矩阵,结合距离均方差值最小化法得到最优变换矩阵,根据最优变换矩阵对叶片场景点云数据进行坐标变换,得到变换后的叶片模型点云数据坐标系下的叶片场景点云数据。
具体地,参见图1和图2,图1是本发明一实施例中基于分块标签描述子的航空发动机叶片点云数据配准方法的流程图,图2是本发明一实施例中基于分块标签描述子的航空发动机叶片点云数据配准方法的整体框图。
基于分块标签描述子的航空发动机叶片点云数据配准方法,首先通过采样得到点云数据,并计算得出点云数据的坐标轴,该点云数据包括叶片场景点云数据和叶片模型点云数据,通过对点云数据提取坐标轴,得到叶片场景点云数据和叶片模型点云数据的坐标轴;接着根据坐标轴分别对点云数据进行分块并附上标签,得到对应的叶片场景点云数据分块标签和叶片模型点云数据分块标签;然后对点云数据分块上的点进行关键点检测,通过分块标签和检测得到的关键点进行分块对应查找和关键点组合配对,得到点云数据分块上的多个关键点对;在每个叶片场景点云分块和对应的叶片模型点云分块中每次选择一个关键点对,组成每个点云数据的四点基,然后针对叶片场景点云数据的四点基和叶片模型点云数据的四点基求解协方差矩阵和变换矩阵,遍历点云数据分块上的多个关键点对,得到多个变换矩阵,结合距离均方差值最小化法得到最优变换矩阵,根据最优变换矩阵对叶片场景点云数据进行坐标变换,得到变换后的叶片模型点云数据坐标系下的叶片场景点云数据。
在一个实施例中,S1具体为:
S11、通过光学测量设备采集叶片多视角场景点云数据,对其进行拼接和采样,得到叶片场景点云数据;
S12、对叶片CAD模型数据进行采样,得到叶片模型点云数据;
S13、采用PCA方法分别对叶片场景点云数据和叶片模型点云数据进行坐标轴提取,得到叶片场景点云数据和叶片模型点云数据的坐标轴。
具体地,由于利用光学测量设备采集到的叶片多视角场景点云数据的数量达到百万级,需要对其进行拼接并采样,通过采样,减少到万级进行点云配准以提高算法效率,针对叶片CAD模型数据进行采样获取叶片模型点云数据,然后采用PCA方法分别计算叶片场景点云数据和叶片模型点云数据的坐标轴。
在一个实施例中,S13具体为:
S131、定义点云数据;
S132、采用PCA方法对点云数据求解协方差矩阵;
S133、提取协方差矩阵的特征值和特征向量,并依据特征值的大小对相应的特征向量进行排列;
具体地,定义叶片场景点云数据为:
定义叶片模型点云数据为:
将特征向量按照其对应的特征值从大到小进行相应的排列(特征值大,对应的特征向量表示的点云数据的对应方向上点的坐标值方差大,分布特征显著),并分别作为叶片场景点云数据和叶片模型点云数据的坐标轴,其中方差最大的特征向量对应于坐标轴的轴,方差最小的特征向量对应于坐标轴的轴,剩下的一个特征向量对应坐标轴的轴,点云数据的质心作为坐标轴的原点。
以计算叶片场景点云的坐标轴为例,采用PCA方法对叶片场景点云求解协方差矩阵,计算该矩阵的特征值,满足,对应的特征向量,表示按照特征值从大到小进行相应排列的三个特征向量,并定义分别表示叶片场景点云数据的轴,轴和轴,点云的质心为坐标轴为原点。叶片模型点云数据的坐标轴求解过程同上,此处不再赘述。
在一个实施例中,S2具体为:
S23、对点云数据分块附标签,并根据夹角的角度对点云数据分块标签赋值,得到点云数据分块标签。
具体地,首先通过点云数据中的坐标轴中的轴和轴对点云数据分块,得到点云数据分块,然后将点云数据的质心与点云数据中的每个点连线形成向量,即针对叶片场景点云,有个向量,针对叶片模型点云,有个向量 ,计算向量和分别与各自的轴和轴形成的夹角的角度;最后对每个点云数据分块均附标签和,根据得到的夹角的角度,对标签赋值,具体赋值方式如下:
若向量与轴夹角成锐角,则对应的点云数据分块的标签被赋值“true”,若向量与轴夹角成钝角,则对应的点云数据分块的标签被赋值“false”,若向量与轴夹角成锐角,则对应的点云数据分块的标签被赋值“true”,若向量与轴夹角成钝角,则对应的点云数据分块的标签被赋值“false”。
在一个实施例中,S3中对点云数据分块上的点进行关键点检测,具体为:
S33、选取点云数据分块上距离最大的点作为检测到的点云数据分块上的关键点。
在一个实施例中,步骤S32中计算点云数据分块上的点到坐标平面和的距离,包括计算叶片场景点云数据分块上的点到叶片场景点云数据的坐标平面和的距离,和计算叶片模型点云数据分块上的点到叶片模型点云数据的坐标平面和的距离,具体公式为:
具体地,参见图3,图3是本发明一实施例中关键点检测过程中的距离求解示意图。
由于光学测量设备所测叶片的叶身在测量过程中被粘贴标记点,故采集的点云数据在这些标记点处表现为孔洞,或被夹具夹持部分的点云数据因被遮挡而缺失,场景点云与模型点云的坐标轴计算存在偏差,故针对于叶片上因粘贴标记点或遮挡导致点云数据缺失从而使直接使用PCA方法配准效果不佳的情况,采用关键点检测的方法进行处理,具体按以下方法进行:
针对叶片场景点云划分的每个分块,提取每个分块上:
针对叶片模型点云划分的每个分块,提取每个分块上:
这样可在叶片场景点云数据和叶片模型点云数据中分别提取8个关键点。
在一个实施例中,S3中通过分块标签描述子对叶片场景点云数据和叶片模型点云数据进行分块对应查找和关键点组合配对,得到点云数据分块上的多个关键点对,具体为:
S34、通过分块标签描述子对叶片场景点云数据和叶片模型点云数据进行分块对应查找,得到对应的点云数据分块;
S35、将叶片场景点云数据分块上的关键点和叶片模型点云数据分块上的关键点按对应分块进行组合配对,得到点云数据分块上的多个关键点对。
具体地,根据关键点的位置特征在点云数据中进行对应查找,首先依据分块标签描述子进行分块对应查找,得到对应的点云数据分块,再在每个对应的点云数据分块中进行关键点对应查找,即将场景点云分块上的关键点和模型点云分块上的关键点任意配对组合,得到多个关键点对。
表1中的“&&”在实例中表示“与”操作,“!”在实例中表示“取反”操作。
这样一来,每个点云数据分块上的关键点都有4种对应情形,每次在每个点云分块中提取一个关键点对应的一种情形,需要遍历所有可能的情形。在每个对应分块中遍历每个关键点对应情形,存在4个关键点对,即在4个对应分块中,每个分块提取一个关键点对,存在种对应情形。
表1
在一个实施例中,S4中从每个点云数据分块中选出一个关键点对,组成点云数据的四点基,通过四点基求解变换矩阵,具体为:
S41、从每个点云数据分块中选出一个关键点对,组成点云数据的四点基;
S42、对叶片场景点云数据的四点基和叶片模型点云数据的四点基求解协方差矩阵;
S43、对协方差矩阵采用奇异值分解法求解,得出点云数据的变换矩阵。
具体地,在每个点云数据分块中每次选取一个关键点对,由于点云数据被分成了4块,因此叶片场景点云数据和叶片模型点云数据每次可以得到4个关键点对,这4个关键点对分别构成对应点云数据的4点基,即每次可以得到叶片场景点云数据的四点基和叶片模型点云数据的四点基,对这两个点云数据的四点基求协方差矩阵,并对协方差矩阵采用奇异值分解法求解点云数据的变换矩阵。
在一个实施例中,S4中结合距离均方差值最小化法得到最优变换矩阵,具体为:
S44、通过变换矩阵对叶片场景点云数据进行变换,得到叶片模型点云数据坐标系下的叶片场景点云数据;
S45、查找叶片模型点云数据坐标系下的叶片场景点云数据中的点与叶片模型点云数据中对应的最近点;
S46、计算得到的所有叶片场景点云中的点与叶片模型点云数据中对应的最近点的距离并求解均方差值,均方差值最小时对应的变换矩阵为最优变换矩阵。
具体地,通过变换矩阵对叶片场景点云数据进行变换,得到叶片模型点云数据坐标系下的叶片场景点云坐标,使叶片场景点云数据和叶片模型点云数据在相同的坐标系下,然后针对叶片场景点云数据中的每个点在叶片模型点云数据中找到距离它最近的那个点,计算这两个点之间的距离,并对所有叶片场景点云中的点对应的叶片模型点云中的最近点的距离求解均方差,遍历每一组4点基,均方差最小时对应的变换矩阵为最优变换矩阵。
在求解均方差的过程中,由于场景点云存在部分缺失,因此在实例中以叶片场景点云为基准在叶片模型点云中查找对应点,以避免缺失部分对结果的影响。另外,在求解均方差的过程中,为提高计算效率,以间隔比例(例如1%)在叶片场景点云中取点,并查找这个点在叶片模型点云中的对应点,作为计算均方差的有效点。比如一个叶片场景点云数据中有100万个点,如果要计算每个点对应的距离差会很耗时,所以在这里每次间隔1%个点,也就是说每隔1万个点查找一次它在模型点云中的对应点并且计算它们之间的距离,这样的话计算量就会大大减少。
上述基于分块标签描述子的航空发动机叶片点云数据配准方法,首先通过采样得到点云数据,点云数据包括叶片场景点云数据和叶片模型点云数据,分别提取点云数据的坐标轴;然后根据点云数据的坐标轴对点云数据进行分块,得到点云数据分块,对点云数据分块附上标签并对标签赋值,得到对应的点云数据分块标签;接着对点云数据分块上的点进行关键点检测,通过分块标签和检测得到的关键点进行分块对应查找和关键点组合配对,得到每个分块上的多个关键点对;最后从每个点云数据分块上的多个关键点对中每次选出一个关键点对,组成点云数据的四点基,采用奇异值分解法求解四点基的变换矩阵,遍历每个点云分块上的所有关键点对,得到多个变换矩阵,结合距离均方差值最小化法从多个点变换矩阵中得到最优变换矩阵,根据最优变换矩阵对叶片场景点云数据进行坐标变换,得到变换后的叶片模型点云数据坐标系下的叶片场景点云数据。该方法解决了现有配准方法中需要预先调整参数、针对自由曲面现有配准方法的配准效果不佳,以及需要手动选择对应点等问题,能够为精配准提供良好的初始值。
以上对本发明所提供的基于分块标签描述子的航空发动机叶片点云数据配准方法进行了详细介绍。本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的核心思想。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以对本发明进行若干改进和修饰,这些改进和修饰也落入本发明权利要求的保护范围内。
Claims (9)
1.基于分块标签描述子的航空发动机叶片点云数据配准方法,其特征在于,所述方法包括以下:
S1、通过采样得到点云数据,并计算得出所述点云数据的坐标轴,所述点云数据包括叶片场景点云数据和叶片模型点云数据;
S2、根据所述点云数据的坐标轴对所述点云数据进行分块,得到点云数据分块,对所述点云数据分块附标签,得到对应的点云数据分块标签;所述S2具体为:
S23、对所述点云数据分块附标签,并根据所述夹角的角度对所述点云数据分块标签赋值,得到点云数据分块标签;
S3、对所述点云数据分块上的点进行关键点检测,采用所述点云数据分块标签对检测得到的关键点进行描述,得到分块标签描述子,通过所述分块标签描述子对所述叶片场景点云数据和叶片模型点云数据进行分块对应查找和关键点组合配对,得到点云数据分块上的多个关键点对;
S4、从每个所述点云数据分块中选出一个关键点对,组成所述点云数据的四点基,通过四点基求解变换矩阵,遍历所述点云数据分块上的多个关键点对,得到多个变换矩阵,结合距离均方差值最小化法得到最优变换矩阵,根据所述最优变换矩阵对所述叶片场景点云数据进行坐标变换,得到变换后的叶片模型点云数据坐标系下的叶片场景点云数据。
2.如权利要求1所述的基于分块标签描述子的航空发动机叶片点云数据配准方法,其特征在于,所述S1具体为:
S11、通过光学测量设备采集叶片多视角场景点云数据,对其进行拼接和采样,得到叶片场景点云数据;
S12、对叶片CAD模型数据进行采样,得到叶片模型点云数据;
S13、采用PCA方法分别对所述叶片场景点云数据和所述叶片模型点云数据进行坐标轴提取,得到所述叶片场景点云数据和所述叶片模型点云数据的坐标轴。
6.如权利要求5所述的基于分块标签描述子的航空发动机叶片点云数据配准方法,其特征在于,所述步骤S32中计算所述点云数据分块上的点到所述坐标平面和的距离,包括计算叶片场景点云数据分块上的点到叶片场景点云数据的坐标平面和的距离,和计算叶片模型点云数据分块上的点到叶片模型点云数据的坐标平面和的距离,具体公式为:
7.如权利要求6所述的基于分块标签描述子的航空发动机叶片点云数据配准方法,其特征在于,所述S3中通过所述分块标签描述子对所述叶片场景点云数据和叶片模型点云数据进行分块对应查找和关键点组合配对,得到点云数据分块上的多个关键点对,具体为:
S34、通过所述分块标签描述子对所述叶片场景点云数据和叶片模型点云数据进行分块对应查找,得到对应的点云数据分块;
S35、将叶片场景点云数据分块上的关键点和叶片模型点云数据分块上的关键点按对应分块进行组合配对,得到点云数据分块上的多个关键点对。
8.如权利要求7所述的基于分块标签描述子的航空发动机叶片点云数据配准方法,其特征在于,所述S4中从每个所述点云数据分块中选出一个关键点对,组成所述点云数据的四点基,通过四点基求解变换矩阵,具体为:
S41、从每个所述点云数据分块中选出一个关键点对,组成所述点云数据的四点基;
S42、对叶片场景点云数据的四点基和叶片模型点云数据的四点基求解协方差矩阵;
S43、对所述协方差矩阵采用奇异值分解法求解,得出点云数据的变换矩阵。
9.如权利要求8所述的基于分块标签描述子的航空发动机叶片点云数据配准方法,其特征在于,所述S4中结合距离均方差值最小化法得到最优变换矩阵,具体为:
S44、通过所述变换矩阵对所述叶片场景点云数据进行变换,得到叶片模型点云数据坐标系下的叶片场景点云数据;
S45、查找所述叶片模型点云数据坐标系下的叶片场景点云数据中的点与所述叶片模型点云数据中对应的最近点;
S46、计算得到的所有叶片场景点云中的点与叶片模型点云数据中对应的最近点的距离并求解均方差值,均方差值最小时对应的变换矩阵为最优变换矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211498139.4A CN115526922B (zh) | 2022-11-28 | 2022-11-28 | 基于分块标签描述子的航空发动机叶片点云数据配准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211498139.4A CN115526922B (zh) | 2022-11-28 | 2022-11-28 | 基于分块标签描述子的航空发动机叶片点云数据配准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115526922A CN115526922A (zh) | 2022-12-27 |
CN115526922B true CN115526922B (zh) | 2023-04-28 |
Family
ID=84704667
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211498139.4A Active CN115526922B (zh) | 2022-11-28 | 2022-11-28 | 基于分块标签描述子的航空发动机叶片点云数据配准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115526922B (zh) |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105118059A (zh) * | 2015-08-19 | 2015-12-02 | 哈尔滨工程大学 | 一种多尺度坐标轴夹角特征的点云快速配准方法 |
GB2550567A (en) * | 2016-05-20 | 2017-11-29 | Nokia Technologies Oy | Point Cloud Matching Method |
CN109147030B (zh) * | 2018-07-05 | 2020-06-30 | 厦门大学 | 基于线特征的室内外场景联合建模方法 |
CN113643270B (zh) * | 2021-08-24 | 2024-04-26 | 凌云光技术股份有限公司 | 一种基于点云数据的图像配准方法及装置 |
-
2022
- 2022-11-28 CN CN202211498139.4A patent/CN115526922B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN115526922A (zh) | 2022-12-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110426051B (zh) | 一种车道线绘制方法、装置及存储介质 | |
US20140267614A1 (en) | 2D/3D Localization and Pose Estimation of Harness Cables Using A Configurable Structure Representation for Robot Operations | |
CN108876852B (zh) | 一种基于3d视觉的在线实时物体识别定位方法 | |
Jiang et al. | Dual-camera-based method for identification and location of scattered self-plugging rivets for robot grasping | |
CN110930442B (zh) | 基于标定块的机器人手眼标定中关键点位置确定方法与装置 | |
CN111598172B (zh) | 基于异构深度网络融合的动态目标抓取姿态快速检测方法 | |
CN110146017B (zh) | 工业机器人重复定位精度测量方法 | |
CN110084842A (zh) | 一种机器人云台伺服二次对准方法及装置 | |
CN109241901A (zh) | 一种对具有孔洞的三维点云的检测与识别方法 | |
CN103729631A (zh) | 一种基于视觉的连接器表面特征自动识别方法 | |
Pi et al. | Stereo visual SLAM system in underwater environment | |
CN109255815B (zh) | 一种基于零阶球谐函数的物体检测与识别方法 | |
CN112561989B (zh) | 一种建造场景下吊装对象的识别方法 | |
CN115526922B (zh) | 基于分块标签描述子的航空发动机叶片点云数据配准方法 | |
CN111340884A (zh) | 一种双目异构相机与rfid的双重目标定位与身份辨识方法 | |
CN111475016A (zh) | 一种基于计算机视觉的装配过程几何参数自适应测量系统及方法 | |
IL123566A (en) | Detection of prominent contours in a pair of stereoscopic figures | |
JPWO2018131163A1 (ja) | 情報処理装置、データベース生成装置、方法、プログラム、及び記憶媒体 | |
CN110517323A (zh) | 基于机械手单相机多目视觉的三维定位系统及方法 | |
Hashimoto et al. | Current status and future trends on robot vision technology | |
CN117994484A (zh) | 一种大型建筑预制构件三维图像拼接装置及方法 | |
Tunca et al. | Determining the pointer positions of aircraft analog indicators using deep learning | |
CN113192134B (zh) | 一种航空发动机卡箍定位及偏移检测方法 | |
Jiang et al. | A multi-view structured light measurement method based on pose estimation using deep learning | |
CN113920049B (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 |