毕业论文(图像匹配)

毕业设计（论文）PAGE35毕业论文(图像匹配)摘要数字图像配准是指将从同一场景拍摄的具有重叠区域的图像通过特征匹配方法，找出图像之间的对应关系。目前，图像配准技术广泛应用于医学、生物、信息处理和其它领域，它已成为图像处理应用中不可或缺的技术。本文主要主要介绍了灰度相关的配准方法，灰度相关的图像配准算法是图像配准算法中比较经典的一种，很多配准技术都以它为基础进行延伸和扩展。它是从待拼接图像的灰度值出发，对待配准图像中一块区域与参考图像中的相同尺寸的区域使用最小二乘法或者其它数学方法计算其灰度值的差异，对此差异比较后来判断待拼接图像重叠区域的相似程度，由此得到待拼接图像重叠区域的范围和位置，从而实现图像配准。对位移量比较大的图像，可以先校正图像的旋转，然后建立两幅图像之间的映射关系。当以两块区域像素点灰度值的差别作为判别标准时，最简单的一种方法是直接把各点灰度的差值累计起来。另一种方法是计算两块区域的对应像素点灰度值的相关系数，相关系数越大，则两块图像的匹配程度越高。该方法的配准效果要更好，配准成功率有所提高。基于灰度相关的配准方法的优点是原理简单，直观性好，计算量较小；缺点是只适用于具有水平或者垂直位移方向上的两幅图像的配准，根据所选择的模板的不同，基于灰度相关的配准方法可以分为：线模板匹配法、比值模板匹配法和块模板匹配法。本文分别对3种匹配方法的算法原理和具体实现进行深入研究。通过研究发现，线模板匹配法精确直观，但是它只能处理简单的平移变换下的图像配准，容易受光照的影响，不能实现图像旋转和缩放情况下的配准。比值模板匹配法算法清晰简单，容易理解，实现起来比较方便。在匹配计算的时候，计算量小，速度快；但是，不能处理旋转缩放情况下的配准。块模板匹配法能准确找到匹配点，能在小的旋转和缩放变形的情况下实现配准；但计算量较大。关键词图像配准，灰度相关，线模板匹配法，比值模板匹配，块模板匹配AbstractDigitalimageregistrationisthesamesceneshotfromtheoverlappingregionoftheimagebyafeaturematchingmethodtofindthecorrespondingrelationshipbetweentheimages.Currently,imageregistrationtechniquesarewidelyusedinmedicine,biology,informationprocessingandotherarea.Ithasbecomeanintegralpartofimageprocessingapplicationsoftechnology.Thispaperintroducesregistrationmethodbasedgray-correlation.Themothodismoreclassicone,manyregistrationtechniquesarewereextendedandexpandedbasedonit.Itisamosaicimagefromthegrayvaluetobeproceedingtowardsregistrationwiththereferenceimageinaregionofthesamesizeintheimageareausingtheleastsquaresorothermathematicalmethodstocalculatethegrayvaluedifferences,thiscomparisonwastodeterminedifferencesoverlaptheimagemosaictobethedegreeofsimilarity,theresultantmosaicimagetobetheextentandlocationofoverlap,inordertoachieveimageregistration.Largerthantheimageofthedisplacementcanbecorrectedimagerotation,thenestablishthemappingbetweentwoimages.Whentworegionaldifferencesingrayvalueofpixelsasthecriterion,themostsimplewayistopointdirectlytothecumulativedifferencebetweenthegrayup.Anothermethodistocalculatethetworegionsofgrayvaluepixelcorrespondingtothecorrelationcoefficient,correlationcoefficientislarger,thehigherthedegreeofmatchingtwoimages.Effectofthemethodbybetterregistration,registrationhasincreasedthesuccessrate.Basedongraycorrelationregistrationmethodhastheadvantageofsimpleinprinciple,straightforward,good,lesscomputation;disadvantageisonlyapplicabletothedirectionofahorizontalorverticaldisplacementoftwoimagesontheregistration,accordingtothedifferentchoiceoftemplatesbasedongraycorrelationregistrationmethodcanbedividedinto:linetemplatematching,templatematchingmethodandtheblockratiotemplatematchingmethod.Thispaperonthelinetemplatematching,templatematchingmethodandtheblockratiotemplatematchingalgorithmprincipleandconcreteachieve-depthstudy.Ithasfoundthatlinetemplatematchingmethodhasprecisevisual,butitcanhandleonlysimpletranslationtransformationoftheimageregistrationisvolunerablecannotbeachievedincaseofregistration.Templatematchingalgorithmistheratioofclearandsimple,moreconvenienttoimplement.Whenthematchingcalculation,computationspeed;however,cannothandletherotationscalingincaseofregistration.Blocktemplatematchingmethodcanaccuratelyfindthematchingpointsthatcanrotateandzoominthesmalldeformationregistrationunderthesituation.KeywordsImageregistration,intensity-related,linetemplatematchingmethod,theratiooftemplatematching,blockmatching,rotatingaroundtheopticalaxis

