CN111563920A - 一种全局优化与多约束条件迭代的3d彩色点云配准方法 - Google Patents

一种全局优化与多约束条件迭代的3d彩色点云配准方法 Download PDF

Info

Publication number
CN111563920A
CN111563920A CN202010295862.7A CN202010295862A CN111563920A CN 111563920 A CN111563920 A CN 111563920A CN 202010295862 A CN202010295862 A CN 202010295862A CN 111563920 A CN111563920 A CN 111563920A
Authority
CN
China
Prior art keywords
point cloud
formula
point
registration
calculating
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
CN202010295862.7A
Other languages
English (en)
Other versions
CN111563920B (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.)
Xian Polytechnic University
Original Assignee
Xian Polytechnic University
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 Xian Polytechnic University filed Critical Xian Polytechnic University
Priority to CN202010295862.7A priority Critical patent/CN111563920B/zh
Publication of CN111563920A publication Critical patent/CN111563920A/zh
Application granted granted Critical
Publication of CN111563920B publication Critical patent/CN111563920B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/337—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods involving reference images or patches
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/23—Clustering techniques
    • G06F18/232—Non-hierarchical techniques
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/24—Classification techniques
    • G06F18/243—Classification techniques relating to the number of classes
    • G06F18/2433—Single-class perspective, e.g. one-against-all classification; Novelty detection; Outlier detection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • 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/10024—Color image
    • 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/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
    • 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
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30—Computing systems specially adapted for manufacturing

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种全局优化与多约束条件迭代的3D彩色点云配准方法。首先，采用基于标记的多状态离群算法有效剔除噪声、离群点及离群点簇后；其次，利用PSO粒子群优化算法对其进行初步粗配准处理；再次，使用蚁群优化算法对其进行全局优化处理；最后，采用全局多约束条件迭代最近点精确配准算法对其进行精配准处理。在粗配准过程中，利用蚁群优化算法对其结果进行全局优化，不仅减少粗配准的错误配对，提高了粗配准精度，还为后续的精配准处理初始化一个精度较高的初始值，进而缩短整个配准的时间；在精配准过程中，采用对应点之间欧氏距离和曲率双约束条件ICP迭代精配准处理，使得配准的精度得到一定的提高。

Description

