CN111080684B - 一种点邻域尺度差异描述的点云配准方法 - Google Patents

一种点邻域尺度差异描述的点云配准方法 Download PDF

Info

Publication number
CN111080684B
CN111080684B CN201911276127.5A CN201911276127A CN111080684B CN 111080684 B CN111080684 B CN 111080684B CN 201911276127 A CN201911276127 A CN 201911276127A CN 111080684 B CN111080684 B CN 111080684B
Authority
CN
China
Prior art keywords
point
point cloud
key
neighborhood
points
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
Application number
CN201911276127.5A
Other languages
English (en)
Other versions
CN111080684A (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.)
Harbin Engineering University
Original Assignee
Harbin Engineering 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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201911276127.5A priority Critical patent/CN111080684B/zh
Publication of CN111080684A publication Critical patent/CN111080684A/zh
Application granted granted Critical
Publication of CN111080684B publication Critical patent/CN111080684B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • 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
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)

Abstract

本发明提供一种点邻域尺度差异描述的点云配准方法,通过对源点云和目标点云中离散点进行最小二乘曲面拟合得到局部曲面,求出曲面的形状指数SI,即该离散点的形状指数,选取形状指数在邻域内最大或最小且满足阈值的点作为点云的关键点;进行特征描述符构造,计算关键点在不同邻域半径下的特征归一化向量差值和法向量夹角差值组合成点领域尺度差异描述符;根据特征描述符的相似程度找出对应点并使用二重筛选和基于全局距离的最优查找算法分别滤除错误点对和估计对应关系。本发明得到的关键点具有很好的代表性和区别性,对点云分布密度差异较大或存在噪声点的情况效果明显,计算简单,提高点云配准速度和精度,具有很好的抗干扰能力。

Description