目录TOC\o"1-2"\h\z\u第1章 绪论 11.1 研究背景及意义 11.2 图像配准方法概述 21.3 研究现状 31.4 研究问题及内容 4第2章 图像配准基本理论 12.1 图像配准的基本介绍 12.2 图像配准的相关概念 12.3 灰度相关的配准方法 5第3章 线匹配法 73.1 线匹配法基本介绍及原理 73.2 线性变换图像匹配模型 73.3 线匹配法具体的算法实现 73.4 实验结果和分析 9第4章 比值匹配法 114.1 比值匹配法基本介绍及原理 114.2 比值匹配法具体的算法实现 124.3 实验结果和分析 12第5章 块匹配法 155.1 块匹配法基本介绍及原理 155.2 块匹配法具体的算法实现 165.3 实验结果和分析 16第6章 总结与展望 206.1 总结 206.2 展望 21参考文献 23致谢 24绪论研究背景及意义数字图像，又称数码图像或数位图像，是二维图像用有限数字数值像素的表示。通常，像素在计算机中保存为二维整数数组的光栅图像，这些值经常用压缩格式进行传输和储存。数字图像可以由许多不同的输入设备和技术生成，例如数码相机、扫描仪、坐标测量机等，也可以从任意的非图像数据合成得到，例如数学函数或者三维几何模型，三维几何模型是计算机图形学的一个主要分支。数字图像处理领域就是研究它们的变换算法。数字图像处理(DigitalImageProcessing)是通过计算机对图像进行去除噪声、增强、复原、分割、提取特征等处理的方法和技术。数字图像处理的产生和迅速发展主要受三个因素的影响:一是计算机的发展;二是数学的发展(特别是离散数学理论的创立和完善);三是广泛的农牧业、林业、环境、军事、工业和医学等方面的应用需求的增长。图像配准(Imageregistration)就是将不同时间、不同传感器(成像设备)或不同条件下(天候、照度、摄像位置和角度等)获取的两幅或多幅图像进行匹配、叠加的过程，它已经被广泛地应用于遥感数据分析、计算机视觉、图像处理等领域。图像配准的方法迄今为止，在国内外的图像处理研究领域，已经报道了相当多的图像配准研究工作，产生了不少图像配准方法。总的来说，各种方法都是面向一定范围的应用领域，也具有各自的特点。比如计算机视觉中的景物匹配和飞行器定位系统中的地图匹配，依据其完成的主要功能而被称为目标检测与定位，根据其所采用的算法称之为图像相关等等。基于灰度信息的图像配准方法一般不需要对图像进行复杂的预先处理，而是利用图像本身具有灰度的一些统计信息来度量图像的相似程度。主要特点是实现简单，但应用范围较窄，不能直接用于校正图像的非线性形变，在最优变换的搜索过程中往往需要巨大的运算量。经过几十年的发展，人们提出了许多基于灰度信息的图像配准方法，大致可以分为三类：互相关法（也称模板匹配法）、序贯相似度检测匹配法、交互信息法。目前主要图像配准方法有基于互信息的配准方法，基于相关性的配准方法和基于梯度的配准方法。其中基于梯度的方法基本很少单独使用，而作为一个辅助性的测度与其它方法相结合起来使用。基于灰度的配准算法是医学图像配准研究的发展方向，也是目前研究的热点之一。基于灰度的配准方法与基于特征的配准方法的区别在于前者没有提取图像特征的步骤，直接对图像中的灰度进行处理。基于灰度的配准方法计算复杂度高、对图像的灰度、旋转、形变以及遮挡都比较敏感。灰度相关的配准方法是从待拼接图像的灰度值出发,图像拼接故而成为灰度相关的配准算法的一个基础。图像拼接(imagemosaic)技术是将一组相互间重叠部分的图像序列进行空间匹配对准,经重采样合成后形成一幅包含各图像序列信息的宽视角场景的、完整的、高清晰的新图像的技术。图像拼接在摄影测量学、计算机视觉、遥感图像处理、医学图像分析、计算机图形学等领域有着广泛的应用价值。灰度相关的图像配准算法在图像处理技术中起着十分关键的作用，它是图像处理技术得以发展的一个重要基础。它推动着图像处理技术在医学、生物、信息处理和其他很多高科技领域内的应用，它已渐渐发展成社会生活中不可分离的一种技术，对于图像处理技术发展及应用具有重要意义。图像配准方法概述配准技术的流程如下：首先对两幅图像进行特征提取得到特征点；通过进行相似性度量找到匹配的特征点对；然后通过匹配的特征点对得到图像空间坐标变换参数：最后由坐标变换参数进行图像配准。而特征提取是配准技术中的关键，准确的特征提取为特征匹配的成功进行提供了保障。因此，寻求具有良好不变性和准确性的特征提取方法，对于匹配精度至关重要。图像配准的方式可以概括为相对配准和绝对配准两种：相对配准是指选择多图像中的一张图像作为参考图像，将其它的相关图像与之配准，其坐标系统是任意的。绝对配准是指先定义一个控制网格，所有的图像相对于这个网格来进行配准，也就是分别完成各分量图像的几何校正来实现坐标系的统一。本文主要研究大幅面多图像的相对配准，因此如何确定多图像之间的配准函数映射关系是图像配准的关键。通常通过一个适当的多项式来拟合两图像之间的平移、旋转和仿射变换，由此将图像配准函数映射关系转化为如何确定多项式的系数，最终转化为如何确定配准控制。目前，根据如何确定的方法和图像配准中利用的图像信息区别可将图像配准方法分为三个主要类别：基于灰度信息法、变换域法和基于特征法，其中基于特征法又可以根据所用的特征属性的不同而细分为若干类别。以下将根据这一分类原则来讨论目前已经报道的各种图像配准方法和原理。研究现状国外从20世纪60年代就开始在图像配准领域进行研究，但直到1980年代才开始引起学者们的关注。到上世纪末，单模图像配准问题已基本解决，但多模图像配准由于涉及模式和领域的复杂性，仍需密切关注。国际上对图像配准技术曾做过调查，其结论是1990年代初技术就明显增加。而国内从1990年代初才开始涉足此领域。与灰度相关的图像配准算法是图像配准算法中比较经典的一种，很多配准技术都以它为基础进行延伸和扩展。朱近，司美玲具体阐述了用局部灰度极值方法进行多光谱图像配准算法的研究过程。针对多光谱遥感图像，提出了一种基于局部灰度极值的配准方法：通过在基准图像和待配准图像中同步寻找含有灰度极值的小区域，再用多项式对极值区域进行曲面拟合，最后，分别计算小区域的极值点作为特征点进行配准。并用真实和模拟多光谱图像进行了试验结果显示该课题提出具有算法简单和配准精度高的特点。这是与灰度相关图像配准算法有关的一个扩展应用。林相波,邱天爽提出一种新的灰度和形状信息相结合的全自动同模态医学图像非刚性配准——分割算法[1]，将欧氏距离表示的形状信息融入基于灰度的配准算法中，构造出新的代价函数。该算法在医学图像多目标分割的应用中，能够较好地完成灰度相近、边缘模糊、间距较小的不同结构的分割，结果表明，该算法优于基于灰度信息的图像配准算法。张密,吴效明阐述图像配准在放疗中应用的关键问题[2]，对基于灰度的3种配准方法的性能做深入研究，包括均方测度、归一化相关测度以及互信息测度。方法：分析各配准要素的算法原理后，基于加以实现，提出使用综合配准误差来评价不同配准算法的性能，并与传统目标配准误差的评价结果作对比。结果：3种测度都能对近模态的图像实施准确的配准，其中互信息测度驱动的配准在配准精度和速度上表现更为稳定，综合配准误差仅为另外两个测度的一半左右。它的结论得出利用综合配准误差得到的评价结果更为客观，互信息测度是放疗中实施配准的较理想测度。在医学、生物、信息处理和其他很多高科技领域内图像配准技术越来越显示出它的重要性，将会越来越受到人们的关注。在现今的科技水平上，图像配准未来的发展趋势也会更好。本文主要主要介绍了灰度相关的配准方法是从待拼接图像[3]的灰度值出发，对待配准图像中一块区域与参考图像中的相同尺寸的区域使用最小二乘法或者其它数学方法计算其灰度值的差异，对此差异比较后来判断待拼接图像重叠区域的相似程度，由此得到待拼接图像重叠区域的范围和位置，从而实现图像拼接。也可以通过变换将图像由时域变换到频域，然后再进行配准。对位移量比较大的图像，可以先校正图像的旋转，然后建立两幅图像之间的映射关系。当以两块区域像素点灰度值的差别作为判别标准时，最简单的一种方法是直接把各点灰度的差值累计起来。这种办法效果不是很好，常常由于亮度、对比度的变化及其它原因导致拼接失败。另一种方法是计算两块区域的对应像素点灰度值的相关系数，相关系数越大，则两块图像的匹配程度越高。该方法的拼接效果要好一些，成功率有所提高。根据所选择的模板的不同，基于灰度相关的配准方法可以分为：线匹配法、比值匹配法和块匹配法。研究问题及内容本文在分析了灰度相关的图像配准算法中的线匹配法、比值匹配法和块匹配法，利用这三种方法分别实现两幅图像在水平垂直位移上的配准，而本课题研究的内容是提出一种基于灰度相关的算法，不仅能实现两幅图在水平和垂直位移的配准，同时也能实现在绕光轴旋转情况下的图像配准。这里提出了一种方法，多尺度模块匹配法。在这三种匹配的环境下，它能实现水平垂直位移上的匹配、缩放以及旋转。同时通过在Matlab编程环境下编程实现相关算法，通过实际图像的配准试验，利用这些结论最终得到精确地配准结果。图像配准基本理论图像配准的基本介绍图像配准的描述图像配准是对取自不同时间，不同传感器或不同视角的同一场景的两幅图像或者多幅图像匹配的过程。图像配准广泛用于多模态图像分析，是医学图像处理的一个重要分支，也是遥感图像处理，目标识别，图像重建，机器人视觉等领域中的关键技术之一，也是图像融合中要预处理的问题，待融合图像之间往往存在偏移、旋转、比例等空间变换关系，图像配准就是将这些图像变换到同一坐标系下，以供融合使用。图像配准的定义对于二维图像配准可定义为两幅图像在空间和灰度上的映射[4]，如果给定尺寸的二维矩阵和代表两幅图像，和分别表示相应位置上的灰度值，则图像间的映射可表示为：,式中表示一个二维空间坐标变换，即，且是一维灰度变换。图像配准的步骤图像配准的基本过程可以分为三个步骤：第一步是为每一个图像信息模式各定义一个坐标系，然后再定义这些参考特征之间的失调或相似函数；第二步是分割出图像的参考特征，再定义这些参数特征之间的失调或相似函数；第三步是应用优化算法，使第二步中失调（相似）函数达到全局最小（最大）值，达到两幅图像的配准。其中参考特征和对应优化算法的选择是配准的核心，也是不同配准算法的差异所在。图像配准的相关概念配准基准通常，图像配准中根据配准基准的特性，可分为基于外部基准的配准和基于内部基准的配准[5]，外部基准是指强加于待配准对象的各种人造标记，这些标记必须在各种配准模式中都清晰可见且可准确检测到。内部基准是指由图像本身得到的位置相对固定且图像特征明晰的各种配准标识。映射变换与配准区域设和表示两幅待匹配的图像，和分别表示两幅图像的密度函数，其中和分别表示在图像和中的像素坐标。图像匹配就是要找到一个把图像映射到图像的变换,使得变换后的图像和具有几何对应性。这种映射变换有刚体变换、仿射变换、投影变换以及曲线变换等。配准时的变换区域根据实际需要又分为局部配准和全局配准。局部变换一般很少直接使用，因为它会破坏图像的局部连续性，且变换的双映射性会影响图像的再采样。从近期关于图像配准方面的文章看，一般刚性和仿射多用于全局变换，而曲线变换多用于局部变换。配准的交互性与优化根据人的参与程序配准又可分为全自动式，交互式和半自动式三种。全自动式中使用者仅需给相应算法提供图像数据以及图像获取的一些可能信息；交互式中使用者必须亲自进行配准，软件仅给目前变换提供一个可视的或数字的感官印象以及初始变换的一个可能参数；半自动式中，交互式有两种方式：一种是使用者须初始化算法，如分割数据，另一种是指导算法，如拒绝或接受配准假设。配准变换的参数可以是直接计算出的，也可以是搜索计算出的。直接计算的最优化方法一般已完全由实例决定，所能研究的工作也仅限于如何使用非常少的信息把此计算方法应用于实际。搜索计算的最优化方法大多都可以用待优化的变换参数的一个标准数学函数来表达配准实例，此函数力图使图像在某一变换时两幅图像可达到最大相似。这些函数通常在单模配准中能简单一些，因为此时图像的相似性更能容易直接定义。我们可以通过使用一个标准的、合适的最优化方法使相似函数达到最优。目前应用比较广泛的方法有Powell的方法、DownhillSimplex方法、Brent的方法以及一系列一维搜索算法、Levenberg-Marquardt最优化算法、Newton-Raphson迭代算法、stochastic搜索算法、梯度下降法(gradientdescentmethods)、遗传算法(geneticmethods)、模拟退火法(simulatedannealing)，粒子群算法(particesworm)，蚁群算法(ant)，几何散列法(geometrichashing)。多分辨率(如金字塔)和多尺度方法可以加速最优化的收敛速度。许多实际应用中使用了不止一种最优化方法，一般是先使用一种粗糙但快速的算法，然后再接着使用一种准确但运算速度慢的算法。图像成像模式与配准方法的分类有单模式和多模式等，单模(monomodality)图像配准是指待配准的两幅图像是同一种成像设备获取的。多模(multimodality)图像配准[6]是指待配准的两幅图像来源于不同的成像设备。基于灰度信息的图像配准方法一般不需要对图像进行复杂的预先处理，而是利用图像本身具有灰度的一些统计信息来度量图像的相似程度。主要特点是实现简单，但应用范围较窄，不能直接用于校正图像的非线性形变，在最优变换的搜索过程中往往需要巨大的运算量。经过几十年的发展，人们提出了许多基于灰度信息的图像配准方法，大致可以分为三类：互相关法（也称模板匹配法）、最大互信息法和基于小波变换的图像配准法。互相关法对于同一物体由于各种图像获取条件的差异或物体自身发生的空间位置的改变而产生的单模图像配准问题常常应用互相关法。在互相关法中互相关值的大小反映了配准的效果。互相关法的思路是找出使各图像之间相关性最大的空间变换参数来实现图像的配准。该方法通过优化两幅图像间的相似性测度来估计空间变换参数（刚体的平移和旋转参数），采用的相似性测度可以是多种多样的，例如相关系数，差值的平方和及相关函数等。其中最经典的相似性测度是归一化的相关系数(correlationcoefficient,CC)，即：,(2.1)式中，为模板图像，，为图像的灰度；为与有相同大小的目标图像，为图像的灰度；和分别为图像和灰度的均方值。由于要对每种变换参数可能的取值都要计算一次相似性测度，互相关法的计算量比较庞大，因此近年来发展了快速搜索算法，例如，用相位相关傅立叶法估算平移和旋转参数；用遗传算法和模拟退火技术减少搜索时间和克服局部极值问题。尤其注意的是互相关法受到不同模态成像特点的影响，例如同一物体在不同的模态图像中表现出纹理和密度的非线性差异，使相关性计算无意义，故互相关性法主要局限于单模图像配准[7]。对于条件不好或曲线不完全闭合的图像配准，Kaneko等提出了一种选择性相关系数法(selectivecorrelationcoefficient，即)，实际上是的扩展，在每次为其计算时间仅仅依靠两幅图像灰度的比较过程，故其代价非常小甚至可以忽略不计。最大互信息法(MaximizationofMutualInformation)互信息是信息论的一个基本概念，是两个随机变量统计相关性的测度。最大互信息法几乎可以用在任何不同模式图像的配准，特别是当其中一个图像的数据部分缺损时也能得到很好的配准效果。当含有相同内容的两幅图像通过几何变换在空间对齐时，它们所包含的灰度值的互信息量最大。因此最大化的互信息量可以作为图像配准准则。基于最大互信息[8]的配准过程实质上是搜索最佳的几何变换参数，使两幅图像的互信息达到最大。该方法采用整幅图像的所有像素共同组成特征空间，再根据特征空间确定一种空间变换，使一幅图像经过该变换后和另一幅图像的互信息最大，最终实现配准。互信息(MutualInformation,MI)是信息论中的一个测度，用于描述两个变量间的统计相关性，或一个变量中包含的另一个变量中的信息的多少，表示两个随机变量之间的依赖程度，一般用熵来表示。熵表达的是一个系数的复杂性和不确定性。变量A的熵定义为：,(2.2),(2.3)将待配准的两幅医学图像定义为浮动图像A和参考图像B，它们是关于图像灰度的两个随机变量集。设它们的边缘概率分布分别为、，联合概率分布，则它们的互信息为：,(2.4)当两幅图像的空间位置达到一致时，其中一幅图像表达另一幅图像的信息，即其互信息应为最大。继互信息测度提出后，学者们对基于Shannon熵的方法做了进一步的研究，相继提出了比互信息更为稳定的，其它一些形式的熵测度，称为归一化的互信息，例如Studholme提出了归一化互信息测度(nrimalizedmutualinformation,):,（2.5）基于小波的图像配准方法近年来图像配准的重要发展之一是采用小波变换进行图像局部特征提取，该方法的关键技术是二维离散小波分解。设在平面内的二维图像,基于二维离散小波变换的图像分解是将该原始图像在某一尺度上分别在方向上进行小波分解，每次分解后的低频部分用表示，高频部分用表示。在某一尺度上，图像可以经过方向和方向的离散小波变换后分解为4个子图像，在方向和方向都是高频子图像,在方向是低频，在方向是高频子图像和在方向是高频方向是低频的子图像。低频子图像给出了原图像的概貌，高频子图像给出了原图像的细貌。对于二维正交小波变换有其快速算法-Mallat算法，它把小波变换的计算问题转化为小波变换后系数的计算问题：在实际操作中，给出尺度层上的离散采样值数据，要计算尺度层上的小波变换系数，即分解算法的问题。基于小波变换的图像配准方法有多分辨率分析的优势，能够提高配准的速度。灰度相关的配准方法基于灰度信息的图像配准方法一般不需要对图像进行复杂的预先处理，而是利用图像本身具有灰度的一些统计信息来度量图像的相似程度。主要特点是实现简单，但应用范围较窄，不能直接用于校正图像的非线性形变，在最优变换的搜索过程中往往需要巨大的运算量。经过几十年的发展，人们提出了许多基于灰度信息的图像配准方法，大致可以分为三类：互相关法（也称模板匹配法）、序贯相似度检测匹配法、交互信息法。互相关法是最基本的基于灰度统计的图像配准[9]方法，通常被用于进行模板匹配和模式识别。它是一种匹配度量，通过计算模板图像和搜索窗口之间的互相关值，来确定匹配的程度，互相关值最大时的搜索窗口位置决定了模板图像在待配准图像中的位置。序贯相似度检测匹配法(SequentialSimilarityDetectionAlgorithms，)是由Barnea等人提出来的。方法的最主要的特点是处理速度快。该方法先选择一个简单的固定门限，若在某点上计算两幅图像残差和的过程中，残差和大于该固定门限，就认为当前点不是匹配点，从而终止当前的残差和的计算，转向别的点去计算残差和，最后认为残差和增长最慢的点就是匹配点。这种方法的基本思想是基于对误差的积累进行分析。所以对于大部分非匹配点来说，只需计算模板中的前几个像素点，而只有匹配点附近的点才需要计算整个模板。这样平均起来每一点的运算次数将远远小于实测图像的点数，从而达到减少整个匹配过程计算量的目的。交互信息法最初是Viola等人于1995年把交互信息引入到图像配准的领域的，它是基于信息理论的交互信息相似性准则。初衷是为了解决多模态医学图像的配准问题。交互信息用来比较两幅图像的统计依赖性。首先将图像的灰度视作具有独立样本的空间均匀随机过程，相关的随机场可以采用高斯—马尔科夫随机场模型建立，用统计特征[10]概率密度函数来描述图像的统计性质。交互信息是两个随机变量和之间统计相关性的量度，或是一个变量包含另一个变量的信息量的量度。交互信息图像配准方法一经提出，有不少基于此类的研究，尤其在医学图像的配准问题上。比如将交互信息和梯度结合起来改善其极值性能的算法、多分辨率图像金字塔法等等。但交互信息是建立在概率密度估计的基础上的，有时需要建立参数化的概率密度模型，它要求的计算量很大，并且要求图像之间有很大的重叠区域，由此函数可能出现病态，且有大量的局部极值。本文接下来将要陈述灰度相关的几种配准方法，大致可分为：线匹配法、比值匹配法和块匹配法。比值匹配法是指将待配准图像的一定间隔的行或列的像素的比值作为模板；块匹配法是指将待配准图像的正方形区域的像素的集合作为模板。比值匹配法是从一幅图像的重叠区域中部分相邻的两列上取出部分像素，然后以它们的比值作为模板，在另一幅图像中搜索最佳匹配，这种算法计算量较小，但精度低；块匹配法是以一幅图像重叠区域中的一块作为模板，在另一幅图像中搜索与此模板最相似的匹配块，这种算法精度较高，但计算量过大。而设计一种基于灰度相关的算法，既能实现水平、垂直位移上的配准，同时也能实现绕光轴旋转这种情况是本文的核心。线匹配法线匹配法基本介绍及原理针对图像的平移、旋转、尺度变换等整体匹配，提出了一种图像线性变换的匹配算法[11]。首先定义待匹配图像匹配点间的线性变换模型，以对应像素灰度差平方和作为图像匹配误差函数，然后借助最小化误差函数确定参数迭代增量，由迭代法求得最佳线性变换参数。为减少计算量与提高收敛速度提出了三种改进策略：增加权函数、图像网格点采样和增加加速运动量。实验显示对于小范围平移、旋转及尺度变换的图像能进行准确快速的整体匹配，而改进策略能有效提高匹配速度。线性变换图像匹配模型令待匹配的两幅二维图像为和，两幅图像中任意一对匹配点的坐标满足线性变换关系。令图像中某一坐标为，它与图像中的坐标对应。上式中上标表示转置运算。坐标与之间存在一个偏移量，定义如下：(3.1)(3.2)(3.3)因此坐标变换可以表示为,即像素与对应。线匹配法具体的算法实现如果给定条件，两幅图像内容整体间存在某种线性变换，则通过求解变换系数能实现配准。这对很多问题是一种合理假设。基于这一思想，为解决图像整体匹配问题，本文针对图像的平移、旋转、尺度变换等整体匹配，提出一种带6个参数的坐标线性变换图像匹配新模型，该模型将匹配误差定义为图像整体误差，通过最小化误差函数确定参数迭代增量，由迭代法求得最佳线性变换参数。并且针对收敛速度提出了增加权函数、图像网格点采样和增加加速动量项三种策略及具体实现方法。（1）误差函数及参数求解算法图像匹配的误差函数定义为：，(3.4)该公式中领域与图像的分辨率一致。这种领域设置与光流模型的领域设置存在显著不同，它对噪声更不敏感。通过最小化误差函数可以求出变换参数，它决定各像素的偏移量为建立迭代求解模型，引入增量，代入式(3.1)得到，新的模型如下：(3.5)(3.6)为求解增量，上述函数需要对求偏导。一种可行方法是将函数一阶泰勒展开：(3.7)上式中点号表示积。将(3.7)代入(3.6)，对求偏导，并令偏导为零，则得到下式：(3.8)上式中表示，它用图像F中坐标X+T处得灰度梯度近似。表示六维列向量。求只需对方程(3.8)乘以的逆矩阵即可。（2）迭代方法求解变换参数求解图像与图像间的匹配，可以用增量迭代法计算变换参数向量，算法如下：①初始化向量为零向量；②按式(3.8)计算参数增量；③更新参数；④如果的二范数小于某一阈值，则表示迭代收敛，程序退出；若大于该阈值则转到步骤②继续迭代；如果迭代次数超过最多限定次数则程序退出，提示不收敛。判断收敛的条件也可设置为最近次迭代的二范数之和小于某一阈值，的典型取值为3。实验结果和分析首先利用配准图像库中的图像，生成具有不同相对位移的子图像；平移配准参数可以在程序中进行调整，图3.1及图3.2为配准参数为时生成一对尺寸为400×400的两幅图像。图3.1Lena参考图像图3.2Lena待配准图像图3.3及图3.4为配准参数为时生成尺寸为288×288的两幅图像。图3.3Lena参考图像图3.4Lena待配准图像图3.5及图3.6为配准参数为生成尺寸为456×456的两幅图像。图3.5Lena参考图像图3.6Lena待配准图像图3.7为采用线模板匹配方法的实验结果。黑色线段表示选取的线模板，连线表示线段两端点为匹配点对。图3.7Lena线匹配效果图由实验可知，线模板匹配是在参考图像中选定一行上取出部分像素的灰度值，在搜索图中搜索最佳匹配。线匹配法的优点是精确直观，缺点是操作复杂，重复步骤较多，实用性不强。它只能处理简单的平移变换下的图像配准，容易受光照的影响，不能实现图像旋转和缩放情况下的配准，故而存在一定的局限性，适用性不强。

