CN109215129B - 一种基于三维点云的局部特征描述方法 - Google Patents
一种基于三维点云的局部特征描述方法 Download PDFInfo
- Publication number
- CN109215129B CN109215129B CN201710540090.7A CN201710540090A CN109215129B CN 109215129 B CN109215129 B CN 109215129B CN 201710540090 A CN201710540090 A CN 201710540090A CN 109215129 B CN109215129 B CN 109215129B
- Authority
- CN
- China
- Prior art keywords
- point
- dimensional
- feature
- points
- coordinate system
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/30—Polynomial surface description
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及一种基于三维点云的局部特征描述方法,本发明包括以下步骤:从场景点云中提取若干特征点,以每个特征点为中心,根据其球邻域内的点建立三维局部坐标系;将特征点球邻域内的点变换到相应的局部坐标系下,并沿径向对球邻域进行空间区域划分。针对每个被划分的空间区域,计算位于其中的每个点与坐标系x轴和z轴的夹角余弦值α和β,并将它们分别映射到两个独立的一维直方图。将所有区域的一维直方图串联起来,并除以特征点球邻域内点的总数进行归一化进而得到最终的三维局部特征描述。本发明提出的三维局部特征描述符具有区分性好、鲁棒性强和计算效率高的特点,提高了场景点云之间相似部位的正确匹配率。
Description
技术领域
本发明属于计算机视觉领域,具体的说涉及一种基于三维点云的局部特征描述方法。
背景技术
在两个不同的点云场景中寻找相似点是很多应用的基础,例如三维场景重建、三维目标识别以及三维目标检索等。寻找相似点的过程称为特征匹配。而对三维场景进行特征描述又是特征匹配的前提。由于传感器自身精度的限制,加上采集视角的不同,获取的场景点云数据可能会出现分辨率不均匀,夹带噪声,某些区域甚至出现孔洞以及感兴趣目标被遮挡的情形,这些因素都会加大特征描述的难度。如何设计一种区分性好、鲁棒性强的三维局部特征正是本发明要解决的问题。
自旋图特征和Signature of Histograms of OrienTations(SHOT)特征是目前研究领域效果较好的局部特征描述符。自旋图特征首先以特征点的法向量为主轴建立局部圆柱坐标系,然后将特征点球邻域内每个点在该坐标系下的二元表示映射到一个二维直方图形成对特征点局部的描述。该方法利用了直方图的统计特性,具有较好的鲁棒性。但是该特征描述没有考虑邻域点的三维信息,故而区分性不强。SHOT特征首先根据特征点的邻域信息建立一个三维局部坐标系,然后以该坐标系为参考轴,将特征点的球邻域沿径向、经线方向以及纬线方向共划分为32个子区域。对于每个子区域,然后计算其中每个点的法向量与特征点法向量的夹角,并将其映射到一个长度为11的一维直方图中,最终得到一个352维的局部特征。该方法通过建立局部坐标系以及区域划分,使得邻域点的三维空间信息得到了很好的保留。但是需要估计每个点的法向量,而法向量的估计容易受到噪声以及分辨率不均匀的影响,所以最终得到的特征描述符也会对这些因素敏感。
发明内容
针对现有技术中存在的不足,本发明直接对特征点球邻域点的三维位置信息进行编码,避免使用点的法向量信息。提出了一种基于三维点云的局部特征描述方法。
本发明为实现上述目的所采用的技术方案是:一种基于三维点云的局部特征描述方法,包括以下步骤:
从场景点云中提取若干特征点,以每个特征点为中心,根据其球邻域内的点建立三维局部坐标系;
将特征点球邻域内的点变换到该特征点的局部坐标系下,并沿径向对球邻域进行空间区域划分;
针对每个被划分的空间区域,计算位于其中的每个点与坐标系x轴和z轴的夹角余弦值α和β,并将它们分别映射到两个独立的一维直方图;
将所有区域的一维直方图串联起来,并除以该特征点球邻域内点的总数进行归一化进而得到最终的三维局部特征描述。
所述根据其球邻域内的点建立三维局部坐标系包括以下步骤:
给定一个特征点p及其半径为R的球邻域点集合S(p)={q1,q2,...,qn},满足||qi-p||2≤R,其中n代表球邻域内点的数量;
首先,以特征点p为局部坐标系原点,按如下公式计算加权协方差矩阵C:
其中,p=(x,y,z)T,qi=(xi,yi,zi)T,di=||qi-p||2,(x,y,z)T为特征点p的三维坐标,(xi,yi,zi)T表示特征点p的第i个邻域点的三维坐标;
然后对矩阵C进行特征值分解:
CV=EV, (10)
其中,E=diag{λ1,λ2,λ3},λk,k=1,2,3是矩阵C的3个特征值,矩阵V={v1,v2,v3},vk,k=1,2,3分别是每个特征值对应的特征向量;
按照特征值的绝对值大小进行降序排列,记最小特征值对应的特征向量为vz;则局部坐标系的z轴由如下公式确定:
记最大特征值对应的特征向量为vx;局部坐标系的x轴由下面的公式确定:
最后,局部坐标系的y轴为z×x。
将特征点球邻域内的点变换到该特征点的局部坐标系下包括以下步骤:
q′i=[x,y,z](qi-p),i=1,2,...,n (13)
以特征点为坐标原点的邻域点集为Q(p)={q′1,q′2,...,q′n};然后沿径向将球邻域均匀划分为Nr份。
所述计算位于其中的每个点与坐标系x轴和z轴的夹角余弦值α和β,并将它们分别映射到两个独立的一维直方图包括以下步骤:
某被划分的空间区域中的每一个点q′i,它与该特征点的局部坐标系原点连线形成的方向向量与x轴、z轴的夹角余弦值α和β的计算公式如下:
α=q′i·x/||q′i||2=q′ix/||q′i||2 (14)
β=q′i·z/||q′i||2=q′iz/||q′i||2 (15)
其中,q′ix与q′iz分别是点q′i的x坐标和z坐标;
夹角余弦值的范围为[-1,1],将该空间区域均匀划分为Nh份,然后将α和β分别映射到两个独立的直方图中;从而每个被划分的空间区域得到两个一维统计直方图hα和hβ。
所述最终的局部特征描述符为:
一种基于三维点云的局部特征描述方法,用于实现两个点云场景间相似部位的匹配。
本发明具有以下优点及有益效果:
1.将特征点的球邻域点变换到以特征点为原点建立的三维局部坐标系下使得最终的特征描述符具有旋转和平移不变性;
2.结合空间区域划分以及直方图统计的优点,最终建立的局部特征描述符具有较强的区分性和鲁棒性。
3.最终的特征描述符维度为2NrNh,研究发现取Nr=6,Nh=11时,有利于节省存储空间以及加速后续的特征匹配效率。
附图说明
图1为基于点云的三维局部特征生成示意图;
具体实施方式
下面结合附图及实施例对本发明做进一步的详细说明。
本发明面向实际工程应用,提出了一种基于三维点云的局部特征描述方法,可以实现两个点云场景间相似部位的匹配,能够适用于复杂多变的背景条件。本发明包括以下四个步骤:1、从场景点云中提取若干特征点,以每个特征点为中心,根据其球邻域内的点建立三维局部坐标系;2、将特征点球邻域内的点变换到相应的局部坐标系下,并沿径向对球邻域进行空间区域划分;3、针对每个被划分的空间区域,计算位于其中的每个点与坐标系x轴和z轴的夹角余弦值α和β,并将它们分别映射到两个独立的一维直方图。4、将所有区域的一维直方图串联起来,并除以特征点球邻域内点的总数进行归一化进而得到最终的三维局部特征描述。本发明提出的三维局部特征描述符具有区分性好、鲁棒性强和计算效率高的特点,提高了场景点云之间相似部位的正确匹配率。
一种基于三维点云的局部特征描述方法,包括以下步骤:
从场景点云中提取若干特征点,以每个特征点为中心,根据其球邻域内的点建立三维局部坐标系;
将特征点球邻域内的点变换到相应的局部坐标系下,并沿径向对球邻域进行空间区域划分;
针对每个被划分的空间区域,计算位于其中的每个点与坐标系x轴和z轴的夹角余弦值α和β,并将它们分别映射到两个独立的一维直方图;
将所有区域的一维直方图串联起来,并除以特征点球邻域内点的总数进行归一化进而得到最终的三维局部特征描述。
给定一个特征点p及其半径为R的球邻域点集合S(p)={q1,q2,...,qn}满足||qi-p||2≤R,其中n代表球邻域内点的数量。所述三维局部坐标系的建立步骤如下:
首先,以特征点p为局部坐标系原点,按如下公式计算加权协方差矩阵C:
其中,p=(x,y,z)T,qi=(xi,yi,zi)T,di=||qi-p||2。然后对矩阵C进行特征值分解:
CV=EV, (18)
其中,E=diag{λ1,λ2,λ3},λk,k=1,2,3是矩阵C的3个特征值,矩阵V={v1,v2,v3},vk,k=1,2,3分别是每个特征值对应的特征向量。按照特征值的绝对值大小进行降序排列,记最小特征值对应的特征向量为vz。则局部坐标系的z轴由如下公式确定:
最后,局部坐标系的y轴为z×x。
所述球邻域内的点变换到局部坐标系的公式如下:
q′i=[x,y,z](qi-p),i=1,2,...,n (21)
则以特征点为坐标原点的邻域点集为Q(p)={q′1,q′2,...,q′n}。然后沿径向将球邻域进均匀划分为6份。
所述每个子区域中的每一个点q′i,它与原点连线形成的方向向量与x轴、z轴的夹角余弦值α和β的计算公式如下:
α=q′i·x/||q′i||2=q′ix/||q′i||2 (22)
β=q′i·z/||q′i||2=q′iz/||q′i||2 (23)
其中,q′ix与q′iz分别是点q′i的x坐标和z坐标。由于夹角余弦值的范围为[-1,1],将该区间均匀划分为11个箱子,将α和β分别映射到两个独立的直方图中。从而每个子区域可以得到两个一维统计直方图hα和hβ。
所述最终的局部特征描述符的生成方式如下式所示:
最后将每个特征点的三维局部特征描述存入特征库,构成对整个点云场景的特征描述,并用于后续的特征匹配。
参阅图1,一种基于点云的局部特征描述方法,包含以下步骤:
1.对于某个特征点p,根据K近邻算法获取与之距离不超过R的邻域点集合S(p)={q1,q2,...,qn},||qi-p||2≤R。
2.以特征点p为中心,根据下式计算加权协方差矩阵C:
对矩阵C进行特征值分解,得到三个特征值λ1,λ2,λ3以及对应的三个特征向量v1,v2,v3。对特征值按照其绝对值的大小进行降序排列,记最小的那个特征值对应的特征向量为vz。然后由下式确定局部坐标系的z轴:
同理,记最大的那个特征值对应的特征向量为vx。接着由下式确定局部坐标系的x轴:
最后y轴由z轴叉乘x轴得到,即z×x。
3.根据第二步中建立好的局部坐标系,将特征点p的所有球邻域点变换到该局部坐标系下,得到以点p为坐标原点的邻域点集合Q(p)={q′1,q′2,...,q′n}。其中,q′i=[x,y,z](qi-p),i=1,2,...,n。
4.将点p的球邻域按照径向等间隔划分为6个球壳空间(图1为了显示清晰只划分为了3个)。针对每个球壳区间,找出位于其中的点,并计算点与原点连线所形成的向量pq′分别于x轴、z轴的夹角和θ。为了提高计算效率,直接计算与它们等效的夹角余弦值α和β。然后将α和β分别映射到两个独立的长度为11的直方图中。该直方图的箱子尺寸为2/11,是将余弦值域区间[-1,1]等间隔划分为11份确定的。
5.根据第4步的计算,每个球壳区域都可以生成两个长度为11的直方图。则6个球壳区域共产生了12个一维直方图。将这12个一维直方图依次串联形成一个长度为132的直方图,然后对直方图的每个元素都除以特征点p的邻域点总数n进行归一化,从而得到最终的局部特征描述符。
Claims (6)
1.一种基于三维点云的局部特征描述方法,其特征在于,包括以下步骤:
从场景点云中提取若干特征点,以每个特征点为中心,根据其球邻域内的点建立三维局部坐标系;
将特征点球邻域内的点变换到该特征点的局部坐标系下,并沿径向对球邻域进行空间区域划分;
针对每个被划分的空间区域,计算位于其中的每个点与坐标系x轴和z轴的夹角余弦值α和β,并将它们分别映射到两个独立的一维直方图;
将所有区域的一维直方图串联起来,并除以该特征点球邻域内点的总数进行归一化进而得到最终的三维局部特征描述。
2.根据权利要求1所述的一种基于三维点云的局部特征描述方法,其特征在于所述根据其球邻域内的点建立三维局部坐标系包括以下步骤:
给定一个特征点p及其半径为R的球邻域点集合S(p)={q1,q2,...,qn},满足||qi-p||2≤R,其中n代表球邻域内点的数量;
首先,以特征点p为局部坐标系原点,按如下公式计算加权协方差矩阵C:
其中,p=(x,y,z)T,qi=(xi,yi,zi)T,di=||qi-p||2,(x,y,z)T为特征点p的三维坐标,(xi,yi,zi)T表示特征点p的第i个邻域点的三维坐标;
然后对矩阵C进行特征值分解:
CV=EV, (2)
其中,E=diag{λ1,λ2,λ3},λk,k=1,2,3是矩阵C的3个特征值,矩阵V={v1,v2,v3},vk,k=1,2,3分别是每个特征值对应的特征向量;
按照特征值的绝对值大小进行降序排列,记最小特征值对应的特征向量为vz;则局部坐标系的z轴由如下公式确定:
记最大特征值对应的特征向量为vx;局部坐标系的x轴由下面的公式确定:
最后,局部坐标系的y轴为z×x。
3.根据权利要求1所述的一种基于三维点云的局部特征描述方法,其特征在于,所述将特征点球邻域内的点变换到该特征点的局部坐标系下包括以下步骤:
q′i=[x,y,z](qi-p),i=1,2,...,n (5)
以特征点为坐标原点的邻域点集为Q(p)={q′1,q′2,...,q′n};然后沿径向将球邻域均匀划分为Nr份。
4.根据权利要求1所述的一种基于三维点云的局部特征描述方法,其特征在于,所述计算位于其中的每个点与坐标系x轴和z轴的夹角余弦值α和β,并将它们分别映射到两个独立的一维直方图包括以下步骤:
某被划分的空间区域中的每一个点q′i,它与该特征点的局部坐标系原点连线形成的方向向量与x轴、z轴的夹角余弦值α和β的计算公式如下:
α=q′i·x/||q′i||2=q′ix/||q′i||2 (6)
β=q′i·z/||q′i||2=q′iz/||q′i||2 (7)
其中,q′ix与q′iz分别是点q′i的x坐标和z坐标;
夹角余弦值的范围为[-1,1],将该空间区域均匀划分为Nh份,然后将α和β分别映射到两个独立的直方图中;从而每个被划分的空间区域得到两个一维统计直方图hα和hβ。
6.根据权利要求1所述的一种基于三维点云的局部特征描述方法,其特征在于用于实现两个点云场景间相似部位的匹配。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710540090.7A CN109215129B (zh) | 2017-07-05 | 2017-07-05 | 一种基于三维点云的局部特征描述方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710540090.7A CN109215129B (zh) | 2017-07-05 | 2017-07-05 | 一种基于三维点云的局部特征描述方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109215129A CN109215129A (zh) | 2019-01-15 |
CN109215129B true CN109215129B (zh) | 2022-10-04 |
Family
ID=64992581
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710540090.7A Active CN109215129B (zh) | 2017-07-05 | 2017-07-05 | 一种基于三维点云的局部特征描述方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109215129B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111062266B (zh) * | 2019-11-28 | 2022-07-15 | 东华理工大学 | 基于圆柱坐标的人脸点云关键点定位方法 |
US11625454B2 (en) | 2019-12-09 | 2023-04-11 | Shenzhen University | Method and device for 3D shape matching based on local reference frame |
US20230015645A1 (en) * | 2019-12-09 | 2023-01-19 | Shenzhen University | 3d shape matching method and device based on 3d local feature description using sghs |
CN113033270B (zh) * | 2019-12-27 | 2023-03-17 | 深圳大学 | 采用辅助轴的3d物体局部表面描述方法、装置及存储介质 |
CN111553343B (zh) * | 2020-04-01 | 2023-04-25 | 青岛联合创智科技有限公司 | 一种激光点云特征的提取方法 |
CN111931604A (zh) * | 2020-07-22 | 2020-11-13 | 哈尔滨工业大学(深圳) | 基于二值描述符的3d特征提取方法、系统及存储介质 |
CN112001261A (zh) * | 2020-07-28 | 2020-11-27 | 哈尔滨工业大学(深圳) | 基于toldi描述符的3d特征提取方法、系统及介质 |
CN112330804B (zh) * | 2020-11-12 | 2021-10-19 | 中国人民解放军火箭军工程大学 | 一种局部可变形三维模型接触检测方法 |
CN113111741B (zh) * | 2021-03-27 | 2024-05-07 | 西北工业大学 | 一种基于三维特征点的装配状态识别方法 |
CN114972459B (zh) * | 2022-05-31 | 2023-09-08 | 哈尔滨理工大学 | 一种基于低维点云局部特征描述符的点云配准方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104134216A (zh) * | 2014-07-29 | 2014-11-05 | 武汉大学 | 基于16维特征描述的激光点云自动配准方法及系统 |
CN105160344A (zh) * | 2015-06-18 | 2015-12-16 | 北京大学深圳研究生院 | 一种三维点云的局部特征提取方法及装置 |
WO2016201671A1 (zh) * | 2015-06-18 | 2016-12-22 | 北京大学深圳研究生院 | 一种三维点云的局部特征提取方法及装置 |
-
2017
- 2017-07-05 CN CN201710540090.7A patent/CN109215129B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104134216A (zh) * | 2014-07-29 | 2014-11-05 | 武汉大学 | 基于16维特征描述的激光点云自动配准方法及系统 |
CN105160344A (zh) * | 2015-06-18 | 2015-12-16 | 北京大学深圳研究生院 | 一种三维点云的局部特征提取方法及装置 |
WO2016201671A1 (zh) * | 2015-06-18 | 2016-12-22 | 北京大学深圳研究生院 | 一种三维点云的局部特征提取方法及装置 |
Non-Patent Citations (2)
Title |
---|
复杂场景下基于C-SHOT特征的3D物体识别与位姿估计;张凯霖等;《计算机辅助设计与图形学学报》;20170515(第05期);全文 * |
用于三维点云表示的扩展点特征直方图算法;庄祉昀等;《国防科技大学学报》;20161228(第06期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109215129A (zh) | 2019-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109215129B (zh) | 一种基于三维点云的局部特征描述方法 | |
Guo et al. | A novel local surface feature for 3D object recognition under clutter and occlusion | |
CN109559340B (zh) | 一种并行的三维点云数据自动化配准方法 | |
JP2014081347A (ja) | 3dシーンにおける3d物体の認識および姿勢決定のための方法 | |
CN111340862B (zh) | 一种基于多特征融合的点云配准方法、装置及存储介质 | |
KR20120094353A (ko) | 시각적 형상 디스크립터 생성 방법 및 그 기록매체 | |
CN113592030B (zh) | 基于复值奇异谱分析的图像检索方法和系统 | |
CN112651408B (zh) | 一种基于点对变换特征的三维局部表面描述方法和系统 | |
CN110633749A (zh) | 基于改进视点特征直方图的三维点云识别方法 | |
CN114332172A (zh) | 一种基于协方差矩阵改进的激光点云配准方法 | |
CN114358166A (zh) | 一种基于自适应k均值聚类的多目标定位方法 | |
CN111401385B (zh) | 一种图像局部拓扑结构特征描述符的相似度计算方法 | |
Bhagyalakshmi et al. | A survey on content based image retrieval using various operators | |
CN111009004B (zh) | 一种加速图像匹配的硬件优化方法 | |
CN115830340B (zh) | 一种点云目标识别方法和装置、存储介质 | |
CN108932726A (zh) | 一种目标跟踪方法及装置 | |
Moon et al. | Fast image-matching technique robust to rotation in spherical images | |
CN105913068A (zh) | 一种用于描述图像特征的多维方向梯度表示方法 | |
Yang et al. | Feature points extraction based on improved ORB-SLAM | |
Wang et al. | Feature extraction algorithm based on improved ORB with adaptive threshold | |
CN110134693A (zh) | 基于哈希和pca的时空索引建立方法 | |
Ni et al. | Research on image segmentation algorithm based on fuzzy clustering and spatial pyramid | |
Ma et al. | Improving KAZE feature matching algorithm with alternative image gray method | |
CN112750144B (zh) | 点云匹配中一种基于特征直方图的点云特征提取方法 | |
Fan et al. | A Semantic-Based Loop Closure Detection of 3D Point Cloud |
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 |