一种基于SIFT的图像配准方法、系统及存储介质

技术邻域

本发明涉及图像处理技术邻域，更具体的，涉及一种基于SIFT的图像配准 方法、系统及存储介质。

背景技术

SIFT算法具有良好的稳定性，在基于特征的图像匹配算法中拥有十分广泛 的应用，SIFT描述符对在原有算法的基础上提高了算法的图像尺度不变性和旋 转不变性，且特征都具有较高的唯一性，能在不同视角、不同拍摄距离、不同 仿射变换、不同光照条件下都保持很好的匹配鲁棒性。

但传统的SIFT算法仍然存在计算复杂程度高、匹配耗时较长、无法实时计 算以及无法满足人们日益增长的对匹配精度的需求、抗噪性仍有待加强等缺陷， 因此从降低算法的复杂程度、提高算法的运算速度、提高算法的匹配精度和抗 噪性等方面对SIFT算法进行改进优化成为了当下主流的研究方向。

在现有理论的基础上针对SIFT算法中存在的计算复杂程度高、匹配耗时较 长、无法实时计算等问题，拟提出一种改进的算法，通过减少需要匹配的特征 点数量和需要匹配的区域以及描述符向量维度，提高匹配效率和算法的速度。 并在现有理论的基础上针对SIFT算法中存在的精度较低、抗噪性仍有待加强等 缺陷，提出相应的改进方案，使其更加贴合当今社会的需求。

发明内容

为了解决上述技术问题，本发明提出了一种基于SIFT的图像配准方法、系 统及存储介质。

本发明第一方面提供了一种基于SIFT的图像配准方法，包括：

获取待配准图像及参考图像，将待配准图像进行预处理，构建高斯差分尺 度空间；

获取所述高斯差分尺度空间中不同尺度空间中的极值点，根据所述极值点 选取特征点，并为所述特征点生成特征描述符；

计算特征点领域的局部信息熵，根据所述局部信息熵与特征描述符确定特 征点的最终特征描述符；

根据遗传算法优化特征点匹配，生成匹配点对，并对匹配点对进行校正， 完成待配准图像的图像配准。

本方案中，所述的获取所述高斯差分尺度空间中不同尺度空间中的极值点， 根据所述极值点选取特征点，具体为：

对高斯差分尺度空间进行精简，去除高斯差分尺度空间中第一组的第一层 级尺度空间，通过精简后的高斯差分尺度空间进行极值点的检测；

通过将待检测点与邻近像素点的位置与尺度信息进行比较，选取极值点， 并获取极值点的对比度；

预设第一对比度阈值，将极值点的对比度与所述第一对比度阈值对比分析， 将对比度大于第一对比度阈值的极值点作为待选特征点；

同时，预设第二对比度阈值，所述第二对比度阈值大于第一对比度阈值， 将对比度大于第二对比度阈值的极值点继续保存为待选特征点；

获取所述待选特征点的Hessian矩阵，去除不稳定的边缘响应点，根据所 述待选特征点生成特征点集。

本方案中，为所述特征点生成特征描述符，具体为：

计算特征点集中特征点的尺度信息获取高斯平滑图像，获取特征点圆形邻 域区域内高斯平滑图像的图像梯度的幅值及幅角；

通过直方图统计特征点邻域内像素的幅值和梯度方向，获取直方图的峰值 作为特征点的主方向信息；

根据特征点的主方向信息将特征点邻域进行旋转，将特征点邻域划分为若 干子区域，在各子区域中进行直方图统计生成特征向量，将所述特征向量进行 降维获取特征描述符。

本方案中，所述的计算特征点领域的局部信息熵，根据所述局部信息熵与 特征描述符确定特征点的最终特征描述符，具体为：

获取特征点领域中各子区域的局部信息熵，预设局部信息熵阈值，将各子 区域的局部信息熵与所述局部信息熵阈值进行对比判断；

若小于，则将对应子区域中的特征点舍弃，若大于，计算对应子区域的局 部信息熵与整个尺度空间的信息熵总和的比值信息；

根据所述比值信息生成权重信息，将所述权重信息与特征描述符结合生成 特征点的最终特征描述符。

