基于稀疏字典学习和剪切波的图像融合方法与流程

本发明涉及图像融合
技术领域：
，更具体的说是涉及一种应用于军事侦察、医疗诊断和遥感等领域的基于稀疏字典学习和剪切波的图像融合方法。
背景技术：
：传感器技术的发展呈现出成像机理多样化、工作环境复杂化、工作波段多元化和功能模块工作的协同化等特点，极大地拓展了人类认识世界的深度和广度。前视红外、激光成像雷达、合成孔径雷达、光学相干断层成像、光学多孔径成像等新型传感器实现了对自然世界的多元模式表达，被广泛应用于军事和民用领域。图像融合通过挖掘反映同一场景不同角度或不同频率波段的有用信息，利用特定的数据处理方法提取其互补分量进行综合表示，实现对场景内容更为全面、精确的表达。多源图像融合的优点体现在：融合图像比源图像分辨率更高、信息集聚度更强，有利于后续的检测识别或目标跟踪；融合图像能够抑制噪声影响，减少决策过程中的不确定因素；融合图像实现了多源信息更高效的表达，能提供更丰富的有用信息；单一传感器性能提升代价高昂，多源图像融合技术则方便经济，融合系统具有良好的应用前景。在图像融合文献中，多尺度变换和稀疏表示是两种最广泛使用的信号/图像表示理论。多尺度变换分解工具包括金字塔变换、小波变换、曲波变换、轮廓波变换、非下采样轮廓波变换，然而上述方法的稀疏表示能力较差，分解后的图像方向子带信息有限。而稀疏字典学习融合方法可以精确的拟合数据，并且分解系数越稀疏越能反映信号内在特征，可以通过稀疏表示来提高低频子带稀疏度，从而提高融合效果，但稀疏表示不能精确表示高频信息，即不能多尺度多方向的分析数据，如果高频采用稀疏表示进行融合，则会存在高频信息丢失的问题。因此，如何提供一种能够准确捕获图像的边缘和纹理信息以提高图像融合质量的基于稀疏字典学习和剪切波的图像融合方法是本领域技术人员亟需解决的问题。技术实现要素：有鉴于此，本发明采用非下采样剪切波变换将待融合图像分解成低频子带系数和高频子带系数；低频子带系数反映了图像的基础轮廓，同时对稀疏度较差的低频子带系数进行字典学习，高频子带系数反映了图像的边缘和纹理信息，采用同一尺度下方向子带绝对值和取最大的方法融合以捕获源图像中的显著特征，最终提高融合效果。为了实现上述目的，本发明采用如下技术方案：一种基于稀疏字典学习和剪切波的图像融合方法，包括如下步骤：步骤1，对已经配准好的两幅源图像{ia,ib}进行剪切波变换，以获得对应的低频子带系数{la,lb}和高频子带系数{ha,hb}，图像大小为m×n；步骤2：对低频子带系数进行稀疏字典学习，并进行融合得到低频融合系数lf；步骤3，对高频子带系数采用同一尺度下方向子带绝对值和取最大值的方法进行融合，得到高频融合系数hf；步骤4，对所述低频融合系数lf和高频融合系数hf进行剪切波逆变换，得到融合图像。剪切波具有更加优异稀疏表示能力，能够准确地捕获图像中的细节信息、同时去除了伪影现象、有良好的重构效果。图像融合中，子带系数融合规则包括绝对值取大和加权平均，绝对值取大是选择对应子带系数大的作为融合系数，可以提高融合图像的对比度，加权平均根据不同权重因子将子带系数加权平均，可以较好保留待融合图像有用信息、轮廓和形状。优选的，所述步骤2具体包括：(1)滑窗取块，按照从左上到右下的顺序，将所述低频子带系数{la,lb}以步长l，大小为的滑动窗口进行分块提取，得到大小为的分块其中，对于每个位置i，再将所有分块转置排列为列向量然后通过将每个所述列向量的平均值归一化为零以获得构成低频子带稀疏字典学习样本集；其中，表示全一值n×1向量，和分别为和的平均值；(2)学习低频子带稀疏字典以构成稀疏表示的基原子；(3)低频子带系数融合，的系数表示为用最大范数l1规则合并和以获得融合的稀疏矢量和的融合结果由下式计算其中，为合并的平均值，由下式得出再将融合向量转置即得到融合分块将插入到lf的原始位置；当分块之间有重叠的部分时，对应lf中的每个像素上的值取平均。。优选的，首先从所述低频子带稀疏字典学习样本集中随机选取p个样本构成训练集，然后随机初始化一个大小为n×m的字典d0，其中n＜m，n为字典原子的大小，m为每个字典的原子数；使用正交匹配追踪算法来计算稀疏系数向量，即并不断迭代更新低频子带稀疏字典d，其中，d是学习的低频子带稀疏字典，αa、αb分别为对应的全部系数，||αa||0、||αb||0分别为αa、αb中所包含的非零分量个数，ε为容错误差，并保存所述低频子带稀疏字典d。优选的，所述源图像包括但不限于配准好的多聚焦图像、中红外和远红外图片。优选的，所述字典d0中，n＝64，m＝256。经由上述的技术方案可知，本发明提供一种基于稀疏字典学习和剪切波的图像融合方法，具体有益效果如下：一、本发明采用剪切波变换进行图像多尺度分解，相比传统小波变换、曲波变换和轮廓波变换，可以将图像分解为更多的方向子带，为后续图像融合处理提供更多方向子带信息，能够准确捕获图像中的边缘信息，且不会引入伪吉布斯效应。二、本发明融合方法低频子带采用稀疏字典学习方法进行融合，提高了低频子带稀疏度，高频子带采用同一尺度下方向子带绝对值最大方法融合，最终提高了图像融合质量。附图说明为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。图1附图为本发明基于稀疏字典学习和剪切波的图像融合方法的原理框图；图2附图为本发明两幅源图像的左右聚焦图像；图3附图为非下采样轮廓波变换和稀疏字典学习融合结果与本发明剪切波和稀疏字典学习融合结果的对比图。具体实施方式下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。参照说明书附图1，本发明的基于稀疏字典学习和剪切波的图像融合方法包括如下步骤：步骤1，参见说明书附图2为两幅源图像，图2(a)为左聚焦图像，图2(b)为右聚焦图片，对已经配准好的两幅源图像{ia,ib}进行剪切波变换，以获得对应的低频子带系数{la,lb}和高频子带系数{ha,hb}，图像大小为m×n；步骤2：对低频子带系数进行稀疏字典学习，并进行融合得到低频融合系数lf；(1)滑窗取块，按照从左上到右下的顺序，将低频子带系数{la,lb}以步长l，大小为的滑动窗口进行分块提取，得到大小为的分块其中，对于每个位置i，再将所有分块转置排列为列向量然后通过将每个列向量的平均值归一化为零以获得构成低频子带稀疏字典学习样本集；其中，表示全一值n×1向量，和分别为和的平均值；(2)学习低频子带稀疏字典以构成稀疏表示的基原子；首先从低频子带稀疏字典学习样本集中随机选取p个样本构成训练集，然后随机初始化一个大小为n×m的字典d0，其中n＜m，n为字典原子的大小，m为每个字典的原子数；字典学习是将高维信号用较少的字典原子进行表示，对于有限维空间，如果列向量ai(i＝1,2,3...,n)，a＝{a1,a2,...,an}是此有限维空间的一组基，那么n维列向量b可以表示为其中，x＝(x1,x2,...,x3)t为向量b的表示系数，如果选取a为a＝{a1,a2,...,am}m＜n，那么求得x的解通常是无穷个的。在解空间中，所求的系数向量x中非零分量的个数最少，就是寻找一个最稀疏的表示，记号||x||0表示为向量x中所含非零分量的个数，要使得||x||0最小，数学表示为s.tb＝ax，s.t表示在约束条件下。而在实际中，约束条件通常为近似，即为s.t||b-ax||≤ε。在稀疏字典学习中，即为其中，yi为信号，为信号所对应的表示系数，ε为容错误差，d为字典。使用正交匹配追踪算法来计算稀疏系数向量，即并不断迭代更新低频子带稀疏字典d，其中，d是学习的低频子带稀疏字典，αa、αb分别为对应的全部系数，||αa||0、||αb||0分别为αa、αb中所包含的非零分量个数，ε为容错误差，并保存低频子带稀疏字典d。在一个实施例中，字典d0的n＝64，m＝256。(3)低频子带系数融合，的表示系数为用最大范数l1规则合并和以获得融合的稀疏矢量和的融合结果由下式计算其中，为合并的平均值，由下式得出再将融合向量转置即得到融合分块将插入到lf的原始位置；当分块之间有重叠的部分时，对应lf中的每个像素上的值取平均。因为相邻的分块步长为1，分块边长为所以之间有重叠的部分，所以在重构时，融合后的分块之间每重叠一次应该在原来基础上做一次平均。步骤3，对高频子带系数采用同一尺度下方向子带绝对值和取最大值的方法进行融合，得到高频融合系数hf；高频子带系数{ha,hb}，ha类似低频子带系数，与原图像等大的矩阵构成，绝对值取大是指a，b对应矩阵同一位置的值取绝对值然后比较大小，绝对值大的值作为融合后的高频值，即：得到高频融合系数hf；步骤4，对低频融合系数lf和高频融合系数hf进行剪切波逆变换，得到融合图像。参见说明书附图3所示，3(a)为非下采样轮廓波变换和稀疏字典学习融合方法结果，3(b)为剪切波和稀疏字典学习融合方法结果。为了定量评价不同融合方法用于多聚焦图片融合性能，本发明采用信息熵(en)、平均梯度(avg)、边缘强度(ei)以及互信息(mi)等评价参数，用于对比的融合方法采用非下采样轮廓波变换和稀疏字典学习融合方法，以及剪切波和稀疏字典学习融合方法记过如下表1所示，参数越大表示效果越好，图像边缘纹理信息越丰富、融合图像质量越好。表1非下采样轮廓波变换和稀疏字典学习融合结果与本发明剪切波和稀疏字典学习融合结果参数比对表nsct_srnsst_sren(信息熵)7.3292282377.34036016avg(平均梯度)5.6769979665.683472982ei(边缘强度)60.4243593760.77937653fd(图像清晰度)6.5904721776.641491989mi(互信息)5.3197326975.323883589如表1中数据所示，第一列为评价参数，第二列为非下采样轮廓波变换和稀疏字典学习融合方法，第三列为本发明剪切波和稀疏字典学习融合方法。可以看出信息熵、边缘强度、图像清晰度等评价参数有较大的提升，表明本发明方法能够增强图像的细节表现能力，从源图像提取了更多子带信息并加以融合，所以具有更好的融合效果。为了进一步优化上述技术方案，源图像包括但不限于配准好的多聚焦图像、中红外和远红外图片。本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。当前第1页12