一种点邻域尺度差异描述的点云配准方法
技术领域
本发明属于计算机视觉领域,具体涉及一种点邻域尺度差异描述的点云配准方法。
背景技术
随着众多新型传感技术的发展,点云数据的获取越来越高效。三维点云数据处理技术也得以快速发展,并且在计算机图形学、CAD造型设计、文物古迹保护、汽车制造业和模具制造业等众多领域都有广泛应用。然而由于外界环境、物体外形大小和扫描设备的自身的局限性等因素制约,不能通过单次扫描的点云数据得到物体的完整面形。因此,需要高效的点云配准算法对从不同视角获取的点云数据进行拼接。在点云配准的众多方法中,最经典的是迭代最近点(Iterative Closest Point,ICP)算法。ICP算法对于待配准点云初始位置要求很高,即需要目标点云与源点云大致重合时,才能得到高精度的配准结果,否则迭代易陷入局部最优,收敛得到错误的配准结果。为了得到更好的配准效果,国内外许多学者提出了其他方法,其中基于特征描述的点云配准方法应用最为广泛。该方法主要是为了在初始条件未知的情况下快速估算出一个点云配准矩阵。主要配准思想是通过计算点的特征(比如FPFH,Fast Point Feature Histogram,LBP,SHOT)来寻找对应点,从而获得多视角点云之间的旋转平移关系。FPFH特征需要计算任意点以及它的k近邻点之间两两相互连接的角度特征,计算量比较大,计算效率较低;LBP算子特征信息的描述过于单一;SHOT算子纬度高计算复杂。点邻域尺度差异描述的配准算法提出了点邻域尺度差异描述符来表示点特征,其具有计算简单,辨识能力强的优点;并且在对应关系查找方面,为了有效地选取并剔除潜在的错误对应点对设计了二重筛选和全局最优查找方法,提高了配准的精度和效率。
发明内容
本发明的目的在于提供了一种点邻域尺度差异描述的点云配准方法,对点云分布密度差异较大或存在噪声点的情况效果明显,计算简单,提高配准效率。
本发明的目的是这样实现的:
一种点邻域尺度差异描述的点云配准方法,具体的实现步骤为:
步骤1.关键点查找,对源点云和目标点云中离散点进行最小二乘曲面拟合得到局部曲面,求出曲面的形状指数SI,即该离散点的形状指数,选取形状指数在邻域内最大或最小且满足阈值的点作为点云的关键点;
步骤2.特征描述符计算,进行特征描述符构造,计算关键点在不同邻域半径下的特征归一化向量差值和法向量夹角差值组合成点领域尺度差异描述符;
步骤3.对应关系估计,根据特征描述符的相似程度找出对应点并使用二重筛选和基于全局距离的最优查找算法分别滤除错误点对和估计对应关系。
所述步骤1的具体步骤为:
步骤1.1.对点云每一点使用最小二乘法进行曲面拟合,得到离散点邻域的局部曲面z=r(u,v);
步骤1.2.在得到拟合的二次参数曲面后,由二次参数曲面的第一类基本量E、F、G和第二类基本量L、M、N计算得到曲面的高斯曲率K和平均曲率H
由高斯曲率和平均曲率计算得到曲面的最大主曲率k1和最小主曲率k2
其中E=ru·ru,F=ru·rv,G=rv·rv,L=ruu·n,N=ryy·n,M=rxy·n,其中ru、rv、ruu、rvv、ruv是拟合曲面的偏微分,n是曲面的单位法向量,曲面P的形状指数定义为
步骤1.3.计算每点对应曲面Pi的形状指数SI(Pi);
步骤1.4.选择点云的关键点,源点云pi邻域内所有点的形状指数均值为μ(pi),满足下面两个条件中任意一个的都选作关键点
且SI(pi)≥(1+a)μ(pi)
且SI(pi)≤(1-b)μ(pi)
其中a、b是用来调整关键点数量的两个标量参数,将选取后的关键点集记为
所述步骤2的具体步骤为:
步骤2.1.对每个关键点p,取l个大小不同的邻域半径,相邻邻域半径之差为△r,有r1<r2<....<rl,对于邻域半径rl的邻域,建立加权协方差矩阵Cl
式中pi为该邻域中的点;
步骤2.2.对加权协方差矩阵Cl进行特征值分解得特征值{λl1,λl2,λl3}(λl1>λl2>λl3)及对应的特征向量{vl1,vl2,vl3},由此建立采样点的局部邻域参考坐标系,三个特征向量构成坐标轴的三个方向;
步骤2.3.对于协方差矩阵Cl构建特征值归一化向量用hl表示为
根据不同邻域半径下的hl差异计算Δhk=Δhk+1-Δhk(k=1,2,…,l-1),随后将△h组合后将得到一个3l-3维向量用H表示为H=(Δh1 Δh2 … Δhl-1);
步骤2.4.对采样点选取了l个大小不同的邻域半径后,能够采样点不同尺度下的l个法向量n1,n2,…,nl,而通过法向量点乘的方法计算得到两个法向量之间的夹角余弦,连续不同尺度下的法向量的夹角余弦αi
计算出连续不同尺度下的法向量的夹角余弦αi,随后将这些αi组合得到一个l-1维的向量N表示为
N=(α1,α2,…αl-1);
步骤2.5.对每个关键点通过选取不同尺度的邻域空间并计算上述局部特征差异用H和N表示,将H和N组合作为该关键点的局部特征描述符得到一个4l-4维的向量F表示为F=[Δh1 Δh2 … Δhl-1 α1 α2 … αl-1];
步骤2.6.对于源点云和目标点云的关键点集和/>中的每一个点计算其特征描述符,将得到源点云关键点特征描述子集/>和目标点云关键点特征描述符集/>
所述步骤3的具体步骤为:
步骤3.1.源点云P和目标点云Q的关键点特征描述符集记为对于/>中的每个/>在/>中寻找与其特征描述符欧式距离最近的将/>记为候选对应点对;
步骤3.2.搜索的k近邻特征描述符集合/>定义阈值η,当时,认为pi和qi1是一组对应点,记为对应点对ei={pi,qi1};若不满足,将剔除,遍历关键点集/>最终得到对应点对集记为E={e1,e2,…en};
步骤3.3.从对应点对集中随机抽选n(n≥3)组对应点对并计算变换矩阵Ti,根据Ti将源点云P变换到新的坐标系下为计算全局距离/>
其中P={pi|i=1,2...,N},Q={qj|j=1,2...,M}分别表示源点云和目标点云,表示在点云Q中的点/>与其相对应的欧式距离/>小于阈值Dε的点云P中的点,/>这样的点的总数,ξ为可调参数;
步骤3.4.设定迭代停止阈值Dε,当时停止迭代,Ti被选取为最优变换,如果不满足/>则再次抽选新的对应点对子集,迭代重复计算/>直至满足/>或达到最大迭代次数Nmax,当达到最大迭代次数时,选取/>最小的变换矩阵为最优变换矩阵。
本发明的有益效果在于:提出了基于形状指数的关键点查找方法,得到的关键点具有很好的代表性和区别性;提出描述关键点不同邻域半径特征信息差异的局部特征描述符,对点云分布密度差异较大或存在噪声点的情况效果明显,计算简单,提高配准效率;提出了二重筛选方法及全局最优查找算法,提高点云配准速度和精度;提出的点云配准方法计算简单,配准准确率高,具有很好的抗干扰能力。
附图说明
图1为点云数据初始位置。
图2为点云关键点分布。
图3为匹配点连线。
图4为配准结果。
图5为本发明的流程图。
图6为配准过程及结果参数。
具体实施方式
下面结合附图对本发明做进一步的描述:
实施例1
本发明使用的点云是斯坦福大学点云模型库和和自研的结构光测量仪获得。使用Bunny、Dragon和Armadillo、Toy四组点云数据给出实施方式,并结合附图加以说明。图1是点云数据初始位置,同一组点云数据中含有两个视角的点云数据。
步骤1.关键点查找,对点云每一点使用最小二乘法进行曲面拟合,得到离散点邻域的局部曲面z=r(u,v)。在得到拟合的二次参数曲面后,由二次参数曲面的第一类基本量E,F,G和第二类基本量L,M,N可以计算得到曲面的高斯曲率K和平均曲率H:
由高斯曲率和平均曲率还可得到曲面的最大和最小主曲率k1,k2
式中:E=ru·ru;F=ru·rv;G=rv·rv;L=ruu·n;N=ryy·n;M=rxy·n;其中ru、rv、ruu、rvv、ruv是拟合曲面的偏微分,n是曲面的单位法向量。形状指数定义为:
其中,P为曲面;k1和k2分别表示曲面的最大和最小主曲率。
计算出每点对应曲面Pi的形状指数SI(Pi)。源点云pi邻域内所有点的形状指数均值为μ(pi),满足下面两个条件中任意一个的都选作关键点:
且SI(pi)≥(1+a)μ(pi);
且SI(pi)≤(1-b)μ(pi).
其中a,b是用来调整关键点数量的两个标量参数。将选取后的关键点集记为
步骤2.特征描述符计算,对每个关键点p,取l个大小不同的邻域半径,相邻邻域半径之差为△r,有r1<r2<....<rl,对于邻域半径rl的邻域,建立加权协方差矩阵Cl
式中:pi为该邻域中的点。对加权协方差矩阵Cl进行特征值分解可得特征值{λl1l2l3}(λl1>λl2>λl3)及对应的特征向量{vl1,vl2,vl3},由此建立采样点的局部邻域参考坐标系,三个特征向量构成坐标轴的三个方向。对于协方差矩阵Cl构建特征值归一化向量用hl表示为:
根据不同邻域半径下的hl差异计算Δhk=Δhk+1-Δhk(k=1,2,…,l-1),随后将△h组合后将得到一个3l-3维向量用H表示为:
H=(Δh1 Δh2 … Δhl-1)
对采样点选取了l个大小不同的邻域半径后,能够采样点不同尺度下的l个法向量n1,n2,…,nl,而通过法向量点乘的方法可得到两个法向量之间的夹角余弦,连续不同尺度下的法向量的夹角余弦αi
计算出连续不同尺度下的法向量的夹角余弦αi,随后将这些αi组合得到一个l-1维的向量N表示为
N=(α1,α2,…αl-1)
对每个关键点通过选取不同尺度的邻域空间并计算上述局部特征差异用H和N表示,将H和N组合作为该关键点的局部特征描述符可得到一个4l-4维的向量F表示为:
F=[Δh1 Δh2 … Δhl-1 α1 α2 … αl-1]
对于源点云和目标点云的关键点集和/>中的每一个点计算其特征描述符,将得到源点云关键点特征描述子集/>和目标点云关键点特征描述符集
步骤3.对应关系估计,源点云P和目标点云Q的关键点特征描述符集记为和/>对于/>中的每个/>在/>中寻找与其特征描述符欧式距离最近的/>将/>记为候选对应点对。搜索/>的k近邻特征描述符集合定义阈值η,当/>时,认为pi和qi1是一组对应点,记为对应点对ei={pi,qi1};若不满足,将/>剔除,遍历关键点集/>最终得到对应点对集记为E={e1,e2,...en}。
从对应点对集中随机抽选n(n≥3)组对应点对并计算变换矩阵Ti,根据Ti将源点云P变换到新的坐标系下为计算全局距离/>
其中P={pi|i=1,2...,N},Q={qj|j=1,2…,M}分别表示源点云和目标点云,表示能够在点云Q中找到一点/>与其相对应的欧式距离/>小于阈值Dε的点云P中的点,/>是/>这样的点的总数,ξ为可调参数。
设定迭代停止阈值Dε,当时停止迭代,Ti被选取为最优变换,如果不满足/>则再次抽选新的对应点对子集,迭代重复计算/>直至满足或达到最大迭代次数Nmax,当达到最大迭代次数时,选取/>最小的变换矩阵为最优变换矩阵。
图1为源点云和目标点云的初始位置,图2为源点云与目标点云的关键点分布图,可以看出关键点分布均匀,位于曲面曲率变化较大的位置,对于源点云和目标点云,关键点的位置对应也较明显。图3为对应点连线,可以看出候选对应点对中夹杂着许多不正确的对应点对。从对比图中可以发现,二重筛选算法确定的对应点对将错误对应关系基本去除了。对应点对的分布也比较均匀,很好的避免了配准的局部最优错误。图4为配准结果图,可以看出配准精度较高。
图6为配准过程及结果参数,可以看出配准精度良好,且配准消耗的时间短,充分说明了算法配准性能的优越性。
综上分析,基于点邻域尺度差异描述的点云配准算法配准速度快,配准精度高,抗干扰能力强,是一种性能良好的点云配准算法。
实施例2
发明的目的是公开一种点邻域尺度差异描述的点云配准方法。首先对源点云和目标点云中离散点进行最小二乘曲面拟合得到局部曲面,求出曲面的形状指数SI,即该离散点的形状指数。选取形状指数在邻域内最大或最小且满足阈值的点作为点云的关键点。其次,进行特征描述符构造,计算关键点在不同邻域半径下的特征归一化向量差值和法向量夹角差值组合成点领域尺度差异描述符。最后,根据特征描述符的相似程度找出对应点并使用二重筛选和基于全局距离的最优查找算法分别滤除错误点对和估计对应关系。此方法设计的特征描述符计算简单,辨识性强且对应关系查找算法提高了配准的准确率。
本发明的具体实施步骤是:
步骤1.关键点查找:由于点云点数较多,且大部分点位于平面等局部特征较弱的位置,为了加快配准速度,先查找点云的关键点为配准做准备。利用形状指数理论,点云中某点的局部拟合曲面的形状指数相比周围其他点形状指数差异越大,说明该点周围表面变化越剧烈,则该点很可能是被测物体表面中的尖锐点,这样的点具有很好的代表性和区别性,符合作为关键点的要求,提出了基于形状指数的关键点选取方法。首先对点云中点使用最小二乘法进行曲面拟合得到离散点邻域的局部曲面z=r(u,v)。然后计算出曲面的最大和最小主曲率k1,k2作为该点的曲率。再由k1,k2计算出每点p的形状指数SI(p)。p邻域内所有点的形状指数均值为μ(pi),满足下面两个条件中任意一个的都选作关键点:
且SI(pi)≥(1+a)μ(pi)
且SI(pi)≤(1-b)μ(pi)
其中a,b是用来调整关键点数量的两个标量参数。将选取后的关键点集记为
步骤2.特征描述符计算:对于传统的局部特征描述符计算方法,当采样点附近点云分布密度差异较大或存在噪声点时,使局部邻域半径选取困难,且提取特征匹配率低。提出一种描述关键点不同邻域半径特征信息差异的局部特征描述符,提高配准精度。对每个关键点p,取l个大小不同的邻域半径,相邻邻域半径之差为△r,有r1<r2<....<rl,对于邻域半径rl的邻域,建立加权协方差矩阵Cl。对加权协方差矩阵Cl进行特征值分解可得特征值{λl1,λl2,λl3}(λl1>λl2>λl3)及对应的特征向量{vl1,vl2,vl3},由此建立采样点的局部邻域参考坐标系,三个特征向量构成坐标轴的三个方向。根据特征值计算特征值归一化向量hl,计算不同邻域半径下hl的差异△hk(k=1,..,l-1),将△h组合后将得到一个3l-3维向量H=(Δh1 Δh2 … Δhl-1)。对不同尺度下的l个法向量n1,n2,…,nl,计算连续不同尺度下的法向量的夹角余弦αi,随后将这些αi组合得到一个l-1维的向量N=(α1,α2,…αl-1)。将H和N组合得到一个4l-4维的向量F,即关键点P对应的特征描述符。遍历源点云和目标点云的所有关键点得到源点云和目标点云特征描述符集。
步骤3.对应关系估计,在计算初始对应关系中,常用的算法是计算源点云关键点特征描述子在目标点云关键点特征描述子中的最近点及其距离,若距离小于一定阈值则将此点对作为对应点对,若源点云中关键点特征描述子计算出现误差,极易出现误匹配问题,加大配准难度。为了提高初始对应关系准确率,提出了二重筛选方法与全局最优查找算法。首先找出目标点云中特征描述符之间的欧式距离与源点云中关键点p最近的点q,搜索q的k近邻特征描述符集合,计算集合中每点与源点云中关键点p的特征描述符欧式距离,并设定阈值筛选,满足阈值的邻域点加入到对应点对集E。从对应点对集中随机抽选n(n≥3)组对应点对并计算变换矩阵Ti,根据Ti将源点云P变换到新的坐标系下为目标点云为Q,计算全局距离/>设定迭代停止阈值Dε,当/>时停止迭代,Ti被选取为最优变换,如果不满足/>则再次抽选新的对应点对子集,迭代重复计算/>直至满足/>或达到最大迭代次数Nmax,当达到最大迭代次数时,选取/>最小的变换矩阵为最优变换矩阵。
本发明的目的是公开一种点邻域尺度差异描述的点云配准方法。首先,对输入点云中离散点进行最小二乘曲面拟合得到局部曲面,求出曲面的最大和最小主曲率,进而求出曲面的形状指数,选取形状指数在邻域内最大或最小且满足阈值的点作为点云的关键点。然后,对每个关键点在l个邻域半径下建立局部坐标系,计算得到3l个特征值和l个法向量,由特征值计算出特征归一化向量,将不同邻域半径下的特征归一化向量差值和法向量夹角差值组合成点邻域尺度差异描述符。最后,找出目标点云中与源点云关键点的特征描述符欧氏距离最近的点,搜索该点的k近邻特征描述符集合,计算集合中每点与源点云中相关关键点的特征描述符欧式距离,并设定阈值筛选,满足阈值的邻域点加入到对应点对集E。从对应点对集E中随机选取n组对应点对求出初始对应关系,再通过基于全局距离的最优查找算法优化对应关系。实验结果表明,本发明基于形状指数选取的关键点特征明显,分布均匀,设计的特征描述符辨识性强且对应关系查找算法提高了配准的准确率。