本方案中，所述的根据遗传算法优化特征点匹配，生成匹配点对，并对匹 配点对进行校正，完成待配准图像的图像配准，具体为

根据特征点的最终特征描述符进行染色体编码表示特征点的特征向量，并 初始化设置种群参数；

通过特征点间的欧式距离作为适应度函数，进行迭代训练，获取染色体的 适应度信息，选择、交叉变异产生下一代染色体；

选取适应度信息最小的染色体进行解码获取匹配点的信息，同时根据 RANSAC算法进行校正，消除误匹配点，进行待配准图像的图像配准。

本发明第二方面还提供了一种基于SIFT的图像配准系统，该系统包括：存 储器、处理器，所述存储器中包括一种基于SIFT的图像配准方法程序，所述一 种基于SIFT的图像配准方法程序被所述处理器执行时实现如下步骤：

获取待配准图像及参考图像，将待配准图像进行预处理，构建高斯差分尺 度空间；

获取所述高斯差分尺度空间中不同尺度空间中的极值点，根据所述极值点 选取特征点，并为所述特征点生成特征描述符；

计算特征点领域的局部信息熵，根据所述局部信息熵与特征描述符确定特 征点的最终特征描述符；

根据遗传算法优化特征点匹配，生成匹配点对，并对匹配点对进行校正， 完成待配准图像的图像配准。

本方案中，所述的获取所述高斯差分尺度空间中不同尺度空间中的极值点， 根据所述极值点选取特征点，具体为：

对高斯差分尺度空间进行精简，去除高斯差分尺度空间中第一组的第一层 级尺度空间，通过精简后的高斯差分尺度空间进行极值点的检测；

通过将待检测点与邻近像素点的位置与尺度信息进行比较，选取极值点， 并获取极值点的对比度；

预设第一对比度阈值，将极值点的对比度与所述第一对比度阈值对比分析， 将对比度大于第一对比度阈值的极值点作为待选特征点；

同时，预设第二对比度阈值，所述第二对比度阈值大于第一对比度阈值， 将对比度大于第二对比度阈值的极值点继续保存为待选特征点；

获取所述待选特征点的Hessian矩阵，去除不稳定的边缘响应点，根据所 述待选特征点生成特征点集。

本方案中，为所述特征点生成特征描述符，具体为：

计算特征点集中特征点的尺度信息获取高斯平滑图像，获取特征点圆形邻 域区域内高斯平滑图像的图像梯度的幅值及幅角；

通过直方图统计特征点邻域内像素的幅值和梯度方向，获取直方图的峰值 作为特征点的主方向信息；

根据特征点的主方向信息将特征点邻域进行旋转，将特征点邻域划分为若 干子区域，在各子区域中进行直方图统计生成特征向量，将所述特征向量进行 降维获取特征描述符。

本方案中，所述的计算特征点领域的局部信息熵，根据所述局部信息熵与 特征描述符确定特征点的最终特征描述符，具体为：

获取特征点领域中各子区域的局部信息熵，预设局部信息熵阈值，将各子 区域的局部信息熵与所述局部信息熵阈值进行对比判断；

若小于，则将对应子区域中的特征点舍弃，若大于，计算对应子区域的局 部信息熵与整个尺度空间的信息熵总和的比值信息；

根据所述比值信息生成权重信息，将所述权重信息与特征描述符结合生成 特征点的最终特征描述符。

本方案中，所述的根据遗传算法优化特征点匹配，生成匹配点对，并对匹 配点对进行校正，完成待配准图像的图像配准，具体为

根据特征点的最终特征描述符进行染色体编码表示特征点的特征向量，并 初始化设置种群参数；

通过特征点间的欧式距离作为适应度函数，进行迭代训练，获取染色体的 适应度信息，选择、交叉变异产生下一代染色体；

选取适应度信息最小的染色体进行解码获取匹配点的信息，同时根据 RANSAC算法进行校正，消除误匹配点，进行待配准图像的图像配准。

本发明第三方面还提供一种计算机可读存储介质，所述计算机可读存储介 质中包括一种基于SIFT的图像配准方法程序，所述一种基于SIFT的图像配准 方法程序被处理器执行时，实现如上述任一项所述的一种基于SIFT的图像配准 方法的步骤。

