CN109447100A - 一种基于b样条曲面相似性检测的三维点云识别方法 - Google Patents
一种基于b样条曲面相似性检测的三维点云识别方法 Download PDFInfo
- Publication number
- CN109447100A CN109447100A CN201811001394.7A CN201811001394A CN109447100A CN 109447100 A CN109447100 A CN 109447100A CN 201811001394 A CN201811001394 A CN 201811001394A CN 109447100 A CN109447100 A CN 109447100A
- Authority
- CN
- China
- Prior art keywords
- point
- similitude
- equidistant
- point cloud
- feature
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- 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
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V40/00—Recognition of biometric, human-related or animal-related patterns in image or video data
- G06V40/10—Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
- G06V40/16—Human faces, e.g. facial parts, sketches or expressions
- G06V40/172—Classification, e.g. identification
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Probability & Statistics with Applications (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Human Computer Interaction (AREA)
- Multimedia (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及一种基于B样条曲面相似性检测的三维点云识别方法,其技术特点在于:包括以下步骤:步骤1、利用B样条曲面方程对三维点云进行建模和拟合;步骤2、采用不同参数采样方法对拟合后的B样条曲面进行采样;步骤3、通过定义一种描述点的几何局部特征计算得出相似点,并进行点对匹配;步骤4、对计算出的点对的特征,进行等距计算;步骤5、对点对内两点的等距变换类型进行分类;步骤6、对等距分类后的点对之间进行等距比较;步骤7、应用谱聚类算法,得出物体之间局部对应相似部分,完成相似区域的检测;步骤8、通过统计相似点云和采样点云的比例是否大于设定阈值的判定来完成整体点云的识别。本发明实现了高效准确的物体识别。
Description
技术领域
本发明属于图像处理和图形变换技术领域,涉及图形处理和相似性检测 的三维点云识别方法,尤其是一种基于B样条曲面相似性检测的三维点云识 别方法。
背景技术
目前,基于特征提取技术探究高效鲁棒性的物体形状的判别算法是机器 视觉领域解决物体识别问题的核心思路。目前,基于识别的特征提取方法主 要有两大类:基于全局特征的方法和基于局部特征的方法。在这些方法中, 往往识别精度高的算法,往往鲁棒性较差,应用场景往往受限。在点云数据 预处理中,步骤繁琐,导致识别效率的降低。还有很多识别方法虽然都能够 解决各自的识别问题,但是并不具备完成识别的前提下,对物体表面的几何 形状精准描述的特点。
经检索,未发现与本发明相同或相似的已公开的专利文献。
发明内容
本发明的目的在于克服现有技术的不足,提供一种设计合理、识别速度 快且识别精度高的基于B样条曲面相似性检测的三维点云识别方法。
本发明解决其技术问题是采取以下技术方案实现的:
一种基于B样条曲面相似性检测的三维点云识别算法,包括以下步骤:
步骤1、利用B样条曲面方程对三维点云进行建模和拟合;
步骤2、采用不同方法对步骤1拟合后的B样条曲面进行采样和采样分 析,评估不同采样方法对点对特征计算结果的影响;
步骤3、利用步骤2采样得到的点云,通过定义了一种点云的局部几何 特征计算得出相似点,并进行点对匹配;
步骤4、对步骤3计算出的点对,对其进行等距计算;
步骤5、对步骤4得到计算的等距进行分类;
步骤6、对步骤5中等距分类后的点对之间进行等距比较的计算;
步骤7、对步骤6得到的相同等距距离的点对组,应用谱聚类算法,得 出物体之间局部对应相似部分,完成相似区域检测;
步骤8、根据步骤7得到的不同等距分类下的相似点云,对给出的不同 参数采样方法,通过统计相似点云和采样点云的比例是否大于设定阈值的判 定来完成整体点云的识别
而且,所述步骤1的具体步骤包括:
(1)利用B样条曲面方程对三维点云进行建模;
对K次B样条参数曲面给定控制点Cij,(i=0,1,2...m;j=0,1,2....n),基函数递推公式得到:
规定了0/0=0,节点矢量U=Ui(i=0,1,2......m+k+1),V=Vj(j=0,1,2......n+k+1) 可以得到曲面方程:
(2)对生成的B样条曲面进行拟合:
反求控制点,令
①固定j,对pi,j(j=0,1,.........n)沿u方向应用曲线反算,分别求出n+1条等参数曲线(截面曲线)的控制顶点;
②固定i,对bi,j(i=0,1,........m)沿v方向应用曲线反算,分别求出m+1条插值曲线的控制顶点bi,j(i=0,1,......m;n=0,1,.......j),即为双三次B样条插值曲面的控制顶点。
而且,所述步骤2的采用不同方法对拟合后的B样条曲面进行采样的具 体方法为:在3D对象的表面上采样一组点,通过统一取样两个参数u和v, 按照参数区间均匀化采样取得的样本或采用分割参数化的采样取得的样本。
而且,所述步骤3的具体步骤包括:
(1)定义参数曲面的第一个基本形式:
ds2=Edu2+2Fdudv+Gdv2
其中,E=Pu·Pu,F=Pu·Pv,G=Pv·Pv.
第一个基本形式表明,对于一个给定点P,偏导数Pu和Pv产生一个与原点 P表面相切的平面;
(2)得出统一的法向量是:
该统一的法向量与非归一化向量Pu,Pv形成一个局部向量坐标系;
(3)定义参数曲面的第二个基本形式:
κcosφds2=Ldu2+2Mdudv+Ndv2
其中,L=Puu·n,M=Puv·n,N=Pvv·n,;
式中,Puu,Puv,Pvv都是定点P的二阶偏导;
参数曲面的第二个基本形式公式意味着对于给定方向du/dv在u,v平面和 给定的角度φ,通过第一和第二基本形式,计算曲面上描绘的曲线的曲率K, 也就是切线指向;
(4)引入了两个对称矩阵:
因为Pu和Pv是线性独立的,所以矩阵F1 -1F2也是对称的,因此可以解出 F1 -1F2的实特征值和正交特征向量;两个特征值K1和K2就是两个主曲率和两 个特征向量t1=(ζ1,η1)T,t2=(ζ2,η2)T定义了两个主方向:
t1=ζ1Pu+η1Pv;t2=ζ2Pu+η2Pv
通过以上对法向量n和主方向t1,t2的定义和计算,对任意给定点P的特征 描述表示为向量形状,P—(n,t1,t2)
在匹配过程中由于一点的主曲率决定了两个主方向,因此通过匹配主曲率的 方式来配对点对,这里不考虑脐点(两个主曲率相等的点)提高了匹配的质 量和效率。对于待测两点pi和pj的主曲率k1,k2和k1’,k2’的匹配过程中, 满足下列判定:
(k1-k1'<r1)&&(k2-k2'<r2)
其中r1,r2分别为满足两个主曲率近似相等的最大差值。
而且,所述步骤4的具体计算方法为:计算等距就是配对的(pi,pj) 计算pi的局部特征向量到pj的局部特征向量的转换;
用Eij表示变换的旋转,tij表示相应的平移变换,计算如下:
tij=pj-Eij*pi
因此存在Tij:pi(ni,t1 i,t2 i)→pj(nj,t1 j,t2 j),等距变换下的空间集合,而Tij就是一次等距变换。
而且,所述步骤5的具体方法为:
给定存在一个唯一的变换其中g是一个具有非空集合的 等距变换,这里是在固定点A下的平移,如果T=g,并且d=0,那么至少有一 个固定点;如果将T规定为直接等距;那么把T规定 为间接等距。
而且,所述步骤6的等距比较具体计算方法是:
(1)直接等距的点对等距距离:
dist(T,T’)=(1-|cos(D·D’)|)+|θ-θ’/2π|+ω1dist(D,D’)+ω2|(||t||-||t’||)
(2)间接等距的点对等距距离:
而且,所述步骤7的具体方法为:令dist(T,T’)与两个点对(pi,pj)和(pi’,pj’)对应的同一类别中的两个等距T和T'之间的距离;dist(T,T’)以及相似度度量将 根据等距的类别而变化;令S={(pi,pj)}{i,j=1,2,.....Nl}∈Γl,l∈[|1,5|]为第l个等 距空间中的Nl个点对的集合,令k为聚类数。
而且,所述步骤7的具体步骤包括:
(1)输入样本相似矩阵S和聚类数K;
(2)选择高斯相似度函数来构造权重矩阵A∈RNl×k和度矩阵D:
(3)构造标准化的拉普拉斯矩阵解出L矩阵的K个最大特征 向量构造矩阵X=[X1,X2,......Xk]∈RNl×k;
(4)通过X矩阵的按行标准化构造成特征Y矩阵,用K-means算法将它 们分成K个聚类,当且仅当矩阵Y的第i行被分配给类θ时,将原始点对(pi,pj) 分配给类别θ;
(5)确定高斯相似度参数σ
高斯相似度函数中定义的描述相似度公式取决于参数σ,定义了点对 (pi,pj)之间的变换距离的阈值;为了设置它,需要考虑点的均匀分布;定义各 向分布的元素被包括在大小为Dmax=maxdist(T,T’)的边界框中,通过将这个盒子 分成Nl个相同的体积,两个点之间的(均匀的)距离Dun:
Dun=max(pi,pj)≠(pi',pj')d(T,T')/Nl 1/m
其中,m是等距空间的维度,根据等距性质而变化;如果一个聚类存在, 那么就有一些点被小于Dun距离的点分开;类似地,高斯参数σ被用作距离 的一部分σ=Dun/2;因此,将m维空间中的点集成在一起,并导出了待定点 的阈值封闭;
(6)确定最优聚类数K
通过利用高斯亲和矩阵A并基于Frobeniu范数的比例定义一个质量度量; 设αk是要识别的簇的数量的界限;对于一个值k'∈[|2,αk|],亲和矩阵是每 个群集的索引;然后定义块矩阵:非对角块表示簇之间的亲和度,对角块表 示簇内的亲和度;从这个结构中,可以评估所有非对角线区块与Frobenius范数 的对角块之间的平均比率,称为rk';由此,在K'属于[|2,αk|]的值中,比率rk'的最小值定义了最优类数K:
K=arg min rk'k'∈[|2,αk|]
这个最小值是在集群内的高斯亲和度最高情况下,对应于集群之间的亲 和力最低的值。
而且,所述步骤8的具体步骤包括:
(1)通过谱聚类算法可以得到等距分类中物体之间相似部分点对的数量 Ni(i∈[1,5]),其中计算后得到的点云对数Ni(i∈[1,5])对应着被划分为5个等距 空间中的对应相似点云,采样到的点云对数为Nu;
(2)对于整体参数的统一采样,定义相似点对数与采样点云对数比值λ:
如果λ≥实验阈值η,说明验证了两个点云是相似的;认为是相同的对象, 否则它们不相似,就是不同的对象;
(3)局部分割式的参数采样,对应分割的不完整点云多次进行判定,定 义比例阈值ε:
ε=Ni/Ns
其中,Ni是相似的局部点云数量,Ns是完整点云数量;若ε≥η阈值, 那么就认为是相同,即完成目标之间的识别。
本发明的优点和积极效果是:
1、本发明提供一种基于B样条曲面相似性检测的三维点云识别方法,对 复杂场景中物体的识别问题,实现三维点云表示的物体高精度识别。该方法 的整体思路是通过对B样条参数曲面物体表面的几何形状的精准描述,根据 匹配后的点对,对其特征进行等距变换的计算,将代表点对内两点特征之间 的等距变换进行分类,比较在同类等距变换下,点对之间特征的变换的差异, 利用谱聚类算法对具有相同近似等距变换的点对进行相似度度量,得出相似 点对群。最后统计在不同等距分类下的相似点对的点云数量,判断相似点云 与采样点云数的比值是否大于设定的阈值来完成高效率的识别。
2、本发明步骤1至3,通过对B样条参数曲面的拟合,优化了点云模型, 更好地表现了物体表面的几何形状,定义了一种对局部形状描述的特征,特 征由该点的主曲率和主方向和法向量组成,通过计算其特征并匹配点对,有 效筛选出高质量的特征匹配。
3、本发明步骤4至步骤7通过计算和分类点对之间的等距变换,验证了 目标之间相似部分的转换,最后统计在不同等距分类下的相似点对的点云数 量,判断相似点云与采样点云比值是否大于设定阈值来完成识别,提高了识 别的效率和准确率。
3、通过统计不同等距分类下的相似点云数,进行相似点云与采样点云比 值的分析和判断来完成识别。
4、本发明的步骤7采用无监督的谱聚类方法来改进计算相似度的结果, 相较于传统的采用欧式距离计算相似度的方法,精度更高。
附图说明
图1是本发明的处理流程图;
图2是本发明的步骤3的处理流程图;
图3是本发明的步骤3的两个相似点pi和pj的局部向量组成的特征示意 图;
图4(a)是本发明的步骤5的第一种等距变换分类示意图;
图4(b)是本发明的步骤5的第二种等距变换分类示意图;
图4(c)是本发明的步骤5的第三种等距变换分类示意图;
图5是本发明的人脸库数据样本示意图;
图6是本发明的步骤8中在数据集测试中的识别效果示意图。
具体实施方式
以下结合附图对本发明实施例作进一步详述:
一种基于B样条曲面相似性检测的三维点云识别算法,包括以下步骤:
步骤1、利用B样条曲面方程对三维点云进行建模和拟合;
所述步骤1的具体步骤包括:
(1)利用B样条曲面方程对三维点云进行建模;
对K次B样条参数曲面给定控制点Cij,(i=0,1,2...m;j=0,1,2....n),基函数递推公式得到:
规定了0/0=0,节点矢量U=Ui(i=0,1,2......m+k+1),V=Vj(j=0,1,2......n+k+1)
可以得到曲面方程:
其中均匀节点矢量和准均匀节点矢量与曲线取法相同,非均匀依赖于数 据点布,对于它们两端节点的重复度均取为r=4.准均匀双三次B样条插值曲 面两个参数方向的节点矢量与数据点的分布情况无关,其取法与曲线相同.非 均匀双三次B样条插值曲面两个参数方向的节点矢量要由数据点确定.曲面 沿同一参数方向要求有一个共同的节点矢量,但该方向各行(列)数据点的分 布一般却不相同,因而由各行(列)数据点确定的节点矢量也就不同.一个共同 的节点矢量通常采用取算术平均值的法确定某方向一个共同的节点参数。
(2)对生成的B样条曲面进行拟合:
反求控制点,令
①固定j,对pi,j(j=0,1,.........n)沿u方向应用曲线反算,分别求出n+1条等参数曲线(截面曲线)的控制顶点;
②固定i,对bi,j(i=0,1,........m)沿v方向应用曲线反算,分别求出m+1条插值曲线的控制顶点bi,j(i=0,1,......m;n=0,1,.......j),即为双三次B样条插值曲面的控制顶点。
步骤2、采用不同方法对步骤1拟合后的B样条曲面进行采样和采样分 析,评估不同采样方法对点对特征计算结果的影响;
在本实施例中,所述步骤2的采用不同方法对拟合后的B样条曲面进行 采样的具体方法为:在3D对象的表面上采样一组点,通过统一取样两个参数 u和v,按照参数区间均匀化采样取得的样本或采用分割参数化的采样取得 的样本。
所述步骤2的工作原理是:
参数曲面表面上的每个点都有以局部差分特性为特征的特点,可能会被 抽样用于以后的计算,通过统一取样两个参数u和v,由于模型中曲面之间 的u和v参数变化往往不同,按照参数变化统一的规则选择来决定有限数量的 样本,按照参数区间均匀化采样取得的样本能够比较全面描述初始样本。分 割参数化的采样,有利于计算精度。不同采样方式直接影响相似性计算的步 骤和结果,因此,对于参数一致的方向,要根据实际情况来选择。
步骤3、利用步骤2采样得到的点云,通过定义了一种点云的几何局部 特征计算得出相似点,并进行点对匹配;
所述步骤3的具体步骤包括:
计算特征的过程中,每个点的特征由两个主曲率和一个原点在该点的法 向量组成,单位向量是法向量和两个主方向(即与主曲率相关的切向量).在 B样条曲面上特定点的特征计算是基于局部微分特性,因此,可以从第一个 和第二个基本形式进行评估。
(1)定义参数曲面的第一个基本形式:
ds2=Edu2+2Fdudv+Gdv2 (5)
其中E=Pu·Pu,F=Pu·Pv,G=Pv·Pv.
第一个基本形式表明,对于一个给定点P,偏导数Pu和Pv产生一个与原点 P表面相切的平面;
(2)得出统一的法向量是:
该统一的法向量与非归一化向量Pu,Pv形成一个局部向量坐标系;
(3)定义参数曲面的第二个基本形式:
κcosφds2=Ldu2+2Mdudv+Ndv2 (7)
L=Puu·n,M=Puv·n,N=Pvv·n,Puu,Puv,Pvv都是定点P的二阶偏导。
公式(7)意味着对于给定方向du/dv在u,v平面和给定的角度φ,通过第 一和第二基本形式,我们计算曲面上描绘的曲线的曲率K,也就是切线指向。
(4)引入了两个对称矩阵:
因为Pu和Pv是线性独立的,所以矩阵F1 -1F2也是对称的,因此可以解出 F1 -1F2的实特征值和正交特征向量;两个特征值K1和K2就是两个主曲率和两 个特征向量t1=(ζ1,η1)T,t2=(ζ2,η2)T定义了两个主方向:
t1=ζ1Pu+η1Pv;t2=ζ2Pu+η2Pv (8)
这里不考虑脐点(即两个主曲率相等的点)。
通过以上对法向量n和主方向t1,t2的定义和计算,对任意给定点P的 特征描述表示为向量形状,P—(n,t1,t2)
在匹配过程中由于一点的主曲率决定了两个主方向,因此通过匹配主曲 率的方式来配对点对,这里不考虑脐点(两个主曲率相等的点)提高了匹配 的质量和效率。对于待测两点pi和pj的主曲率k1,k2和k1’,k2’的匹配过 程中,满足下列判定:
(k1-k1'<r1)&&(k2-k2'<r2)
其中r1,r2分别为满足两个主曲率近似相等的最大差值。
所述步骤3的工作原理是:在3D对象的表面上采样一组点,这里不同物 体采样的点云数量保证是相等的。由于点位置不足以确定一般变换的特征, 所以需要计算每个样本的几何特征(主曲率和由主方向和法向矢量组成的局 部坐标系)。特征是两主曲率,如果它们有相同的特征,将它们配对成相似 点对,每一对点都与转换空间中的投票对应的变换相关联。给定两个点pi和pj和它们的局部向量组(标准正交)由两个切线和一个法线组成的(图1),将变 换Tij进行计算,使pi和它的局部坐标系被映射到pi和pj的向量坐标系中。这 个变换被转换成转换空间的投票。
如图3所示,本发明给出了两个相似点pi和pj的局部特征示意图,通过 步骤3的特征计算,说明两个点是否相似和是否匹配。
步骤4、对步骤3计算出的点对,进行点对内的等距计算;
在本实施例中,所述步骤4的点对的等距计算具体方法是:在变换空间 中,每个等距变换表示两个相似点特征之间的变换。因此,可以识别类似转换 的点云,因为它们可以描述对象的两个相似部分。
计算等距就是配对的(pi,pj)计算pi的局部特征向量到pj的局部特征 向量的转换。用Eij表示变换的旋转,tij表示相应的平移变换,计算如下:
tij=pj-Eij*pi (10)
因此存在Tij:pi(ni,t1 i,t2 i)→pj(nj,t1 j,t2 j),等距变换下的空间集合,而Tij就是一次等距变换。
步骤5、对步骤4得到计算的等距进行分类;
在本实施例中,所述步骤5的具体方法为:
给定存在一个唯一的变换其中g是一个具有非空集合的 等距变换,这里是在固定点A下的平移,。如果T=g,并且d=0,那么至少有 一个固定点。如果将T规定为直接等距;那么把T 规定为间接等距。
在本实施例中,等距分类按如下规则划分:
直接等距:1)存在一个固定点的轴线D:等距T就是一个围绕轴线D 的旋转;2)没有固定点:等距变换T就是d的非零平移或由向量a和非零 平移共线所指向的轴线D。
如图4(a)、图4(b)和图4(c)所示,间接等距:1)一个特点的固 定点A:由一个围绕一个轴线D的旋转构成,该轴线垂直平面G。2)一个固 定点的平面G,a=常数C,向量a=0:等距T是一个相对于由向量n定义 的平面G的对称。3)没有固定点:等距T由相对于平面G的对称性构成, 平面G的法向量为n,平行于该平面的非零平移。
步骤6、对步骤5中等距分类后的点对之间进行等距比较的计算;
对等距变换的划分表达了一个适当的维度转换,明确等距的分类;
在本实施例中,所述步骤6的等距比较具体计算方法是:
通过将等距变换的分类,我们最多有五个不同的等距变换,并且每个等距 转换空间都将应用一个聚类算法。聚类需要在每个空间中都有一定的距离, 可以得到两个相同类型的等距的距离。
(1)直接等距的点对等距距离:
dist(T,T’)=(1-|cos(D·D’)|)+|θ-θ’/2π|+ω1dist(D,D’)+ω2|(||t||-||t’||)| (11)
(3)间接等距的点对等距距离:
对于直接等距,是由旋转轴(D)和角度θ,平移组成。由于旋转轴和平 移方向相同,因此平移矢量d和轴上的点P就足够了。为了比较旋转,我们 使用角度和两个轴之间的距离,以及角度的差异;对于平移,我们仍然比较平 移向量的长度。对于间接等距,除了对称平面G之外,分析与直接等距是相 同的。平面之间的比较包括比较这些平面的法线并计算中点和平面之间的距 离。dist(T,T’)对应两个点对(pi,pj)和(pi’,pj’)的两个等距T和T'之间的距 离;Mij,Mi'j'是(pi,pj)和(pi',pj')的中点;dist(p,G)表示从点,线或平面到另一个 点的距离
步骤7、对步骤6得到的相同等距距离的点对组,应用谱聚类算法,得 出物体之间局部对应相似部分,完成相似区域检测。
在本实施例中,所述步骤7的具体方法为:令dist(T,T’)与两个点对(pi,pj)和(pi’,pj’)对应的同一类别中的两个等距T和T'之间的距离。dist(T,T’)以及相似度 度量将根据等距的类别而变化。令S={(pi,pj)}{i,j=1,2,.....Nl}∈Γl,l∈[|1,5|]为 第l个等距空间中的Nl个点对的集合,令k为集群数量。
在本实施例中,所述步骤7的具体步骤包括:
(1)输入样本相似矩阵S和聚类数K;
(2)选择高斯相似度函数来构造权重矩阵A∈RNl×k和度矩阵D:
(3)构造标准化的拉普拉斯矩阵解出L矩阵的K个最大特征 向量构造矩阵
(4)通过X矩阵的按行标准化构造成特征Y矩阵,用K-means算法将它 们分成K个聚类,当且仅当矩阵Y的第i行被分配给类θ时,将原始点对(pi,pj) 分配给类别θ。
在整个算法的过程中,仅依赖于高斯相似度参数σ和类别数K两个参数, 为了使这个方法完全不受监督,我们要确定这两个参数。
(5)确定高斯相似度参数σ
高斯相似度函数中定义的描述相似度公式取决于参数σ,定义了点对 (pi,pj)之间的变换距离的阈值。为了设置它,需要考虑点的均匀分布。定义各 向分布的元素被包括在大小为Dmax=maxdist(T,T’)的边界框中,通过将这个盒子 分成Nl个相同的体积,两个点之间的(均匀的)距离Dun:
Dun=max(pi,pj)≠(pi',pj')d(T,T')/Nl 1/m (14)
其中,m是等距空间的维度,根据等距性质而变化。我们可以认为,如果 一个聚类存在,那么就有一些点被小于Dun距离的点分开。类似地,高斯参 数σ被用作距离的一部分σ=Dun/2;因此,将m维空间中的点集成在一起, 并导出了待定点的阈值封闭。
(6)确定最优聚类数K
簇的数量的选择是所有无监督聚类算法的一个普遍问题。为了确定聚类 数K,通过利用高斯亲和矩阵A并基于Frobeniu范数的比例定义一个质量度量。 设αk是要识别的簇的数量的界限。对于一个值k'∈[|2,αk|],亲和矩阵是每 个群集的索引。然后定义块矩阵:非对角块表示簇之间的亲和度,对角块表 示簇内的亲和度。从这个结构中,我们可以评估所有非对角线区块与Frobenius 范数的对角块之间的平均比率,称为rk'。由此,在K'属于[|2,αk|]的值中,比 率rk'的最小值定义了最优类数K:
K=arg min rk'k'∈[|2,αk|] (15)
这个最小值是在集群内的高斯亲和度最高情况下,对应于集群之间的亲 和力最低的值。
步骤8、根据步骤7得到的相似点云,对步骤2中给出的不同参数采样 方法,通过对比例阈值的识别判定来完成整体点云的识别,并进一步给出了 体现算法特点的识别效果。
在本实施例中,所述步骤8的具体步骤包括:
(1)通过谱聚类算法可以得到等距分类中物体之间相似部分点对的数量 Ni(i∈[1,5]),其中计算后得到的点云对数Ni(i∈[1,5])对应着被划分为5个等距 空间中的对应相似点云,采样到的点云对数为Nu。
(2)对于整体参数的统一采样,定义相似点对数与采样点云对数比值λ:
如果λ≥实验阈值η,说明验证了两个点云是相似的。认为是相同的对 象,否则它们不相似,就是不同的对象。最终可以识别采样到的局部点云是 否相同,由于采样的不同方式会影响识别速度和精度,这里统一按照B样条 的两个参数方向均匀等长采样,可以得到相对完整的点云轮廓,提高了识别 速度。
(3)局部分割式的参数采样,对应分割的不完整点云多次进行判定,定 义比例阈值ε:
ε=Ni/Ns (17)
其中,Ni是相似的局部点云数量,Ns是完整点云数量。ε≥η阈值,那么就 认为是相同,即完成目标之间的识别。实验中阈值因测试数据集的不同,取 不同待选阈值下识别率最高的为最优阈值。
本发明的工作原理是:
利用自由参数曲面模型提出一种基于B样条曲面相似性检测的三维点云 识别算法,应用投票变换空间的思想,在相似性检测中计算物体相似部分的 转换,该投票变换空间与霍夫变换有相似之处:具有相似特征的模型上的点 是成对的。点对对应于两点之间的转换及其特征;这些转变投向转型空间,形 成投票空间,每一个点对之间的变换相当于投票空间中的一点。这些投票的 点对群是在模型中定义相似部分的候选人,因为它可以检测潜在转换,并且 能够识别三维物体中的相似部分和表征转换,并定义了一种对局部形状描述 的特征,特征由该点的主曲率和主方向和法向量组成,通过计算其特征并匹 配点对,有效筛选出高质量的特征匹配。根据匹配后的点对,对其特征进行 等距变换的计算,将代表点对内两点特征之间的等距变换进行分类,比较在 同类等距变换下,点对之间特征的变换的差异,利用谱聚类算法对具有相同 近似等距变换的点对进行相似度度量,得出相似点对群。最后统计在不同等 距分类下的相似点对群点云数,通过相似点云与采样点云比值是否大于设定 阈值来完成识别,提高了识别的准确率和效率。
识别和相似度计算流程如图1和图2所示,这里把步骤8中针对两种情 况统计的相似点云数记为Nr。在实验中,我们将相似点对的主曲率的初始阈 值设为0~0.1,然后逐渐调试。如图5和图6所示,在公开的TOSCA数据集 下,对人脸库中样本进行测试,结果表明对于很多物体三维点云识别方法往 往会识别出同样的形状,而我们的方法可以通过调节阈值范围来提高识别精 度,从而达到同一张人脸样本能区分出不同的区域,也意味着区分了不同的 表情,这体现了该算法精确描述形状的特点。
本发明通过在不同的公开数据集下进行同类模型样本和人脸样本的测 试,其测试结果进一步说明了本发明的算法具备能够精确描述物体形状的特 点,本发明在工业产品的检测和分类以及人脸识别等领域也有很好的应用前 景。
需要强调的是,本发明所述的实施例是说明性的,而不是限定性的,因 此本发明包括并不限于具体实施方式中所述的实施例,凡是由本领域技术人 员根据本发明的技术方案得出的其他实施方式,同样属于本发明保护的范围。
Claims (10)
1.一种基于B样条曲面相似性检测的三维点云识别方法,其特征在于:包括以下步骤:
步骤1、利用B样条曲面方程对三维点云进行建模和拟合;
步骤2、采用不同方法对步骤1拟合后的B样条曲面进行采样和采样分析,评估不同采样方法对点对特征计算结果的影响;
步骤3、利用步骤2采样得到的点云,通过定义了一种描述点的几何局部特征计算得出相似点,并进行点对匹配;
步骤4、对步骤3计算出的点对,进行点对内两点之间的等距计算;
步骤5、对步骤4根据计算后点对内的等距变换进行分类;
步骤6、对步骤5中等距分类后的点对之间进行等距比较的计算;
步骤7、对步骤6得到的相同等距距离的点对组,应用谱聚类算法,得出物体之间对应相似部分,完成物体相似区域检测;
步骤8、根据步骤7得到的不同等距分类下的相似点云,对给出的不同参数采样方法,通过统计相似点云和采样点云的比例是否大于设定阈值的判定来完成整体点云的识别。
2.根据权利要求1所述的一种基于B样条曲面相似性检测的三维点云识别方法,其特征在于:所述步骤1的具体步骤包括:
(1)利用B样条曲面方程对三维点云进行建模;
对K次B样条参数曲面给定控制点Cij,(i=0,1,2...m;j=0,1,2....n),基函数递推公式得到:
规定了0/0=0,节点矢量U=Ui(i=0,1,2......m+k+1),V=Vj(j=0,1,2......n+k+1)可以得到曲面方程:
(2)对生成的B样条曲面进行拟合:
反求控制点,令
①固定j,对pi,j(j=0,1,.........n)沿u方向应用曲线反算,分别求出n+1条等参数曲线(截面曲线)的控制顶点;
②固定i,对bi,j(i=0,1,........m)沿v方向应用曲线反算,分别求出m+1条插值曲线的控制顶点bi,j(i=0,1,......m;n=0,1,.......j),即为双三次B样条插值曲面的控制顶点。
3.根据权利要求1所述的一种基于B样条曲面相似性检测的三维点云识别方法,其特征在于:所述步骤2的采用不同方法对拟合后的B样条曲面进行采样的具体方法为:在3D对象的表面上采样一组点,通过统一取样两个参数u和v,按照参数区间均匀化采样取得的样本或采用分割参数化的采样取得的样本。
4.根据权利要求1所述的一种基于B样条曲面相似性检测的三维点云识别方法,其特征在于:所述步骤3的具体步骤包括:
(1)定义参数曲面的第一个基本形式:
ds2=Edu2+2Fdudv+Gdv2
其中,E=Pu·Pu,F=Pu·Pv,G=Pv·Pv.
第一个基本形式表明,对于一个给定点P,偏导数Pu和Pv产生一个与原点P表面相切的平面;
(2)得出统一的法向量是:
该统一的法向量与非归一化向量Pu,Pv形成一个局部向量坐标系;
(3)定义参数曲面的第二个基本形式:
κcosφds2=Ldu2+2Mdudv+Ndv2
其中,L=Puu·n,M=Puv·n,N=Pvv·n,;
式中,Puu,Puv,Pvv都是定点P的二阶偏导;
参数曲面的第二个基本形式公式意味着对于给定方向du/dv在u,v平面和给定的角度φ,通过第一和第二基本形式,计算曲面上描绘的曲线的曲率K,也就是切线指向;
(4)引入了两个对称矩阵:
因为Pu和Pv是线性独立的,所以矩阵F1 -1F2也是对称的,因此可以解出F1 -1F2的实特征值和正交特征向量;两个特征值K1和K2就是两个主曲率和两个特征向量t1=(ζ1,η1)T,t2=(ζ2,η2)T定义了两个主方向:
t1=ζ1Pu+η1Pv;t2=ζ2Pu+η2Pv通过以上对法向量n和主方向t1,t2的定义和计算,对任意给定点P的特征描述表示为向量形状,P—(n,t1,t2)在匹配过程中由于一点的主曲率决定了两个主方向,因此通过匹配主曲率的方式来配对点对,这里不考虑脐点(两个主曲率相等的点)提高了匹配的质量和效率。对于待测两点pi和pj的主曲率k1,k2和k1’,k2’的匹配过程中,满足下列判定:
(k1-k1'<r1)&&(k2-k2'<r2)
其中r1,r2分别为满足两个主曲率近似相等的最大差值。
5.根据权利要求1所述的一种基于B样条曲面相似性检测的三维点云识别方法,其特征在于:所述步骤4的具体计算方法为:计算等距就是配对的(pi,pj)计算pi的局部特征向量到pj的局部特征向量的转换;
用Eij表示变换的旋转,tij表示相应的平移变换,计算如下:
tij=pj-Eij*pi
因此存在Tij:pi(ni,t1 i,t2 i)→pj(nj,t1 j,t2 j),等距变换下的空间集合,而Tij就是一次等距变换。
6.根据权利要求1所述的一种基于B样条曲面相似性检测的三维点云识别方法,其特征在于:所述步骤5的具体方法为:
给定存在一个唯一的变换其中g是一个具有非空集合的等距变换,这里是在固定点A下的平移,如果T=g,并且d=0,那么至少有一个固定点;如果将T规定为直接等距;那么把T规定为间接等距。
7.根据权利要求1所述的一种基于B样条曲面相似性检测的三维点云识别方法,其特征在于:所述步骤6的等距比较具体计算方法是:
(1)直接等距的点对等距距离:
dist(T,T')=(1-|cos(D·D')|)+|(θ-θ')/2π|+ω1dist(D,D')+ω2|(||t||-||t'||)|
(2)间接等距的点对等距距离:
8.根据权利要求1所述的一种基于B样条曲面相似性检测的三维点云识别方法,其特征在于:所述步骤7的具体方法为:令dist(T,T’)与两个点对(pi,pj)和(pi’,pj’)对应的同一类别中的两个等距T和T'之间的距离;dist(T,T’)以及相似度度量将根据等距的类别而变化;令S={(pi,pj)}{i,j=1,2,.....Nl}∈Γl,l∈[|1,5|]为第l个等距空间中的Nl个点对的集合,令k为聚类数。
9.根据权利要求8所述的一种基于B样条曲面相似性检测的三维点云识别方法,其特征在于:所述步骤7的具体步骤包括:
(1)输入样本相似矩阵S和聚类数K;
(2)选择高斯相似度函数来构造权重矩阵A∈RNl×k和度矩阵D:
(3)构造标准化的拉普拉斯矩阵解出L矩阵的K个最大特征向量构造矩阵X=[X1,X2,....Xk]∈RNl×k;
(4)通过X矩阵的按行标准化构造成特征Y矩阵,用K-means算法将它们分成K个聚类,当且仅当矩阵Y的第i行被分配给类θ时,将原始点对(pi,pj)分配给类别θ;
(5)确定高斯相似度参数σ
高斯相似度函数中定义的描述相似度公式取决于参数σ,定义了点对(pi,pj) 之间的变换距离的阈值;为了设置它,需要考虑点的均匀分布;定义各向分布的元素被包括在大小为Dmax=maxdist(T,T')的边界框中,通过将这个盒子分成Nl个相同的体积,两个点之间的(均匀的)距离Dun:
Dun=max(pi,pj)≠(pi',pj')d(T,T')/Nl 1/m:
其中,m是等距空间的维度,根据等距性质而变化;如果一个聚类存在,那么就有一些点被小于Dun距离的点分开;类似地,高斯参数σ被用作距离的一部分σ=Dun/2;因此,将m维空间中的点集成在一起,并导出了待定点的阈值封闭;
(6)确定最优聚类数K
通过利用高斯亲和矩阵A并基于Frobenius范数的比例定义一个质量度量;设αk是要识别的簇的数量的界限;对于一个值k'∈[|2,αk|],亲和矩阵是每个群集的索引;然后定义块矩阵:非对角块表示簇之间的亲和度,对角块表示簇内的亲和度;从这个结构中,可以评估所有非对角线区块与Frobenius范数的对角块之间的平均比率,称为rk';由此,在K'属于[|2,αk|]的值中,比率rk'的最小值定义了最优类数K:
k=argminrk'k'∈[|2,αk|]
这个最小值是在集群内的高斯亲和度最高情况下,对应于集群之间的亲和力最低的值。
10.根据权利要求1所述的一种基于B样条曲面相似性检测的三维点云识别方法,其特征在于:所述步骤8的具体步骤包括:
(1)通过谱聚类算法可以得到等距分类中物体之间相似部分点对的数量Ni(i∈[1,5]),其中计算后得到的点云对数Ni(i∈[1,5])对应着被划分为5个等距空间中的对应相似点云,采样到的点云对数为Nu;
(2)对于整体参数的统一采样,定义相似点对数与采样点云对数比值λ:
如果λ≥实验阈值η,说明验证了两个点云是相似的;认为是相同的对象,否则它们不相似,就是不同的对象;
(3)局部分割式的参数采样,对应分割的不完整局部多次进行判定,定义比例阈值ε:
ε=Ni/Ns
其中,Ni是相似的局部点云数量,Ns是完整点云数量;若ε≥η阈值,那么就认为是相同,即完成目标之间的识别。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811001394.7A CN109447100B (zh) | 2018-08-30 | 2018-08-30 | 一种基于b样条曲面相似性检测的三维点云识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811001394.7A CN109447100B (zh) | 2018-08-30 | 2018-08-30 | 一种基于b样条曲面相似性检测的三维点云识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109447100A true CN109447100A (zh) | 2019-03-08 |
CN109447100B CN109447100B (zh) | 2020-12-01 |
Family
ID=65530156
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811001394.7A Expired - Fee Related CN109447100B (zh) | 2018-08-30 | 2018-08-30 | 一种基于b样条曲面相似性检测的三维点云识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109447100B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110061869A (zh) * | 2019-04-09 | 2019-07-26 | 中南民族大学 | 一种基于关键词的网络轨迹分类方法及装置 |
CN110246225A (zh) * | 2019-04-15 | 2019-09-17 | 上海交通大学医学院附属第九人民医院 | 舌癌术后缺损舌体重建方法、系统、可读存储介质及设备 |
CN110427956A (zh) * | 2019-04-18 | 2019-11-08 | 中国林业科学研究院资源信息研究所 | 一种基于谱聚类算法的LiDAR点云数据单木提取方法 |
CN110728327A (zh) * | 2019-10-18 | 2020-01-24 | 中国科学技术大学 | 一种具有可解释性的直推学习方法及系统 |
CN111369609A (zh) * | 2020-03-04 | 2020-07-03 | 山东交通学院 | 一种基于点云曲面特征约束的建筑物局部变形分析方法 |
CN111402393A (zh) * | 2019-12-06 | 2020-07-10 | 温州大学 | 参数曲面模拟点云生成方法 |
CN112488127A (zh) * | 2020-12-09 | 2021-03-12 | 南京邮电大学 | 基于层次b样条的多尺度点云特征点检测与匹配 |
CN113076870A (zh) * | 2021-03-31 | 2021-07-06 | 广东电网有限责任公司 | 一种基于激光点云的杆塔识别方法及系统 |
CN117828379A (zh) * | 2024-03-05 | 2024-04-05 | 山东省地质测绘院 | 基于多源数据融合的地下资源探测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101216956A (zh) * | 2007-12-29 | 2008-07-09 | 浙江工业大学 | 基于nurbs的心脏三维表示方法 |
CN103810271A (zh) * | 2014-01-29 | 2014-05-21 | 辽宁师范大学 | 基于路径跟随的三维点云物体形状特征匹配方法 |
US20150339541A1 (en) * | 2014-05-22 | 2015-11-26 | Nokia Technologies Oy | Point cloud matching method |
CN106932271A (zh) * | 2017-03-10 | 2017-07-07 | 厦门大学 | 一种基于逆向工程的球压试验压痕尺寸测量方法 |
US20170374342A1 (en) * | 2016-06-24 | 2017-12-28 | Isee, Inc. | Laser-enhanced visual simultaneous localization and mapping (slam) for mobile devices |
-
2018
- 2018-08-30 CN CN201811001394.7A patent/CN109447100B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101216956A (zh) * | 2007-12-29 | 2008-07-09 | 浙江工业大学 | 基于nurbs的心脏三维表示方法 |
CN103810271A (zh) * | 2014-01-29 | 2014-05-21 | 辽宁师范大学 | 基于路径跟随的三维点云物体形状特征匹配方法 |
US20150339541A1 (en) * | 2014-05-22 | 2015-11-26 | Nokia Technologies Oy | Point cloud matching method |
US20170374342A1 (en) * | 2016-06-24 | 2017-12-28 | Isee, Inc. | Laser-enhanced visual simultaneous localization and mapping (slam) for mobile devices |
CN106932271A (zh) * | 2017-03-10 | 2017-07-07 | 厦门大学 | 一种基于逆向工程的球压试验压痕尺寸测量方法 |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110061869B (zh) * | 2019-04-09 | 2022-04-15 | 中南民族大学 | 一种基于关键词的网络轨迹分类方法及装置 |
CN110061869A (zh) * | 2019-04-09 | 2019-07-26 | 中南民族大学 | 一种基于关键词的网络轨迹分类方法及装置 |
CN110246225A (zh) * | 2019-04-15 | 2019-09-17 | 上海交通大学医学院附属第九人民医院 | 舌癌术后缺损舌体重建方法、系统、可读存储介质及设备 |
CN110246225B (zh) * | 2019-04-15 | 2023-04-21 | 上海交通大学医学院附属第九人民医院 | 舌癌术后缺损舌体重建方法、系统、可读存储介质及设备 |
CN110427956A (zh) * | 2019-04-18 | 2019-11-08 | 中国林业科学研究院资源信息研究所 | 一种基于谱聚类算法的LiDAR点云数据单木提取方法 |
CN110728327A (zh) * | 2019-10-18 | 2020-01-24 | 中国科学技术大学 | 一种具有可解释性的直推学习方法及系统 |
CN110728327B (zh) * | 2019-10-18 | 2021-11-23 | 中国科学技术大学 | 一种具有可解释性的直推学习方法及系统 |
CN111402393A (zh) * | 2019-12-06 | 2020-07-10 | 温州大学 | 参数曲面模拟点云生成方法 |
CN111369609A (zh) * | 2020-03-04 | 2020-07-03 | 山东交通学院 | 一种基于点云曲面特征约束的建筑物局部变形分析方法 |
CN111369609B (zh) * | 2020-03-04 | 2023-06-30 | 山东交通学院 | 一种基于点云曲面特征约束的建筑物局部变形分析方法 |
CN112488127A (zh) * | 2020-12-09 | 2021-03-12 | 南京邮电大学 | 基于层次b样条的多尺度点云特征点检测与匹配 |
CN112488127B (zh) * | 2020-12-09 | 2022-08-09 | 南京邮电大学 | 基于层次b样条的多尺度点云特征点检测与匹配 |
CN113076870A (zh) * | 2021-03-31 | 2021-07-06 | 广东电网有限责任公司 | 一种基于激光点云的杆塔识别方法及系统 |
CN117828379A (zh) * | 2024-03-05 | 2024-04-05 | 山东省地质测绘院 | 基于多源数据融合的地下资源探测方法 |
CN117828379B (zh) * | 2024-03-05 | 2024-05-07 | 山东省地质测绘院 | 基于多源数据融合的地下资源探测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109447100B (zh) | 2020-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109447100A (zh) | 一种基于b样条曲面相似性检测的三维点云识别方法 | |
Guo et al. | Efficient center voting for object detection and 6D pose estimation in 3D point cloud | |
Vesanto | Data exploration process based on the self-organizing map | |
CN109359684A (zh) | 基于弱监督定位和子类别相似性度量的细粒度车型识别方法 | |
Shamir | Segmentation and shape extraction of 3D boundary meshes. | |
Demisse et al. | Deformation based curved shape representation | |
CN106228539B (zh) | 一种三维点云中多种几何基元自动识别方法 | |
Zhang et al. | KDD: A kernel density based descriptor for 3D point clouds | |
CN110222642A (zh) | 一种基于全局图聚类的平面建筑构件点云轮廓提取方法 | |
Liu et al. | A novel rock-mass point cloud registration method based on feature line extraction and feature point matching | |
CN109993213A (zh) | 一种用于服装部件图的自动识别方法 | |
CN105930859B (zh) | 基于线性流形聚类的雷达信号分选方法 | |
CN106682575A (zh) | Elm算法的人眼点云特征定位 | |
Chen et al. | Mesh convolution: a novel feature extraction method for 3d nonrigid object classification | |
CN114648445B (zh) | 一种基于特征点提取及精配准优化的多视角高分辨率点云拼接方法 | |
Akgül et al. | Density-based 3D shape descriptors | |
Tu et al. | Feature extraction using multidimensional spectral regression whitening for hyperspectral image classification | |
US20210150078A1 (en) | Reconstructing an object | |
Xin et al. | Accurate and complete line segment extraction for large-scale point clouds | |
CN106951873B (zh) | 一种遥感图像目标识别方法 | |
Liu et al. | Neighbor reweighted local centroid for geometric feature identification | |
CN103559488B (zh) | 基于定性空间关系的卡通人物的脸部特征提取方法 | |
Amsaleg et al. | Intrinsic dimensionality estimation within tight localities: A theoretical and experimental analysis | |
Xiao et al. | Hierarchical aggregation for efficient shape extraction | |
Ju et al. | A robust parity test for extracting parallel vectors in 3D |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201201 Termination date: 20210830 |
|
CF01 | Termination of patent right due to non-payment of annual fee |