专利名称：：基于混合多分辨率分解的图像融合方法
技术领域：
：本方法涉及一种图像融合方法，特别涉及一种基于混合多分辨率分解的图像融合方法。
背景技术：
：图像融合指按照一定的规则，把同一目标或同一场景的多个传感器的成像或单一传感器的多次成像进行处理，生成一幅新的图像。新图像与原图像相比，信息更全面、精确和稳定。它综合了传感器技术、信号处理、图像处理和人工智能等新兴技术，在军事领域和非军事领域如遥感图像、医学图像、机器视觉上得到了广泛地应用。图像融合一般可分为像素级图像融合、特征级图像融合和决策级图像融合。其中像素级图像融合是最基本的融合，是特征级和决策级图像融合的基础。图像融合最简单的方法是加权平均法，即将原图像对应像素的灰度值进行加权平均，生成新的图像。这种方法虽然简单但效果比较差。基于Laplacian金字塔、比率低通金字塔、梯度金字塔等图像多分辨率分解的融合方法是图像融合最常用的方法。但是塔形分解是冗余分解，而且不能够很好的表达图像的方向信息。90年代小波变换在图像处理中得到广泛应用，因为其具有的方向性、非冗余性以及具有快速算法等特点，在图像融合中得到了广泛的应用，可获得良好的融合效果。小波变换在分析点状瞬态特性的奇异性时是最优的，但在表示图像结构的直线和曲线奇异性时却不是最优的。
发明内容为了解决基于小波变换的图像融合方法存在的上述技术问题，本发明提供一种基于混合多分辨率分解的图像融合方法。本发明方法将基于非下采样的轮廓波变换和静态小波变换的融合方法结合起来，可有效提高融合结果图像的质量，达到比较理想的融合效果。本发明解决上述技术问题的技术方案包括以下歩骤将两幅源图像进行非下采样轮廓波分解，得到分解后的子带轮廓波高频和低频系数；根据高频于带系数计算其能量值，选择能量值最大的系数作为融合的高频轮廓波变换系数；对子带轮廓波分解系数屮的低频于带系数进--步进行静态小波变换，得到对应的静态小波变换系数，对静态小波变换系数中的高频分量，采用选择小波系数能量值最大的方法进行融合，对静态小波变换系数中的低频分量通过求平均的方法进行融合，得到融合图像对应的静态小波变换系数，对静态小波变换系数进行逆变换，得到融合的低频轮廓波变换系数；对融合的轮廓波变换系数进行逆变换，得到融合的图像。上述的基于混合多分辨率分解的图像融合方法中，所述能量值的计算包括以下步骤采用基于窗口的加权平均法，先取每个像素对应的所有高频系数的平方和，再对其进行基"？窗口的加权平均。本发明的技术效果在于本发明在图像融合的过程中对每-《昌源图像分别进行非下采样轮廓波变换，得到它们的轮廓波分解系数，对其高频分量，釆用选择能量值最大的方法进行融合；对低频分量采用基于静态小波变换的图像融合方法进行融合；对融合的轮廓波变换系数进行逆变换，得到最终的融合图像。本发明充分结合了静态小波变换的能够保留源图像角点、纹理等细节信息的特点以及非下采样轮廓波变换具有的方向性和各向异性，能够很好的表达源图像中边缘信息的优点。经实验对比，本发明的融合效果较好，能够使融合后的图像的质量得到较大的提高。下面结合附图对本发明f1进-一步的说明。图l为本发明的流程图。图2为本发明中多聚焦图像融合实例(a)聚焦在右边的图像；(b)聚焦在左边的源图像；(c)基-丁轮廓波变换方法得到的融合结果；(d)基于静态小波变换得到的结果；(e)本发明方法得到的融合结果。图3为医学图像融合实例(a)CT图像；(b)MRJ图像；(c)基于轮廓波变换方法得到的融合结果；(d)基于静态小波变换得到的结果；(e)本发明方法得到的融合结果。图4为红外可见光图像融合实例(a)红外图像；(b)可见光图像；(C)基于轮廓波变换方法得到的融合结果；(d)基于静态小波变换得到的结果；(e)本发明方法得到的融合结果。具体实施例方式参见图l，图l为木发明的流程图。其具体实施歩骤如下1.将输入的已经配准的源图像分别进行非下采样的轮廓波(NSCT)分解t下釆样轮廓波变换将多尺度分解和方向分解分开进行，首先使用非卜釆样塔式滤波器组(NonsubsampledPyramidFilterBank,NSPFB)对源图像进行多尺度分解，然后利用非下采样方向滤波器组(NonsubsampledDirectionalFilterBank,NSDFB)对多尺度分解得到的子带系数进行方向分解，从而得到不同尺度、不同方向的子带系数。非下采样塔式分解是一个在拉普拉斯分解过程中去掉下采样，对滤波器进行上采样的过程。当分解层数为J时，其冗余度为J+1，并且满足完全重构的Bezout恒等式//0(梶2)+"1(1)其中/AU，G。(z)为塔式分解的低通分解滤波器和合成滤波器，q("为塔式分解的卨通分解滤波器和合成滤波器。非下采样的方向滤波是一种在方向滤波中无下采样的滤波过程，通过atrous小波算法实现。^trous算法是将图像分解为不同频带上的细节信息和最低频带的近似信息，该细节信恩称为小波面，其大小与原图像大小相同。对于图像/(x，力，可逐级得到以下图像序列》,，聰x，力],,力=丄2[你,力]②乂:、(二>0=丄/^—其中，./;(x,力为尺度r下的近似图像，/。"力是/oc,力，^为低通滤波。it42,….,V。相邻尺度的近似图像间差异构成小波变换的高频系数即小波面AO,>o=y;(x，_y)—y;—,(x,力(3)其中《/(0c，.y)为第A层小波面。本发明对源图像A、B依次进行非下采样塔式滤波和非下采样方向滤波变换得到一系列的子带系数，记为WC7;"和A^C"^1，其中5=0,1,...,&"=1,2,...,2"，S为分解的子带数，2"为各子带的方向总数，AWC7;"和MSC7f"为低频子带，其它为方向高频子带。2.基于能量最大规则的高频子带融合6在非F采样轮廓波分解中，能量值较高的系数对应于图像中对比度变化较大的边缘特征，因此高频子带融合的关键是选取能量较高的部分。本发明采用基于窗口的能量值最大的融合规则，取对应像素能量值较大者作为高频子带融合结果。^先定义-一个人小为5x5的巻积模板ff:ff2564444441616164416641644161616444444妖后利用此模板对所有高频子带系数进行巻积，再对巻积结果在对应位置上的各个方向和尺度卜.求和，并将求和结果作为中心位置的能量。能量的整体讣算公式为'(,，/)=11过,+3，"+3)l服7^)(,++")1(4),—2〃=_2,、'=2ts代表尺度总数，d代表方向频率总数，/,y表示当前点在图像中的位置，/为横坐标，乂为纵坐标，w,为变量，其取值由巻积模板^的大小决定。最后高频系数的融合规则为<formula>formulaseeoriginaldocumentpage7</formula>(5)WSC"、"'和」V义T》、力分别为源图像a和b经过nsct分解后在尺度s方向d上的子带系数，」VSC7;^为融合后的对应系数。式(5)表示当a图像在(;c,力处能量大于b图像在此处的能量时，则在此处的各个尺度和方向的高频系数选取a图像在此处的系数，否则选取b在此处的高频系数。3.基于静态小波的低频系数融合低频子带系数w双T，'1和wser，'采用基于静态小波变换的图像融合方法进行融合。静态小波分解不对滤波结果进行下采样操作，而是对滤波器进行上采样。由于静态小波分解的滤波器是可分离的，因而可在图像的行和列分别进行滤波。以F是由j层得到j+1层的分解参数<formula>formulaseeoriginaldocumentpage8</formula>式中，4"是原始图像在尺度2'下的低频近似，『。和吋"分别是j层图像在水平，垂直和对角方向的高频细节信息，a,b表示当前点在图像中的位置，a为横坐标，b为纵坐标，k,l为变量，其取值由巻积模板A和g的大小决定。基P静态小波的低频系数融合实施方式如下1)对AWCTT1和A^C7f"分别做3层的静态小波分解，得到它们对应的静态小波金字塔系数，这些系数分为"频系数和低频系数。2)对丁-分解后的低频系数，采用加权平均算子进行融合F〃=—,+(1—w)xS〃(10)w为加权系数，在本发明中w取的取值为0.5，A,和L分别为MSCTf11禾口WlTf卩经静态小波分解后的低频子带系数。3)对于高频变换系数，仍然采用能量最大的方法进行融合。在各个尺度内的所有卨频于带上，首先基于5x5的窗口W获得对某--像素的区域特性描述(该像素居于正方形窗口的中心处)，然后合并该尺度上所有子带在对应々间位置处的区域特性来衡量该系数的重要程度，采用系数的能量来表示系数的軍要程度，其计算公式为<formula>formulaseeoriginaldocumentpage8</formula>窗口W为W=2564444441616164416641644161616444444，L代表尺度总数，《代表方向总数<比较两源图像在各个尺度内不同空间位置上变换系数的重要程度，通过选择最大重要程度的方法完成变换域内的系数融合，可用如下公式表示'/《、《(12),WTf"和'S,7;)"'分别是对WSC7f和MC7f"进行静态小波分解后在尺度/方向it上的子带系数，融合后的对应系数。4)对融合得到的多分辨率系数作3层静态小波逆变换(ff『T)，得到低频融合4.采用轮廓波逆变换得到最终融合结果对所得到的轮廓波系数W^^"1作非下采样的轮廓波逆变换，即非下采样的方向滤波重构和非下釆样的塔式重构，得到最终的融合结果图像。将本发明所提供的方法所得到的融合结果与单独使用静态小波变换方法和单独使用非下采样轮廓波变换方法所得到的结果进行了比较。图2、图3和图4给出了实验结果。表l列出了客观评价的结果。这里用了常用的信息熵(EN)，标准差(STD)以及QAB/F作为评价标准。对于三种评价标准，数值大表示融合结果越好。从实验结果我们可以看出本发明提出的方法优于传统的基于小波变换的方法。9表1.图像融合结果指标评价图2<table>tableseeoriginaldocumentpage10</column></row><table>权利要求1.一种基于混合多分辨率分解的图像融合方法，包括以下步骤将两幅源图像进行非下采样轮廓波分解，得到分解后的子带轮廓波高频和低频系数；根据高频子带系数计算其能量值，选择能量值最大的系数作为融合的高频轮廓波变换系数；对子带轮廓波分解系数中的低频子带系数进一步进行静态小波变换，得到对应的静态小波变换系数，对静态小波变换系数中的高频分量，采用选择小波系数能量值最大的方法进行融合，对静态小波变换系数中的低频分量通过求平均的方法进行融合，得到融合图像对应的静态小波变换系数，对静态小波变换系数进行逆变换，得到融合的低频轮廓波变换系数；对融合的轮廓波变换系数进行逆变换，得到融合的图像。2.根据权利要求1所述的基于混合多分辨率分解的图像融合方法，所述能量值的计算包括以下步骤采用基T窗口的加权平均法，先取每一个像素对应的所有高频系数的平方和，再对其进行基于窗口的加权平均。全文摘要本发明公开了一种基于混合多分辨率分解的图像融合方法。它包括以下步骤将两幅源图像进行非下采样轮廓波分解，得到高频和低频子带系数；根据高频子带系数计算其能量值，选择能量值最大的系数作为融合的高频轮廓波变换系数；对低频子带系数进一步进行静态小波变换，对其高频分量，采用选择小波系数能量值最大的方法进行融合，对其低频分量通过求平均的方法进行融合，得到静态小波变换系数，对静态小波变换系数进行逆变换，得到低频轮廓波变换系数；对轮廓波变换系数进行逆变换，得到融合的图像。本发明结合了静态小波变换和轮廓波变换能表达不同图像特性的优点，可有效提高融合结果图像的质量，达到比较理想的融合效果。文档编号G06T5/50GK101504766SQ20091004295公开日2009年8月12日申请日期2009年3月25日优先权日2009年3月25日发明者李树涛,斌杨申请人:湖南大学