耦合点线特征的多模态亚米/米级卫星影像配准技术

技术领域

本发明涉及遥感图像处理领域，更具体的说是涉及耦合点线特征的多模态亚米/米级卫星影像配准技术。

背景技术

随着航天技术、计算机技术、传感器技术的快速发展，高空间分辨率以及广覆盖的卫星遥感数据层出不穷。利用多模态(多源多尺度多时相等)高空间分辨率卫星影像实现协同对地观测的研究已成为当前遥感学界的研究热点。在实际应用中，由于天气状况的多变性、地形的复杂性以及应用的多样性，单一传感器类型的遥感数据往往不能提供足够的信息以满足应用的需求，综合利用不同源遥感卫星获取的数据，可以较为全面地刻画观测对象的几何和物理特性，因此需要将多种传感器获取的多角度、多尺度、多时相的卫星影像融合在一起，得到更加丰富完整的信息，提取更为精确的信息指数，从而快速地制定出更为合理的决策方案。

遥感数据应用是遥感的最终目的，为了提高对海量遥感数据的处理速度和精度，研究和开发高效的遥感影像处理技术显得尤为重要。由于遥感卫星搭载的各种传感器之间的物理特性和成像方式不同，而且成像位置和成像角度可能存在差异，因此，其拍摄的影像存在灰度、分辨率、位移和旋转角度等的不同。然而，多模态影像综合利用是发展趋势，如多传感器信息融合、变化检测、影像镶嵌及信息验证分析等，应用中需要对多模态影像的位置进行对齐，因此，影像配准是其中必不可少的关键技术。但是，随着影像空间分辨率的提高，影像配准的干扰因素增多，尤其是亚米/米级影像的全自动配准成为行业的难题。

针对多模态高分辨率卫星影像间非线性辐射及尺度差异较大造成配准精度较低的问题，本发明公开一种耦合点线特征的多模态亚米/米级卫星影像高精度配准技术。首先，该技术采用Bhattacharyya距离和特征斜率约束优化SIFT算法以获取更多优质同名点；其次，在多尺度空间下使用CannyLines算法提取多尺度线特征；然后，通过每条线特征与由SIFT同名点构建的Delaunay三角网的空间几何信息构建一种可以抵抗多源影像间非线性辐射差异和尺度差异的全局线特征描述符，并且利用设计的成本函数获取正确的线匹配对；最后，使用多种约束条件获取定位准确的线特征交点，并采用耦合SIFT同名点和线特征交点的方式迭代求解高精度的配准模型参数。

发明内容

有鉴于此，本发明提出一种耦合点线特征的多模态亚米/米级卫星影像配准技术，图1为本发明多模态亚米/米级卫星影像配准技术的总体流程，包括基于改进SIFT的点特征获取、线特征提取与描述、特征耦合和影像配准四部分。

为了实现上述目的，本发明采用如下技术方案：

耦合点线特征的多模态亚米/米级卫星影像配准技术，包括：基于改进SIFT的点特征获取、线特征提取与描述、特征耦合和影像配准四部分；

1.基于改进SIFT的点特征获取步骤如下：

S1：首先采用双线性插值的方式对原始影像进行降采样以构建影像金字塔，利用SIFT在不同尺度的影像上提取稳定的点特征，然后对这些稳定的点特征进行描述；

S2：以Bhattacharyya距离为相似性度量，使用FLANN方法进行特征粗匹配，获取初始匹配点集；

S3：计算初始匹配结果的斜率以求得置信区间，并保留位于该置信区间内的匹配点，等待下一步检验；

S4：利用以单应矩阵为几何约束模型的PROSAC方法对上一步得到的匹配点集进行检验，进一步剔除外点，从而获取可靠度较高的优质匹配点集；

2.线特征提取与描述步骤如下：

S5：利用CannyLines检测器提取多尺度线特征；

S6：基于空间几何信息构建线特征描述符；

S7：基于定义的匹配代价函数进行线特征匹配，得到同名线特征；

3.特征耦合步骤如下：

S8：利用得到的同名线特征进行交点拟合；

S9：将拟合得到的同名交点与SIFT同名点特征合并为一个同名特征集；

4.影像配准步骤如下：

S10：利用同名特征集迭代求解影像的配准模型参数，得到配准结果；

S11：结束。

本发明技术先进、科学，可以很好地抵抗多模态高分卫星影像间较大的非线性辐射及尺度差异，且能够实现多模态亚米/米级卫星影像的高精度自动配准。通过实验表明，该技术具有很好的普适性，计算速度很快，效率很高，配准误差优于1个像素，具有较好的使用价值。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施示例或现有技术描述中所需要使用的附图作简单介绍，显而易见地，下面描述中的附图仅仅是本发明的示意图，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。

图1为本发明提供的耦合点线特征的多模态亚米/米级卫星影像配准技术的流程图；

图2为本发明提出的一种基于空间几何信息的线特征描述符；

图3为本发明提出的线特征交点拟合方法；

图4为本发明提供的实验所用的多模态亚米/米级高分卫星影像数据；

图5为本发明提供的配准结果图；

表1为本发明提供的实验所用影像数据说明。

表2为本发明提供的实验结果。

表1

表2

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

以下步骤为基于改进SIFT算法的点特征获取部分：

Step 1：首先采用双线性插值的方式对原始影像进行降采样以构建影像金字塔，利用SIFT在不同尺度的影像上提取稳定的点特征，然后对这些稳定的点特征进行描述；

Step 2：以Bhattacharyya距离为相似性度量，使用FLANN方法进行特征粗匹配，获取初始匹配点集。Bhattacharyya距离定义如下：

其中，X和Y为两个L1归一化的特征向量。

Step 3：计算初始匹配结果的斜率以求得置信区间，并保留位于该置信区间内的匹配点，等待下一步检验。特征斜率的定义如下：

其中，(x

Step 4：利用以单应矩阵为几何约束模型的PROSAC方法对上一步得到的匹配点集进行检验，进一步剔除外点，从而获取可靠度较高的优质SIFT匹配点集。

以下步骤为线特征提取和描述部分：

Step 5：利用CannyLines检测器提取多尺度线特征。为了确保提取稳定的多尺度线特征，通过不同尺度下线特征的梯度方向是否一致进行稳定性判定，线特征上每个像素的梯度方向的定义如下：

其中，L(i，j，σ)为像素(i，j)在图像高斯尺度空间的取值。

Step 6：基于空间几何信息构建线特征描述符。线要素l

DLO(i，j)＝{D(i，j)，L(i，j)，O(i，j)}

其中，d

Step 7：基于定义的匹配代价函数进行线特征匹配，得到同名线特征。匹配代价函数计算公式如下：

式中，其中，h

以下步骤为特征耦合部分：

Step 8：通过多种约束条件，利用得到的同名线特征进行交点拟合。约束条件包括：①区域约束：交点必须位于线特征的搜索区域内；②夹角约束：两线特征的夹角必须大于θ

Step 9：将拟合得到的同名交点与SIFT同名点特征合并为一个同名特征集；

以下步骤为影像配准部分：

Step 10：利用同名特征集迭代求解影像的配准模型参数，得到配准结果。配准模型定义如下：

其中，其中(x,y)和(x′＝u/w′,y′＝v/w′)分别表示两个影像上同名点的坐标，a

Step 11：结束。

对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，该方案所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于该方案所示的这些实施例，而是要符合与本发明所公开的原理和新颖特点相一致的最宽的范围。