本发明公开了一种基于SIFT的图像配准方法、系统及存储介质，包括：获 取待配准图像及参考图像，将待配准图像进行预处理，构建高斯差分尺度空间； 获取所述高斯差分尺度空间中不同尺度空间中的极值点，根据所述极值点选取 特征点，并为所述特征点生成特征描述符；计算特征点领域的局部信息熵，根 据所述局部信息熵与特征描述符确定特征点的最终特征描述符；根据遗传算法 优化特征点匹配，生成匹配点对，并对匹配点对进行校正，完成待配准图像的 图像配准。本发明通过精简需要匹配的特征点数量和需要匹配的区域，优化了 匹配效率和算法的速度，提高了SIFT算法搜索效率及配准精度。

附图说明

图1示出了本发明一种基于SIFT的图像配准方法的流程图；

图2示出了本发明根据遗传算法优化特征点匹配方法的流程图；

图3示出了本发明一种基于SIFT的图像配准系统的框图。

具体实施方式

为了能够更清楚地理解本发明的上述目的、特征和优点，下面结合附图和 具体实施方式对本发明进行进一步的详细描述。需要说明的是，在不冲突的情 况下，本申请的实施例及实施例中的特征可以相互组合。

在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发 明还可以采用其他不同于在此描述的其他方式来实施，因此，本发明的保护范 围并不受下面公开的具体实施例的限制。

图1示出了本发明一种基于SIFT的图像配准方法的流程图。

如图1所示，本发明第一方面提供了一种基于SIFT的图像配准方法，包括：

S102，获取待配准图像及参考图像，将待配准图像进行预处理，构建高斯 差分尺度空间；

S104，获取所述高斯差分尺度空间中不同尺度空间中的极值点，根据所述 极值点选取特征点，并为所述特征点生成特征描述符；

S106，计算特征点领域的局部信息熵，根据所述局部信息熵与特征描述符 确定特征点的最终特征描述符；

S108，根据遗传算法优化特征点匹配，生成匹配点对，并对匹配点对进行 校正，完成待配准图像的图像配准。

需要说明的是，由于拍摄设备和环境等因素的影响，待配准图像通常质量 不高，存在噪声等额外因素影响，在进行特征提取前先对参考图像和待配准图 像通过灰度变换等滤波去噪方式进行预处理，为了提高检测的效率及准确度， 将生成无用特征点最多的高斯差分尺度空间的第一组第一层跳过，对高斯差分 尺度空间进行精简，去除高斯差分尺度空间中第一组的第一层级尺度空间，通 过精简后的高斯差分尺度空间进行极值点的检测；通过将待检测点与邻近像素 点的位置与尺度信息进行比较，选取极值点，并获取极值点的对比度；预设第 一对比度阈值，将极值点的对比度与所述第一对比度阈值对比分析，将对比度 大于第一对比度阈值的极值点作为待选特征点；同时，预设第二对比度阈值， 所述第二对比度阈值大于第一对比度阈值，将对比度大于第二对比度阈值的极 值点继续保存为待选特征点；获取所述待选特征点的Hessian矩阵，去除不稳 定的边缘响应点，根据所述待选特征点生成特征点集。

需要说明的是，为所述特征点生成特征描述符，具体为：计算特征点集中 特征点的尺度信息获取高斯平滑图像，获取特征点圆形邻域区域内高斯平滑图 像的图像梯度的幅值及幅角；通过直方图统计特征点邻域内像素的幅值和梯度 方向，获取直方图的峰值作为特征点的主方向信息；直方图的横轴是梯度方向 角，纵轴是梯度方向角对应的梯度幅值累加，直方图将0-360度的范围分为36 个柱，每10度一个柱；根据特征点的主方向信息将特征点邻域进行旋转，将特 征点邻域划分为若干子区域，在各子区域中进行直方图统计生成特征向量，将 所述特征向量进行降维获取特征描述符。

需要说明的是，所述的计算特征点领域的局部信息熵，根据所述局部信息 熵与特征描述符确定特征点的最终特征描述符，具体为：获取特征点领域中各 子区域的局部信息熵，预设局部信息熵阈值，将各子区域的局部信息熵与所述 局部信息熵阈值进行对比判断；若小于，则将对应子区域中的特征点舍弃，若 大于，计算对应子区域的局部信息熵与整个尺度空间的信息熵总和的比值信息； 根据所述比值信息生成权重信息，将所述权重信息与特征描述符结合生成特征 点的最终特征描述符。通过局部信息熵来表征特征点邻域内像素灰度值分布的 混乱程度，其中局部信息熵的计算公式为：