一种全局优化与多约束条件迭代的3D彩色点云配准方法
技术领域
本发明属于3D点云数据处理技术领域，具体涉及一种全局优化与多约束条件迭代的3D彩色点云配准方法。
背景技术
随着人工智能的发展，图像配准技术逐渐从二维图像过度到三维图像，而3D点云作为三维图像的典型代表之一，已逐渐被广泛应用。3D点云配准是首先利用计算机对3D点云进行预处理，然后采用合适的算法对于处理后的数据进行粗配准，最后再选取合适的方法对粗配准后的结果进行精配准。其点云数据配准的整个过程中，滤波器的选取、粗配准准则及精配准约束条件的选取直接决定点云数据最终的配准结果。
近些年，3D点云数据配准技术应用越来越广泛，引起了一些学者的关注。这种配准方法的基本思想是：首先先去一定的滤波器对原点云数据进行预处理，其次采取一定的算法对预处理结果进行粗配准，最后对其结果进行精配准。目前，点云数据的配准方法虽然由很多，可以较好的进项点云数据的配准处理，但仍旧存在一些缺陷：第一，滤波器选取得好坏直接影响噪声、离群点和离群点簇剔除的结果，从而决定其粗配准的精度高低；第二，在粗配准的过程中，往往会陷入局部最优解，使得粗配准的误差较大；第三，在精配准的过程中，约束条件的单一化选择导致精配准后的误差较大，配准精度较低。
发明内容
本发明的目的是提供一种全局优化与多约束条件迭代的3D彩色点云配准方法，解决了现有3D彩色点云配准方法中配准精度较低，配准时间较长的问题。
本发明所采用的技术方案是，一种全局优化与多约束条件迭代的3D彩色点云配准方法，具体按照以下步骤实施：
步骤1，采用基于标记的多状态离群点剔除原点云数据中的噪声、离群点及离群点簇；
步骤2，经过步骤1后，采用PSO粒子群优化算法对滤除后源点云数据和目标点云数据进行粗配准处理；
步骤3，经过步骤2后，采用蚁群优化算法对其进行全局优化处理，得到更加精准的粗配准结果；
步骤4，经过步骤3后，采用改进的多约束条件迭代最近点精确配准算法对其进行精配准处理。
本发明的特点还在于，
步骤1中，具体按照以下步骤实施：
步骤1.1，用k-d树进行点云的邻域搜索；
步骤1.2，经步骤1.1后，用移动最小二乘法局部表面拟合方法对点云法向量进行计算；
步骤1.3，经步骤1.2后，对不同点云数据设置不同的阈值，以点法线差值的正交分量的比值与设置阈值进行比较，以作为孤立点的判断条件，标记离群点；
步骤1.4，经步骤1.3后，以改进的DBSCAN聚类方法对标记离群点进行聚类；
步骤1.5，经步骤1.4后，统计每个聚类中标记离群点的数量占该聚类中点云总数的比值，当该比值大于设定的比例阈值时，判定并删除离群点簇；
步骤1.6，重复步骤1.5，直到删除所有的离群点簇，同时删除剩余的标记离群点；得到滤除后源点云数据和目标点云数据。
步骤2中，具体按照以下步骤实施：
步骤2.1，初始化PSO粒子群算法的迭代次数M＝80，随机生成平移变量tx,ty,tz以及旋转变量rx,ry,rz，初始化学习因子c1＝2、c2＝2；
步骤2.2，经过步骤2.1后，采用单位四元数法求出源点云P和目标点云Q的旋转变换矩阵R和平移矩阵T，并计算目标点云Q旋转、平移后的矩阵Q′；具体按照以下步骤实施：
步骤2.21，计算源点云P与目标点云Q的重心，如式(1)所示；
Figure BDA0002452166470000031
式(1)中，pi和qi是源点云与目标点云的点，i为其下标，表示第i个点，n为每个点云的总数，μp和μq分别为源点云P与目标点云Q的重心；
步骤2.22，利用公式(2)计算其协方差矩阵M；
Figure BDA0002452166470000032
步骤2.23，根据协方差矩阵M，利用公式(3)构造4阶矩阵S；
Figure BDA0002452166470000041
步骤2.24，计算矩阵S的特征值和特征向量，并将最大特征值对应的特征向量作为所求的旋转变量R；
步骤2.25，利用公式(4)计算平移矩阵T；
T＝μP-RμQ(4)；
步骤2.26，经过步骤2.24和2.25后，获得旋转矩阵R和平移矩阵T，并利用R和T及公式(5)计算目标点云旋转、平移后的矩阵Q′；
Q′＝RQ+T (5)；
步骤2.3，利用公式(6)计算适应度函数fitness；
fitness＝|RT(Q-Q′)|，i＝1,2,3,...,m(6)；
式(6)中，变换矩阵RT包含6个参数，平移变量tx,ty,tz以及旋转变量rx,ry,rz，Q’为变换后的点云集；
步骤2.4，利用公式(7)计算源点云和目标点云对应点对的局部最小距离fpbest，表示粒子当前最佳位置；
f(pi，qi)＝|p′i-q′i| (7)；
式(7)中，p′i和q′i表示待配准点对pi和qi的k邻域的重心；
步骤2.5，令fp＝fpbest，fitnessbest＝fitness，RT＝RT(R0,T0)，fp表示粒子历史最佳位置，fitnessbest表示所有粒子的最佳位置函数，RT表示最佳变换矩阵，RT(R0,T0)表示初始化的变换矩阵，实现源点云和目标点云粒子群粗配准；具体按照以下步骤实施：
步骤2.51，设置r1i＝rand()，r2i＝rand()，其中r1i和r2i都是伪随机数，rand()是随机函数，在(0,1)的均匀分布中生成随机数；
步骤2.52，以空间矩阵RT代替速度V，并利用公式(8)更新RT，对目标点云Q进行空间矩阵变换；
Figure BDA0002452166470000051
式(8)中，D为维数，V为粒子i的速度，pi为粒子i在一次迭代中的最佳位置，c1和c2为加速器，设置为2；r1和r2是伪随机数，设置为0-1；
步骤2.53，利用公式(6)计算此次变换的适应度函数，并利用公式(7)计算源点云P与目标点云Q的各点对的局部距离函数，并确定其中最小值fpbest；
步骤2.54，比较fp和fpbest，如果fpbest小于fp，则继续比较fitnessbest和fitness，如果fitness小于fitnessbest，则fitnessbest＝fitness，fp＝fpbest；否则，fitnessbest不变；如果fpbest大于fp，则fp不变；
步骤2.55，确定满足实验数据的最佳迭代次数M，如果迭代M次后，标记此时的最近点对集合；否则，重复步骤2.51至步骤2.55。
步骤3中，具体按照以下步骤实施：
步骤3.1，对目标点集中的点与点之间赋予初始信息素值C＝2，并在目标点云上随机选择n个点作为“蚂蚁”的初始位置；
步骤3.2，利用公式(9)计算每个“蚂蚁”点在其邻域内的状态转移概率并进行状态转移；
Figure BDA0002452166470000052
式(9)中，pij k表示编号为k的蚂蚁从i位置向j位置的“转移”的概率；τij(t)表示t时刻位置i与位置j的信息素值；ηij通常设为i位置和j位置之间距离的倒数；Jk(i)表示i位置上的蚂蚁k未访问的邻近位置点的集合；α,β为启发式因子，其中α＝0.9,β＝1.06，分别表示确定信息素浓度τij和距离间的评价参数；
步骤3.3，对所有“蚂蚁”点搜索到一条路径后，遍历目标点集中所有点；
步骤3.4，利用公式(10)进行信息素迭代计算；
Figure BDA0002452166470000061
式(10)中，ρ为信息素蒸发参数，ρ＝0.1，0＜ρ＜1；Δτk ij表示蚂蚁k在i位置和j位置之间信息素的浓度，并利用公式(11)计算获得；
Figure BDA0002452166470000062
式(11)中，
Figure BDA0002452166470000063
和
Figure BDA0002452166470000064
是选择的对应点对的法向量，且两法向量相似度高时，
Figure BDA0002452166470000065
k表示第k个粒子的搜索过程；λ＝8；
步骤3.5，确定最终比较的阈值，若路径存在更新的信息素值τij大于预先设定的阈值，则搜索结束，该路径(点集合)与源点云组成的点对作为配准的点对；否则重复步骤3.2至步骤3.4；
步骤3.6，利用单位四元数法求解最终选择的点对之间的空间变换矩阵，完成粗配准。
步骤4中，在多约束条件的ICP精配准算法中，对源点云P与目标点云Q预先设置最大迭代次数M＝200，目标函数阈值δ＝10-6m，利用公式(12)表示迭代配准中的目标函数；
Figure BDA0002452166470000071
其中，NS为点云数量，R和T为两点云配准中的旋转矩阵和平移矩阵，pi和qi分别为源点云与目标点云中的点；
具体按照以下步骤实施：
步骤4.1，搜索源点云P和目标点云Q，确定对应点对；具体按照以下步骤实施：
步骤4.11，搜索源点云P和目标点云Q，使得欧式距离‖pi-qi‖＝min最小的点对为初始对应点对；
步骤4.12，利用公式(13)计算所有点对中每点的曲率，并将这些点对存储到集合TS{(pi,qi),i＝1,2,3,...,n′}中，该集合包含初始点对及这些点对的曲率信息；
Figure BDA0002452166470000072
式(13)中，k1(pi)、k2(pi)分别是点曲率的最大值和最小值，也称为点的主曲率；
步骤4.13，利用公式(22)计算每对点对的曲率相似度ψ，并判断误差是否在取值区间之内，若在，保留该点对作为待配准的对应点对；若不在，搜索目标点qi的k邻域，并利用公式(13)计算邻域内各点的曲率信息；
ψ＝|S(pi)-S(qi)|(22)；
步骤4.14，计算点对中点pi和qi邻域各点的曲率相似度及欧式距离，将满足欧式距离最小且曲率相似度符合误差区间[0,0.08]的邻域点作为新的qi；若找不到满足以上两条约束的邻域点，则将该点对从集合TS中删除；
步骤4.15，将最终确定的点对在源点云与目标点云集合进行标记；
步骤4.2，目标函数最小化，并利用奇异值分解法计算此次迭代过程中的旋转矩阵Rk、平移矩阵Tk，利用公式(23)计算此次误差函数d；
Figure BDA0002452166470000081
其中，计算旋转矩阵Rk、平移矩阵Tk的具体步骤为：
步骤4.21，利用公式(24)表示空间变换矩阵函数；
Figure BDA0002452166470000082
步骤4.22，利用公式(25)计算源点云与目标点云的重心坐标，并利用公式(26)确定去重心化后的两个点集p′和q′，则利用公式(27)表示新的变换矩阵函数；
Figure BDA0002452166470000083
p′＝p-μP,q′＝q-μQ(26)；
Figure BDA0002452166470000084
步骤4.23，对目标函数的总误差进行化简，化简后的函数由公式(28)表示，并对其进一步表示为公式(29)；
Figure BDA0002452166470000085
Figure BDA0002452166470000086
步骤4.24，要使目标函数变为最小，即使得Tranc(RH)最小化；线性代数中，M是一个对称矩阵且B是正交矩阵，则有Tranc(M)≥Tranc(BM)，对H进行奇异值分解，得到公式(30)；
H＝U∑VT(30)
式(30)中，∑为3阶非负元素对角矩阵；
步骤4.25，令旋转矩阵R＝VUT时，RH为对称正定矩阵；此时的R为最优旋转矩阵，并计算平移矩阵T；
步骤4.3，利用公式(31)更新目标点集；
Q＝RkP+Tk(31)；
步骤4.4，再次计算新的对应点对使之满足欧式距离最小，并判断点对是否被标记；若标记，则直接根据标记索引，在集合TS中找到该点对；否则按照步骤4.12计算点对的曲率信息并对这些点对进行归类或者删除；
步骤4.5，重复迭代并计算两次误差的差值|d-d’|，若满足|d-d’|＜δ，则算法结束；否则继续迭代，直至达到最大迭代次数。
任意一点p的主曲率，具体实现步骤为：
设二次平面的一般方程，如公式(14)所示；
z＝f(x,y)＝a0x2+a1y2+a2x+a2y+a4(14)；
在目标点pi的k邻域内，根据最小二乘原理利用公式(15)该曲面方程的最小值；
Figure BDA0002452166470000091
由公式(14)和公式(15)求出该曲面的拟合曲面方程，拟合曲面方程由公式(16)表示；
r(x,y)＝(x,y,a0+a1x+a2y+a3x2+a4y2)(16)；
计算拟合曲面r(x,y)的偏导rx，ry，rxx，ryy，rxy，并利用公式(17)计算曲面的单位法矢量；
Figure BDA0002452166470000101
利用公式(18)和(19)分别计算平均曲率和高斯曲率；
Figure BDA0002452166470000102
Figure BDA0002452166470000103
其中，E＝rx·rx,F＝rx·ry,G＝ry·ry,
Figure BDA0002452166470000104
Figure BDA0002452166470000105
利用公式(20)和公式(21)计算点p的主曲率；
Figure BDA0002452166470000106
Figure BDA0002452166470000107
本发明的有益效果是，
该方法采用基于标记的多状态离群算法有效剔除噪声、离群点及离群点簇后，首先利用PSO粒子群优化算法对其进行初步粗配准处理；其次，使用蚁群优化算法对其进行全局优化处理；最后，采用全局多约束条件迭代最近点精确配准算法对其进行精配准处理。在粗配准过程中，由于PSO粒子群优化算法的处理过程收敛速度快，容易陷入局部最优解，导致错误配准，再利用蚁群优化算法对其结果进行全局优化，不仅减少粗配准的错误配对，还为后续的精配准处理初始化一个精度较高的初始值，进而缩短整个配准的时间；在精配准过程中，采用了对应点之间欧氏距离和曲率双约束条件ICP迭代精配准处理，使得配准的精度得到一定的提高。
附图说明
图1是本发明一种全局优化与多约束条件迭代的3D彩色点云配准方法的流程图；
图2a是人工加噪的斯坦福bunny3D点云数据图；
图2b是人工加噪的斯坦福bunny3D点云数据滤波后的结果图；
图2c是使用kinect相机采集得点云数据图；
图2d是使用kinect相机采集得点云数据滤波后的结果图；
图3a是本实施例中待配准源点云和目标点云数据图；
图3b是现有基于SAC-IA粗配准(正面)结果图；
图3c是现有基于PSO粒子群的粗配准(正面)结果图；
图3d是本实施例中进行粗配准(正面)结果图；
图3e是现有基于SAC-IA粗配准(侧面)结果图；
图3f是现有基于PSO粒子群的粗配准(侧面)结果图；
图3g是本实施例中进行粗配准(侧面)结果图；
图4a是现有3DNDT算法精配准(正面)结果图；
图4b是现有ICP算法精配准(正面)结果图；
图4c是本实施例中进行精配准(正面)结果图；
图4d是现有3DNDT算法精配准(侧面)结果图；
图4e是现有ICP算法精配准(侧面)结果图；
图4f是本实施例中进行精配准(侧面)结果图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明一种全局优化与多约束条件迭代的3D彩色点云配准方法，如图1所示，具体按照以下步骤实施：
步骤1，采用基于标记的多状态离群点剔除原点云数据中的噪声、离群点及离群点簇；具体按照以下步骤实施：
步骤1.1，用k-d树进行点云的邻域搜索；
步骤1.2，经步骤1.1后，用移动最小二乘法局部表面拟合方法对点云法向量进行计算；
步骤1.3，经步骤1.2后，对不同点云数据设置不同的阈值，以点法线差值的正交分量的比值与设置阈值进行比较，以作为孤立点的判断条件，标记离群点；本发明的两个实例的阈值是bunny模型取值0.25，椅子模型取0.3；
步骤1.4，经步骤1.3后，以现有改进的DBSCAN聚类方法对标记离群点进行聚类；
步骤1.5，经步骤1.4后，统计每个聚类中标记离群点的数量占该聚类中点云总数的比值，当该比值大于经多次实验调整设定的比例阈值时，判定并删除离群点簇；本发明的两个实例的阈值是bunny模型取值0.2，椅子模型取0.33；
步骤1.6，重复步骤1.5，直到删除所有的离群点簇，同时删除剩余的标记离群点；得到滤除后源点云数据和目标点云数据；
图2a是人工加噪的斯坦福bunny3D点云数据，图2b是滤波后的结果，图2c是使用kinect相机采集得点云数据，图2d是滤波后的结果，由处理结果可以看出，该滤波方法不仅能够滤除人工噪音，而且还能有效滤除相机采集得3D点云数据。
步骤2，经过步骤1后，采用PSO粒子群优化算法对滤除后源点云数据和目标点云数据进行粗配准处理；具体按照以下步骤实施：
步骤2.1，初始化PSO粒子群算法的迭代次数M＝80，随机生成平移变量tx,ty,tz以及旋转变量rx,ry,rz，初始化学习因子c1＝2、c2＝2；
步骤2.2，经过步骤2.1后，采用单位四元数法求出源点云P和目标点云Q的旋转变换矩阵R和平移矩阵T，并计算目标点云Q旋转、平移后的矩阵Q′；
具体按照以下步骤实施：
步骤2.21，计算每个源点云P与目标点云Q的重心，如式(1)所示；
Figure BDA0002452166470000131
式(1)中，pi和qi是源点云与目标点云的点，i为其下标，表示第i个点，n为每个点云的总数，μp和μq分别为源点云P与目标点云Q的重心；
步骤2.22，利用公式(2)计算每个源点云P和目标点云Q的协方差矩阵M；
Figure BDA0002452166470000132
步骤2.23，根据协方差矩阵M，利用公式(3)构造4阶矩阵S；
Figure BDA0002452166470000141
步骤2.24，计算矩阵S的特征值和特征向量，并将最大特征值对应的特征向量作为所求的旋转变量R；
步骤2.25，利用公式(4)计算源点云P和目标点云Q的平移矩阵T；
T＝μP-RμQ(4)；
步骤2.26，经过步骤2.24和2.25后，获得旋转矩阵R和平移矩阵T，并利用R和T及公式(5)计算目标点云旋转、平移后的矩阵Q′；
Q′＝RQ+T (5)；
步骤2.3，利用公式(6)计算适应度函数fitness；
fitness＝|RT(Q-Q′)|，i＝1,2,3,...,m(6)；
式(6)中，变换矩阵RT包含6个参数，平移变量tx,ty,tz以及旋转变量rx,ry,rz，Q’为变换后的点云集；
步骤2.4，利用公式(7)计算源点云和目标点云对应点对的局部最小距离fpbest，表示粒子当前最佳位置；
f(pi，qi)＝|p′i-q′i| (7)；
式(7)中，p′i和q′i表示待配准点对pi和qi的k邻域的重心；
步骤2.5，令fp＝fpbest，fitnessbest＝fitness，RT＝RT(R0,T0)，fp表示粒子历史最佳位置，fitnessbest表示所有粒子的最佳位置函数，RT表示最佳变换矩阵，RT(R0,T0)表示初始化的变换矩阵，实现源点云和目标点云粒子群粗配准；
具体按照以下步骤实施：
步骤2.51，设置r1i＝rand()，r2i＝rand()，其中r1i和r2i都是伪随机数，rand()是随机函数，在(0,1)的均匀分布中生成随机数；
步骤2.52，以空间矩阵RT代替速度V，并利用公式(8)更新RT，对目标点云Q进行空间矩阵变换；
Figure BDA0002452166470000151
式(8)中，D为维数，V为粒子i的速度，pi为粒子i在一次迭代中的最佳位置，c1和c2为加速器，设置为2；r1和r2是伪随机数，设置为0-1；
步骤2.53，利用公式(6)计算此次变换的适应度函数，并利用公式(7)计算源点云P与目标点云Q的各点对的局部距离函数，并确定其中最小值fpbest；
步骤2.54，比较fp和fpbest，如果fpbest小于fp，则继续比较fitnebest和fitness，如果fitness小于fitnessbest，则fitnessbest＝fitness，fp＝fpbest；否则，fitnessbest不变；如果fpbest大于fp，则fp不变；
步骤2.55，经过多次实验调整，确定满足实验数据的最佳迭代次数M，如果迭代M次后，标记此时的最近点对集合；否则，重复步骤2.51至步骤2.55；
步骤3，经过步骤2后，得到初配准的初始配准点对，再采用蚁群优化算法对其进行全局优化处理，得到更加精准的粗配准结果；具体按照以下步骤实施：
步骤3.1，对目标点集中的点与点之间赋予初始信息素值C＝2，并在目标点云上随机选择n个点作为“蚂蚁”的初始位置；不同的点云数据，n的取值不同，本发明的两个实例n的值分别为：bunny模型n＝20，椅子n＝38；
步骤3.2，利用公式(9)计算每个“蚂蚁”点在其邻域内的状态转移概率并进行状态转移；
Figure BDA0002452166470000161
式(9)中，pij k表示编号为k的蚂蚁从i位置向j位置的“转移”的概率；τij(t)表示t时刻位置i与位置j的信息素值；ηij通常设为i位置和j位置之间距离的倒数；Jk(i)表示i位置上的蚂蚁k未访问的邻近位置点的集合；α,β为启发式因子，其中α＝0.9,β＝1.06，分别表示确定信息素浓度τij和距离间的评价参数；
步骤3.3，对所有“蚂蚁”点搜索到一条路径后，遍历目标点集中所有点；
步骤3.4，利用公式(10)进行信息素迭代计算；
Figure BDA0002452166470000162
式(10)中，ρ为信息素蒸发参数，ρ＝0.1，0＜ρ＜1；Δτk ij表示蚂蚁k在i位置和j位置之间信息素的浓度，并利用公式(11)计算获得；
Figure BDA0002452166470000163
式(11)中，
Figure BDA0002452166470000164
和
Figure BDA0002452166470000165
是选择的对应点对的法向量，且两法向量相似度高时，
Figure BDA0002452166470000166
k表示第k个粒子的搜索过程；λ＝8为一个常数，是为保证信息素浓度函数的取值符合算法要求，否则信息素浓度函数会永远存在于区间[0,1]之间，无法进行后续算法实施；
步骤3.5，根据本算法中信息素浓度更新函数的具体实验结果及具体数据模型的多次实验结果进行调整，并确定最终比较的阈值；若路径存在更新的信息素值τij大于预先设定的阈值，则搜索结束，该路径(点集合)与源点云组成的点对作为配准的点对；否则重复步骤3.2至步骤3.4；不同的点云数据，其阈值设置不同，本发明的两个实例阈值分别设置为：bunny模型，阈值＝470，椅子模型，阈值＝760；
步骤3.6，利用单位四元数法求解最终选择的点对之间的空间变换矩阵，完成粗配准；配准结果如图3所示；
图3a是源点云和目标点云滤除后待配准的数据，图3b和图3e是采用基于SAC-IA粗配准算法结果的正面图和侧面图，由该图可以看出，在模型坐垫部分融合比较杂乱，有较大程度的分离；图3c和图3f是采用基于PSO粒子群的粗配准结果的正面图和侧面图，由该图可以看出，采集的数据模型在正面上平移误差较SAC-IA粗配准算法结果大，但椅背处该方法的配准效果明显优于SAC-IA粗配准算法结果，但从侧面的配准结果显示，该算法旋转误差更小，总的分析，PSO粒子群算法粗配准在该模型效果优于基于SAC-IA的粗配准结果；图3d和图3g是采用本发明的粗配准方法结果的正面图和侧面图，由该图可以看出，在模型坐垫部分层次分明且整齐，分离程度更小，源点云与目标点云的位姿调整的更合理，更有利于精确配准。
步骤4，经过步骤3后，得到优化后的粗配准结果，再采用改进的多约束条件迭代最近点精确配准算法对其进行精配准处理；
在多约束条件的ICP精配准算法中，对源点云P与目标点云Q预先设置最大迭代次数M＝200，目标函数阈值δ＝10-6m，利用公式(12)表示迭代配准中的目标函数；
Figure BDA0002452166470000181
其中，NS为点云数量，R和T为两点云配准中的旋转矩阵和平移矩阵，pi和qi分别为源点云与目标点云中的点。
具体按照以下步骤实施：
步骤4.1，搜索源点云P和目标点云Q，确定对应点对；具体按照以下步骤实施：
步骤4.11，搜索源点云P和目标点云Q，使得欧式距离‖pi-qi‖＝min最小的点对为初始对应点对；
步骤4.12，利用公式(13)计算所有点对中每点的曲率，并将这些点对存储到集合TS{(pi,qi),i＝1,2,3,...,n′}中，该集合包含初始点对及这些点对的曲率信息；
Figure BDA0002452166470000182
式(13)中，k1(pi)、k2(pi)分别是点曲率的最大值和最小值，也称为点的主曲率；
任意一点p的主曲率具体实现步骤为：
设二次平面的一般方程，如公式(14)所示；
z＝f(x,y)＝a0x2+a1y2+a2x+a3y+a4 (14)；
1.在目标点pi的k邻域内，根据最小二乘原理利用公式(15)该曲面方程的最小值；
Figure BDA0002452166470000191
2.由公式(14)和公式(15)求出该曲面的拟合曲面方程，拟合曲面方程由公式(16)表示；
r(x,y)＝(x,y,a0+a1x+a2y+a2x2+a4y2)(16)；
3.计算拟合曲面r(x,y)的偏导rx，ry，rxx，ryy，rxy，并利用公式(17)计算曲面的单位法矢量；
Figure BDA0002452166470000192
4.利用公式(18)和(19)分别计算平均曲率和高斯曲率；
Figure BDA0002452166470000193
Figure BDA0002452166470000194
其中，E＝rx·rx,F＝rx·ry,G＝ry·ry,
Figure BDA0002452166470000195
Figure BDA0002452166470000196
5.利用公式(20)和公式(21)计算点p的主曲率；
Figure BDA0002452166470000197
Figure BDA0002452166470000198
步骤4.13，利用公式(22)计算每对点对的曲率相似度ψ，并判断误差是否在取值区间之内，若在，保留该点对作为待配准的对应点对；若不在，搜索目标点qi的k邻域，并利用公式(13)计算邻域内各点的曲率信息；
ψ＝|S(pi)-S(qi)|(22)；
步骤4.14，计算点对中点pi和qi邻域各点的曲率相似度及欧式距离，将满足欧式距离最小且曲率相似度符合误差区间[0,0.08]的邻域点作为新的qi；若找不到满足以上两条约束的邻域点，则将该点对从集合TS中删除；
步骤4.15，将最终确定的点对在源点云与目标点云集合进行标记；
步骤4.2，目标函数最小化，并利用奇异值分解法计算此次迭代过程中的旋转矩阵Rk、平移矩阵Tk，利用公式(23)计算此次误差函数d；
Figure BDA0002452166470000201
其中，计算旋转矩阵Rk、平移矩阵Tk的具体步骤为：
步骤4.21，利用公式(24)表示空间变换矩阵函数；
Figure BDA0002452166470000202
步骤4.22，利用公式(25)计算源点云与目标点云的重心坐标，并利用公式(26)确定去重心化后的两个点集p′和q′，则利用公式(27)表示新的变换矩阵函数；
Figure BDA0002452166470000203
p′＝p-μP,q′＝q-μQ (26)；
Figure BDA0002452166470000204
步骤4.23，对目标函数的总误差进行化简，化简后的函数由公式(28)表示，并对其进一步表示为公式(29)；
Figure BDA0002452166470000205
Figure BDA0002452166470000206
步骤4.24，要使目标函数变为最小，即使得Tranc(RH)最小化；线性代数中，M是一个对称矩阵且B是正交矩阵，则有Tranc(M)≥Tranc(BM)，对H进行奇异值分解，得到公式(30)；
H＝U∑VT (30)
式(30)中，∑为3阶非负元素对角矩阵；
步骤4.25，令旋转矩阵R＝VUT时，RH为对称正定矩阵。此时的R为最优旋转矩阵，再利用公式(4)求出平移矩阵T；
步骤4.3，利用公式(31)更新目标点集；
Q＝RkP+tk (31)；
步骤4.4，再次计算新的对应点对使之满足欧式距离最小，并判断点对是否被标记；若标记，则直接根据标记索引，在集合TS中找到该点对；否则按照步骤4.12计算点对的曲率信息并对这些点对进行归类或者删除；
步骤4.5，重复迭代并计算两次误差的差值|d-d’|，若满足|d-d’|＜δ，则算法结束；否则继续迭代，直至达到最大迭代次数；精配准执行的结果如图4所示；
图4a和图4d分别是采用3DNDT算法精配准结果的正面图和侧面图，由该图可以看出，源点云与目标点云在采集数据的椅背处有明显的分离现象，扶手和坐垫处重合不完全，该算法精确配准误差较大；图4b和图4e是采用ICP算法精配准结果的正面图和侧面图，由该图可以看出，该算法最终的配准结果误差比3DNDT算法误差小，源点云与目标点云在采集数据的扶手、椅背顶端及坐垫处的垂直方向有分离的情况，该算法精确配准误差较小；图4c和图4f是采用本发明的精配准方法执行结果的正面图，由该图可以看出，该发明最终的配准结果误差最小，源点云与目标点云在采集数据的精确配准上，所有位置均得到更好的配准效果，源点云与目标点云的重合度最高。
将现有的粗配准方法和精配准方法和本发明的方法在多个性能方面进行了对比，对比结果如表1、表2所示。表1是将本发明的粗配准方法与现有的粗配准方法在运行时间和配准精度方面进行了对比。由表1可知，本发明粗配准算法比基于局部特征描述子FPFH结合SAC-IA算法的配准时间更短，说明本发明的计算效率高于传统粗配准算法的配准效率。而且，从配准误差结果分析，可以确定本发明配准的精度比传统粗配准有较大提高，能够为精确配准提供更准确的初值。表2是将本发明的精配准方法与现有的精配准方法进行了对比，结果表明，本发明的配准误差都是最小的；在配准时间上，由于3DNDT算法的机制，该算法耗时更短，该算法耗时更长，但与原始ICP算法比较，本算法用时更短。综合分析，本发明3D彩色点云配准方法更优异。
表1各种粗配准算法的性价对比
Figure BDA0002452166470000221
表2各种精配准算法性价对比
Figure BDA0002452166470000231

