基于傅里叶变换红外热波图像拼接.docVIP

基于傅里叶变换红外热波图像拼接 　　摘要：红外热波无损检测中，红外热波空间序列图像间的变化主要是位移和小幅度旋转。针对该特点，在可能的旋转角度范围内，利用傅里叶变换的相位相关技术判定旋转的角度和平移量，实现图像的配准，然后实现图像的拼接。实验结果表明，该方法能够对重叠面积较小的图像序列实现快速拼接。 　　关键词：傅里叶变换;红外热波;图像拼接 　　中图法分类号：TP391文献标识码：A 　　文章编号：1001－3695(2007)01－0227－02 　　 　　红外热成像技术是借助于物体的热辐射以得到其热图的一项技术。红外热波无损检测是一种新兴的无损检测技术，针对被检物的材质、结构、缺陷类型以及特定的检测条件，设计不同特性的热源，对检测材料进行主动加热，利用红外热成像技术对材料表面的辐射分布进行数据采集，生成热波图像。将不均匀材料内部热传导的不连续性反映到热像图上，从而进行材料检测。该技术具有快速、高效、直观的优点，特别适合复合材料的检测[1，2]。 　　红外热波无损检测技术在实际应用中，每次检测面积有限，当被测物较大时，不能一次获得整个被测物的图像，这就需要划分区域，分别获得各个部分的图像，然后将各次所成的热像拼接成一幅全图进行分析。?? 　　图像拼接的步骤可以分为图像配准和图像融合。图像配准的方法可以分为基于图像灰度信息的方法、基于图像特征的方法和基于变换域的方法。基于图像灰度的方法一般不需要对图像进行复杂的预先处理，而是利用图像本身具有的灰度统计信息来度量图像的相似程度，较容易实现。但是该方法对图像的变形、噪声、对比度差异敏感和光照条件敏感，且运算量很大。基于图像特征的方法和基于变换域的方法首先对待配准图像进行预处理，提取特征，通过特征的匹配关系建立图像之间的配准映射关系。该方法的局限性在于待配准的图像特征要易于提取，且要易于建立两幅图像之间同名点的关联。最主要的基于变换域的方法就是傅里叶变换方法。图像的平移、旋转、镜像和缩放等变换在傅氏变换域都有相应体现，适合光源变化采集的图像。但是傅里叶变换方法只能用来配准灰度属性有线性正相关的图像，图像之间也必须严格满足定义好的变换关系，实际使用中要求有较大范围的图像重叠区域[3]。?? 　　红外热波无损检测主要针对各种复合材料，热波图像本身的特征提取困难，空间相邻的两幅图像间存在局部小面积重叠，由于工作人员是手持工具，还会造成上下错位和小幅度旋转。本文针对以上特点，对待配准图像在可能的旋转角度范围内，先按一定角度旋转，然后利用傅里叶相位相关法计算平移量和脉冲幅值，通过比较这些脉冲幅值的大小，确定实际旋转的角度和偏移量。?? 　　 　　1理论[?C4~6]?? 　　 　　1．1两幅图像仅有平移变化?? 　　 　　1．2两幅图像存在小幅度旋转和平移?? 　　式（1）进行反变换所得冲击函数δ(x-x??0,y-y??0)的幅值大小，同时反应了两幅图像的相关程度，两幅图像越相关，幅值越大；当两幅图像内容完全一样，只存在平移时，则幅值为1。当两幅图像存在小幅度旋转时，旋转幅度越小，冲击函数的幅值越大。在可能的旋转角度范围内，按照一定的步长，将待配准的图像旋转，根据1.1节中的算法分别计算出一系列脉冲幅值，比较其大小，取出最大值，与之相对应的旋转角度和平移量即为待配准图像的旋转角度和位移量。?? 　　 　　2实验?? 　　 　　为验证本算法的适用性和精度，本文将飞机表皮铝合金材料的热波图像分别裁剪为如图1所示的两张图。图1（a)位于原图左侧，图1（b）位于右侧且相对于原图有逆时针3.6°的旋转，两幅图像重叠部分约占图像宽度的20%。假设在实际操作中，待配准图像即第二幅图像相对偏转角度不超过±10°。具体配准步骤如下: ?? 　　(1)对待配准图像在[-10,10]的角度范围内作旋转，每次旋转2°，共旋转11次。?? 　　(2)对旋转后的图像与图1（a)按照1.1节中的方法，分别求傅里叶变换，进而求得互功率谱，对互功率谱求傅里叶反变换，取反变换中的最大值，即脉冲幅值。各旋转角度图像对应的脉冲幅值如图2所示。?? 　　(3)比较求得的所有脉冲值，找出最大者，与之相对应的旋转角度和位移量即为所求。?? 　　如图2所示，可得旋转角度为-4°，由于步长为2°，进一步在 [-5,-3]角度范围内求精，每次旋转0.2°，分别求得脉冲值如图3所示。 　　 　　旋转误差在一个像素范围之内，完全能够满足要求。将待配准图像顺时针旋转3.6°，按照1.1节中的方法求得位移量，然后将两幅图像拼合在一起，如图4所示。 　　 　　3结论?? 　　 　　实验表明:本算法对于两幅相邻的热波图像间存在局部小范围重叠，小幅度旋转时，利用傅里叶变换的相位相关法（其中的脉