其中，s

图2示出了本发明根据遗传算法优化特征点匹配方法的流程图。

根据本发明实施例，所述的根据遗传算法优化特征点匹配，生成匹配点对， 并对匹配点对进行校正，完成待配准图像的图像配准，具体为

S202，根据特征点的最终特征描述符进行染色体编码表示特征点的特征向 量，并初始化设置种群参数；

S204，通过特征点间的欧式距离作为适应度函数，进行迭代训练，获取染 色体的适应度信息，选择、交叉变异产生下一代染色体；

S206，选取适应度信息最小的染色体进行解码获取匹配点的信息，同时根 据RANSAC算法进行校正，消除误匹配点，进行待配准图像的图像配准。

需要说明的是，特征点的匹配通过欧氏距离计算相似性进行判断，选择某 个特征点在待配准图像中进行搜索，获取欧氏距离符合预设条件的特征点，获 取匹配点对，根据遗传算法优化特征点匹配，通过实体编码染色体长度表示特 征点的特征向量，初始化参数，设置种群规模大小为20，最大迭代次数为200， 根据特征点的欧氏距离计算公式作为适应度函数，计算种群中染色体的适应度， 适应度适应度大的染色体被移除，并在迭代训练中染色体之间进行交叉杂交、 变异，以至于快速找到最小欧式距离的特征点。

图3示出了本发明一种基于SIFT的图像配准系统的框图。

本发明第二方面还提供了一种基于SIFT的图像配准系统3，该系统包括： 存储器31、处理器32，所述存储器中包括一种基于SIFT的图像配准方法程序， 所述一种基于SIFT的图像配准方法程序被所述处理器执行时实现如下步骤：

获取待配准图像及参考图像，将待配准图像进行预处理，构建高斯差分尺 度空间；

获取所述高斯差分尺度空间中不同尺度空间中的极值点，根据所述极值点 选取特征点，并为所述特征点生成特征描述符；

计算特征点领域的局部信息熵，根据所述局部信息熵与特征描述符确定特 征点的最终特征描述符；

根据遗传算法优化特征点匹配，生成匹配点对，并对匹配点对进行校正， 完成待配准图像的图像配准。

需要说明的是，由于拍摄设备和环境等因素的影响，待配准图像通常质量 不高，存在噪声等额外因素影响，在进行特征提取前先对参考图像和待配准图 像通过灰度变换等滤波去噪方式进行预处理，为了提高检测的效率及准确度， 将生成无用特征点最多的高斯差分尺度空间的第一组第一层跳过，对高斯差分 尺度空间进行精简，去除高斯差分尺度空间中第一组的第一层级尺度空间，通 过精简后的高斯差分尺度空间进行极值点的检测；通过将待检测点与邻近像素 点的位置与尺度信息进行比较，选取极值点，并获取极值点的对比度；预设第 一对比度阈值，将极值点的对比度与所述第一对比度阈值对比分析，将对比度 大于第一对比度阈值的极值点作为待选特征点；同时，预设第二对比度阈值， 所述第二对比度阈值大于第一对比度阈值，将对比度大于第二对比度阈值的极 值点继续保存为待选特征点；获取所述待选特征点的Hessian矩阵，去除不稳 定的边缘响应点，根据所述待选特征点生成特征点集。

需要说明的是，为所述特征点生成特征描述符，具体为：计算特征点集中 特征点的尺度信息获取高斯平滑图像，获取特征点圆形邻域区域内高斯平滑图 像的图像梯度的幅值及幅角；通过直方图统计特征点邻域内像素的幅值和梯度 方向，获取直方图的峰值作为特征点的主方向信息；直方图的横轴是梯度方向 角，纵轴是梯度方向角对应的梯度幅值累加，直方图将0-360度的范围分为36 个柱，每10度一个柱；根据特征点的主方向信息将特征点邻域进行旋转，将特 征点邻域划分为若干子区域，在各子区域中进行直方图统计生成特征向量，将 所述特征向量进行降维获取特征描述符。

