CN106250898B - 一种基于尺度预测的图像局部区域特征提取方法 - Google Patents

Abstract

一种基于尺度预测的图像局部区域特征提取方法，涉及行星着陆中图像处理技术领域，具体涉及一种基于尺度预测的特征提取方法。为了解决SURF特征提取方法不具有仿射不变性的问题。本发明首先对当前所拍摄图像在目标天体表面的位置进行初步估计，选定全局特征库中的搜索范围；然后根据探测器拍照时的位姿信息以及全局特征库中相应特征点的特征尺度和方向信息，预测该特征在拍摄图像中的特征尺度和旋转角度；并将拍摄的图像按照所预测的旋转角度旋转获得旋转图像，根据二次型求极大值方法以及预测的特征尺度生成特征检测模板，将特征检测模板与旋转图像卷积并在位置空间进行非极大值抑制从而实现特征提取。本发明适用于行星着陆自主导航的图像处理。

Description

一种基于尺度预测的图像局部区域特征提取方法

技术领域

本发明涉及行星着陆自主导航中图像处理技术领域，具体涉及一种基于尺度预测的特征提取方法。

背景技术

在行星着陆自主导航任务设计中，需要通过前期在轨观测器观测图像及行星车收集的高清图像中建立一个标定有行星固联坐标系的二维图像数据库，在着陆过程中，探测器通过从星载相机对地拍摄的图像中提取视觉信息并与数据库进行匹配来获得自身的绝对位姿信息，因此，高效、准确的识别图像中的视觉信息是行星自主导航任务的重要前提。由于组建的数据库与行星着陆所拍摄的图像存在光照、拍摄角度、尺度等差异，从图像中提炼的视觉信息也应具有相应的不变性以保证有效的特征匹配；与此同时，在星载计算机处理能力有限情况下，视觉特征的提取与匹配还应具备实时性。

图像局部特征的识别隶属计算机视觉领域，在近年来受到了学者们的广泛关注。其中最有代表性的算法是D.Lowe提出的尺度不变特征描述子，又称SIFT。秉承SIFT的构思，亦有学者提出了更高效的算法，如SURF，PCA-SIFT以及具有仿射不变性的Affine-SIFT算法等。概括而言，图像局部特亦可以代表为图像中局部纹理变化较为明显的区域。虽然SIFT描述子在大部分情况下能够正确匹配，但仍有数处误匹配存在。

SURF算法与SIFT算法均可以归类为尺度不变特征家族当中。具体而言，该类算法都是将图像放在尺度空间中，并在多尺度描述中寻找一个类似于“高斯斑点”结构的最佳相应以及对应的特征尺度。为了找到上述结构，SURF与SIFT都使用了“多尺度特征检测子”来分析图像的尺度空间表征。秉承SIFT的设计思想，SURF算法做出了更简洁的近似，进一步降低了计算量。不同于DoG算子，SURF算法使用了Hessian矩阵与积分图像来近似二阶高斯微分，并用盒子滤波器对积分图像进行操作，从而减少了计算量，增加了计算速度。但是由于拍摄条件的不同，在轨拍摄和下降过程中的在线拍摄对同一特征目标的成像结果是不同的，SURF特征提取方法不具有仿射不变性，而且计算量仍然较大。

发明内容

本发明为了解决现有的SURF特征提取方法不具有仿射不变性的问题。

一种基于尺度预测的图像局部区域特征提取方法，包括以下步骤：

步骤一：根据行星着陆过程中探测器相对于预定着陆点的位置、探测器体坐标系相对于在轨拍摄行星表面图像是的姿态及相机焦距、视场角等信息，对当前所拍摄图像在目标天体表面的位置进行初步估计，选定全局特征库中的搜索范围；

步骤二：根据探测器拍照时的位姿信息以及全局特征库中相应特征点的特征尺度，预测该特征在拍摄图像中的特征尺度；

所述步骤二的具体实现过程如下：

建立全局特征库拍摄图像时探测器与该特征的距离为d₁，建立全局特征库拍摄图像时所用相机的焦距为f₁，建立全局特征库拍摄图像时相机光心与特征点之间连线和光轴的夹角为α₁，设全局特征库中的某特征的特征尺度为σ₁，则该特征在下降段图像中会产生放缩变换，其尺度变换为

