[发明专利]一种基于二维EMD和遗传算法的多聚焦数字图像融合方法

权利要求书

1.一种基于二维EMD和遗传算法的多聚焦数字图像融合方法，其特征在于步骤如下：

步骤1：采用序贯相似度检测匹配法对两幅聚焦不同的源图像进行图像配准，并利用线性变换方法将两幅源图像的灰度方位映射到一个一致的灰度区间，得到预处理后的两幅图像A和B；

步骤2：对预处理后的图像A进行二维EMD分解得到IMF分量对预处理后的图像B进行二维EMD分解得到IMF分量

步骤3：对分量进行T检验得到低频分量和高频分量对分量进行T检验得到低频分量和高频分量

步骤4：按低频分量融合规则和高频分量融合规则对和进行融合，得到融合后的低频分量IMF_F,L(m,n)和高频分量IMF_F,H(m,n)：

所述低频分量融合规则：

按IMFF,L(m,n)=IMFA,L(m,n)Entropy(IMFA,Lwindow(m,n))≥Entropy(IMFB,Lwindow(m,n))IMFB,L(m,n)Entropy(IMFA,Lwindow(m,n))<Entropy(IMFB,Lwindow(m,n))]]>对图像A和图像B的IMF分量进行融合，得到融合后的低频分量；其中：IMF_F,L(m,n)表示融合后的低频分量；是图像A的低频分量IMF_A,L(m,n)中以m行n列像素为中心的区域窗口内的信息熵；是图像B的低频分量IMF_B,L(m,n)中以m行n列像素为中心的区域窗口内的信息熵，

所述高频分量融合规则：

利用相关系数

corr(IMFA,Hwindow(m,n),IMFB,Hwindow(m,n))=Σm,n[(IMFA,Hwindow(m,n)-IMFA,Hwindow(m,n)‾)×((IMFB,Hwindow(m,n)-IMFB,Hwindow(m,n)‾)]Σm,n[(IMFA,Hwindow(m,n)-IMFA,Hwindow(m,n)‾)2]Σi,j[(IMFB,Hwindow(m,n)-IMFB,Hwindow(m,n)‾)2]]]>

判别图像A高频分量IMF_A,H(m,n)和图像B高频分量IMF_B,H(m,n)的相应高频子带窗口邻域内的相关性；当时表明两个区域的相关性差，其中T为阈值；并按IMFF,H(m,n)=IMFA,H(m,n)SFA,Hwindow(m,n)≥SFB,Hwindow(m,n)IMFB,H(m,n)SFA,Hwindow(m,n)<SFB,Hwindow(m,n)]]>对图像A和图像B的高频分量进行融合处理，得到融合后的高频分量IMF_F,H(m,n)，其中是处理后的图像A的高频分量IMF_A,H中以m行n列像素为中心的区域窗口内的空间频率，是处理后的图像B的高频分量IMF_B,H中以m行n列像素为中心的区域窗口内的空间频率；若时，根据区域的空间频率大小确定加权系数分别为：w=0.5+0.5×1-corr(IMFA,Hwindow(m,n),IMFB,Hwindow(m,n))1-TSFA,Hwindow(m,n)>SFB,Hwindow(m,n)w=0.5-0.5×1-corr(IMFA,Hwindow(m,n),IMFB,Hwindow(m,n))1-TSFA,Hwindow(m,n)≤SFB,Hwindow(m,n)]]>

故融和图像的高频信息为IMF_F,H(m,n)＝w×IMF_A,H(m,n)+(1-w)×IMF_B,H(m,n)；

步骤5：对步骤4得到的融合后的IMF分量IMF_F,L(m,n)和IMF_F,H(m,n)进行EMD逆变换，得到融合图像。