需要说明的是，所述的计算特征点领域的局部信息熵，根据所述局部信息 熵与特征描述符确定特征点的最终特征描述符，具体为：获取特征点领域中各 子区域的局部信息熵，预设局部信息熵阈值，将各子区域的局部信息熵与所述 局部信息熵阈值进行对比判断；若小于，则将对应子区域中的特征点舍弃，若 大于，计算对应子区域的局部信息熵与整个尺度空间的信息熵总和的比值信息； 根据所述比值信息生成权重信息，将所述权重信息与特征描述符结合生成特征 点的最终特征描述符。通过局部信息熵来表征特征点邻域内像素灰度值分布的 混乱程度，其中局部信息熵的计算公式为：

其中，s

根据本发明实施例，所述的根据遗传算法优化特征点匹配，生成匹配点对， 并对匹配点对进行校正，完成待配准图像的图像配准，具体为

根据特征点的最终特征描述符进行染色体编码表示特征点的特征向量，并 初始化设置种群参数；

通过特征点间的欧式距离作为适应度函数，进行迭代训练，获取染色体的 适应度信息，选择、交叉变异产生下一代染色体；

选取适应度信息最小的染色体进行解码获取匹配点的信息，同时根据 RANSAC算法进行校正，消除误匹配点，进行待配准图像的图像配准。

需要说明的是，特征点的匹配通过欧氏距离计算相似性进行判断，选择某 个特征点在待配准图像中进行搜索，获取欧氏距离符合预设条件的特征点，获 取匹配点对，根据遗传算法优化特征点匹配，通过实体编码染色体长度表示特 征点的特征向量，初始化参数，设置种群规模大小为20，最大迭代次数为200， 根据特征点的欧氏距离计算公式作为适应度函数，计算种群中染色体的适应度， 适应度适应度大的染色体被移除，并在迭代训练中染色体之间进行交叉杂交、 变异，以至于快速找到最小欧式距离的特征点。

本发明第三方面还提供一种计算机可读存储介质，所述计算机可读存储介 质中包括一种基于SIFT的图像配准方法程序，所述一种基于SIFT的图像配准 方法程序被处理器执行时，实现如上述任一项所述的一种基于SIFT的图像配准 方法的步骤。

在本申请所提供的几个实施例中，应该理解到，所揭露的设备和方法，可 以通过其它的方式实现。以上所描述的设备实施例仅仅是示意性的，例如，所 述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式， 如：多个单元或组件可以结合，或可以集成到另一个系统，或一些特征可以忽 略，或不执行。另外，所显示或讨论的各组成部分相互之间的耦合、或直接耦 合、或通信连接可以是通过一些接口，设备或单元的间接耦合或通信连接，可 以是电性的、机械的或其它形式的。

上述作为分离部件说明的单元可以是、或也可以不是物理上分开的，作为 单元显示的部件可以是、或也可以不是物理单元；既可以位于一个地方，也可 以分布到多个网络单元上；可以根据实际的需要选择其中的部分或全部单元来 实现本实施例方案的目的。

另外，在本发明各实施例中的各功能单元可以全部集成在一个处理单元中， 也可以是各单元分别单独作为一个单元，也可以两个或两个以上单元集成在一 个单元中；上述集成的单元既可以采用硬件的形式实现，也可以采用硬件加软 件功能单元的形式实现。

本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可 以通过程序指令相关的硬件来完成，前述的程序可以存储于计算机可读取存储 介质中，该程序在执行时，执行包括上述方法实施例的步骤；而前述的存储介 质包括：移动存储设备、只读存储器(ROM，Read-Only Memory)、随机存取存 储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码 的介质。

或者，本发明上述集成的单元如果以软件功能模块的形式实现并作为独立 的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。基于这样 的理解，本发明实施例的技术方案本质上或者说对现有技术做出贡献的部分可 以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包 括若干指令用以使得一台计算机设备(可以是个人计算机、服务器、或者网络 设备等)执行本发明各个实施例所述方法的全部或部分。而前述的存储介质包 括：移动存储设备、ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。

以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于 此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到 变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应 以所述权利要求的保护范围为准。