Claims (3)

1.一种点邻域尺度差异描述的点云配准方法,其特征在于,具体的实现步骤为:
步骤1:关键点查找,对源点云和目标点云中离散点进行最小二乘曲面拟合得到局部曲面,求出曲面的形状指数SI,即该离散点的形状指数,选取形状指数在邻域内最大或最小且满足阈值的点作为点云的关键点;
步骤2:特征描述符计算,进行特征描述符构造,计算关键点在不同邻域半径下的特征归一化向量差值和法向量夹角差值,组合成点邻域尺度差异描述符;
步骤2.1:对每个关键点p,取l个大小不同的邻域半径,相邻邻域半径之差为△r,有r1<r2<…rl;对于半径为rl的邻域,建立加权协方差矩阵Cl
其中,pi为关键点p的半径为rl的邻域中的点;
步骤2.2:对加权协方差矩阵Cl进行特征值分解得到特征值{λl1,λl2,λl3}及对应的特征向量{vl1,vl2,vl3},由此建立采样点的局部邻域参考坐标系,三个特征向量构成坐标轴的三个方向;其中,λl1>λl2>λl3
步骤2.3:对于协方差矩阵Cl,构建特征值归一化向量hl
根据不同邻域半径下的hl差异计算Δhk=Δhk+1-Δhk,将Δhk组合后得到向量H,H=(Δh1Δh2…Δhl-1);其中,k=1,2,...,l-1;
步骤2.4:对关键点p选取了l个大小不同的邻域半径后,得到关键点p不同尺度下的l个法向量n1,n2,...,nl,通过法向量点乘的方法计算得到两个法向量之间的夹角余弦,连续不同尺度下的法向量的夹角余弦αk为:
计算出连续不同尺度下的法向量的夹角余弦αk,将αk组合后得到向量N,N=(α1,α2,…αl-1);
步骤2.5:对每个关键点,将H和N组合为向量f,作为该关键点的局部特征描述符;
f=[Δh1Δh2…Δhl-1α1α2…αl-1]
步骤2.6:对于源点云和目标点云的关键点集和/>中的每一个点计算其特征描述符,得到源点云关键点特征描述子集/>和目标点云关键点特征描述符集
步骤3:对应关系估计,根据特征描述符的相似程度找出对应点并使用二重筛选和基于全局距离的最优查找算法分别滤除错误点对和估计对应关系。
2.根据权利要求1所述的一种点邻域尺度差异描述的点云配准方法,其特征在于,所述步骤1的具体步骤为:
步骤1.1:对点云每一点使用最小二乘法进行曲面拟合,得到离散点邻域的局部曲面z=r(u,v);
步骤1.2:在得到拟合的二次参数曲面后,由二次参数曲面的第一类基本量E、F、G和第二类基本量L、M、N计算得到曲面的高斯曲率K和平均曲率H
由高斯曲率和平均曲率计算得到曲面的最大主曲率k1和最小主曲率k2
其中E=ru·ru,F=ru·rv,G=rv·rv,L=ruu·n,N=rvv·n,M=ruv·n,其中ru、rv、ruu、rvv、ruv是拟合曲面的偏微分,n是曲面的单位法向量,曲面P的形状指数定义为
步骤1.3:计算每点对应曲面Pi的形状指数SI(Pi);
步骤1.4:选择点云的关键点,源点云pi邻域内所有点的形状指数均值为μ(pi),满足下面两个条件中任意一个的都选作关键点
且SI(pi)≥(1+a)μ(pi)
且SI(pi)≤(1-b)μ(pi)
其中a、b是用来调整关键点数量的两个标量参数,将选取后的关键点集记为
3.根据权利要求1或2所述的一种点邻域尺度差异描述的点云配准方法,其特征在于,所述步骤3的具体步骤为:
步骤3.1:源点云P和目标点云Q的关键点特征描述符集记为对于/>中的每个/>在/>中寻找与其特征描述符欧式距离最近的将/>记为候选对应点对;
步骤3.2:搜索的k近邻特征描述符集合/>定义阈值η,当/>时,认为pi和qi1是一组对应点,记为对应点对ei={pi,qi1};若不满足,将/>剔除,遍历关键点集/>最终得到对应点对集记为E={e1,e2,...en};
步骤3.3:从对应点对集中随机抽选β组对应点对并计算变换矩阵Ti,根据Ti将源点云P变换到新的坐标系下为计算全局距离/>
其中,β≥3;P={pi|i=1,2...,N},Q={qj|j=1,2...,M}分别表示源点云和目标点云,表示在点云Q中的点/>与其相对应的欧式距离/>小于阈值Dε的点云P中的点,/>这样的点的总数,ξ为可调参数;
步骤3.4:设定迭代停止阈值Dε,当时停止迭代,Ti被选取为最优变换,如果不满足/>则再次抽选新的对应点对子集,迭代重复计算/>直至满足或达到最大迭代次数Nmax,当达到最大迭代次数时,选取/>最小的变换矩阵为最优变换矩阵。
CN201911276127.5A 2019-12-12 2019-12-12 一种点邻域尺度差异描述的点云配准方法 Active CN111080684B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911276127.5A CN111080684B (zh) 2019-12-12 2019-12-12 一种点邻域尺度差异描述的点云配准方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911276127.5A CN111080684B (zh) 2019-12-12 2019-12-12 一种点邻域尺度差异描述的点云配准方法