比值匹配法比值匹配法基本介绍及原理比值匹配法[12]算法思路是利用图像中两列上的部分像素的比值作为模板，即在参考图像的重叠区域中分别在两列上取出部分像素，用它们的比值作为模板，然后在搜索图中搜索最佳的匹配。匹配的过程是在搜索图中，由左至右依次从间距相同的两列上取出部分像素，并逐一计算其对应像素值比值;然后将这些比值依次与模板进行比较，其最小差值对应的列就是最佳匹配。这样在比较中只利用了一组数据，而这组数据利用了两列像素及其所包含的信息。图4.1为图像模板选取示意图，其中，图1为像素的图像，图1为像素的图像，和可相等，也可不等。图1和图1为左右重叠关系，图1在图1的左边。jj+spanjj+span像素的图像像素的图像图4.1图像模板选取示意图在图4.1的重叠区域选取间隔为span的2列像素（第列和第列），计算其对应像素比值，即为模板。(4.1)式中，，为选定的列。在图4.1中从第1列开始依次取间隔为span的2列，计算其对应像素的比值，即为模板。(4.2)计算a模板与b模板差值，即为c模板，函数表达式为：（4.3）式中，，。为二维数组，对对应的列向量求得到。（4.4）的大小就反映图1所示图像选定像素对应列的差异，的最小值对应的列坐标即为最佳匹配。比值匹配法具体的算法实现该算法的具体实现步骤如下：(1)在参考图像中间隔为个像素的距离上的两列像素中，各取个像素，计算这个像素的比值，将个比值存入数组中，将其作为比较的模板。(2)从搜索图中在同样相隔个像素的距离上的两列，各取出个像素，计算其比值，将个比值存入数组。假定垂直错开距离不超过个像素，多取的个像素则可以解决图像垂直方向上的交错问题。(3)利用参考图像中的比值模板在搜索图中寻找相应的匹配。首先进行垂直方向上的比较，即记录下搜索图中每个比值数组内的最佳匹配。再将每个数组的组内最佳匹配进行比较，即进行水平方向的比较，得到的最小值就认为是全局最佳匹配。此时全局最佳匹配即为图像间在水平方向上的偏移距离，该全局最佳匹配队应的组内最佳匹配即为图像间垂直方向上的偏移距离。本算法的思路是在第1幅图像的重叠区域中分别在2列取出部分像素，用它们的比值作为模板；在第2幅图像中由左至右依次从间距相同2列取出部分像素，并逐一计算其对应像素值比值；然后将这些比值依次与模板进行比较，其最小差值所对应的列就是最佳匹配。实验结果和分析以下是光照效果由弱到强的比值匹配仿真实验结果图。比值匹配和线匹配相比不同之处是通过不同亮度的调节可检测到配准参数，由此可以看出它们的区别。比值匹配是利用两条线段的比值进行匹配。灰度比值匹配法较线匹配法多了一项光照对平移参数的影响，下面我们以平移参数为主，着重研究输入图像存在亮度差异时，仿真图像配准的实验效果。图4.2为待配准图像与参考图像之间存在亮度差异为-20%时的配准效果图。图4.2Lena灰度比值匹配图（亮度差异-20%）图4.3为待配准图像与参考图像之间存在亮度差异为+5%时的配准效果图。图4.3Lena灰度比值匹配图（亮度差异+5%）图4.4为待配准图像与参考图像之间存在亮度差异为+20%时的配准效果图。图4.4Lena灰度比值匹配图（亮度差异+20%）由图4.2-图4.4可以看出灰度比值匹配的优缺点如下：(1)比值匹配法的优点:a.算法思路清晰简单，容易理解，实现起来比较方便。b.匹配计算的时候，计算量小，速度快。