Claims (6)

1.一种全局优化与多约束条件迭代的3D彩色点云配准方法，其特征在于，具体按照以下步骤实施：
步骤1，采用基于标记的多状态离群点剔除原点云数据中的噪声、离群点及离群点簇；
步骤2，经过步骤1后，采用PSO粒子群优化算法对滤除后源点云数据和目标点云数据进行粗配准处理；
步骤3，经过步骤2后，采用蚁群优化算法对其进行全局优化处理，得到更加精准的粗配准结果；
步骤4，经过步骤3后，采用改进的多约束条件迭代最近点精确配准算法对其进行精配准处理。
2.根据权利要求1所述的一种全局优化与多约束条件迭代最近点算法的3D彩色点云配准方法，其特征在于，所述步骤1中，具体按照以下步骤实施：
步骤1.1，用k-d树进行点云的邻域搜索；
步骤1.2，经步骤1.1后，用移动最小二乘法局部表面拟合方法对点云法向量进行计算；
步骤1.3，经步骤1.2后，对不同点云数据设置不同的阈值，以点法线差值的正交分量的比值与设置阈值进行比较，以作为孤立点的判断条件，标记离群点；
步骤1.4，经步骤1.3后，以改进的DBSCAN聚类方法对标记离群点进行聚类；
步骤1.5，经步骤1.4后，统计每个聚类中标记离群点的数量占该聚类中点云总数的比值，当该比值大于设定的比例阈值时，判定并删除离群点簇；
步骤1.6，重复步骤1.5，直到删除所有的离群点簇，同时删除剩余的标记离群点；得到滤除后源点云数据和目标点云数据。
3.根据权利要求1所述的所述的一种全局优化与多约束条件迭代的3D彩色点云配准方法，其特征在于，所述步骤2中，具体按照以下步骤实施：
步骤2.1，初始化PSO粒子群算法的迭代次数M＝80，随机生成平移变量tx,ty,tz以及旋转变量rx,ry,rz，初始化学习因子c1＝2、c2＝2；
步骤2.2，经过步骤2.1后，采用单位四元数法求出源点云P和目标点云Q的旋转变换矩阵R和平移矩阵T，并计算目标点云Q旋转、平移后的矩阵Q′；具体按照以下步骤实施：
步骤2.21，计算源点云P与目标点云Q的重心，如式(1)所示；
Figure FDA0002452166460000021
式(1)中，pi和qi是源点云与目标点云的点，i为其下标，表示第i个点，n为每个点云的总数，μp和μq分别为源点云P与目标点云Q的重心；
步骤2.22，利用公式(2)计算其协方差矩阵M；
Figure FDA0002452166460000022
步骤2.23，根据协方差矩阵M，利用公式(3)构造4阶矩阵S；
Figure FDA0002452166460000031
步骤2.24，计算矩阵S的特征值和特征向量，并将最大特征值对应的特征向量作为所求的旋转变量R；
步骤2.25，利用公式(4)计算平移矩阵T；
T＝μP-RμQ (4)；
步骤2.26，经过步骤2.24和2.25后，获得旋转矩阵R和平移矩阵T，并利用R和T及公式(5)计算目标点云旋转、平移后的矩阵Q′；
Q′＝RQ+T (5)；
步骤2.3，利用公式(6)计算适应度函数fitness；
fitness＝|RT(Q-Q′)|，i＝1,2,3,...,m (6)；
式(6)中，变换矩阵RT包含6个参数，平移变量tx,ty,tz以及旋转变量rx,ry,rz，Q’为变换后的点云集；
步骤2.4，利用公式(7)计算源点云和目标点云对应点对的局部最小距离fpbest，表示粒子当前最佳位置；
f(pi，qi)＝|pi′-qi′| (7)；
式(7)中，p′i和q′i表示待配准点对pi和qi的k邻域的重心；
步骤2.5，令fp＝fpbest，fitnessbest＝fitness，RT＝RT(R0,T0)，fp表示粒子历史最佳位置，fitnessbest表示所有粒子的最佳位置函数，RT表示最佳变换矩阵，RT(R0,T0)表示初始化的变换矩阵，实现源点云和目标点云粒子群粗配准；具体按照以下步骤实施：
步骤2.51，设置r1i＝rand()，r2i＝rand()，其中r1i和r2i都是伪随机数，rand()是随机函数，在(0,1)的均匀分布中生成随机数；
步骤2.52，以空间矩阵RT代替速度V，并利用公式(8)更新RT，对目标点云Q进行空间矩阵变换；
Figure FDA0002452166460000041
式(8)中，D为维数，V为粒子i的速度，pi为粒子i在一次迭代中的最佳位置，c1和c2为加速器，设置为2；r1和r2是伪随机数，设置为0-1；
步骤2.53，利用公式(6)计算此次变换的适应度函数，并利用公式(7)计算源点云P与目标点云Q的各点对的局部距离函数，并确定其中最小值fpbest；
步骤2.54，比较fp和fpbest，如果fpbest小于fp，则继续比较fitnebest和fitness，如果fitness小于fitnessbest，则fitnebest＝fitness，fp＝fpbest；否则，fitnessbest不变；如果fpbest大于fp，则fp不变；
步骤2.55，确定满足实验数据的最佳迭代次数M，如果迭代M次后，标记此时的最近点对集合；否则，重复步骤2.51至步骤2.55。
4.根据权利要求1所述的一种全局优化与多约束条件迭代的3D彩色点云配准方法，其特征在于，所述步骤3中，具体按照以下步骤实施：
步骤3.1，对目标点集中的点与点之间赋予初始信息素值C＝2，并在目标点云上随机选择n个点作为“蚂蚁”的初始位置；
步骤3.2，利用公式(9)计算每个“蚂蚁”点在其邻域内的状态转移概率并进行状态转移；
Figure FDA0002452166460000051
式(9)中，pij k表示编号为k的蚂蚁从i位置向j位置的“转移”的概率；τij(t)表示t时刻位置i与位置j的信息素值；ηij通常设为i位置和j位置之间距离的倒数；Jk(i)表示i位置上的蚂蚁k未访问的邻近位置点的集合；α,β为启发式因子，其中α＝0.9,β＝1.06，分别表示确定信息素浓度τij和距离间的评价参数；
步骤3.3，对所有“蚂蚁”点搜索到一条路径后，遍历目标点集中所有点；
步骤3.4，利用公式(10)进行信息素迭代计算；
Figure FDA0002452166460000052
式(10)中，ρ为信息素蒸发参数，ρ＝0.1，0＜ρ＜1；Δτk ij表示蚂蚁k在i位置和j位置之间信息素的浓度，并利用公式(11)计算获得；
Figure FDA0002452166460000053
式(11)中，
Figure FDA0002452166460000054
和
Figure FDA0002452166460000055
是选择的对应点对的法向量，且两法向量相似度高时，
Figure FDA0002452166460000056
k表示第k个粒子的搜索过程；λ＝8；
步骤3.5，确定最终比较的阈值，若路径存在更新的信息素值τij大于预先设定的阈值，则搜索结束，该路径(点集合)与源点云组成的点对作为配准的点对；否则重复步骤3.2至步骤3.4；
步骤3.6，利用单位四元数法求解最终选择的点对之间的空间变换矩阵，完成粗配准。
5.根据权利要求1所述的一种全局优化与多约束条件迭代的3D彩色点云配准方法，其特征在于，所述步骤4中，在多约束条件的ICP精配准算法中，对源点云P与目标点云Q预先设置最大迭代次数M＝200，目标函数阈值δ＝10-6m，利用公式(12)表示迭代配准中的目标函数；
Figure FDA0002452166460000061
其中，NS为点云数量，R和T为两点云配准中的旋转矩阵和平移矩阵，pi和qi分别为源点云与目标点云中的点；
具体按照以下步骤实施：
步骤4.1，搜索源点云P和目标点云Q，确定对应点对；具体按照以下步骤实施：
步骤4.11，搜索源点云P和目标点云Q，使得欧式距离‖pi-qi‖＝min最小的点对为初始对应点对；
步骤4.12，利用公式(13)计算所有点对中每点的曲率，并将这些点对存储到集合TS{(pi,qi),i＝1,2,3,...,n′}中，该集合包含初始点对及这些点对的曲率信息；
Figure FDA0002452166460000062
式(13)中，k1(pi)、k2(pi)分别是点曲率的最大值和最小值，也称为点的主曲率；
步骤4.13，利用公式(22)计算每对点对的曲率相似度ψ，并判断误差是否在取值区间之内，若在，保留该点对作为待配准的对应点对；若不在，搜索目标点qi的k邻域，并利用公式(13)计算邻域内各点的曲率信息；
ψ＝|S(pi)-S(qi)| (22)；
步骤4.14，计算点对中点pi和qi邻域各点的曲率相似度及欧式距离，将满足欧式距离最小且曲率相似度符合误差区间[0,0.08]的邻域点作为新的qi；若找不到满足以上两条约束的邻域点，则将该点对从集合TS中删除；
步骤4.15，将最终确定的点对在源点云与目标点云集合进行标记；
步骤4.2，目标函数最小化，并利用奇异值分解法计算此次迭代过程中的旋转矩阵Rk、平移矩阵Tk，利用公式(23)计算此次误差函数d；
Figure FDA0002452166460000071
其中，计算旋转矩阵Rk、平移矩阵Tk的具体步骤为：
步骤4.21，利用公式(24)表示空间变换矩阵函数；
Figure FDA0002452166460000072
步骤4.22，利用公式(25)计算源点云与目标点云的重心坐标，并利用公式(26)确定去重心化后的两个点集p′和q′，则利用公式(27)表示新的变换矩阵函数；
Figure FDA0002452166460000073
p′＝p-μP,q′＝q-μQ (26)；
Figure FDA0002452166460000074
步骤4.23，对目标函数的总误差进行化简，化简后的函数由公式(28)表示，并对其进一步表示为公式(29)；
Figure FDA0002452166460000081
Figure FDA0002452166460000082
步骤4.24，要使目标函数变为最小，即使得Tranc(RH)最小化；线性代数中，M是一个对称矩阵且B是正交矩阵，则有Tranc(M)≥Tranc(BM)，对H进行奇异值分解，得到公式(30)；
H＝U∑VT (30)
式(30)中，∑为3阶非负元素对角矩阵；
步骤4.25，令旋转矩阵R＝VUT时，RH为对称正定矩阵；此时的R为最优旋转矩阵，并计算平移矩阵T；
步骤4.3，利用公式(31)更新目标点集；
Q＝RkP+Tk (31)；
步骤4.4，再次计算新的对应点对使之满足欧式距离最小，并判断点对是否被标记；若标记，则直接根据标记索引，在集合TS中找到该点对；否则按照步骤4.12计算点对的曲率信息并对这些点对进行归类或者删除；
步骤4.5，重复迭代并计算两次误差的差值|d-d’|，若满足|d-d’|＜δ，则算法结束；否则继续迭代，直至达到最大迭代次数。
6.根据权利要求5所述的一种全局优化与多约束条件迭代的3D彩色点云配准方法，其特征在于，所述任意一点p的主曲率，具体实现步骤为：
设二次平面的一般方程，如公式(14)所示；
z＝f(x,y)＝a0x2+a1y2+a2x+a3y+a4 (14)；
在目标点pi的k邻域内，根据最小二乘原理利用公式(15)该曲面方程的最小值；
Figure FDA0002452166460000091
由公式(14)和公式(15)求出该曲面的拟合曲面方程，拟合曲面方程由公式(16)表示；
r(x,y)＝(x,y,a0+a1x+a2y+a3x2+a4y2) (16)；
计算拟合曲面r(x,y)的偏导rx，ry，rxx，ryy，rxy，并利用公式(17)计算曲面的单位法矢量；
Figure FDA0002452166460000092
利用公式(18)和(19)分别计算平均曲率和高斯曲率；
Figure FDA0002452166460000093
Figure FDA0002452166460000094
其中，E＝rx·rx,F＝rx·ry,G＝ry·ry,
Figure FDA0002452166460000095
Figure FDA0002452166460000096
利用公式(20)和公式(21)计算点p的主曲率；
Figure FDA0002452166460000097
Figure FDA0002452166460000098
CN202010295862.7A 2020-04-15 2020-04-15 一种全局优化与多约束条件迭代的3d彩色点云配准方法 Active CN111563920B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010295862.7A CN111563920B (zh) 2020-04-15 2020-04-15 一种全局优化与多约束条件迭代的3d彩色点云配准方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010295862.7A CN111563920B (zh) 2020-04-15 2020-04-15 一种全局优化与多约束条件迭代的3d彩色点云配准方法

