多聚焦图像融合算法研究

西北工业大学明德学院本科毕业设计论文本科毕业设计论文题 目 多聚焦图像融合算法研究 专业名称学生姓名指导教师毕业时间西北工业大学明德学院本科毕业设计论文毕业 任务书一、题目多聚焦图像融合算法研究二、指导思想和目的要求本题目来源于科研，主要研究多聚焦图像的概念，学习多聚焦图像的常用融合算法，进而实现相关算法。希望通过该毕业设计，学生能达到：1利用已有的专业知识，培养学生解决实际工程问题的能力；2锻炼学生的科研工作能力和培养学生团队合作及攻关能力。三、主要技术指标1学习多聚焦图像的特点；2研究多聚焦图像的融合算法；3实现多聚焦图像的融合。四、进度和要求第 01 周-第 02 周： 参考翻译英文文献；第 03 周-第 04 周： 学习多聚焦图像的特点；第 05 周-第 08 周： 研究多聚焦图像的融合算法； 第 09 周-第 14 周： 编写多聚焦图像的融合程序；第 15 周-第 16 周： 撰写毕业设计论文，论文答辩。5、主要参考书及参考资料1.张德丰.MATLAB 数字图像处理M.北京：机械工业出版社，2012.2. 敬忠良. 图像融合理论与应用M.北京：高等教育出版社，2010.3. 郭雷.图像融合M. 北京：电子工业出版社，2011.4. 孙巍. 孙巍. 像素及多聚焦图像融合算法研究D.长春：吉林大学，2008.5. 马先喜. 多聚焦图像融合算法研究D.无锡：江南大学， 2012.学生 指导教师 系主任 _ _设计论文西北工业大学明德学院本科毕业设计论文I摘 要图像融合是将同一对象的两个或多个图像按一定规则合成为一幅图像。其关键是抽取每幅源图像中的清晰区域，并将这些清晰区域以一定的规则融合起来，从而生成一幅清晰且信息量完整的融合图像。多聚焦图像融合的具体目标在于提高图像的空间分辨率、改善图像的几何精度、增强特征显示能力、改善分类精度、替代或修补图像数据的缺陷等。本文概括了多聚焦图像融合的一些基本概念和相关的基本知识，对 DWT 分解的层数和方向子带的个数对融合结果的影响进行了初步的研究。并就加权平均法、单层 DWT 分解、二层及二层以上 DWT 分解对多聚焦图像的融合进行了算法研究和编程实现，并对这些方法的仿真结果进行了比较分析。仿真结果表明，基于空间域的加权平均法的融合效果非常一般，在图像的细节表现力方面存在很大的不足；而基于变换域中的小波变换的低频取平均、高频取绝对值最大的融合算法在小波分解层数达到三层时，所得融合图像的性能指标，如信息熵、空间频率和清晰度都较为理想，达到了预期目的。可见多层 DWT 分解融合方法具有较高的应用价值，可以将其用于图片的判读分析，如指纹识别、人脸鉴别、不完整图片的复原等。关键词：图像融合，多聚焦图像，加权平均，DWT西北工业大学明德学院本科毕业设计论文IIABSTRACTImage fusion is two or more images of the same object according to certain rules of synthesis of an image. The key is to extract every source image in a clear area, and these rules must be clear and the region to come together. The amount of information to generate a clear and complete fusion image is its ultimate goal.Multi-focus image fusion specific objectives is to improve the spatial resolution of the image, to improve the image geometric accuracy, enhanced features display capabilities, improve classification accuracy, replace or repair the defects and other image data.This paper summarized some basic concepts and knowledge of multi-focus image fusion and did a preliminary study about the impact of DWT decomposition layers and the number of directional subbands related to multi-focus fusion results. And the weighted average method, single DWT decomposition, second floor and the second floor above DWT decomposition of multi-focus image fusion algorithm performed research and programming, and the simulation results of these methods were compared.It was found that the effect of the weighted average method of integration is very general and there is a big shortage in detail expressive aspects of the image.However, based on averaging the low-frequency and getting the largest absolute value of high frequency of wavelet transform fusion algorithm in transform domain performed good when wavelet decomposition level reaches three.The resulting image fusion performed indicators, such as information entropy, spatial frequency and clarity are more ideal to achieve the desired purpose.It is visible that multi-DWT decomposition fusion method has a high value, and its interpretation can be used to image analysis, such as fingerprint recognition, face identification,incomplete pictures recovery and so on.KEY WORDS: image fusion，multi-focus image，weighted average，DWT 西北工业大学明德学院本科毕业设计论文III目 录第一章 绪论 .51.1 课题研究的背景和意义 .51.2 研究现状及存在的问题 .51.2.1 研究现状 .51.2.2 存在的问题 .61.3 多聚焦图像融合技术 .71.3.1 多聚焦图像成像特点 .71.3.2 多聚焦图像融合层次 .81.4 本文章节的安排 .8第二章 多聚焦图像融合的基础理论 .102.1 多聚焦图像成像原理 .102.2 多聚焦图像融合主要算法 .112.2.1 空间域多聚焦图像融合 .122.2.2 变换域多聚焦图像融合 .172.3 多聚焦图像融合质量评价 .192.3.1 主观评价.192.3.2 客观评价.192.4 本章小结 .21第三章 基于小波变换的多聚焦图像融合.223.1 引言 .223.2 小波变换的基础理论 .223.2.1 小波定义.223.2.2 连续小波变换 .233.2.3 离散小波变换 .24西北工业大学明德学院本科毕业设计论文IV3.3 小波的分解.253.3.1 小波基的选择与分解层数 .263.4 算法仿真与结果分析 .273.4.1 算法仿真 .273.4.2 融合结果分析.303.5 本章小结 .30第四章 总结与展望 .314.1 全文总结 .314.2 工作展望 .32参考文献 .33致 谢 .35毕业设计小结 .36附 录 A.37附 录 B.38附 录 C.40附 录 D.43西北工业大学明德学院本科毕业设计论文5第一章 绪论1.1 课题研究的背景和意义自第一台具有光学镜头相机的出现开始，便出现了能生成场景中物体信息的图像。而图像融合 1是将同一对象的两个或多个图像按一定规则合成为一幅图像，以便它比原来的任何一幅更能容易的为人们所理解。这一技术可应用于多频谱图像理解、目标识别、机器视觉及医学图像处理等领域，在这些场合，同一物体部件的图像往往是采用不同的成像机理得到的。传统的图像融合方法主要是在时间域通过算术运算实现融合，具有算法简单直观，融合速度快，适合实时处理等优点，但其没对频率变化进行考虑。而多分辨率图像融合算法则是在频率域实现图像融合。根据分解形式的不同，多分辨率图像融合算法又可分为多分辨率金字塔方法和小波变换方法。近年来，基于小波变换的图像融合越来越受到重视。由于人的视网膜图像是在不同频带上分别以不同算子进行融合，而基于小波分解的图像融合也是以同样方式进行，所以，其可获得与人的视觉特性更为接近的融合效果。多聚焦图像融合是信息融合的一个重要分支，近年来被越来越多的研究机构所重视。由光学系统成像原理可知，物体在镜头的共轭物平面时，成清晰的像，处于共轭物平面以外的物体，将呈现出不同程度的模糊，当模糊程度不超过光学系统的景深时，则物体所成的像仍可以看做清晰像。光学系统的景深指保证在像平面能获得清晰的像的物体在物放空间前后移动的最大距离 2。由于光学镜头的景深有限，使得在拍摄时很难获取一幅所有景物都聚焦清晰的图像 3。因此仅仅依靠成像系统本身很难解决这一问题。但是可以加上一些后续处理，就可方便地获取一个场景清晰的图像，即采用图像融合算法。这一技术定将在现代医疗、军事侦查、产品防伪及文化艺术方面占据愈发重要的地位 3。1.2 研究现状及存在的问题西北工业大学明德学院本科毕业设计论文61.2.1 研究现状就研究层次和研究算法来看，早期的图像融合由于技术不发达，而且研究的人员较少，因此大多是针对像素级的。此阶段的方法有滤波法、加权法、IHS 变换等。早期的算法都是属于空间域上的算法，没有很多计算，复杂性较小。在融合处理阶段，不需要对融合的图像进行分解变换，只是对每幅图片中的单个离散的像素值进行简单的数据运算。80 年代中期金字塔算法被许多学者陆续提出，其中有比率低通金字塔、拉普拉斯金字塔、梯度金字塔等。这些算法都具有较强的实用性，所以有学者开始尝试将这些技术应用到图像处理中来，如多聚焦图像的融合，可见光图像和红外图像的融合等。90 年代初期，又有学者提出变换域的小波方法，随着小波理论的逐渐丰富，小波变换技术为图像融合技术提供了新的方法与思路 8。而且学者及科学家对小波算法进行了众多改进，其主要针对的是小波分解后低频系数和高频系数的选取，使图像融合技术不断上升。随着不断研究，人们发现金字塔算法和小波变换都有其局限性，所以提出了一些改进算法。Contourlet 变换 5是一种灵活的多尺度、局部的、方向性的分析算法，而且可以很好地捕捉图像中的几何结构，克服了小波变换中一个很难解决的问题平移不变性问题，因此陆续有人提出将 Contourlet 变换用于图像融合中。1.2.2 存在的问题随着人们对多聚焦图像融合研究的不断深入，一些问题也随之产生，目前主要有如下几个问题：1.算法研究尚不够完善。常用的算法目前有金字塔变换，小波变换等。但这些算法都存在一些缺点和不足。其中小波变换算法中存在的问题有：(1) 不能最优地表示含有线或面奇异的高维函数 6；(2) 只能捕捉到水平、垂直和对角线方向的信息；(3) 无平移不变性。2. 待融合的多聚焦图像的质量很难得到保证。在多聚焦图像数据方面，图像易受噪声等外界因素的干扰，并且不同传感器所获得的图像差异较大，往往对融合的性能产生较大的影响。3. 图像融合没有统一的评价标准。目前有主观评价和客观评价。主观评价容易因人的情绪、爱好等因素影响；客观评价有信息熵、均方根误差、空间频率、清晰度等众多参数，西北工业大学明德学院本科毕业设计论文7它们只能反映图像质量中的一个信息。一个算法用一个参数来衡量时可能是好的，但用另一个参数来衡量时，可能又是差的。所以需要制定一个统一、完善的融合标准出来。4. 融合模型不统一。在多聚焦图像融合模型方面，目前仍缺乏统一的融合模型，特别是数学模型。1.3 多聚焦图像融合技术多聚焦图像融合的关键是在源图像中寻找每幅图像中的清晰区域，滤除图像中的模糊部分，然后将这些清晰区域以一定的算法进行融合，从而得到一幅清晰且信息量完整的图像。多聚焦图像融合示意图如图 1-1 所示。图像 A 融合图像图像 B 图 1-1 多聚焦图像融合示意图1.3.1 多聚焦图像成像特点多聚焦图像指的是由同一传感器（如相机）在成像条件相同而镜头聚焦目标不同时所得到的多个图像。多聚焦图像的产生主要是由光学镜头的景深所决定的 9。当人或事物处在光学镜头的景深范围外时，就会引起图像模糊；反之，所得图像清晰。所以，多聚焦图像的产生主要是由成像镜头所造成的。可以认为，对于场景中的任一人或物都存在一幅聚焦良好的多聚焦图像。为了便于说明，论文仅对两幅多聚焦图像进行融合，这就要求用于实验的源图像满足以下条件：清晰区域模糊区域清晰区域清晰区域模糊区域西北工业大学明德学院本科毕业设计论文81. 对同一场景成像。2. 对场景中的不同物体聚焦。3. 对于场景中的任一区域，可在其中一幅源图像中找到其对应的清晰区域。在实际应用中，为了得到一个场景中每个人或物都清晰的图像，就需要对多幅多聚焦图像采用一定的算法进行融合操作，提取出每幅图像中的清晰区域，并将这些清晰区域整合起来，生成一幅能够包含完整场景的融合图像。1.3.2 多聚焦图像融合层次按照图像融合所处理的层次的不同，目前主要分为三层：像素级图像融合、特征级图像融合及决策级图像融合 9。这三个层次是从低到高的，具有阶段性和递进性，不仅可以单独进行，而且可以组合使用。实际应用中，针对所给的条件和目标，可以选择及结合不同层次融合的特点，以达到最优效果。像素级图像融合方法简单并且提供的信息丰富、可靠，便于图像的进一步分析、处理与理解 10。其他的两种融合策略，虽然都可以进行实时处理，且都可以进行信息压缩，但是算法的复杂性比像素级的高，而且这两种算法目前还不是研究的重点，还没有实用的算法来实现它们。1.4 本文章节的安排本文共分为五章，具体内容和章节安排如下：第 1 章：绪论本章从总体上简要介绍了课题研究的背景及意义，并指出了目前多聚焦图像融合的研究现状和所存在的一些问题。同时简述了多聚焦图像融合的成像特点和融合层次，为后续章节的展开做铺垫。第 2 章：多聚焦图像融合的基础理论本章首先介绍了多聚焦图像的成像原理，并概要介绍了多聚焦图像融合的基本算法：空间域算法和变换域算法。在空间域算法中主要介绍了基于像素点和基于分块的算法。在变换域算法中主要介绍了小波算法。此外还简述了融合结果的评价标准，并列出了当前常用的一西北工业大学明德学院本科毕业设计论文9些评价指标。第 3 章：基于小波变换的多聚焦图像融合算法本章首先介绍了小波变换的基础理论，包括小波的定义、连续小波变换和离散小波变换，并以图解的方式对二维离散小波分解进行了可视化分析。此外，还简述了小波基及小波分解层数对图像融合质量的影响，并就基于小波变换的低频取平均值、高频取绝对值大的融合算法进行了算法仿真。最后对融合图像进行了主观分析，并以信息熵、空间频率和清晰度三项指标对融合图像的质量进行了客观分析。仿真结果论证了该算法的有效性。第四章：总结与展望本章总结了全文的主要内容，并展望了多聚焦图像融合算法的未来。西北工业大学明德学院本科毕业设计论文10第二章 多聚焦图像融合的基础理论 2.1 多聚焦图像成像原理P 聚焦平面 非聚焦平面fD P 2RU Vs 图 2-1 光学镜头成像原理图图 2-1 表示了相机光学镜头的成像原理 9，其中 U 表示物距，V 表示焦平面的像距，S表示非聚集平面的像距。由成像原理可知，在聚焦平面成像时，图像清晰；在非聚焦平面成像时，图像模糊 4。如上图可，可做以下论证：根据三角形相似定理，有（2-VSRD21）经过变换，可以得到（2-VSR22）光学镜头是凸透镜，由其成像方程可知西北工业大学明德学院本科毕业设计论文11（2-VUf13）将式（2-3 ）代人式（2-2 ）中，可得（2-)1(21SUfDR4）其中，R 表示物体成像时模糊圈的大小。假设 R 一定，改变物距 U 的大小，由（2-1VUf5）可得：（2-fV6） 其中， V等于图中的非聚焦平面的像距 S，由式（2-2）可知 ，即：DRVS2（2-DRV27）时，有V（2-RV28）否则，（2-DV29）将式（2-8 ）和（2-9 ）分别代入式（2-6） ，可得（2-DRUf2)1(max西北工业大学明德学院本科毕业设计论文1210）（2-DRUf2)1(min11）（2-minaxOF12）R 一定时 ， 分别为 U 的最大值和最小值。两者之差 表示镜头景深。maxUin DOF2.2 多聚焦图像融合主要算法多聚焦图像融合算法要解决的主要问题是如何准确提取源图像的清晰区域。其关键在于正确地判断清晰区域的位置。这里所说的清晰是人类视觉系统主观上的感觉，而融合时使用的图像是由像素点组成的，图像的效果由各像素点的灰度值决定。因此，融合算法需要将各像素点的灰度值与人类主观上的“清晰区域”联系起来，寻找“清晰区域”中像素点应具备的特征，这样才能准确找到“清晰区域”并将其提取出来。可见，寻找像素值与“清晰区域”之间的对应关系是多聚焦图像融合算法首先要解决的问题。目前像素级多聚焦图像融合算法主要分为两类：1. 空间域多聚焦图像融合算法。2. 变换域多聚焦图像融合算法。2.2.1 空间域多聚焦图像融合空间域融合算法不需要对配准后的源图像进行任何变换，仅根据源图像中像素点的灰度值进行融合，算法简单，容易执行。根据融合时所取源图像中区域大小的不同，空间域融合算法分为基于像素点的融合算法和基于分块的融合算法。1. 基于像素点的融合算法该类算法是根据各个源图像中对应位置上每个像素点的灰度特征值进行融合的。计算公西北工业大学明德学院本科毕业设计论文13式如下：（2-),(,),(,),( nmLwnmLwnmFBA13）其中，A，B 为源图像， 为图像在 点的灰度值，),(),(m=1， 2， ，M ；n=1， 2，N， 为源图像大小， (m,n)表示图像中像素点的位置，M和 为加权系数，满足 。AwB 1BAw式（2-13 ）中加权系数 和 根据像素点的特征值计算。按照像素点特征值的不同，基于像素点的融合算法可分为基于单个像素点的融合算法和基于区域的融合算法，如图 2- 2所示。基于单个像素点的融合基于区域的融合融合图像 输入图像 图 2-2 基于像素点的融合算法示意图（1）基于单个像素点的融合算法西北工业大学明德学院本科毕业设计论文14这种算法仅根据每个离散的像素点灰度值的大小计算加权系数。该算法下的加权系数是独立的，相邻像素之间没有相关性。常见的灰度值取大法、灰度值取小法、加权平均法 10都属于这种算法。对于式（2-13） ，当 和 中有一个数为 0，另一个数为 1 时，即为灰度值AwB取大（取小）法。当 = =0.5 时，即为加权平均法。这些融合算法计算简单，但是降低了AB图像的对比度，融合图像中的边缘、轮廓等容易出现模糊，纹理不清晰，难以得到令满意的融合结果。这主要是因为单个像素点的灰度值大小不能体现该像素点处的清晰或模糊特征。（2）基于区域的融合算法基于区域的融合算法根据像素点所在区域的特征因子计算加权系数。该算法不仅考虑了单个像素点的灰度值，还考虑了该像素点所在区域内各像素点之间的相关性，在体现像素点所在区域的清晰度方面得到了改善。与基于单个像素点的算法相比，提高了提取源图像清晰区域的准确性。常见的区域特征因子有邻域能量、邻域方差和邻域梯度，计算公式如下：（区域为33）1）邻域能量 （2-snmii nymxfHyxE, 2),(),),(14）（2-08/12/15）2）邻域方差（2-1 2),(),(9/),(mniii yxfnyxfyx16）（2-),(9/11 nyxffmnii17）3）邻域梯度西北工业大学明德学院本科毕业设计论文15（2-),(),),(, nymxfnHyxGsnmii 18）（2-1819）其中，i=A 或 B（A,B 为源图像） ， 为 i 图像在 点的灰度值，S 代表窗口),(yxfi ),(yx区域， H 为掩模矩阵。 加权因子 ， 的求取方法有两种：特征因子取大法和特征因子加权法。AwB特征因子取大法计算公式如下：1）邻域能量取大法（2-),(),(01),( yxEyxwBBAA20）2）邻域方差取大法（2-),(),(01),( yxyxwBAA21）3）邻域梯度取大法 （2-),(),(01AyxGwB22）特征因子加权法计算公式如下：1）邻域能量加权法（2-),(),(),(yxEyxwBAA23）2）邻域方差加权法西北工业大学明德学院本科毕业设计论文16（2-),(),(),(yxyxwBAA24）3）邻域梯度加权法（2-),(),(),(yxGyxwBAA25）特征因子加权法是一种“折中”算法，在提高模糊区域图像质量的同时，降低了清晰区域的图像质量，融合图像不是最佳的。但是，多聚焦图像满足“模糊区域不重叠” ，这使得最佳融合图像的获得存在可能。因此，多聚焦图像融合应该采用特征因子“取大法” ，从而得到最佳的融合图像。2. 基于分块的融合算法基于像素点的融合算法得到的融合图像中相邻像素点之间一致性差，虽然可以对其进行一致性校检，但仍难以避免这一现象。针对这一问题，研究者们提出了分块融合算法。分块融合法是多聚焦图像融合特有的算法，该算法充分利用了多聚焦图像的特点，能够准确地提取源图像的清晰区域，保留源图像的有用信息，提高了像素点之间的相关性。分块融合算法首先将源图像分别分割成若干个子块，通过计算对应子块的特征因子，将源图像中清晰区域的子块放到融合图像对应的位置上。可见，子块特征因子能否准确反映子块的清晰度，直接影响融合图像的质量。融合过程如图 2-3 所示。分割 重构 分割合并子块 融合图像图 2-3 分块融合算法示意图目前常用的子块特征因子有以下几种:源图像 A源图像 B规则西北工业大学明德学院本科毕业设计论文171. 空间频率（SF） 12：反映子块的活性程度，定义如下。（2-2CFRS26）（2-MmNnnff122)1,(),(27）（2-MmNnnffCF212),1(),(28）其中，RF 和 CF 分别代表行频率和列频率，MN 为子块大小。2. 显著度参数（VI） 9：反映子块的均匀度，定义如下（2-kByxkk myxfnmVI),( ),(1)(29）（2-kByxkfn),(,130）（2-)1(kkm31)其中， 为源图像的某个区域，大小为 ， 为 处的灰度值， 为一固定值，kBn),(yxf,在 0.6 到 0.7 之间取值。3. 边界特征 12：反映图像子块的边界数量。对于两幅多聚焦图像，清晰的图像比模糊的图像具备更多的边界。基于空间频率的分块融合算法步骤如下：步骤 1 将源图像 A，B 分割成若干个大小为 的子块，设 ， 分别表示图像 ANMiAiB和 B 的第 i 个子块。步骤 2 分别计算子块 和 的空间频率 和 。ii iASFiB西北工业大学明德学院本科毕业设计论文18步骤 3 比较两幅图像对应子块的空间频率 和 ，得到融合图像的第 i 个子块 。iASFiB iF(2-32)othersTHBAFii iBiiiii ,2/,其中，TH 为阈值。 步骤 4 对所有子块进行上述操作，将融合子块重构得到融合图像。为了进一步提高子块选取的准确性，李树涛 12等通过人工神经网络进行子块选取。这些算法虽然可以准确地选取清晰子块，但是子块尺寸的选择难以统一：子块尺寸过大将导致清晰区域、模糊区域位于同一子块；子块尺寸过小则不能充分反映子块特征，容易造成子块的误选。由于子块尺寸的选择难以划分源图像的清晰区域与模糊区域，融合图像容易出现块效应，这正是分块融合算法存在的最主要问题。空间域多聚焦图像融合算法根据像素点或子块的灰度值特征因子进行融合。但是，图像的清晰程度主要体现在物体的边缘等细节上，与像素点的灰度值之间不存在直接的对应关系。可见，空间域融合算法缺乏细节表现力，因此难以准确选取源图像的清晰区域。2.2.2 变换域多聚焦图像融合空间域多聚焦图像融合算法缺乏细节表现力，融合图像存在对比度低、模糊以及块效应等问题。人类的视觉系统对事物的局部对比度变化（边界、拐角等）十分敏感。一般来说，图像中物体边缘的对比度变化较大，图像的清晰程度主要体现在物体的边缘等细节上。可见，融合算法的细节表现力决定了算法提取清晰区域的能力。多分辨率分析方法具有良好的细节表现能力，被广泛应用于多聚焦图像的处理，形成了变换域多聚焦图像融合算法。该方法首先对各个源图像进行多分辨率分解，得到源图像在不同分解层上的子带图像。然后对分解得到的各子带图像进行融合，融合时，不同分解层上的子带图像可以采用不同的融合准则。最后，将融合后的子带图像进行逆变换，得到融合图像。根据多分辨率分析方法的不同，变换域多聚焦图像融合算法主要有基于金字塔变换的融合算法、基于小波变换的融合算法和其他多尺度分析的融合算法。 1. 基于小波变换的融合算法西北工业大学明德学院本科毕业设计论文19小波变换是一种类似于金字塔变换的多分辨率分析方法。小波变换将图像分解成低频子带和不同尺度、不同方向的高频子带。其中，低频子带包含图像的平均信息，占整个图像的大部分能量；高频分量包含了图像在不同尺度、不同方向上的细节信息。基于 DWT 的图像融合过程如图 2-4 所示，概括如下：步骤 1 对源图像 A 和 B 进行二维 DWT 分解，得到各源图像的低频分量和高频分量。步骤 2 根据低频和高频分量的特点，按照各自的融合算法进行融合处理。步骤 3 对以上得到的高低频分量，经过小波逆变换得到最终融合图像 F。DWT源图像 A融合规则 融合图像 DWT 融合后的小波分解系数源图像 B 图 2-4 基于 DWT 的图像融合过程在基于多分辨率分析的图像融合算法中，系数融合准则的好坏直间影响到融合图像的质量。因此，大多数基于多分辨率分析的融合算法都致力于系数融合准则的研究。根据融合对象的不同，融合准则分为低频系数融合准则和高频系数融合准则。1. 低频系数融合准则在基于小波变换的融合算法中，低频系数的融合多采用加权平均法，定义如下(2-33) ),(5.0),(5.0),( nmCnnmCBAF其中，A 和 B 为源图像，F 为融合图像， 为图像 i 分解得到的低频子带在i,处的系数。）， nm(加权平均法将清晰区域的系数与模糊区域的系数进行“结合” ，在提高模糊区域质量的同时，造成了清晰区域信息的损失，表现为融合图像的对比度下降。2. 高频系数融合准则在图像的小波分解中，绝对值较大的高频系数对应于图像中对比度变化较大的边缘特征融合算子 F西北工业大学明德学院本科毕业设计论文209，而人眼对于这些特征比较敏感。所以，对于高频系数多采用绝对值取大的融合准则，原理如下:(2-34),(),(,)(S),(, , nmnmSjiBjiAjiBjijijiAjiF若若其中，j=1 ， 2，.，N 代表分解的层数，i=1,2,3 代表每一层高频系数的 3 个方向。2.3 多聚焦图像融合质量评价 采用科学合理的算法对融合图像的质量进行评价，对于在实际应用中选择适当的融合算法、对现有融合算法的改进及研究新的融合算法等都是非常重要的。多聚焦图像融合结果的评价主要分为主观评价和客观评价。2.3.1 主观评价主观评价法就是依靠人眼对融合图像的质量进行主观评估的方法。主观评价法具有简单、直观的优点，对明显的图像信息可以进行快捷、方便的评价。这种算法对图像中一些较为明显的信息，如是否有重影，边缘信息是否损失，纹理和色彩是否丰富等能很快做出判断。但是，图像的视觉质量主要取决于观察者，具有主观性强、片面、可重复性差等缺点。因此，需要与客观评价标准相结合进行综合评价。2.3.2 客观评价常用的判断多聚焦图像融合好坏的评价指标有以下几种：1.均方根误差 6RMSE 表示标准图像与融合图像灰度差的统计。RNSE 越小，说明标准图像与融合图像差异越小，融合效果越好。其定义如下：西北工业大学明德学院本科毕业设计论文21(2-35)NMnmFRRSEmn12),(),(其中，R ，F 分别为标准图像和融合图像，大小为 MN。2. 信息熵 6信息熵用来描述一幅图像所含信息量的大小。一般来说，熵值越大，图像所含的信息量就多，融合效果就越好。其定义如下：(2-36)102lLiipogE其中，L 为图像的灰度级， 为图像中灰度值为 的概率， 为熵值。i iE3. 峰值信噪比 9PSNR 表示图像中所含高频信息量的大小。其值越大，融合图像就越好。定义如下：(2-37)MxNyyxFRPSNR1210 ),(),(,minalog4. 清晰度 12图像在人眼中的清晰程度就是清晰度。一般采用平均梯度作为清晰度指标，平均梯度能敏感地反映出图像中的微小细节反差，平均梯度越大，图像越清晰。其定义如下：(2-38)1( 22),(),()(1NMyxfxfg其中， 为窗口的大小，计算时，一般用 来代替 ，用NM ),1(),jifjifxyf),(来代替 。