配准综述

图像配准是对取自不同时间，不同传感器或者不同视角的同一区域的两幅或者多幅图像进行匹配叠加的过程。随着技术手段的不断发展，人们对图像配准的要求也越来越离。包括配准方法的鲁棒性、算法的难易程度、算法的自主性等都是考量算法的指标。图像配准的方法中，有人提出过被广泛认同的四个巧骤：

a．特征空间的选择

特征空间是指将运用到配准中元素的集合。特征空间包括很多方面，比如像素值，比如点、线，或者平面甚至是曲面。

b．搜索空间

搜索策略是指一系列配准变换操作的集合。搜索空间是建立在几何形变基础上的。而几何形变一般包括全局的和局部的几何形变。典型的全局变形包括平移，缩放，扭曲等变形以及它们的组合。而局部形变，对变换参数要求比较苛刻，因为某一套变换参数只能作用在局部形变区域，而其余区域需采用插值技术补充。然而，另外的变形区域则需要另外的变换参数去表述。在一般的处理中，我们将变换模型看成是一种先验知识。不然则需要考量所有的变换模型。

c．搜索策略

搜索策略是实施变换的依据。它的存在是为了找寻变换模型的最优解。常用的搜索策略有松弛模型法、牛顿法和共扼梯度法。

d．相似性度量

相似性度量是对采用的变换模型的评价。相似性度量主要是采用参考图像和配准图像之间的共有的特征之间的差异比较的方法来评价配准的情况。这其中的差异可以包括点位的误差，或者是灰度值相关的差异，还可以是在变换域之间的差异值等等。最常见的相似性度量是检测特征的欧氏距离。

上述的几个步骤构成了图像配准的基本框架。许许多多的算法充斥其中产生了种类繁多的配准方法。下面主要介绍下几种典型的配准算法。目前典型的配准算法中大致可分成基于灰度、基于特征和基于模型的方法。

基于灰度的配准算法

基于灰度的算法中，重点是对图像配准的过程上，而不是特征集的选择

上。很多时候基于灰度的方法都会采用矩形窗口在图像上遍历。根据窗口的设计，去选取变换的模型，因此如果图像中有比较复杂的变换，则这种方法的适用性就很小了，而且当图像在灰度上的表现比较平滑，那么误配准的情况就会很大。

经典的灰度区域的特征匹配方法主要包括相关系数法、傅里叶变换法、互信息法、最小二乘法以及核线相关法。

（1）相关系数法

相关系数是标准化的协方差函数，协方差函数除以两信号的方差便可求得相关系数。

其具体操作过程为：在参考影像上的某点，为了找到他的同名点，这个

⨯的区域序列，当然，一般来说n的取值为奇数。并点为中心，周围选择n n

且根据所选取的这个点，大致估计出配准影像上点的位置。并在配准影像上取等大的配准窗口，这个窗口称为相关窗口。再用下文将提到的公式计算相关系数。若相关系数接近１或者-１时就表示二者区域的相关程度非常高。也就是被认为是同名点。当相关系数为１是表示完全正相关，当相关系数为-１时，表示二者完全负相关。

其配准过程可以分为几个步骤：

⨯的区域作为目标a．首先在参考图像选取以目标点为中心，大小为m n

区域T1，并且最好可以确保目标区域的中心处于特征点或者特征地物上。然后确定搜索区域S1，使得S1中必须要完整的包含一个模板T1，其位置确定可以是采用大致估计，或者粗加工处理后的坐标相对误差来确定。

b．移动模板，模板遍历搜索区域，取相关系数最大的位置，其中心位置被认为是同名点。其中在找同名点的过程中，可以设定一个阀值，当所求得的相关系数值，不满足阀值，则可表示这个区域内找到的同名点失效。

c．选取下一个区域，照１，２步骤继续选取。选择足够的同名点作为点集。

d．当点集数量足够的时候便可以进行图像配准了。配准的方法有很多，比如可以选择多项式拟合等。

该方法最大的优点就是简明易懂，但是这个方法的计算量相对比较大。

而且对点集的选择有要求。点集不能选择过于集中，最好是可以均匀的分布整个图幅。虽然这个方法能精确的配准具有平行关系的图像。但是不适用于图像的尺度和旋转变化大的情况。而且该方法直接作用在图像强度上，并没有对图像进行结构分析。所以对图像的强度变化、噪声、明暗变化甚至传感器的类型都很敏感。

（2）互信息法

互信息法是最近研究比较多的方法，尤其是在多模态图像配准领域。互信息法以信息论为基础，用信息熵去衡量两个区域的相似程度。这个方法可以避免特征点的选取，它可以根据像素的灰度值直接计算相似性度量函数。或者说是直接用熵来衡量，熵的形式有多种。其中Shannon 熵的相似性度量是目前运用比较广的。在信息论中，熵表示的是不确定性的量度，互信息是基于熵的一个概念，两幅图像A ,B 的熵，以及A ,B 的联合熵分别定义为：

()()lg ()

a H A p a p a =-∑ (5-1)

()()lg ()

b H B p b p b =-∑ (5-2)

,(,)(,)lg (,)a b

H A B p a b p a b =-∑ (5-3) 其中ａ，ｂ表示图像的灰度值，()p a ，()p b 分别表示像素出现的概率。

而互信息是基于摘的概念，它表征两个随机变量的相关性大小。在图像范畴则可以看成是两幅图像之间相互包含的程度。对于两幅图像互信息的定义为：

(,)()()()()MI X Y H Y H Y X H X H X Y =-=- (5-4)

其中，()(lg(()))H Y Ex P X =-代表随机变量的熵，()P X 为X 的分布概率。互信息的重要思想是虽然同一个景物不同传感器成像后，灰度值上会有差异，但是在灰度值的变化率或者灰度值的分布上却是基本一致的。它认为配准图像和参考图像会在配准位置上表现出分布的方差最小。该表达式可表示为：

a ()()()()

a b B A b B A n n a b PIU N a N b σσμμ=+∑∑ (5-5)