Publications (2)

Publication Number Publication Date
CN111563920A true CN111563920A (zh) 2020-08-21
CN111563920B CN111563920B (zh) 2023-04-07

Family

ID=72071770

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010295862.7A Active CN111563920B (zh) 2020-04-15 2020-04-15 一种全局优化与多约束条件迭代的3d彩色点云配准方法

Country Status (1)

Country Link
CN (1) CN111563920B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112017220A (zh) * 2020-08-27 2020-12-01 南京工业大学 一种基于抗差约束最小二乘算法的点云精确配准方法
CN112132876A (zh) * 2020-09-25 2020-12-25 天津大学 2d-3d图像配准中的初始位姿估计方法
CN112396641A (zh) * 2020-11-17 2021-02-23 中国矿业大学(北京) 一种基于全等二基线匹配的点云全局配准方法
CN112488153A (zh) * 2020-11-06 2021-03-12 中国人民解放军国防科技大学 基于聚类趋向性原理的气味源自主搜索及定位方法
CN112669359A (zh) * 2021-01-14 2021-04-16 武汉理工大学 一种三维点云配准方法、装置、设备及存储介质
CN113450395A (zh) * 2021-06-16 2021-09-28 北京理工大学 基于光学跟踪探针划取点的手术导航空间配准方法及装置
CN113793369A (zh) * 2021-09-30 2021-12-14 太原理工大学 一种复杂自由曲面迭代双向法线投影配准方法
CN115797418A (zh) * 2022-09-27 2023-03-14 湖南科技大学 一种基于改进icp的复杂机械零件测量点云配准方法及系统
CN116363178A (zh) * 2023-02-27 2023-06-30 西安电子科技大学 针对点云特征的多任务优化方法
CN116883469A (zh) * 2023-07-20 2023-10-13 中国矿业大学 平面特征约束下基于eiv模型描述的点云配准方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105551015A (zh) * 2015-12-02 2016-05-04 南京邮电大学 一种散乱点云图像配准方法
CN106600639A (zh) * 2016-12-09 2017-04-26 江南大学 遗传算法结合自适应阈值约束的icp位姿定位技术
WO2018072351A1 (zh) * 2016-10-20 2018-04-26 北京工业大学 一种基于粒子群优化算法对支持向量机的优化方法
EP3435029A1 (en) * 2017-07-28 2019-01-30 The Boeing Company Resolution adaptive mesh for performing 3-d metrology of an object

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105551015A (zh) * 2015-12-02 2016-05-04 南京邮电大学 一种散乱点云图像配准方法
WO2018072351A1 (zh) * 2016-10-20 2018-04-26 北京工业大学 一种基于粒子群优化算法对支持向量机的优化方法
CN106600639A (zh) * 2016-12-09 2017-04-26 江南大学 遗传算法结合自适应阈值约束的icp位姿定位技术
EP3435029A1 (en) * 2017-07-28 2019-01-30 The Boeing Company Resolution adaptive mesh for performing 3-d metrology of an object

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杨高朝: "基于特征提取的点云自动配准优化研究", 《计算机工程与应用》 *
段德全等: "基于粒子群优化算法的散乱点云数据配准", 《广西师范大学学报(自然科学版)》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112017220A (zh) * 2020-08-27 2020-12-01 南京工业大学 一种基于抗差约束最小二乘算法的点云精确配准方法
CN112017220B (zh) * 2020-08-27 2023-07-28 南京工业大学 一种基于抗差约束最小二乘算法的点云精确配准方法
CN112132876A (zh) * 2020-09-25 2020-12-25 天津大学 2d-3d图像配准中的初始位姿估计方法
CN112132876B (zh) * 2020-09-25 2022-03-29 天津大学 2d-3d图像配准中的初始位姿估计方法
CN112488153B (zh) * 2020-11-06 2022-07-05 中国人民解放军国防科技大学 基于聚类趋向性原理的气味源自主搜索及定位方法
CN112488153A (zh) * 2020-11-06 2021-03-12 中国人民解放军国防科技大学 基于聚类趋向性原理的气味源自主搜索及定位方法
CN112396641A (zh) * 2020-11-17 2021-02-23 中国矿业大学(北京) 一种基于全等二基线匹配的点云全局配准方法
CN112396641B (zh) * 2020-11-17 2023-07-21 中国矿业大学(北京) 一种基于全等二基线匹配的点云全局配准方法
CN112669359A (zh) * 2021-01-14 2021-04-16 武汉理工大学 一种三维点云配准方法、装置、设备及存储介质
CN112669359B (zh) * 2021-01-14 2023-05-26 武汉理工大学 一种三维点云配准方法、装置、设备及存储介质
CN113450395B (zh) * 2021-06-16 2022-12-02 北京理工大学 基于光学跟踪探针划取点的手术导航空间配准方法及装置
CN113450395A (zh) * 2021-06-16 2021-09-28 北京理工大学 基于光学跟踪探针划取点的手术导航空间配准方法及装置
CN113793369A (zh) * 2021-09-30 2021-12-14 太原理工大学 一种复杂自由曲面迭代双向法线投影配准方法
CN113793369B (zh) * 2021-09-30 2024-03-15 太原理工大学 一种复杂自由曲面迭代双向法线投影配准方法
CN115797418A (zh) * 2022-09-27 2023-03-14 湖南科技大学 一种基于改进icp的复杂机械零件测量点云配准方法及系统
CN116363178A (zh) * 2023-02-27 2023-06-30 西安电子科技大学 针对点云特征的多任务优化方法
CN116363178B (zh) * 2023-02-27 2024-01-30 西安电子科技大学 针对点云特征的多任务优化方法
CN116883469A (zh) * 2023-07-20 2023-10-13 中国矿业大学 平面特征约束下基于eiv模型描述的点云配准方法
CN116883469B (zh) * 2023-07-20 2024-01-19 中国矿业大学 平面特征约束下基于eiv模型描述的点云配准方法