Publications (2)

Publication Number Publication Date
CN111080684A CN111080684A (zh) 2020-04-28
CN111080684B true CN111080684B (zh) 2023-11-17

Family

ID=70314205

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911276127.5A Active CN111080684B (zh) 2019-12-12 2019-12-12 一种点邻域尺度差异描述的点云配准方法

Country Status (1)

Country Link
CN (1) CN111080684B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113850748B (zh) * 2020-06-09 2024-07-26 上海交通大学 点云质量的评估系统及方法
CN111860319B (zh) * 2020-07-20 2024-03-26 阿波罗智能技术(北京)有限公司 车道线的确定方法、定位精度的评测方法、装置、设备
CN112150523B (zh) * 2020-09-24 2022-04-12 中北大学 一种低重叠率的三维点云配准方法
CN112652003B (zh) * 2020-10-29 2024-04-12 西北工业大学 一种基于ransac测度优化的三维点云配准方法
CN112529945B (zh) * 2020-11-17 2023-02-21 西安电子科技大学 一种多视角三维isar散射点集配准方法
CN112489094B (zh) * 2020-11-25 2023-08-04 国网福建省电力有限公司 一种多线机载激光雷达点云数据配准方法
CN112418250B (zh) * 2020-12-01 2024-05-10 怀化学院 一种复杂3d点云的优化匹配方法
CN112488127B (zh) * 2020-12-09 2022-08-09 南京邮电大学 基于层次b样条的多尺度点云特征点检测与匹配
CN114693881A (zh) * 2020-12-30 2022-07-01 合肥美亚光电技术股份有限公司 点云数据的处理方法、配准装置和口内扫描设备
CN112802077B (zh) * 2021-01-12 2023-08-22 广东奥普特科技股份有限公司 一种高精度点云法向计算方法
CN113393504A (zh) * 2021-05-28 2021-09-14 广东电网有限责任公司广州供电局 一种开关柜点云配准方法、系统、计算机设备及存储介质
CN113239500B (zh) * 2021-07-12 2021-09-21 四川大学 基于协方差矩阵的参考点邻域特征匹配方法
CN114742867A (zh) * 2022-04-04 2022-07-12 西北大学 一种文物碎片拼接方法、系统、介质、设备及终端
CN114972459B (zh) * 2022-05-31 2023-09-08 哈尔滨理工大学 一种基于低维点云局部特征描述符的点云配准方法
CN116452648B (zh) * 2023-06-15 2023-09-22 武汉科技大学 一种基于法向量约束纠正的点云配准方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104143210A (zh) * 2014-07-31 2014-11-12 哈尔滨工程大学 一种多尺度法向特征点云配准方法
CN110473239A (zh) * 2019-08-08 2019-11-19 刘秀萍 一种三维激光扫描的高精度点云配准方法
CN110490912A (zh) * 2019-07-17 2019-11-22 哈尔滨工程大学 基于局部灰度顺序模型描述符的3d-rgb点云配准方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104978549B (zh) * 2014-04-03 2019-04-02 北京邮电大学 三维人脸图像特征提取方法和系统
CN108604301B (zh) * 2015-12-04 2022-03-25 欧特克公司 用于大rgb-d扫描的可缩放自动全局配准的基于关键点的点对特征
WO2017113260A1 (zh) * 2015-12-30 2017-07-06 中国科学院深圳先进技术研究院 一种三维点云模型重建方法及装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104143210A (zh) * 2014-07-31 2014-11-12 哈尔滨工程大学 一种多尺度法向特征点云配准方法
CN110490912A (zh) * 2019-07-17 2019-11-22 哈尔滨工程大学 基于局部灰度顺序模型描述符的3d-rgb点云配准方法
CN110473239A (zh) * 2019-08-08 2019-11-19 刘秀萍 一种三维激光扫描的高精度点云配准方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
汪霞 ; 赵银娣 ; 王坚 ; .一种低重叠率激光点云的配准方法.测绘科学.2018,(12),全文. *
陆军 ; 彭仲涛 ; .基于快速点特征直方图的特征点云迭代插值配准算法.国防科技大学学报.2014,(06),全文. *
陆军 ; 方莹 ; 夏桂华 ; .基于曲面拟合系数特征的点云快速拼接算法.光电子・激光.2015,(09),第4-5页. *