其中，下标2对应的参数表示相应物理量在下降段拍摄图像时的值，即：下降段拍摄图像时探测器与该特征的距离为d₂，下降段拍摄图像时所用相机的焦距为f₂，下降段拍摄图像时相机光心与特征点之间连线和光轴的夹角为α₂；

设探测器下降时的俯仰角为θ，偏航角为ψ，滚转角为则所述全局特征库中的某特征在下降段图像中会产生仿射变换，相应的椭圆长、短半轴为

步骤三：根据探测器拍照时的姿态信息以及全局特征库中相应特征点的方向信息，预测该特征在拍摄图像中的旋转角度；

步骤三的具体实现过程如下：

所述全局特征库中的某特征在下降段图像中会产生旋转变换，旋转角度为

步骤四：将拍摄的图像按照所预测的旋转角度旋转获得旋转图像；

步骤五：根据二次型求极大值方法仿照SURF检测原理以及预测的特征尺度生成特征检测模板；

步骤五的具体实现过程如下：

首先根据所预测的特征尺度建立二阶高斯微分盒子滤波器D_xx、D_xy和D_yy，然后根据SURF原理建立求Hessian矩阵行列式最大值问题：

从而将上式转化为一个二次型优化问题：

其中，μ、ν分别为像元像线坐标，D_SURF(μ,ν)为SURF特征检测子在(μ,ν)处的Hessian响应值，I为将图像中与盒子滤波器卷积的对应部分重新排列成的列向量，G为相应的盒子滤波器重新排列成的行向量；

对上式中的二次型，通过进行二次型矩阵的特征值分解求取在仿射变换条件下的特征检测模板。

步骤六：将特征检测模板与旋转图像卷积并在位置空间进行非极大值抑制从而实现特征提取。

本发明对当前所拍摄图像在目标天体表面的位置进行初步估计，选定全局特征数据库中的搜索范围，根据探测器拍照时的位姿信息计算数据库中相应特征在拍摄图像中特征尺度和旋转角度，根据旋转角度旋转图像，根据特征尺度生成特征检测模板并对图像进行检测。本发明与现有技术相比的优点在于：能够实现仿射不变性特征检测，同时计算量小，计算效率高，是着陆过程中进行特征检测的理想方法。

附图说明

图1为本发明的流程图；

图2为三维模拟地景正射影像；

图3-1至图3-3为不同高度、拍摄角度下的下降图像影像示意图；

图4为仿真测试方法示意图；

图5为两种SURF特征检测子仿真比较示例图；其中，图5(a)和图5(c)使用了本发明所提出的仿射不变SURF模板进行SURF特征检测示例，图5(b)和图5(d)为标准的SURF检测算子示例；

图6为两种SURF检测子仿真比较结果；其中，图6(a)为两种算法平均定位误差随姿态角改变量变化的曲线，图6(b)为误匹配率随姿态角该变量变化的曲线。

具体实施方式

具体实施方式一：结合图1说明本实施方式，

一种基于尺度预测的图像局部区域特征提取方法，包括以下步骤：

步骤一：根据行星着陆过程中探测器相对于预定着陆点的位置、探测器体坐标系相对于在轨拍摄行星表面图像是的姿态及相机焦距、视场角等信息，对当前所拍摄图像在目标天体表面的位置进行初步估计，选定全局特征库中的搜索范围；

步骤二：根据探测器拍照时的位姿信息以及全局特征库中相应特征点的特征尺度，预测该特征在拍摄图像中的特征尺度；

步骤三：根据探测器拍照时的姿态信息以及全局特征库中相应特征点的方向信息，预测该特征在拍摄图像中的旋转角度；

步骤四：将拍摄的图像按照所预测的旋转角度旋转获得旋转图像；

步骤五：根据二次型求极大值方法仿照SURF检测原理以及预测的特征尺度生成特征检测模板；

步骤六：将特征检测模板与旋转图像卷积并在位置空间进行非极大值抑制从而实现特征提取。

具体实施方式二：

本实施方式所述步骤二的具体实现过程如下：

