一种基于三维点云的局部特征描述方法与流程

本发明属于计算机视觉领域，具体的说涉及一种基于三维点云的局部特征描述方法。

背景技术：

在两个不同的点云场景中寻找相似点是很多应用的基础，例如三维场景重建、三维目标识别以及三维目标检索等。寻找相似点的过程称为特征匹配。而对三维场景进行特征描述又是特征匹配的前提。由于传感器自身精度的限制，加上采集视角的不同，获取的场景点云数据可能会出现分辨率不均匀，夹带噪声，某些区域甚至出现孔洞以及感兴趣目标被遮挡的情形，这些因素都会加大特征描述的难度。如何设计一种区分性好、鲁棒性强的三维局部特征正是本发明要解决的问题。

自旋图特征和signatureofhistogramsoforientations(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轴由如下公式确定：

其中，当且仅当否则为0；

记最大特征值对应的特征向量为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β。

所述最终的局部特征描述符为：

其中分别是第j个区域的两个直方图。

一种基于三维点云的局部特征描述方法，用于实现两个点云场景间相似部位的匹配。

本发明具有以下优点及有益效果：

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轴由如下公式确定：

其中，当且仅当否则为0。记最大特征值对应的特征向量为vx。类似地，局部坐标系的x轴由下面的公式确定：

最后，局部坐标系的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β。

所述最终的局部特征描述符的生成方式如下式所示：

其中分别是第j个区域的两个直方图。

最后将每个特征点的三维局部特征描述存入特征库，构成对整个点云场景的特征描述，并用于后续的特征匹配。

参阅图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进行归一化，从而得到最终的局部特征描述符。