[发明专利]一种基于NSST与聚焦区域检测的多聚焦图像融合方法

摘要

本发明公开一种基于NSST和聚焦区域检测的多聚焦图像融合方法，首先对待融合的多聚焦图像分别利用NSST及相应的融合规则，得到初始融合图像；将待融合的源图像分别与初始融合图像进行比较，对于源图像中与初始融合图像相似性较大的像素点，我们认为它是来自源图像中的严格聚焦区域的，根据这一理论我们可以确定源图像中的聚焦区域；最后，在已经确定的聚焦区域上利用改进的脉冲耦合神经网络获得最终的融合图像。本发明可以充分的提取各源图像中的有用信息并注入到融合图像中，有效地克服了基于变换域的融合算法不能选取出所有位于聚焦区域内的系数的缺陷，并且获得的融合图像在视觉效果及客观指标上均优于经典的多尺度融合算法。

主权项

一种基于NSST与聚焦区域检测的多聚焦图像融合方法，其特征在于包括下述步骤：步骤1：对已配准的待融合图像A、B分别进行NSST分解，对得到的低频子带系数和高频子带系数分别应用不同的融合规则，从而得到一幅初始融合图像；步骤2：通过均方根误差方法测定源图像与初始融合图像之间的相似性来识别出初始聚焦区域，均方根误差计算公式分别为：${\mathrm{RMSE}}_A(i,j)={\left(\frac{\underset{m=-M}{\overset{M}{\mathrm{\Σ}}}\underset{n=-N}{\overset{N}{\mathrm{\Σ}}}{(I_F(i+m,j+n)-I_A(i+m,j+n))}^2}{(2M+1)\×(2N+1)}\right)}^{1/2}$${\mathrm{RMSE}}_B(i,j)={\left(\frac{\underset{m=-M}{\overset{M}{\mathrm{\Σ}}}\underset{n=-N}{\overset{N}{\mathrm{\Σ}}}{(I_F(i+m,j+n)-I_B(i+m,j+n))}^2}{(2M+1)\×(2N+1)}\right)}^{1/2}$其中，I_A(i,j),I_B(i,j)和I_F(i,j)分别表示图像A、B和初始融合图像F在点(i,j)处的像素值；通过比较RMSE_A和RMSE_B的大小，可以经由下式得到一个二值化图像，$Z(i,j)=\left\{\begin{array}{cc}1& \mathrm{if}{\mathrm{RMSE}}_A(i,j)\≤{\mathrm{RMSE}}_B(i,j)\\ 0& \mathrm{if}{\mathrm{RMSE}}_A(i,j)>{\mathrm{RMSE}}_B(i,j)\end{array}\right.$二值图像中，Z(i,j)＝1表示在图像A中，点(i,j)处的像素位于聚焦区域，否则，图像B中在该点处的像素属于聚焦区域；步骤3：利用小结构元素的形态学开闭运算对步骤2中得到的二值图像进行修正，消除初始聚焦区域中的突起、断点和裂缝；步骤4：通过改进的脉冲耦合神经网络(IPCNN)和步骤3中已修正的聚焦区域获得最终的融合图像，融合规则如下：$I_{\mathrm{FF}}(i,j)=\left\{\begin{array}{cc}{I}_A(i,j)& \mathrm{if}{Z}^{\′}(i,j)=1\mathrm{and\; W}(i,j)=m\×n\\ I_B(i,j)& \mathrm{if}{Z}^{\′}(i,j)=0\mathrm{and\; W}(i,j)=0\\ I_A(i,j)& \mathrm{if}0\<W(i,j)\<m\×\mathrm{n\; and}{T}_A(i,j)\≤T_B(i,j)\\ I_B(i,j)& \mathrm{if}0\<W(i,j)\<m\×\mathrm{n\; and}{T}_A(i,j)>T_B(i,j)\end{array}\right.$其中，I_A(i,j),I_B(i,j)和I_FF(i,j)分别表示图像A、B和最终融合图像F在点(i,j)处的像素值，Z′为修正的二值图像，W(i,j)表示Z′中以点(i,j)为中心、大小为m×n的局部区域像素和，T_A(i,j)和T_B(i,j)分别表示源图像A和B输入到改进的脉冲耦合神经网络后得到的赋时矩阵在点(i,j)处的值；上述的邻域(2M+1)×(2N+1)和m×n的大小分别设为5×5和3×3。