专利名称：多特征多级别的可见光与红外图像高精度配准方法
技术领域：
本发明属于图像处理技术领域，具体涉及遥感图像的配准方法，特别是一种可见光与红外图像的配准方法。本发明可以广泛适用于航天、航空传感器平台获取的遥感影像的配准方法。
背景技术：
红外图像的光谱位于可见光波段之外，包含了很多可见光波段之外的很多特征，这些特征对于准确识别目标具有重要的意义。例如，处于机动状态和休息状态的同一辆坦克，在可见光图像上无法区分，而通过红外图像的反映的温度信息可以正确识别坦克的状态。因此，综合利用不同传感器的互补特性，对可见光与红外 图像进行融合对情报获取和分析具有重要的作用。作为可见光与红外图像融合的前提，可见光与红外图像必须进行自动、高精度配准。可见光与红外图像高精度自动配准技术是制约可见光与红外图像融合应用的关键技术，具有重要的应用价值。然而，由于成像机理的差异，同一物体在可见光与红外图像上的表现差距很大，这种差距给传统的特征匹配方法带来了很大的挑战。实际上，即使让判图专家在可见光与红外图像手工标定快速一些控制点对也是一件很困难的事情。正是因为可见光与红外图像的差异较大，目前可见光与红外图像配准基本上仍处于研究阶段，尚未有专门的可见光-红外图像配准软件。尽管商业软件如Erdas和Envi都提供了配准功能，但这些配准功能都是为可见光图像设计的。在这些软件上对可见光和红外图像进行配准，成功的概率很小。目前，对可见光与红外图像配准的研究大都集中在相似性度量，其中基于互信息的配准方法最受青睐。然而，基于互信息的可见光与红外图像配准方法有致命的缺陷对初始值依赖严重。若没有好的初始值，极易陷入局部极小值或不收敛。该缺陷严重制约了其实用性，因为自动从可见光与红外图像上获取这种初始值信息是很困难的。在这种情况下，必须另辟思路研究新的配准方法，在保证配准精度的同时，尽可能提供其自动化程度和计算速度。
发明内容
(一 )要解决的技术问题本发明所要解决的技术问题是针对可见光图像与红外图像的互补性而提出一种特别适合于可见光与红外图像的自动、高效、高精度的配准方法，以推动遥感识别技术的发展。( 二 )技术方案本发明提出一种可见光图像与红外图像的配准方法，包括如下步骤步骤SI、对可见光图像和红外图像进行多尺度分解，得到低分辨率的可见光图像和红外图像；步骤S2、在所述低分辨率图像对上提取SIFT特征、多尺度角点特征和面点特征，匹配SIFT特征并去除外点，利用匹配的SIFT特征对得到初始变换模型，所谓变换模型是指红外图像和高光谱图像之间的一种几何变换关系；步骤S3、在原始图像对上利用初始变换模型提供的几何约束，进行基于图像块对的多尺度角点特征、面点特征提取，根据初始变换模型和多尺度角点特征、面点特征确定更精确的变换模型；步骤S4、根据所述变换模型对红外图像进行变换，得到变换后的红外图像。(三)有益效果本发明的可见光与红外图像高精度自动配准方法基于多特征多级别的配准，综合利用了模式识别、计算机视觉领域中最新的特征提取和特征匹配技术及并行计算技术，有效地解决了可见光与红外图像的配准问题，弥补了现有配准软件对可见光与红外图像配准方面的缺陷。本发明的配准方法对可见光与红外图像配准具有很好的通用性和实用性，能够大 大推动可见光与红外图像的广泛应用，具有很好的经济效益。
图I是本发明的系统框图；图2是本发明的方法的一个具体实施例的原始可见光图像和生成的低分辨率可见光图像的例图；图3是本发明的方法的一个具体实施例的原始红外图像与生成的低分辨率红外图像的例图；图4SIFT特征提取流程图；图5是本发明的方法的一个具体实施例的SIFT特征提取中DOG构建图例；图6是本发明的方法的一个具体实施例的SIFT特征提取中极值点提取图例；图7是本发明的方法的一个具体实施例的SIFT特征提取中主方向检测图例；图8是本发明的一个具体实施例的为SIFT特征构造描述向量的一个示意图；图9是本发明的一个具体实施例的粗配准流程；图10是本发明的一个具体实施例的⑶BICP核心流程。
具体实施例方式为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明作进一步的详细说明。本发明的多特征多级别的可见光与红外图像的高精度配准方法既可以硬件方式实现，也可以软件方式实现。例如在个人计算机、工控机及服务器上以软件的形式安装并执行，也可本发明的方法做成嵌入式芯片以硬件的形式来体现。下面参照附图对本发明的具体实施方式
进行说明。在下面的描述中，所述的“图像”特指通过遥感设备获得的遥感图像，并且是已进行了数字化的数字图像。然而，本发明并不限于遥感图像，对于其它技术领域中获得的需要配准的可见光图像和红外图像，本发明也可适用。图I给出了一个实施例的系统框图。如图I所示，总的来说，本发明的方法包括如下步骤
步骤SI :多尺度分解。对可见光图像和红外图像进行多尺度分解，得到低分辨率的可见光图像和红外图像；根据本发明的一种具体实施方式
，该实施例中的多尺度分解通过降采样实现。虽然多尺度分解可以由小波金字塔、高斯金字塔等多种方法实现，但计算量比降采样大的多；而在精配准阶段，SIFT特征还要在原始图像上提取。因此，基于降采样的多尺度分解方式在保持配准精度的同时，还大大减少了计算量。步骤S2:粗配准。在低分辨率图像对上提取SIFT(Scale_invariant FeatureTransform,尺度不变特征变换)特征、多尺度角点特征和面点特征,采用最近邻比值法匹配SIFT特征并去除外点，利用匹配的SIFT特征对和GDBICP法得到初始变换模型；根据本发明的一种具体实施方式
，角点特征分为主动匹配角点特征和被动匹配角点特征，面点特征分为主动匹配面点特征和被动匹配面点特征。主动匹配角(面)点特征是参考图像(如可见光图像)上需要匹配的角(面)点特征，被动匹配角(面)点特征是 待配准图像(如红外图像)上可能匹配的角(面)点特征。主动匹配角(面)点特征是被动匹配角(面)点特征的子集。根据本发明的一种具体实施方式
，采用最近邻比值法来匹配SIFT特征，并去除外点，外点指不满足上述变换模型的SIFT特征对；根据本发明的一种具体实施方式
，利用⑶BICP(Generalized Dual BootstrapIterative Closest Point)法得到初始变换模型,所谓变换模型是指红外图像和高光谱图像之间的一种几何变换关系。变换模型包含变换类型和变换参数两层含义。同样，变换模型的确定也包括变换类型选择和变换参数求解。上述变换模型是后续配准变换模型的初始值即初始变换模型，在精配准阶段变换模型的类型及其参数都将不断调整。步骤S3 :精配准。在所述原始图像对上利用初始变换模型提供的几何约束，进行基于图像块对的多尺度角点特征、面点特征提取，根据初始变换模型和多尺度角点特征、面点特征确定更精确的变换模型。在更精确的变换模型下，红外图像和显著波段图像的多尺度角/面点特征之间的拟合误差更小。根据本发明的一种具体实施方式
，利用迭代重加权最小二乘法选择变换类型并求取变换参数。步骤S4 :图像变换。根据所述变换模型对红外图像进行变换，得到变换后的红外图像。根据本发明的一种具体实施方式
，根据变换模型对红外图像进行双三次样条插值，然后根据变换模型及可见光图像和红外图像的尺寸计算参考图像和变换后的红外图像的最大重叠区域的坐标，取出、保存所述最大重叠区域的坐标内的可见光图像和红外图像。为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，对本发明作进一步的详细说明。在该实施例中，步骤SI利用双线性插值方法对可见光图像和红外图像进行降采样，生成低分辨率可见光图像和低分辨率红外图像。图2为原始可见光图像和生成的低分辨率可见光图像的例图，其中上图为原始可见光图像，下图为低分辨率可见光图像。图3为原始红外图像和生成的低分辨率红外图像的例图，其中上图为原始红外图像，下图为低分辨率红外图像。图4是本发明的方法的具体实施例中步骤S2的提取SIFT特征的流程图。在该实施例中，SIFT特征的提取是按照如下流程实现的步骤S21、构建图像的DOG(Difference of Gaussian,高斯差分)金字塔。设图像为1(1,7),贝1]006金字塔中第1^层高斯差分图像为0(1,7,1^0) = L(x, y,
IX2+.V2
(k+1) o )-L(x, y, k o )，其中，L(x, y, o ) = G(x, y, o )*I (x, y) y,a) =---e :,T" *
Incr，
表示卷积运算。图5是DOG构建图例，如图5所示，图像在5个尺度上进行高斯卷积，得到的DOG金字塔中有4个图像。步骤S22、在DOG金字塔的每层上提取极值点，所谓极值点是指在局部邻域内D (x,y，ko )取值最大的点。提取极值点的过程为，选择DOG金字塔上的任一点，若该点在该层及上、下相邻两层的26个邻域内不是极值点，则将该点去除，否则将该点作为极值点。图6是提取极值点的示意图，如图6所示，DOG金字塔中的第k层图像的标“ X ”的点与第k-1层、第k+1层的26个标““的点进行比较，若标“ X ”的点对应的D (x，y，k 0)是这26个邻域点中的最大值，则标“ X ”的点为极值点。步骤S23、对于所提取的极值点，去除局部曲率非常不对称的极值点。在该实施例中，计算差分图像D的局部Hessian矩阵H，去除满足下面条件的极值点tr (H)2/det (H) >10，其中det(H)表示矩阵H的行列式值，tr(H)表示矩阵H的迹；步骤S24、计算SIFT特征的亚像素级别的空间位置、尺度，所述SIFT特征是指保留下来的极值点。假设SIFT特征X= (X, y, 0)，其中乂，7, o分别为极值点X的x、y方向坐标及尺度参数，此时三个方向的坐标都为正数。根据高斯差分图像D(x，y，o)泰勒展开公式
D(X) = + +AH十算出以SIFT特征X为原点相对于X在行、列及层数上的
OA 2 CA -
Cr I) 1 ol)Pi)! I)
偏移量，即i =其中D表示泰勒展开公式的一次项，D、夂和可由SIFT特征
Sr Bx，dX dX 2
X及其邻域的像素根据差分法计算。然后按照如下规则对极值点进行亚像素级插值如果
三个方向上的偏移量都小于0.5像素，则该点就是极值点X +i为所求的亚像素级极值点坐
标；如果某一方向上的偏移量大于等于0. 5像素,如假设X方向的偏移量大于0. 5像素,将
偏移量四舍五入到一整数值a，将a与X相加得到新的SIFT特征X2 = (x+a, y，s)，然后将
极值点X2按上述步骤进行操作，直到三个方向上的偏移量都小于0. 5。步骤S25、确定所述SIFT特征主方向，所谓主方向是指与以SIFT特征为中心的邻域内的梯度方向直方图的峰值对应的梯度方向。在该实施例中，在 以SIFT特征X = (x，y，o)为中心、以1.5 o为半径的邻域窗口内采样，计算高斯平滑图像L(x，y，o)在上述邻域窗口内的每一像素的梯度方向
J ( 少，.+ I ) _ J ( . I)
& = tan^iLix+ ly)-Ux^
直方图。梯度方向直方图是一种关于梯度方向0的统计图，其范围是0 360度，其中每10度一个柱，总共36个柱。该梯度方向直方图的峰值代表了该极值点处邻域梯度的主方向，即作为该SIFT特征的方向。图7是本发明的该实施例的梯度直方图的一个示例图。在该图中显示了采用7个柱为SIFT特征确定主方向的示例。在该梯度方向直方图中，当存在另一个相当于主峰值80%能量的峰值时，则将这个方向认为是该SIFT特征的辅方向。一个SIFT特征可能会被指定具有多个方向(一个主方向，一个以上辅方向)。步骤S26、为SIFT特征构造描述向量，所谓描述向量是指用于刻画该SIFT特征周围图像块统计特征的、由梯度方向直方图构成的向量。在该实施例中，首先将SIFT特征周围图像块的坐标轴旋转为所述SIFT特征的方向，以确保旋转不变性；然后将SIFT特征周围16X16像素的窗口内分成4X4像素的小块，计算每个小块的8个方向的梯度方向直方图，将每个小块的梯度方向直方图连接起来形成128维的特征向量；最后，将其归一化到单位长度。图8是本发明的为SIFT特征构造描述向量的一个不意图。 本发明的步骤S2还涉及到在多个尺度上提取角点特征和面点特征。对图像I，在Crfr 多个尺度上提取角点特征和面点特征。k = 0，…，N，N控制尺度的个数，本实例N = 4。在每个尺度上，角点特征和面点特征提取方法如下步骤S21’、对于图像的每个像素点X计算互相关矩阵M(X)，互相关矩阵是用来描述该像素点与其邻域像素的关系。计算方式如下式，其中VJ表示图像I在水平、竖直方向的梯度图像，《0表示高斯权重函数，其标准差为O，像素点X邻域N(X)的尺寸为30。步骤S22’、根据互相关矩阵的特征值选取角点特征和面点特征，得到角点特征集合和面点特征集合。首先对互相关矩阵M(X)进行特征值分解，即M(X) =E ^lj2Ai(X) Ti(X) Ti(X)1,X1(X)(入2 (X)，X1(X)和X2(X)表示X处M(X)的特征值，r J (x)表示与Xi(X)对应的特征向量，Fi(X)T表示rjx)的转置。然后，根据特征值来选取角点特征和面点特征。若X1(X)7^2(X) >0.1且m (X) = A j (x) + A 2 (X) > tm,则 X 为候选角点特征；若入！ (X) / X 2 (X)彡 0. I 且 m (X)=A1UhA2(X) <tm，x为候选面点特征。对每个像素点进行上述操作，得到角点特征集合C和面点特征集合F。tm是一阈值，表示特征的强度。tm应取一正数，tm越大得到的候选角点特征和面点特征越少，不同的tm取值对后续的特征匹配影响不大，在该实施例中，tffl = I。步骤S23’、去除角点特征集合和面点特征集合中的噪声点。对于角点特征集合C或面点特征集合F中的某特征点X，若m(x) < Um+0. 5om，将X从C或F中去除。其中，Um和O m分别为半径r = 30的邻域内的m(x)的均值和方差。步骤S24’、对于角点特征和面点特征进行非极值点抑制。对于角点特征集合C中的角点特征x，沿特征值X1(X)和A2(X)对应的特征向量rjx)和r2(x)的方向进行非极大值点抑制；对于面点特征集合F中的某面点特征X，沿X2(X)对应的特征向量的方向r2(x)进行非极大值点抑制。具体方法如下对以X为中心、沿rjx)方向的半径r = 30的邻域内，若存在X'满足!!!(X' )> m(x),则X不是好的特征,将其从候选特征集合中去掉。
步骤S25’、根据特征强度的大小来选择被动角点特征和被动面点特征。在该实施例中,对角点特征集合C或面点特征集合F中的特征按特征强度m (X)从大到小排序，选择前n个角点特征和面点特征作为被动匹配角点特征和面点特征的集合。n根据图像大小确定，对于1024X1024像素的图像，n—般取1000。步骤S26’、根据特征强度的大小来选择主动角点特征和主动面点特征。改变特征强度阈值，令tm = 2tm，r = 2r，重复步骤S22’_S25’。对角点特征集合C或面点特征集合F中的特征点按特征强度m (X)从大到小排序,选择前m个角点特征和m个面点特征作为主动匹配角点特征、面点特征集合。m根据图像大小确定，对于1024X1024像素的图像，m—般取500。同一图像的主动匹配角(面)点特征是被动匹配角(面)点特征的子集，在匹配时，参考图像的主动匹配角(面)点特征与待配准图像的被动匹配角(面)点特征(而不是待配准图像的主动匹配角(面)点特征)进行匹配，这样做的好处是可以得到更稠密的角(面)点特征匹配点对，从而有利用提高配准精度。 按照上述步骤在每个尺度O k上提取被动匹配角点特征、面点特征和主动匹配角点特征、面点特征，将各尺度、各类型特征点合并在一起，组成被动匹配角点特征集合Cm、被动匹配面点特征集合Fm、主动匹配角点特征集合Cd、主动匹配面点特征集合Fd。每一个角点特征Ci由三维向量组成，即Ci = (Xi, yi； Si ) , Si等于它所在的尺度参数O k。每一个面点特征fi由四维向量组成，即fi = (xi； yi； Si, Hi), Hi表示方向，Hi对应步骤S26’中在(Xi, Y1)处的 r2(x)。图9是本发明的方法的步骤S3的流程图。该步骤为粗配准步骤。如图9所示，本发明中的步骤S3是按如下流程实现的选择置信度高且未尝试过的SIFT特征对，由该SIFT特征对及被动匹配角点特征集合Cm、被动匹配面点特征集合Fm、主动匹配角点特征集合Cd、主动匹配面点特征集合Fd通过⑶BICP方法进行匹配，若该次匹配成功，粗配准终止；否则，选择下一个置信度高的SIFT特征对通过⑶BICP方法进行匹配。图10是本发明的方法的步骤S3中的⑶BICP方法的流程图。如图10所示，⑶BICP方法包括如下步骤步骤S31、根据选择的SIFT特征对确定初始变换模型(相似变换)，SIFT特征对所在的邻域作为初始匹配区域；设当前选择的SIFT 特征对为 A = (X1, y1； S1, 0 J 和 B = (x2, y2, S2, 0 2),其中 Xi和71是坐标值，^和e2是对应的尺度参数和主方向参数。该步骤为粗配准，变换类型为相
V s cos O s sin 0 t ^
似变换，SIFT特征对之间应满足2 = _sin^ v.cos0 / s、0、tx、ty为所求的变
L少2」L‘I
换参数。根据A和B的尺度参数和主方向参数的对应关系,可求得^ — Y 0 = 0 2- 0 1； Xt
=x2-s cos 0 X1-S sin 0 Y1, yt = y2+s sin 0 x「s cos 9 yp 以(Xi, Yi)为中心、以h+3 Si为半径的区域分别作为A和B的初始匹配区域。步骤S32、根据当前迭代步骤的变换模型(注求变换模型是一个迭代的过程，S4的变换模型是最终的变换模型，当前迭代步骤的变换模型是一个中间结果)和当前匹配区域，分别从前向和后向利用ICP (Iterative Closest Point)方法匹配当前匹配区域内的角点特征和面点特征；前向匹配以可见光图像为参考图像，为红外图像上的每一个角点特征和面点特征寻找其在红外图像上对应的角点特征和面点特征。假设当前变换模型类型为&M (k表示变换类型，k = 0表示相似变换、k = I表示仿射变换、k = 2表示投影变换)，对于红外图像的任一角点特征Ci = (Xi, Yi, Si),它对应可见光图像的坐标记为7在可见光图像上距离/Xe,jf)最近的三个角点特征为么<)，i = 1,2,3.然后计算<中每个元素与Ci的匹配度” =min(</^，^/<)。将与最大的< (丨=1，2，3)对应的<作为Ci的匹配角点特征，记为e”类似的，为红外图像的每一面点特征& = (Xi, Yi, Si, Jli)寻找可见光图像对应的匹配面点特征心面点特征的匹配度计算公式为< =min(04/0|/7fif|，
g!(丨=1，2,3)为可见光图像上距离11/^广)最近的三个面点特征。后向匹配以红外图像为参考图像，为可见光图像上的每一个角点特征和面点特征寻找其在红外图像上对应的角点特征和面点特征，操作过程与前向匹配类似。步骤S33、根据当前匹配的角点特征和面点特征利用迭代重加权最小二乘法重新估计前向和后向变换参数；通过ICP和参数估计的交替迭代，确定最优变换参数和最优的匹配；当前变换模型与当前匹配的角点特征和面点特征紧密相关。为此，在给定变换类型的前提下采用迭代重加权最小二乘法不断调整变换参数及匹配的角点特征和面点特征集合。为方便叙述，在步骤S33中，当变换类型给定后将变换模型0己为0M。迭代重加权
最小二乘法的目标函数为
权利要求
1.一种可见光图像与红外图像的配准方法，其特征在于，包括如下步骤 步骤SI、对可见光图像和红外图像进行多尺度分解，得到低分辨率的可见光图像和红外图像； 步骤S2、在所述低分辨率图像对上提取SIFT特征、多尺度角点特征和面点特征，匹配SIFT特征并去除外点，利用匹配的SIFT特征对得到初始变换模型，所谓变换模型是指红外图像和高光谱图像之间的一种几何变换关系； 步骤S3、在原始图像对上利用初始变换模型提供的几何约束，进行基于图像块对的多尺度角点特征、面点特征提取，根据初始变换模型和多尺度角点特征、面点特征确定更精确的变换模型； 步骤S4、根据所述变换模型对红外图像进行变换，得到变换后的红外图像。
2.如权利要求I所述的可见光图像与红外图像的配准方法，其特征在于，在步骤SI中，通过降采样实现多尺度分解。
3.如权利要求I所述的可见光图像与红外图像的配准方法，其特征在于，在步骤S2中，采用最近邻比值法来匹配SIFT特征，并去除外点，外点指不满足上述变换模型的SIFT特征对。
4.如权利要求I所述的可见光图像与红外图像的配准方法，其特征在于，在步骤S2中，利用GDBICP法得到初始变换模型。
5.如权利要求I所述的可见光图像与红外图像的配准方法，其特征在于，所述步骤S2中的SIFT特征提取步骤包括 步骤S21、构建图像的DOG金字塔； 步骤S22、在DOG金字塔的每层上提取极值点，所谓极值点是指在局部邻域内取值最大的点； 步骤S23、对于所提取的极值点，去除局部曲率不对称的极值点； 步骤S24、步骤S24、计算SIFT特征的亚像素级别的空间位置、尺度，所述SIFT特征是指保留下来的极值点； 步骤S25、确定所述SIFT特征主方向，所谓主方向是指与以SIFT特征为中心的邻域内的梯度方向直方图的峰值对应的梯度方向； 步骤S26、为SIFT特征构造描述向量，所谓描述向量是指用于刻画该SIFT特征周围图像块统计特征的、由梯度方向直方图构成的向量。
6.如权利要求I所述的可见光图像与红外图像的配准方法，其特征在于，所述步骤S2中的中提取角点特征和面点特征的步骤包括 步骤S21’、对于图像的每个像素点计算互相关矩阵，互相关矩阵是用来描述该像素点与其邻域像素的关系； 步骤S22’、根据互相关矩阵的特征值选取角点特征和面点特征，得到角点特征集合和面点特征集合； 步骤S23’、去除角点特征集合和面点特征集合中的噪声点； 步骤S24’、对于角点特征和面点特征进行非极值点抑制； 步骤S25’、根据特征强度的大小来选择被动角点特征和被动面点特征； 步骤S26’、根据特征强度的大小来选择主动角点特征和主动面点特征。
7.如权利要求I所述的可见光图像与红外图像的配准方法，其特征在于，在所述步骤S3中，采用最近邻比值法来匹配SIFT特征，并去除外点，外点指不满足上述变换模型的SIFT特征对。
8.如权利要求I所述的可见光图像与红外图像的配准方法，其特征在于，所述步骤S3包括 步骤S31、根据选择的SIFT特征对确定初始变换模型，SIFT特征对所在的邻域作为初始匹配区域； 步骤S32、根据选择的SIFT特征对确定初始变换模型，SIFT特征对所在的邻域作为初始匹配区域和当前匹配区域，分别从前向和后向利用ICP方法匹配当前匹配区域内的角点特征和面点特征； 步骤S33、根据当前匹配的角点特征和面点特征利用迭代重加权最小二乘法重新估计前向和后向变换参数；通过ICP和参数估计的交替迭代，确定最优变换参数和最优的匹配； 步骤S34、利用当前最优匹配的角点特征和面点特征集合,根据Akaike Information准则进行模型选择； 步骤S35、根据当前变换模型的不确定性进行匹配区域增长，区域增长的速度与当前匹配区域的前向匹配和后向匹配的角点特征和面点特征的匹配误差成反比。
9.如权利要求I所述的可见光图像与红外图像的配准方法，其特征在于，所述步骤S4包括 步骤S41、在原始可见光图像和红外图像上利用步骤S3提供的初始变换模型P进行基于图像块对的多尺度角点特征、面点特征特征提取； 步骤S42、根据初始变换模型和多尺度角点特征、面点特征集合利用迭代重加权最小二乘法选择变换模型并求取变换参数。
10.如权利要求I所述的可见光图像与红外图像的配准方法，其特征在于，所述步骤S5包括根据步骤S4求得的变换模型，对红外图像进行双三次样条插值得到变换后的红外图像。
全文摘要
本发明公开了一种可见光图像与红外图像的配准方法包括对可见光图像和红外图像进行多尺度分解，得到低分辨率的可见光图像和红外图像；在低分辨率图像对上提取SIFT特征、多尺度角点特征和面点特征，匹配SIFT特征并去除外点，利用匹配的SIFT特征对得到初始变换模型，所谓变换模型是指红外图像和高光谱图像之间的一种几何变换关系；在原始图像对上利用初始变换模型提供的几何约束，进行基于图像块对的多尺度角点特征、面点特征提取，根据初始变换模型和多尺度角点特征、面点特征确定更精确的变换模型；根据变换模型对红外图像进行变换，得到变换后的红外图像。
文档编号G06T7/00GK102800097SQ20121025140
公开日2012年11月28日 申请日期2012年7月19日 优先权日2012年7月19日
发明者霍春雷, 郭建恩, 潘春洪, 樊彬, 杜鹃, 常民 申请人:中国科学院自动化研究所, 北京市遥感信息研究所