(2)比值匹配法的缺点:a.利用图像的特征信息太少。只利用了两条竖直的平行特征线段的像素的信息，没有能够充分利用了图像重叠区域的大部分特征信息。虽然算法提到，在搜索图中由左至右依次从间距相同的两列上取出部分像素，计算其对应像素的比值，然后将这些比值依次与模版进行比较，好像是利用了搜索图中的重叠区域的大部分图像信息，但在参考图像中，只是任意选择了两条特征线，没有充分利用到参考图像的重叠区域的特征信息。b.对图片的采集提出了较高的要求。此算法对照片先进行垂直方向上的比较，然后再进行水平方向上的比较，这样可以解决上下较小的错开问题。在采集的时候只能使照相机在水平方向上移动。然而，有时候不可避免的照相机镜头会有小角度的旋转，使得拍摄出来的照片有一定的旋转，在这个算法中是无法解决的。而且对重叠区域无明显特征的图像，比较背景是海洋或者天空，这样在选取特征模版的时候存在很大的问题。由于照片中存在大块纹理相同的部分，所以与模版的差别就不大，这样有很多匹配点，很容易造成误匹配。c.不易对两条特征线以及特征线之间的距离进行确定。算法中在参考图像的重叠区域中取出两列像素上的部分像素，并没有给出选择的限制。然而在利用拼接算法实现自动拼接的时候，如果选取的特征线不是很恰当，那么这样的特征线算出来的模版就失去了作为模版的意义。同时，在确定特征线间距时，选的过大，则不能充分利用重叠区域的图像信息。选择的过小，则计算量太大。

