基于梯度正则化的多聚焦图像融合方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN108399611A(43)申请公布日2018.08.14(21)申请号201810092718.6(22)申请日2018.01.31(71)申请人西北工业大学地址710072陕西省西安市友谊西路127号(72)发明人王健张修飞任萍杨珂(74)专利代理机构西北工业大学专利中心61204代理人刘新琼(51)Int.Cl.G06T5/50(2006.01)G06T5/10(2006.01)权利要求书3页说明书8页附图3页(54)发明名称基于梯度正则化的多聚焦图像融合方法(57)摘要针对多聚焦图像融合后图像细节保持能力有限，配准失调且敏感等问题，本发明提出一种基于梯度正则化多聚焦图像融合方法。首先，将源图像进行两尺度分解得到高频分量和低频分量；然后，将分解得到的高频分量利用梯度正则化的表示得到模型的稀疏系数并通过最大值策略获得融合后图像的高频分量，将两尺度分解得到低频分量利用取最大值融合策略得到低频分量；最后，将得到的融合后图像的高频分量和低频分量相加得到融合图像。相对其他三种融合方法，本发明不论在主观视觉和客观评价指标上还是在计算效率上，可以更好保留源图像的细节等纹理信息和显著性信息。CN108399611ACN108399611A权利要求书1/3页1.一种基于梯度正则化的多聚焦图像融合方法，其特征在于步骤如下：步骤1：两尺度图像分解首先将源图像Ik分解为一个低频分量和一个高频分量所述的低频分量表征多聚焦图像中的对比度和亮度信息，高频分量反映图像中的细节信息；其中，低频分量通过优化方式求解获得：式中：Ik表示输入的源图像，*表示卷积，η表示正则化参数，gx和gy分别表示沿图像行和T列的梯度滤波器，gx＝[-11]，gy＝[-11]；对式(1)进行快速傅里叶变换，其变换后数学表达式如下：利用二次优化可解得：高频分量通过源图像Ik减去的逆傅立叶变化获得，如下式：步骤2：高频分量的多聚焦图像融合方法图像高频部分系数图为ck,m，m∈{1,2,...,M}，通过求解如下正则化的表示模型获得：式中：{dm}表示一组M个字典滤波器集，*表示卷积，{ck,m}表示高频部分系数图，{αm}表示一组l1范数的系数权重，{βm}表示一组l2范数的系数权重，λ和μ表示正则化参数，g0和g1分T别表示沿图像行和列的梯度滤波器，这里g0＝[-11]，g1＝[-11]；将式(5)利用梯度正则化的ADMM算法求解，定义线性算子Dm、Gl使DmCk,m＝dm*ck,m，GlCk,m＝gl*ck,m，式(5)的最后项可相互变换如下：定义矩阵Γl式(6)的最后项进一步可变换如下：2CN108399611A权利要求书2/3页将式(5)进行傅里叶变换，其数学表达式如下：式中：D、α和C分别为块矩阵，并引入一个辅助变量y0，y1，y2，将式(8)变换为：利用对偶变量，引入拉格朗日乘子u0，u1和u2则通过迭代将式(9)的约束优化问题变为非约束优化问题形式：(j+1)(j)(j+1)(j+1)u0＝u0+C-y0(12)(j+1)(j)(j+1)(j+1)u1＝u1+C-y1(13)(j+1)(j)(j+1)(j+1)u2＝u2+C-y2(14)式(11)通过下式给出：Sγ(u)＝sign(u)⊙max(0,|u|-γ)(16)用DFT域中和分别表示Dm、Cm、Γ0、DΓ1、y0、y1、y2、u0、u1、u2和I，式(11)可变换为：对(17)中求偏导数，并令偏导数为0，可解得：式中：矩阵为由M个N×N对角矩阵组成，M为滤波器的维数，N为源图像的维数，为一个MN×MN维的对称矩阵，和分别为对角矩阵；采用谢尔曼-莫里森公式求解式(18)得到：然后对进行逆傅里叶变换，求得ck,m；假如ck,1:M(x,y)表示ck,m在空间域中位置(x,y)处的内容，ck,1:M(x,y)是一个M维向量，利用ck,1:M(x,y)的l1范数作为源图像的活动水平度量，则活动水平图Ak(x,y)通过如下表达式3CN108399611A权利要求书3/3页获得：Ak(x,y)＝||ck,1:M(x,y)||1(20)对Ak(x,y)利用基于窗口的平均策略以获得最终的活动水平图式中：r决定窗口的大小，r的取值越大，这种方法对于图像配准误读会更加稳健，但是同时可能会丢失一些细小的细节，在多聚焦图像中，多个源图像中的物体边缘具有不同的清晰度，使得每个源图像的物体边缘位置不确切的相同，因此，相对较大的r更适合多聚焦图像融合；在多聚焦图像融合取“最大值策略”融合后的系数图为：最后，融合图像的高频部分通过如下式被重构：步骤3：低频分量的多聚焦图像融合方法利用选择“最大融合策略”对源图像进行融合，其低频分量的融合结果表达式如下：步骤4：两