建立全局特征库拍摄图像时探测器与该特征的距离为d₁，建立全局特征库拍摄图像时所用相机的焦距为f₁，建立全局特征库拍摄图像时相机光心与特征点之间连线和光轴的夹角为α₁，设全局特征库中的某特征的特征尺度为σ₁，则该特征在下降段图像中会产生放缩变换，其尺度变换为

其中，下标2对应的参数表示相应物理量在下降段拍摄图像时的值，即：下降段拍摄图像时探测器与该特征的距离为d₂，下降段拍摄图像时所用相机的焦距为f₂，下降段拍摄图像时相机光心与特征点之间连线和光轴的夹角为α₂；

设探测器下降时的俯仰角为θ，偏航角为ψ，滚转角为则所述全局特征库中的某特征在下降段图像中会产生仿射变换，相应的椭圆长、短半轴为

其它步骤及参数与具体实施方式一相同。

具体实施方式三：

本实施方式所述步骤三的具体实现过程如下：

所述全局特征库中的某特征在下降段图像中会产生旋转变换，旋转角度为

其它步骤及参数与具体实施方式一或二相同。

具体实施方式四：

本实施方式所述步骤五的具体实现过程如下：

首先根据所预测的特征尺度建立二阶高斯微分盒子滤波器D_xx、D_xy和D_yy，然后根据SURF原理建立求Hessian矩阵行列式最大值问题：

从而将上式转化为一个二次型优化问题：

其中，μ、ν分别为像元像线坐标，D_SURF(μ,ν)为SURF特征检测子在(μ,ν)处的Hessian响应值，I为将图像中与盒子滤波器卷积的对应部分重新排列成的列向量，G为相应的盒子滤波器重新排列成的行向量；

对上式中的二次型，通过进行二次型矩阵的特征值分解求取在仿射变换条件下的特征检测模板。

其它步骤及参数与具体实施方式一至三之一相同。

实施例

为验证基于仿射不变SURF模板的性能，在此选择标准SURF算法进行比较。采用三维实物地景做为图像采集样本，首先在正射视角下所拍摄的图像作为参照样本来提取SURF特征作为基准，接着变换相机的拍摄角度、高度来模拟行星着陆段相机所拍摄的序列图像，并从中提取SURF特征与样本SURF特征进行匹配。图2和图3所示为行星三维地景以及下降段图像；图2为三维模拟地景及正射影像，图3-1至图3-3为序列图像对应不同高度、拍摄角度下的下降图像的示意图；不难看出，下降图像与正射影像之间存在典型的仿射变换。

此处设计如下比较方案：首先根据多视影像及相机内外方位信息确定三维视觉数据库，并将三维绝对位置标注在正射影像上。在此基础上，根据拍摄下降图像时相机的位姿信息将数据库SURF特征投影至下降图像上，并根据投影特征与同名特征点的距离来比较两种算法的性能，图4给出了该流程示意图；如图4所示，SURF特征点l经多视影像匹配后可以获得全局位置l＝[p_x,p_y,p_z]^{G}(此处将正射视角下相机姿态指向作为全局坐标系指向)，其中[p_x,p_y,p_z]^{G}表示全局坐标系下的坐标。R,T分别为下降段对地拍摄点相机相对正射拍摄视角的姿态矩阵及平移量；假设下降段相机拍摄地点为X_L，则特征点l在F_L相机坐标系下的坐标位置为：

因此，根据上式可以推导出l在下降图像中的二维像素位置m₀：

表示二维坐标；假设在下降图像中特征l对应的同名点为m_L，则定位误差e＝|m_L-m_o|便可以用来衡量特征检测与匹配算法的精度。

图5给出了两种算法与数据库同名特征点的距离示意图，其中相机在拍摄下降图像时，其相对全局系下的姿态角为[0,40,30]deg；图5为两种SURF特征检测子仿真比较示例，其中图5(a)和图5(c)使用了本文所提出的仿射不变SURF模板进行SURF特征检测，图5(b)和图5(d)为标准的SURF检测算子。其中深色实心点(实际是红色实心点，非彩色下显示为深色实心点)为投影特征点，浅色实心点(实际是绿色实心点，非彩色下显示为浅色实心点)为同名特征点。为了使仿真结果更加清晰可辨，在此列举匹配得分前30的特征点。在本示例中，所使用的算法的平均定位误差为4.2个像素，误匹配率为3.7％。使用标准的SURF特征平均定位误差为11.6个像素，误匹配率为3.9％。秉承上述例子的设计思路，通过改变相机的拍摄地点及姿态来模拟下降图像，并统计两种算法的平均定位误差以及误匹配率，得到如下仿真结果：