Also Published As

Publication number Publication date
CN111563920B (zh) 2023-04-07

Similar Documents

Publication Publication Date Title
CN111563920B (zh) 一种全局优化与多约束条件迭代的3d彩色点云配准方法
CN109887015B (zh) 一种基于局部曲面特征直方图的点云自动配准方法
CN107886529B (zh) 一种用于三维重建的点云配准方法
CN109410321B (zh) 基于卷积神经网络的三维重建方法
CN100559398C (zh) 自动的深度图像配准方法
CN107967463B (zh) 一种基于合成图像和深度学习的虚拟人脸识别方法
CN112562081B (zh) 一种用于视觉分层定位的视觉地图构建方法
CN109034065B (zh) 一种基于点云的室内场景物体提取方法
CN108627798B (zh) 基于线性判别分析和梯度提升树的wlan室内定位算法
WO2021248908A1 (zh) 基于高斯过程回归的室外场景三维点云的地面提取方法
CN110599506A (zh) 一种复杂异形曲面机器人三维测量的点云分割方法
CN109493344A (zh) 一种大规模城市三维场景的语义分割方法
CN114200477A (zh) 一种激光三维成像雷达地面目标点云数据处理方法
CN110942476A (zh) 基于二维图像引导的改进三维点云配准方法、系统及可读存储介质
CN111145129A (zh) 基于超体素的点云去噪方法
CN109323697B (zh) 一种针对室内机器人任意点启动时粒子快速收敛的方法
CN112084895B (zh) 一种基于深度学习的行人重识别方法
CN108549905A (zh) 一种严重遮挡情况下的精准目标跟踪方法
CN114358166B (zh) 一种基于自适应k均值聚类的多目标定位方法
CN109934859A (zh) 一种基于特征增强多维权重描述子的icp配准方法
CN111553410B (zh) 基于关键点局部曲面特征直方图和空间关系的点云识别方法
Liu et al. Deep learning of directional truncated signed distance function for robust 3D object recognition
CN116842818A (zh) 基于鹈鹕优化算法的结构面产状分组方法
CN111583098B (zh) 一种基于序列图像的线段聚类及拟合方法和系统
CN112101461A (zh) 一种基于hrtf-pso-fcm的无人机侦察视觉信息可听化方法

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