块匹配法块匹配法基本介绍及原理块匹配法[13]是指将待配准图像的正方形区域的像素的集合作为模板。块匹配法是以一幅图像重叠区域中的一块作为模板，在另一幅图像中搜索与此模板最相似的匹配块，这种算法精度较高，但计算量过大。灰度信息包含了最大量的图像信息，而边缘信息则反映了图像内在的性质，它不易受外界光照条件的影响而产生剧烈的变化。因此相比灰度信息，边缘信息，使得其抗灰度和几何畸变能力强，采用边缘信息构建模板在图像拼接[17]中可获得更加可靠的稳定性。因此我们以图像的边缘为特征寻找基准块。应用块匹配算法，首先要有搜索最佳匹配的标准，这里称之为价值函数：均方误差(MSE)、绝对误差和(SAD)、平均绝对误差(MAD)、方差和(SSE)、绝对变化误差和(SATD)都可以作为价值函数。其中常用的是均方误差(MSE)和平均绝对误差(MAD)，如方程(1)和(2)。其中为块边长像素数（为方便搜索块一般划分为正方形），和分别为当前宏块和参考宏块相应像素的灰度。，（5.1），（5.2）具体步骤首先要将当前和参考帧图像分块，选取的参考帧为当前帧的前一帧，一般宏块为16×16像素大小，然后确定搜索范围大小，一般为以宏块为中心的30×30像素区域，最后在参考图像相应搜索区域中寻找价值函数最小的宏块，运动矢量从参考宏块位置指向当前宏块位置，如图5.1所示。图5.1块匹配原理块匹配法具体的算法实现首先需要对图像进行边缘检测，这里我们使用算子获得原图像的边缘图，然后我们以图像中每个像素点的邻域边缘量来定义此位置的边缘信息大小:,(5.3)其中，为原图像所对应的二值边缘图，A为在点处所取邻域的1/2边长。以值为依据，找出最大值所对应的点作为基准点，然后以此点为中心点选择大小适中的块就可以找到基准特征块[14]。这种算法我们可以称之为传统的块匹配算法。具体步骤首先要将当前和参考帧图像分块，选取的参考帧为当前帧的前一帧，一般宏块为16×16像素大小，然后确定搜索范围大小，一般为以宏块为中心的30×30像素区域，最后在参考图像相应搜索区域中寻找价值函数最小的宏块，运动矢量从参考宏块位置指向当前宏块位置。实验结果和分析通过三、四章的阐述，我们已经得出线匹配法只能处理平移操作下的图像配准，而灰度比值法在线匹配法的基础上多了一个光照变换对平移参数的影响。这两种算法都存在自身的不足，那就是不能实现图像在旋转、缩放情况下的配准。由此我们研究一种更加完善的匹配法——块匹配法，它不仅能实现水平、垂直方向上位移变化的图像配准，同时也能在小的图像旋转缩放变形的情况下的实现配准。图5.2是块匹配的效果图，接下来我们将具体描述平移、旋转和缩放情况下的图像配准。首先利用对图像做小的旋转、缩放，再进行配准，旋转情况下研究旋转对求得平移参数的影响，缩放情况下研究配准参数的误差情况。图5.2Lena块匹配效果图构造平移、旋转、缩放图像利用Lena彩色图像的灰度处理，将生成的灰度图像做小的旋转，缩放，在Matlab上实现图像的微小变化，运行结果如下：图5.2选用Lena原图图5.3Lena灰度处理后图5.4平移（2020）的Lena图像图5.5旋转3°的Lena图像图5.6缩放2%Lena图像块匹配算法下的平移、旋转、缩放配准仿真实验分别对经过平移、旋转、缩放后的图像进行配准，结果如图5.7-图5.9所示。图5.7采用块模板匹配方法的Lena图像平移配准图5.8旋转3°的Lena图像配准图5.9放大2%的Lena图像配准由实验结果可知，块匹配法较以上两种匹配法的优点是它能够实现两幅图像不仅在平移，而且在旋转。缩放的情况下也能实现图像的配准。块匹配法的优点是能准确找到匹配点，打破灰度配准不能实现旋转和缩放情况的结论，缺点是计算量大。在今后的工作中，我们将继续研究，实现该算法的快速匹配。总结与展望总结本文首先介绍了图像配准的研究背景和意义，然后具体阐述了图像配准的方法、特征及分类。接着引出本课题灰度相关的图像配准算法，通过简要的介绍，开始详细的阐述灰度相关的几种方法，分别从线匹配法、比值匹配法和块匹配法三种灰度相关的配准算法进行分析比较，最终得出研究结果，能够设计出一种算法实现平移、旋转及缩放。第1章绪论部分对本课题研究的基础，原理和基本概念做了一下相关的介绍。从现阶段科技发展的需要与社会上相关的各方面的应用，以及图像配准的现状，显示了它的重要意义。第2章就图像配准的基础理论出发，对它进行了基本的介绍，分别从它的相关特征，方法分类等等具体的说明，让我们对图像配准有了进一步的理解。接下来就引出了图像配准算法中比较经典的一种，灰度相关的配准算法，并进行了简要的阐述，引出了灰度相关的几种匹配法，线匹配法、比值匹配法以及块匹配法。第3章具体介绍了线匹配法的基本原理，具体的算法实现以及它的优缺点。此方面源于最早的一种匹配算法，简单直观，精确，当然它也存在自身的不足，如果图片经过旋转或缩放后，该算法就无能为力了。第4章本文就具体的介绍了比值匹配法的基本原理，具体的算法实现，以及它的优缺点。此方法相对来说具有较强的适用性但还有需要改进之处，例如在处理灰度变化较小的图像时，该方法的配准成功率不高，还需要针对这种情况对算法进行改进。比值匹配法与线匹配相比多了一项光照的影响，当输入图像存在亮度差异时，可以看出配准效果的不同。比值匹配法具有较强的实用性，再经过一些改进，可以得到广泛应用。第5章介绍了块匹配法，从它的基本原理出发，阐述此算法的具体实现，分析块匹配法的优缺点。得出了块匹配法较前两种算法更加全面，不仅能实现平移的图像配准，在图像经过旋转和缩放后，也能实现配准，只不过配准在缩放或旋转幅度较大时，误差很大。最后一章是本文研究的重点，如何实现两幅图像不仅在水平、垂直位移上的配准而且能实现空间上绕光轴旋转匹配的这种配准，针对这种想法，可以通过提取对应的灰度特征点，进行特征匹配。找到两幅图像的相对位移，利用图像中的线或者方形领域的灰度比较，来判断位移。我们通过线匹配法、比值匹配法、块匹配法这三种算法具体的程序分别实现两幅图像水平、垂直位移上的配准以及旋转缩放的情况。并且通过大量的实验，结合Matlab程序，以及图像的处理，最终得出结果并作出分析，实验结果表明，我们是可以设计出一种基于灰度相关的算法，来实现两幅图像不仅在水平、垂直位移上，同时在绕光轴旋转这种情况下的配准。也相信随着科技水平的提高，对算法进行更深入的研究，在今后的发展中，会得到较大的应用。展望本文通过灰度相关的基本匹配算法的阐述，提出本文的研究内容，通过线匹配法、比值匹配法、块匹配法三种算法实现两幅图像在水平垂直位移上以及绕光轴旋转的配准。而能实现绕光轴旋转这种配准则需要运用到块匹配法，当然，这里列举的诸多实例阐述了这种情形，实现两幅图像绕光轴旋转的这种情况。当然，各种算法都存在它自身的不足，线匹配法操作复杂，重复步骤较多，实用性不强。它只能处理简单的平移变换下的图像配准，容易受光照的影响，不能实现图像旋转和缩放情况下的配准，故而存在一定的局限性，适用性不强。比值匹配法利用图像的特征信息太少。只利用了两条竖直的平行特征线段的像素的信息，没有能够充分利用了图像重叠区域的大部分特征信息。虽然算法提到，在搜索图中由左至右依次从间距相同的两列上取出部分像素，计算其对应像素的比值，然后将这些比值依次与模版进行比较，好像是利用了搜索图中的重叠区域的大部分图像信息，但在参考图像中，只是任意选择了两条特征线，没有充分利用到参考图像的重叠区域的特征信息。比值匹配对图片的采集提出了较高的要求。此算法对照片先进行垂直方向上的比较，然后再进行水平方向上的比较，这样可以解决上下较小的错开问题。在采集的时候只能使照相机在水平方向上移动。然而，有时候不可避免的照相机镜头会有小角度的旋转，使得拍摄出来的照片有一定的旋转，在这个算法中是无法解决的。而且对重叠区域无明显特征的图像，比较背景是海洋或者天空，这样在选取特征模版的时候存在很大的问题。由于照片中存在大块纹理相同的部分，所以与模版的差别就不大，这样有很多匹配点，很容易造成误匹配。另外，比值匹配不易对两条特征线以及特征线之间的距离进行确定。算法中在参考图像的重叠区域中取出两列像素上的部分像素，并没有给出选择的限制。然而在利用拼接算法实现自动拼接的时候，如果选取的特征线不是很恰当，那么这样的特征线算出来的模版就失去了作为模版的意义。同时，在确定特征线间距时，选的过大，则不能充分利用重叠区域的图像信息。选择的过小，则计算量太大。由实验结果可知，块匹配法较以上两种匹配法的优点是它能够实现两幅图像不仅在平移，而且在旋转。缩放的情况下也能实现图像的配准。块匹配法的优点是能准确找到匹配点，打破灰度配准不能实现旋转和缩放情况的结论，虽然如此它仍然存在着自身的不足，那就是计算量大。在今后的研究中，我们需要克服算法的复杂度，争取研究出一种快速匹配法，能够减少计算量，达到灰度相关的配准的精确性和稳定度。使得配准效果更加快速精确，尽可能减小误差。参考文献[1]Brown,L.G.,1992.Asurveyofimageregistrationtechniques.ACMComput.Surveys24(4),325–376.[2]Hsieh,J.,Liao,H.M.,Fan,K.,Ko,M.,1996.Afastalgorithmforimageregistrationwithoutpredeterminingcorrespondences.In:Proc.IAPRInternat.Conf.onPatternRecognition,