如图6(a)和图6(b)所示，分别相机相对全局坐标系的三轴姿态角改变量及对应的平均定位误差及误匹配率。其中横轴代表了姿态角该变量的均值，在本次仿真试验中设定为5～40deg区间内变化。图6(a)所示为两种算法平均定位误差随姿态角改变量变化的曲线，图6(b)为误匹配率随姿态角该变量变化的曲线。随下降图像与正射影像的仿射角度改变量的增加，两种算法的性能均有所下降，但仿射不变SURF模板在特征检测计算量降低30％的同时，其定位误差与误匹配率均优于标准的SURF算法。

Claims (4)

1.一种基于尺度预测的图像局部区域特征提取方法，其特征在于其包括以下步骤：

步骤一：根据行星着陆过程中探测器相对于预定着陆点的位置、探测器体坐标系相对于在轨拍摄行星表面图像是的姿态及相机焦距、视场角信息，对当前所拍摄图像在目标天体表面的位置进行初步估计，选定全局特征库中的搜索范围；

步骤二：根据探测器拍照时的位姿信息以及全局特征库中相应特征点的特征尺度，预测该特征在拍摄图像中的特征尺度；

步骤三：根据探测器拍照时的姿态信息以及全局特征库中相应特征点的方向信息，预测该特征在拍摄图像中的旋转角度；

步骤四：将拍摄的图像按照所预测的旋转角度旋转获得旋转图像；

步骤五：根据二次型求极大值方法仿照SURF检测原理以及预测的特征尺度生成特征检测模板；

步骤六：将特征检测模板与旋转图像卷积并在位置空间进行非极大值抑制从而实现特征提取。

2.根据权利要求1所述的一种基于尺度预测的图像局部区域特征提取方法，其特征在于所述步骤二的具体实现过程如下：

建立全局特征库拍摄图像时探测器与该特征的距离为d₁，建立全局特征库拍摄图像时所用相机的焦距为f₁，建立全局特征库拍摄图像时相机光心与特征点之间连线和光轴的夹角为α₁，设全局特征库中的某特征的特征尺度为σ₁，则该特征在下降段图像中会产生放缩变换，其尺度变换为

其中，下标2对应的参数表示相应物理量在下降段拍摄图像时的值，即：下降段拍摄图像时探测器与该特征的距离为d₂，下降段拍摄图像时所用相机的焦距为f₂，下降段拍摄图像时相机光心与特征点之间连线和光轴的夹角为α₂；

设探测器下降时的俯仰角为θ，偏航角为ψ，滚转角为则所述全局特征库中的某特征在下降段图像中会产生仿射变换，相应的椭圆长、短半轴为

3.根据权利要求2所述的一种基于尺度预测的图像局部区域特征提取方法，其特征在于所述步骤三的具体实现过程如下：

所述全局特征库中的某特征在下降段图像中会产生旋转变换，旋转角度为

4.根据权利要求2所述的一种基于尺度预测的图像局部区域特征提取方法，其特征在于所述步骤五的具体实现过程如下：

首先根据所预测的特征尺度建立二阶高斯微分盒子滤波器D_xx、D_xy和D_yy，然后根据SURF原理建立求Hessian矩阵行列式最大值问题：

从而将上式转化为一个二次型优化问题：

约束条件：||I||≤1

其中，μ、ν分别为像元像线坐标，D_SURF(μ,ν)为SURF特征检测子在(μ,ν)处的Hessian响应值，I为将图像中与盒子滤波器卷积的对应部分重新排列成的列向量，G为相应的盒子滤波器重新排列成的行向量；

通过进行二次型矩阵的特征值分解求取在仿射变换条件下的特征检测模板。