Also Published As

Publication number Publication date
CN111080684A (zh) 2020-04-28

Similar Documents

Publication Publication Date Title
CN111080684B (zh) 一种点邻域尺度差异描述的点云配准方法
CN109887015B (zh) 一种基于局部曲面特征直方图的点云自动配准方法
CN108376408B (zh) 一种基于曲率特征的三维点云数据快速加权配准方法
CN114170279B (zh) 一种基于激光扫描的点云配准方法
CN110599506B (zh) 一种复杂异形曲面机器人三维测量的点云分割方法
CN111553409B (zh) 一种基于体素形状描述符的点云识别方法
CN109903319B (zh) 一种基于多分辨率的快速迭代最近点配准算法
CN105118059A (zh) 一种多尺度坐标轴夹角特征的点云快速配准方法
CN114972459B (zh) 一种基于低维点云局部特征描述符的点云配准方法
CN105551015A (zh) 一种散乱点云图像配准方法
CN109523582B (zh) 一种顾及法向量和多尺度稀疏特征的点云粗配准方法
CN114677418B (zh) 一种基于点云特征点提取的配准方法
CN113327275B (zh) 一种基于多约束点到局部曲面投影的点云双视角精配准方法
CN105654483A (zh) 三维点云全自动配准方法
CN109299301B (zh) 一种基于形状分布和曲度的三维模型检索方法
CN113344019A (zh) 一种决策值选取初始聚类中心改进的K-means算法
CN115861397A (zh) 一种基于改进fpfh-icp的点云配准方法
CN113961738B (zh) 一种多特征铸件三维模型检索方法及装置
CN115797418B (zh) 一种基于改进icp的复杂机械零件测量点云配准方法及系统
CN108921853B (zh) 基于超像素和免疫稀疏谱聚类的图像分割方法
CN106844620A (zh) 一种基于视图的特征匹配三维模型检索方法
CN111862177A (zh) 一种基于方向直方图签名特征的工件的三维点云配准方法
CN117132630A (zh) 一种基于二阶空间兼容性度量的点云配准方法
CN110503148B (zh) 一种具有尺度不变性的点云对象识别方法
CN117274339A (zh) 一种基于改进的iss-3dsc特征结合icp的点云配准方法

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