基于NSCT的多聚焦图像融合算法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号(10)申请公布号CNCN103632353103632353A(43)申请公布日2014.03.12(21)申请号201210305603.3(22)申请日2012.08.24(71)申请人西安元朔科技有限公司地址710077陕西省西安市高新区锦业路69号创业研发园C区1号瞪羚谷E座503室(72)发明人但春林封长林(51)Int.Cl.G06T5/50(2006.01)权权利要求书2页利要求书2页说明书5页说明书5页附图2页附图2页(54)发明名称基于NSCT的多聚焦图像融合算法(57)摘要本发明公开了一种基于NSCT的图像融合算法，该方法的处理过程是：首先非采样的轮廓波变换(NSCT)对经过均值滤波器滤波(mean)的源图像进行分解，再分别采用平均值和绝对值最大为低频子带和高频方向子带的融合规则，以NSCT的逆运算(inverseNSCT)进行图像重构，获得初始融合图像；其次，应用均方根误差(RMSE)提取初始融合图像的融合区域(fusionregion)；根据融合区域的特性设计融合规则，最后采用inverseNSCT进行图像融合。实验结果表明，本发明方法是非常有效，并且融合后的图像符合人眼视觉效果。CN103632353ACN103625ACN103632353A权利要求书1/2页1.一种基于NSCT的多聚焦图像融合算法，包括如下过程：步骤1：初始融合图像的获得1)采用均值滤波器对两幅多聚焦图像A和B进行滤波处理，得到滤波后的多聚焦图像，分别记为A′和B′；2)应用NSCT变换将图像A′和图像B′进行分解，得到两图像分解后的低频系数为和高频系数为和其中l为尺度分解数，k为方向分解级数；3)采用均值法和绝对值最大法分别作为低频子带和高频方向子带的融合规则，具体见公式(1)和公式(2)；4)根据上述融合策略，经过NSCT逆变换重构图像，即为初始融合图像F。步骤2：融合区域的提取应用公式(3)计算初始融合图像F和多聚焦图像A′、B′之间在像素点(x，y)的均方根误差，分别记为RMSEA′(i，j)和RMSEB′(i，j)，由此应用公式(4)确定标识矩阵Z，提取初始融合图像的融合区域。(3)步骤3：建立融合规则本发明在窗口大小为m1×n1上确定低频子带和高频方向子带的融合规则。具体如下：1)低频子带的融合策略其中，2CN103632353A权利要求书2/2页2)高频方向子带的融合策略其中，步骤4：融合图像的获得根据上述融合策略，经过NSCT逆变换得到图像即为融合图像。3CN103632353A说明书1/5页基于NSCT的多聚焦图像融合算法技术领域[0001]本发明属于图像处理领域，尤其涉及一种基于NSCT的多聚焦图像融合算法。背景技术[0002]图像融合是指将众多源信道所采集到的关于同一目标的图像数据经过图像处理和计算机技术等，最大限度的提取各自信道中的有利信息，最后综合成高质量的图像，以提高图像信息的利用率、改善计算机解译度和可靠性、提升原始图像的空间分辨率和光谱分辨率，利于监测。[0003]多聚焦图像融合算法是最常用的图像融合算法之一，传统的多聚焦图像融合算法主要分为两种。一种是基于空间领域的方法，该方法通过从众多图像中选择较为清晰的像素点或区域进行融合；另一种则是在频域中进行图像融合。这两种方法中，最为简单的算法为基于空间领域的融合算法，但是这类方法产生的不良后果比较多，例如对比度不明显、块状效应等等。另一种方法是在频域中进行图像融合，该方法主要通过多尺度变换见空域变为频域，并通过相应的融合规则，最后通过拟多尺度变换重构图像，这个重构图像即为融合图像。基于多尺度变换的图像融合算法越来越受到广泛关注，研究表明该类方法的融合效果更接近人类视觉感知。常用的多尺度变换包括拉普拉斯金字塔、曲波、轮廓波、非采样轮廓波变换等，通过大量的实验结果表明，基于多尺度变换的图像融合效果较基于空间的图像融合算法的效果更好，这是因为通过多尺度变换可将图像分为系数子带，而不是空间区域中的像素或块，能够较好地保留图像的细节。[0004]基于多尺度图像融合算法需要考虑两种问题，其一是如何选择多尺度变换算法，其二是系数子带的融合规则。大量的研究结果表明：不同的多尺度变换和融合规则直接影响到多聚焦图像融合的效果。[0005]文献[1]表明，在多尺度分析中，轮廓波变换(CT)能够最优化的表示图像轮廓，并且已广泛应用到图像融合中，但是文献[2]研究表明轮廓波变换不具备平移不变性，这将导致伪吉布斯现象出现。为了弥补轮廓波的不足，文献[2]提出了非采样的轮廓波变换(NSCT)，文献[3]将NSCT成功应用到图像融合中，并取得了较好的效果。但文献[3]中以平均值和绝对值最大为融合规则，导致了图