基于高阶奇异值分解和模糊推理的多聚焦图像融合方法
【专利摘要】本发明公开了一种基于高阶奇异值分解和模糊推理的多聚焦图像融合方法。（1）对待融合的两幅多聚焦源图像分别进行滑动窗口分块，将分块后相互对应的源图像子块形成子张量；（2）对子张量进行高阶奇异值分解(HOSVD)，采用基于模糊推理的加权平均融合规则对分解系数进行融合，通过提取分解系数的方差，区域能量，匹配度三个特征设计模糊推理规则；（3）对融合后系数进行HOSVD逆变换得到融合后图像子块；（4）将融合后的图像子块组合生成最终融合图像。本发明克服了传统多聚焦图像融合方法容易造成边缘失真的缺陷，很好地解决了图像融合中源图像对融合图像贡献程度的不确定性问题，使融合后图像质量得到明显提高。
【专利说明】基于高阶奇异值分解和模糊推理的多聚焦图像融合方法
【技术领域】
[0001]本发明涉及图像融合及其应用领域，特别是基于奇异值分解和模糊推理的多聚焦图像融合方法。
【背景技术】
[0002] 图像融合技术是图像处理领域的研究热点，多聚焦图像融合是图像融合领域的一个重要分支。由于光学成像系统景深的限制，很难将场景内的所有目标全部清晰成像。多聚焦图像融合方法能够将成像设备对某一场景焦距改变所得具有不同焦点图像的清晰部分进行整合，有效提取图像中的互补信息，获得场景中所有对象均清晰的融合图像，这为图像的后续处理奠定了良好的基础。
[0003]从目前的研究方法以及多聚焦图像融合分析工具来看主要有两大类方法:一是在空间域完成图像的融合，二是将图像转换变换域内进行融合。系数加权平均是空间域最简单的一种融合方法，方法简单，实时性好。但是融合结果对比度下降，当源图像有较大灰度差异时，融合结果质量欠佳。基于变换域的融合方法将待融合图像执行变换，以变换系数的形式融合图像。变换域的方法利用多分辨率的概念，研究者们提出了基于金字塔，离散小波变换，脊波，轮廓波等变换的融合方法以进行多尺度分解。多分辨率变换域内的系数表示图像的边缘信息，可以检测图像的显著特征，但系数和源图像的像素值没有一一对应的关系，一个系数的改变会影响到几个像素值，在反变换过程中，源图像的信息就会丢失，且由于源图像对融合图像贡献程度的不确定性，多分辨率方法在下采样过程中容易造成融合信息不完全的问题。
[0004]基于张量的图像信息处理方法更适合表示高维数据，提取图像的相关信息。高阶奇异值分解是最有效的张量分解技术之一，相比于其他多分辨率变换方法具有优越的性质。高阶奇异值分解具有完全数据驱动的优良特性，不需要设定阈值和参数。高阶奇异值分解的特征系数提供了更丰富的图像纹理信息。
[0005]因为高阶奇异值分解的分解系数具有纹理表示能力，所以利用高阶奇异值分解的方法为多聚焦图像融合提供了一种新的途径。
【发明内容】
[0006]本发明的目的在于针对上述现有技术中存在的不足，提出一种基于高阶奇异值分解和模糊推理的多聚焦图像融合方法，从而达到提高多聚焦图像融合效果的目的。
[0007]本发明的技术方案包括如下步骤:
[0008]I)准备待融合的两幅多聚焦源图像，对两幅多聚焦图像分别进行滑动窗口分块，将分块后I个相互对应的源图像子块形成I个子张量；
[0009]2)对子张量进行高阶奇异值分解(H0SVD)，得到分解系数；
[0010]3)采用基于模糊推理的加权融合规则(WAFR)对分解系数进行融合；
[0011]3.1)提取HOSVD分解系数的方差区域能量(RE)、匹配度(MDi)三个特征，并将三个特征的值域量化为模糊论域；
[0012]3.2)采用Mamdani极小化原理设计模糊融合规则，根据最大隶属度原则对权值进行模糊决策，得到分解系数的加权值，根据加权值得到融合后系数；
[0013]4)对步骤3)得到的融合系数进行HOSVD逆变换获得融合图像子块；
[0014]5)将步骤4)得到的融合图像子块组合生成最终融合图像。
[0015]所述步骤I)，包括如下步骤:
[0016]a)将两幅的多聚焦源图像MXN分别记为B(I)和B (2)，利用滑动窗口技术把B (I)和B (2)分别分成I个相互对应的Mx//源图像子块,滑动步长大小为r, r< M ,且r < N ,
其中
【权利要求】
1.基于高阶奇异值分解和模糊推理的多聚焦图像融合方法，其特征在于，包括如下步骤: 1)准备待融合的两幅多聚焦源图像，对两幅多聚焦图像分别进行滑动窗口分块，将分块后相互对应的源图像子块形成I个子张量； 2)对子张量进行高阶奇异值分解(HOSVD)，得到分解系数； 3)采用基于模糊推理的加权融合规则(WAFR)对分解系数进行融合: 3.1)提取HOSVD分解系数的方差(加)、区域能量(& )、匹配度(MDi)三个特征，并将三个特征的值域量化为模糊论域； 3.2)采用Mamdani极小化原理设计模糊推理规则,利用最大隶属度原则对各自的权值进行模糊决策得到加权值，根据加权值得到融合后系数； 4)对步骤3)得到的融合后系数进行HOSVD逆变换得到融合后图像子块； 5)将步骤4)得到的融合后图像子块组合生成最终的融合图像。
2.根据权利要求1所述的基于高阶奇异值分解和模糊推理的多聚焦图像融合方法，其特征在于，所述步骤I)包括如下步骤: a)将两幅MXN的多聚焦源图像分别记为B(I)和B(2)，利用滑动窗口技术把B(I)和B (2)分别分成I个相互 对应的Mx及源图像子块，滑动步长大小为r,且元。
3.根据权利要求1所述的基于高阶奇异值分解和模糊推理的多聚焦图像融合方法，其特征在于，所述步骤2)具体为进行HOSVD分解得到Si,其公式SAi =SiX ^X2Vi X具，其中Si为核张量，Ui, Vi, Wi分别为Ai的一维，二维，三维展开的左奇异值向量。为了便于得到融合结果，利用ξ代替Si作为图像块的特征，^ = SfX3巧。对于每一个图像块 Bi (k)，k = I，2，用公式表示为 (/f) = C/, X 瓦(:,:,f)x If.>
4.根据权利要求1所述的基于高阶奇异值分解和模糊推理的多聚焦图像融合方法，其特征在于，所述步骤3.1)包括如下步骤:a)对于瓦，以(x，y)为中心开tXt的邻域窗口，计算每个窗口的为了便于表示，记4=瓦k=l，2，表示第i块中第k个图像子块的系数，其计算公式为:
5.根据权利要求1所述的基于高阶奇异值分解和模糊推理的多聚焦图像融合方法，其特征在于，所述步骤3.2)包括如下步骤: a)采用Mamdani极小化原理设计8种模糊推理规则，即: ①如果Cil的方差占优，Cil的区域能量和占优，且Cil和Ci2不匹配，则规则R1定义为: R1:如果
6.根据权利要求1所述的基于高阶奇异值分解和模糊推理的多聚焦图像融合方法，其特征在于，所述步骤4)具体为:根据得到的融合后系数Di，融合后的图像子块Fi通过下面的公式得到:
Fi=UiXDiXVi1O
7.根据权利要求1所述的基于高阶奇异值分解和模糊推理的多聚焦图像融合方法，其特征在于，所述步骤5)具体为:按照源图像重叠分块的顺序对融合后图像子块进行拼接重构，生成最终的融合图像，并且将该融合图像输出显示。
【文档编号】G06T5/50GK103985104SQ201410057924
【公开日】2014年8月13日 申请日期:2014年2月20日 优先权日:2014年2月20日
【发明者】罗晓清, 张战成, 张翠英, 吴小俊, 李丽兵, 吴兆明, 高瑞超 申请人:江南大学