【智能系统】多小波和NSDFB组合域递归滤波多聚焦图像融合编辑部

第11卷第2期 智能系统学报 Vol.11No.2 2016年4月 CAAI Transactions on Intelligent Systems Apr.2016 D0I:10.11992/is.201509017 网络出版地址：http://www.cnki.net/kcms/detail/23.1538.TP.20160314.1431.004.html 多小波和NSDFB组合域递归滤波多聚焦图像融合 任晓霞，孙秀明，耿鹏2，苏醒2 (1.张家口学院理学系，河北张家口075000：2.石家庄铁道大学信息科学与技术学院，河北石家庄050043) 摘要：多小波能同时满足正交、紧支、对称等对信号处理十分重要的特性，结合多小波变换的多尺度特点和非子采 样方向滤波器组变换的多方向性，提出了一种新的基于多小波和非子采样方向滤波器组的多尺度多方向变换。对 于高频系数首先计算其修改空间频率，然后利用域变换递归滤波进行滤波的融合规则：低频系数采用了修改拉普拉 斯能量和的(SML)融合规则。通过与其他融合方法进行实验对比，实验结果表明：本文提出的融合方法能够更加有 效地选择源图像中的聚焦良好区域，并且引入的伪影信息较少：此外，与其他融合方法相比本文方法的客观评价结 果也是最好的。 关键词：图像处理；图像融合：递归滤波；改进空间频率：多小波 中图分类号：TP391文献标志码：A文章编号：1673-4785(2016)02-0241-08 中文引用格式：任晓霞，孙秀明，耿鹏，等.多小波和NSDFB组合域递归滤波多聚焦图像融合[J].智能系统学报，2016,11(2)： 241-248. 英文引用格式：REN Xiaoxia,SUN Xiuming,GENG Peng,etal.Multifocus image fusion using a recursive filter in the combined domain of multiwavelets and NSDFB[].CAAI transactions on intelligent systems,2016,11(2):241-248. Multifocus image fusion using a recursive filter in the combined domain of multiwavelets and NSDFB REN Xiaoxia',SUN Xiuming',GENG Peng2,SU Xing? (1.Science department,Zhangjiakou University,Zhangjiakou 075000,China;2.School of Information Science and Technology,Shi- jiazhuang Tiedao University,Shijiazhuang 050043,China) Abstract:The multiwavelet transform has properties of orthogonality,tight frame,and symmetry,which are vital to signal processing.In this study,a new transform,called as MNSDFB,is proposed by combining the multi-wavelet and nonsubsampled directional filter banks.The domain transform recursive filter is adopted to fuse the filters after the spatial frequency of the high frequency coefficient is calculated.The modified sum-modified-Laplacian (SML) is employed in the low pass subbands as a focus measurement to select fused coefficients.The presented fusion rule in the high pass subband can distinguish the focused regions from the blurred regions.The proposed fusion method was compared with three other fusion methods.The experimental results demonstrate that the proposed fusion meth- od can select the focused regions while introducing few artifacts into the final merged image.Furthermore,its objec- tive criteria,such as MI and QAB/F,are better than those of the other three methods. Keywords:image processing;image fusion;recursive filter;modified spatial frequency;multiwavelett 图像融合是信息融合的一个重要分支，其目的 一幅新的具有更高可信度清晰度以及可理解性的图 是将不同传感器获取的同一目标的互补信息合并为 像[山。因此，经过图像融合技术得到的合成图像可 以更全面、更精确地描述所研究的对象，为进一步图 收稿日期：2015-09-09.网络出版日期：2016-03-14. 像处理和分析提供高质量数据。因此，图像融合技 基金项目：河北省自然科学基金项目(F2013210094,F2013210109). 通信作者：耿鹏.E-mail:Gengpeng(@stdu.cdu.cn. 术在军事、医学、遥感、计算机视觉等领域得到了广

第 １１ 卷第 ２ 期 智 能 系 统 学 报 Ｖｏｌ．１１№．２ ２０１６ 年 ４ 月 ＣＡＡＩ Ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ Ｉｎｔｅｌｌｉｇｅｎｔ Ｓｙｓｔｅｍｓ Ａｐｒ． ２０１６ ＤＯＩ：１０．１１９９２ ／ ｔｉｓ．２０１５０９０１７ 网络出版地址：ｈｔｔｐ： ／ ／ ｗｗｗ．ｃｎｋｉ．ｎｅｔ ／ ｋｃｍｓ／ ｄｅｔａｉｌ ／ ２３．１５３８．ＴＰ．２０１６０３１４．１４３１．００４．ｈｔｍｌ 多小波和 ＮＳＤＦＢ 组合域递归滤波多聚焦图像融合 任晓霞１ ，孙秀明１ ，耿鹏２ ，苏醒２ （１．张家口学院 理学系， 河北 张家口 ０７５０００； ２．石家庄铁道大学 信息科学与技术学院， 河北 石家庄 ０５００４３） 摘 要：多小波能同时满足正交、紧支、对称等对信号处理十分重要的特性，结合多小波变换的多尺度特点和非子采 样方向滤波器组变换的多方向性，提出了一种新的基于多小波和非子采样方向滤波器组的多尺度多方向变换。 对 于高频系数首先计算其修改空间频率，然后利用域变换递归滤波进行滤波的融合规则；低频系数采用了修改拉普拉 斯能量和的（ＳＭＬ）融合规则。 通过与其他融合方法进行实验对比，实验结果表明：本文提出的融合方法能够更加有 效地选择源图像中的聚焦良好区域，并且引入的伪影信息较少；此外，与其他融合方法相比本文方法的客观评价结 果也是最好的。 关键词：图像处理；图像融合；递归滤波；改进空间频率；多小波 中图分类号： ＴＰ３９１ 文献标志码：Ａ 文章编号：１６７３⁃４７８５（２０１６）０２⁃０２４１⁃０８ 中文引用格式：任晓霞，孙秀明，耿鹏，等． 多小波和 ＮＳＤＦＢ 组合域递归滤波多聚焦图像融合［ Ｊ］． 智能系统学报， ２０１６， １１（ ２）： ２４１⁃２４８． 英文引用格式：ＲＥＮ Ｘｉａｏｘｉａ， ＳＵＮ Ｘｉｕｍｉｎｇ， ＧＥＮＧ Ｐｅｎｇ， ｅｔ ａｌ． Ｍｕｌｔｉｆｏｃｕｓ ｉｍａｇｅ ｆｕｓｉｏｎ ｕｓｉｎｇ ａ ｒｅｃｕｒｓｉｖｅ ｆｉｌｔｅｒ ｉｎ ｔｈｅ ｃｏｍｂｉｎｅｄ ｄｏｍａｉｎ ｏｆ ｍｕｌｔｉｗａｖｅｌｅｔｓ ａｎｄ ＮＳＤＦＢ［Ｊ］． ＣＡＡＩ ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ ｉｎｔｅｌｌｉｇｅｎｔ ｓｙｓｔｅｍｓ， ２０１６， １１（２）： ２４１⁃２４８． Ｍｕｌｔｉｆｏｃｕｓ ｉｍａｇｅ ｆｕｓｉｏｎ ｕｓｉｎｇ ａ ｒｅｃｕｒｓｉｖｅ ｆｉｌｔｅｒ ｉｎ ｔｈｅ ｃｏｍｂｉｎｅｄ ｄｏｍａｉｎ ｏｆ ｍｕｌｔｉｗａｖｅｌｅｔｓ ａｎｄ ＮＳＤＦＢ ＲＥＮ Ｘｉａｏｘｉａ １ ， ＳＵＮ Ｘｉｕｍｉｎｇ １ ， ＧＥＮＧ Ｐｅｎｇ ２ ， ＳＵ Ｘｉｎｇ ２ （１． Ｓｃｉｅｎｃｅ ｄｅｐａｒｔｍｅｎｔ ， Ｚｈａｎｇｊｉａｋｏｕ Ｕｎｉｖｅｒｓｉｔｙ， Ｚｈａｎｇｊｉａｋｏｕ ０７５０００， Ｃｈｉｎａ； ２． Ｓｃｈｏｏｌ ｏｆ Ｉｎｆｏｒｍａｔｉｏｎ Ｓｃｉｅｎｃｅ ａｎｄ Ｔｅｃｈｎｏｌｏｇｙ， Ｓｈｉ⁃ ｊｉａｚｈｕａｎｇ Ｔｉｅｄａｏ Ｕｎｉｖｅｒｓｉｔｙ， Ｓｈｉｊｉａｚｈｕａｎｇ ０５００４３， Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｔｈｅ ｍｕｌｔｉｗａｖｅｌｅｔ ｔｒａｎｓｆｏｒｍ ｈａｓ ｐｒｏｐｅｒｔｉｅｓ ｏｆ ｏｒｔｈｏｇｏｎａｌｉｔｙ， ｔｉｇｈｔ ｆｒａｍｅ， ａｎｄ ｓｙｍｍｅｔｒｙ， ｗｈｉｃｈ ａｒｅ ｖｉｔａｌ ｔｏ ｓｉｇｎａｌ ｐｒｏｃｅｓｓｉｎｇ． Ｉｎ ｔｈｉｓ ｓｔｕｄｙ， ａ ｎｅｗ ｔｒａｎｓｆｏｒｍ， ｃａｌｌｅｄ ａｓ ＭＮＳＤＦＢ， ｉｓ ｐｒｏｐｏｓｅｄ ｂｙ ｃｏｍｂｉｎｉｎｇ ｔｈｅ ｍｕｌｔｉ⁃ｗａｖｅｌｅｔ ａｎｄ ｎｏｎｓｕｂｓａｍｐｌｅｄ ｄｉｒｅｃｔｉｏｎａｌ ｆｉｌｔｅｒ ｂａｎｋｓ． Ｔｈｅ ｄｏｍａｉｎ ｔｒａｎｓｆｏｒｍ ｒｅｃｕｒｓｉｖｅ ｆｉｌｔｅｒ ｉｓ ａｄｏｐｔｅｄ ｔｏ ｆｕｓｅ ｔｈｅ ｆｉｌｔｅｒｓ ａｆｔｅｒ ｔｈｅ ｓｐａｔｉａｌ ｆｒｅｑｕｅｎｃｙ ｏｆ ｔｈｅ ｈｉｇｈ ｆｒｅｑｕｅｎｃｙ ｃｏｅｆｆｉｃｉｅｎｔ ｉｓ ｃａｌｃｕｌａｔｅｄ． Ｔｈｅ ｍｏｄｉｆｉｅｄ ｓｕｍ⁃ｍｏｄｉｆｉｅｄ⁃Ｌａｐｌａｃｉａｎ （ＳＭＬ） ｉｓ ｅｍｐｌｏｙｅｄ ｉｎ ｔｈｅ ｌｏｗ ｐａｓｓ ｓｕｂｂａｎｄｓ ａｓ ａ ｆｏｃｕｓ ｍｅａｓｕｒｅｍｅｎｔ ｔｏ ｓｅｌｅｃｔ ｆｕｓｅｄ ｃｏｅｆｆｉｃｉｅｎｔｓ． Ｔｈｅ ｐｒｅｓｅｎｔｅｄ ｆｕｓｉｏｎ ｒｕｌｅ ｉｎ ｔｈｅ ｈｉｇｈ ｐａｓｓ ｓｕｂｂａｎｄ ｃａｎ ｄｉｓｔｉｎｇｕｉｓｈ ｔｈｅ ｆｏｃｕｓｅｄ ｒｅｇｉｏｎｓ ｆｒｏｍ ｔｈｅ ｂｌｕｒｒｅｄ ｒｅｇｉｏｎｓ． Ｔｈｅ ｐｒｏｐｏｓｅｄ ｆｕｓｉｏｎ ｍｅｔｈｏｄ ｗａｓ ｃｏｍｐａｒｅｄ ｗｉｔｈ ｔｈｒｅｅ ｏｔｈｅｒ ｆｕｓｉｏｎ ｍｅｔｈｏｄｓ． Ｔｈｅ ｅｘｐｅｒｉｍｅｎｔａｌ ｒｅｓｕｌｔｓ ｄｅｍｏｎｓｔｒａｔｅ ｔｈａｔ ｔｈｅ ｐｒｏｐｏｓｅｄ ｆｕｓｉｏｎ ｍｅｔｈ⁃ ｏｄ ｃａｎ ｓｅｌｅｃｔ ｔｈｅ ｆｏｃｕｓｅｄ ｒｅｇｉｏｎｓ ｗｈｉｌｅ ｉｎｔｒｏｄｕｃｉｎｇ ｆｅｗ ａｒｔｉｆａｃｔｓ ｉｎｔｏ ｔｈｅ ｆｉｎａｌ ｍｅｒｇｅｄ ｉｍａｇｅ． Ｆｕｒｔｈｅｒｍｏｒｅ， ｉｔｓ ｏｂｊｅｃ⁃ ｔｉｖｅ ｃｒｉｔｅｒｉａ， ｓｕｃｈ ａｓ ＭＩ ａｎｄ ＱＡＢ ／ Ｆ， ａｒｅ ｂｅｔｔｅｒ ｔｈａｎ ｔｈｏｓｅ ｏｆ ｔｈｅ ｏｔｈｅｒ ｔｈｒｅｅ ｍｅｔｈｏｄｓ． Ｋｅｙｗｏｒｄｓ： ｉｍａｇｅ ｐｒｏｃｅｓｓｉｎｇ； ｉｍａｇｅ ｆｕｓｉｏｎ； ｒｅｃｕｒｓｉｖｅ ｆｉｌｔｅｒ； ｍｏｄｉｆｉｅｄ ｓｐａｔｉａｌ ｆｒｅｑｕｅｎｃｙ； ｍｕｌｔｉｗａｖｅｌｅｔｔ 收稿日期：２０１５⁃０９⁃０９． 网络出版日期：２０１６⁃０３⁃１４． 基金项目：河北省自然科学基金项目（Ｆ２０１３２１００９４， Ｆ２０１３２１０１０９）． 通信作者：耿鹏．Ｅ⁃ｍａｉｌ：Ｇｅｎｇｐｅｎｇ＠ ｓｔｄｕ．ｅｄｕ．ｃｎ． 图像融合是信息融合的一个重要分支，其目的 是将不同传感器获取的同一目标的互补信息合并为 一幅新的具有更高可信度清晰度以及可理解性的图 像［１］ 。 因此，经过图像融合技术得到的合成图像可 以更全面、更精确地描述所研究的对象，为进一步图 像处理和分析提供高质量数据。 因此，图像融合技 术在军事、医学、遥感、计算机视觉等领域得到了广

.242. 智能系统学报 第11卷 泛的应用)。一般来说，图像融合方法可以分为两 通过等距变换映射到五维空间中，并保持曲线上 类：空域融合方法和变换域融合方法。空域方法根 每个点之间的距离不变，通过自适应地改变输入 据图像的清晰判别准则选择源图像中的像素，从而 信号，将滤波运算降到了一维上，使得域变换滤 形成一幅融合后的图像。空域方法一般有空间频率 波计算复杂度不受相关参数影响，因此，DTRF在 (SF)、拉普拉斯能量(EOL)、修改拉普拉斯能量和 保持图像边缘的同时，具有实时性的优点。本文 (SML)等方法[)。根据采用不同的变换，变换域方 在提出的MNSDFB变换基础上，提出了结合修改 法可以分为小波变换、Contourlet变换、非子采样 空间频率的域变换递归滤波的多聚焦图像融合 Contourlet变换以及Shearlet变换[a等方法。 方法。 苗启广[)提出了一种小波域基于区域局部能量 1 域变换递归滤波 的不同聚焦点图像融合方法，取得了很好的融合效 果。传统的小波变换虽能高效处理和分析一维分段 Gastal!山在域变换基础上提出了边缘保持平滑 连续信号，但由一维小波通过张量积而形成的二维小 递归滤波器(DTRF),它的基本思路可以理解为先 波基不满足各向异性的尺度关系，小波变换无法精确 降低输入信号的维数再进行低维滤波。在五维空间 地表述图像边缘信息，基于小波变换融合的图像容易 中，一个边缘保持的平滑滤波可以定义成一个空间 产生细节成分模糊现象[6)。多小波能同时满足正交、 不变的核H,它的响应会随着像素点间距离的变化 紧支、对称等对信号处理十分重要的特性，被广泛应 而变化，如果在某个低维空间中这些距离保持不变 用于图像融合领域[)]。王迎春)在多小波变换域对 那么核的响应也将保持不变，因而，达到了保留图像 低频和高频小波系数采用不同的融合方法，对低频系 边缘的目的)。DTRF首先将输入信号I变换到变 数采用取平均的方法，而对高频系数采用边缘梯度对 换域2.，，随后，变换域信号由递归滤波器(recursive 比的方法。该方法能够很好保存图像的边缘和细节 filter,RF)进行处理如下： 信息。但多小波是对传统的小波的扩展，仍不能解决 J[m]=(1-a)1[m]+aJ[m-1](1) 类似于小波变换无法精确地表述图像边缘方向信息 式中：J[m]是滤波后的结果，反馈系数a= 的问题。Cunhat9提出的非子采样Contourlet变换 exp(-√2/8，)∈[0,1]，δ是空间参数； (nonsubsampled contourlet transform,NSCT),由于同 I[m]=[xm]为输入的离散信号。b是变换域 时具有灵活的多分辨率和多方向性表示的优点，相对 2w相邻采样点xm和xm-1之间的距离，可以由 小波变换来说，能够更好地表示图像的边缘信息。米 b=ct(xm)一ct(xm-1)计算得到。对于变换域信 德伶[1O对NSCT分解后的低频子带和高频方向子带 号I(x),ct(u)可以由式(2)得到 分别以邻域梯度及合成邻域模值作为清晰度指标，采 (2） 用自适应选择法实现对多聚焦图像的融合处理。然 c()=小1+分lr(x)lk,ueD 0 而，NSCT变换中的大量的冗余使得这些融合方法具 式中：I'(x)是I(x)的导数，6和分别是DTRF的 有较高的时间复杂度。为了同时兼顾较小的分解冗 空间参数和范围参数。由式(2)可知：随着b的增 余和分辨率多方向的优点，本文充分利用NSCT中非 大，a°趋近于零，这逐渐使得式(2)的递归过程收 子采样方向滤波器组(nonsubsampled direction filter 敛，使得滤波后输出结果中同一侧图像边缘的像素 bank,NSDFB)能够更好捕捉图像的边缘的优点和 会取得相近的值，不同侧图像边缘的像素会有很大 Multiwavelet能同时满足正交、紧支、对称的特点，将 差别，从而图像的边缘可以更好地保留下来。 Multiwavelet和NSDFB结合构建MNSDFB(multiwave- let nonsubsampled direction filter bank)变换，并利用该 2融合方法 变换进行图像融合。 2.1组合多小波与NSDFB变换 除了变换方法外，融合规则对融合的效果也 非子采样Contourlet变换中首先对图像进行非 会产生重要影响。近年来，许多学者提出了各种 子采样塔式分解，然后进行非子采样方向分解(NS 边缘保持滤波方法：双边滤波、引导滤波、域变换 DFB),从而建立了图像的多尺度、多方向的表示。 滤波、Cost-volume filtering!]等，这些边缘保持滤 根据NSCT的想法，本文利用多小波将图像分解为 波能够很好地保持图像的边缘结构。其中 一个低频系数和3个高频系数，然后利用NSDFB将 Gastal]等提出的域变换递归滤波(domain trans- 4个高频系数进行NSDFB方向分解。低频系数可 form recursive filter,DTRF)将图像中的二维曲线 以再次进行多小波分解，并再次对产生的高频系数

泛的应用［２］ 。 一般来说，图像融合方法可以分为两 类：空域融合方法和变换域融合方法。 空域方法根 据图像的清晰判别准则选择源图像中的像素，从而 形成一幅融合后的图像。 空域方法一般有空间频率 （ＳＦ）、拉普拉斯能量（ＥＯＬ）、修改拉普拉斯能量和 （ＳＭＬ）等方法［３］ 。 根据采用不同的变换，变换域方 法可以分为小波变换、Ｃｏｎｔｏｕｒｌｅｔ 变换、非子采样 Ｃｏｎｔｏｕｒｌｅｔ 变换以及 Ｓｈｅａｒｌｅｔ 变换［４］等方法。 苗启广［５］提出了一种小波域基于区域局部能量 的不同聚焦点图像融合方法，取得了很好的融合效 果。 传统的小波变换虽能高效处理和分析一维分段 连续信号，但由一维小波通过张量积而形成的二维小 波基不满足各向异性的尺度关系，小波变换无法精确 地表述图像边缘信息，基于小波变换融合的图像容易 产生细节成分模糊现象［６］ 。 多小波能同时满足正交、 紧支、对称等对信号处理十分重要的特性，被广泛应 用于图像融合领域［７］ 。 王迎春［８］ 在多小波变换域对 低频和高频小波系数采用不同的融合方法，对低频系 数采用取平均的方法，而对高频系数采用边缘梯度对 比的方法。 该方法能够很好保存图像的边缘和细节 信息。 但多小波是对传统的小波的扩展，仍不能解决 类似于小波变换无法精确地表述图像边缘方向信息 的问题。 Ｃｕｎｈａ ［９］ 提出的非子采样 Ｃｏｎｔｏｕｒｌｅｔ 变换 （ｎｏｎｓｕｂｓａｍｐｌｅｄ ｃｏｎｔｏｕｒｌｅｔ ｔｒａｎｓｆｏｒｍ， ＮＳＣＴ），由于同 时具有灵活的多分辨率和多方向性表示的优点，相对 小波变换来说，能够更好地表示图像的边缘信息。 米 德伶［１０］对 ＮＳＣＴ 分解后的低频子带和高频方向子带 分别以邻域梯度及合成邻域模值作为清晰度指标，采 用自适应选择法实现对多聚焦图像的融合处理。 然 而，ＮＳＣＴ 变换中的大量的冗余使得这些融合方法具 有较高的时间复杂度。 为了同时兼顾较小的分解冗 余和分辨率多方向的优点，本文充分利用 ＮＳＣＴ 中非 子采样方向滤波器组（ｎｏｎｓｕｂｓａｍｐｌｅｄ ｄｉｒｅｃｔｉｏｎ ｆｉｌｔｅｒ ｂａｎｋ， ＮＳＤＦＢ）能够更好捕捉图像的边缘的优点和 Ｍｕｌｔｉｗａｖｅｌｅｔ 能同时满足正交、紧支、对称的特点，将 Ｍｕｌｔｉｗａｖｅｌｅｔ 和 ＮＳＤＦＢ 结合构建 ＭＮＳＤＦＢ（ｍｕｌｔｉｗａｖｅ⁃ ｌｅｔ ｎｏｎｓｕｂｓａｍｐｌｅｄ ｄｉｒｅｃｔｉｏｎ ｆｉｌｔｅｒ ｂａｎｋ）变换，并利用该 变换进行图像融合。 除了变换方法外，融合规则对融合的效果也 会产生重要影响。 近年来，许多学者提出了各种 边缘保持滤波方法：双边滤波、引导滤波、域变换 滤波、Ｃｏｓｔ⁃ｖｏｌｕｍｅ ｆｉｌｔｅｒｉｎｇ ［ １１］ 等，这些边缘保持滤 波 能 够 很 好 地 保 持 图 像 的 边 缘 结 构。 其 中 Ｇａｓｔａｌ ［ １１］等提出的域变换递归滤波（ ｄｏｍａｉｎ ｔｒａｎｓ⁃ ｆｏｒｍ ｒｅｃｕｒｓｉｖｅ ｆｉｌｔｅｒ， ＤＴＲＦ） 将图像中的二维曲线 通过等距变换映射到五维空间中，并保持曲线上 每个点之间的距离不变，通过自适应地改变输入 信号，将滤波运算降到了一维上，使得域变换滤 波计算复杂度不受相关参数影响，因此，ＤＴＲＦ 在 保持图像边缘的同时，具有实时性的优点。 本文 在提出的 ＭＮＳＤＦＢ 变换基础上，提出了结合修改 空间频率的域变换递归滤波的多聚焦图像融合 方法。 １ 域变换递归滤波 Ｇａｓｔａｌ ［１１］在域变换基础上提出了边缘保持平滑 递归滤波器（ＤＴＲＦ），它的基本思路可以理解为先 降低输入信号的维数再进行低维滤波。 在五维空间 中，一个边缘保持的平滑滤波可以定义成一个空间 不变的核 Ｈ ，它的响应会随着像素点间距离的变化 而变化，如果在某个低维空间中这些距离保持不变， 那么核的响应也将保持不变，因而，达到了保留图像 边缘的目的［１１］ 。 ＤＴＲＦ 首先将输入信号 Ｉ 变换到变 换域 Ωω ，随后，变换域信号由递归滤波器（ｒｅｃｕｒｓｉｖｅ ｆｉｌｔｅｒ， ＲＦ）进行处理如下： Ｊ［ｍ］ ＝ （１ － ａ ｂ ）Ｉ［ｍ］ ＋ ａ ｂ Ｊ［ｍ － １］ （１） 式中： Ｊ［ ｍ］ 是 滤 波 后 的 结 果， 反 馈 系 数 ａ ＝ ｅｘｐ（ － ２ ／ δｓ） ∈ ［ ０，１］ ， δｓ 是 空 间 参 数； Ｉ［ ｍ］ ＝ Ｉ［ ｘｍ ］ 为输入的离散信号。 ｂ 是变换域 ΩＷ 相邻采 样 点 ｘｍ 和 ｘｍ － １ 之 间 的 距 离，可 以 由 ｂ ＝ ｃｔ（ ｘｍ ） － ｃｔ（ ｘｍ － １ ） 计算得到。 对于变换域信 号 Ｉ（ ｘ） ， ｃｔ（ ｕ） 可以由式（ ２）得到 ｃｔ（ｕ） ＝ ∫ ｕ ０ １ ＋ δｓ δｒ Ｉ′（ｘ） ｄｘ，ｕ ∈ Ωｗ （２） 式中： Ｉ′（ｘ） 是 Ｉ（ｘ） 的导数， δｓ 和分别是 ＤＴＲＦ 的 空间参数和范围参数。 由式（２） 可知：随着 ｂ 的增 大， ａ ｂ 趋近于零，这逐渐使得式（２） 的递归过程收 敛，使得滤波后输出结果中同一侧图像边缘的像素 会取得相近的值，不同侧图像边缘的像素会有很大 差别，从而图像的边缘可以更好地保留下来。 ２ 融合方法 ２．１ 组合多小波与 ＮＳＤＦＢ 变换 非子采样 Ｃｏｎｔｏｕｒｌｅｔ 变换中首先对图像进行非 子采样塔式分解，然后进行非子采样方向分解（ＮＳ⁃ ＤＦＢ），从而建立了图像的多尺度、多方向的表示。 根据 ＮＳＣＴ 的想法，本文利用多小波将图像分解为 一个低频系数和 ３ 个高频系数，然后利用 ＮＳＤＦＢ 将 ４ 个高频系数进行 ＮＳＤＦＢ 方向分解。 低频系数可 以再次进行多小波分解，并再次对产生的高频系数， ·２４２· 智 能 系 统 学 报 第 １１ 卷

第2期 任晓霞，等：多小波和NSDFB组合域递归滤波多聚焦图像融合 ·243. 并再次对高频系数进行NSDFB方向分解，从而建立 图像进行较好地稀疏表示，同时能够更完美地表示 图像的MNSDFB变换。该变换不仅有正交、紧支撑 图像的边缘和细节信息。该变换的2尺度2方向的 的特点，同时还具有多尺度和多方向的特点，可以对 分解如图1所示。 NSDFB 高频 HH 方向子带 NSDFB HL 高频 LL- NSDEB 高频 方向子带 HH 方向子带 图像 多小波 NSDFB 高频 LL- NSDFB 高频 LH 方向子带 HL 方向子带 LL- NSDFB 高频 LH 方向子带 多小波 Lowpass coefficient 图1 MNSDFB分解过程 Fig.1 The Decomposition of MNSDFB 2.2低频系数融合规则 的图像特征细节信息，高频信息中绝对值较大的系 改进拉普拉斯能量和(sum-modified-Lapla- 数对应着一些突变，如图像的边缘、纹理等重要特征 cian,SML)是空间域内典型的清晰度评价指 信息。由于图像的特征信息不是由单一的像素所表 标【)，反映了图像的边缘特征，能恰当地表征图 征的，而是由这一区域的多个像素来表征和体 像的聚焦特性和清晰度。由于SML反映的是邻 现)。空间频率就是利用了相邻像素的特征来表 域窗口内像素的能量和，SML规则也能够体现局 示的行频率和列频率，来区分图像的清晰程度。然 部区域内多个像素的综合特征。因此，本文在低 而，空间频率没有利用相邻像素间的方向信息，不能 频系数中采用SML规则。其中，拉普拉斯(modi- 有效地区分图像的清晰区域和模糊区域。Das提出 fied-Laplacian,ML)定义为 的MSF14)综合利用了图像的方向信息，同时结合列 ML'(i,j)=|2C(i,》-C(i-1,)-C(i+1,)|+ 频率和行频率来作为图像的显著性特征准则。相对 |2c(i,)-C(i,j+1)-C(i,j+1)|(3) 于空间频率方法，MSF能更加有效地区分图像的清 式中：C(i,j)表示经过MNSDFB变换后，在l尺度 晰区域和模糊区域。图像空间频率的计算方法为 分解低频系数在(i,j)处的像素值。ML'(i,)表示 SF(m,n)= 其对应的ML值。改进的拉普拉斯能量和(SML)定 义为 2Au.-1-u.-1 sML(i,)=∑∑[ML'(i+oj+p)]2(4) (7) 0=-0p=-P 式中：Imn是图像I中m行n列像素值。图像的方 式中：0和p表示大小为(20+1)×(2P+1)的窗口的 向信息(direction information,DF)可以按照式(8)进 中心位置。设C(i,)和两幅源图像经过在1尺度 行计算： MNSDFB分解后得到的低频分解系数，C(i,j)为 DF(m,n)= 融合后的系数。由SML'(i,)得到的决策图如下： 1, SML(i,j)≥SML6(i,j) M(i,j)= (5) /日Σ∑[0a-1-a-i)2+0.-1-1-)2 0, SML(i,j)<SMLe(i,j） (8) 因此，可以根据决策图选择融合后的MNSDFB 因此，MSF可以表示为 低频系数C(i,j): MSF(m,n)=DF (m,n)2+SF (m,n)2 (9) (CA(ij), M(i,j)=1 Cr(i,j)= 由于DTRF可以忽略空间域图像的纹理信息， (6) \C(ij), M'(i,j）=0 使得图像中具有相似特征的像素连接成为相似像素 2.3高频系数融合规则 值区域，可以更好分聚焦区域和模糊。因此，在图像 MNSDFB变换中的高频细节信息中包含着丰富 的多尺度变换域可以利用变DTRF滤波器对MSF

并再次对高频系数进行 ＮＳＤＦＢ 方向分解，从而建立 图像的 ＭＮＳＤＦＢ 变换。 该变换不仅有正交、紧支撑 的特点，同时还具有多尺度和多方向的特点，可以对 图像进行较好地稀疏表示，同时能够更完美地表示 图像的边缘和细节信息。 该变换的 ２ 尺度 ２ 方向的 分解如图 １ 所示。 图 １ ＭＮＳＤＦＢ 分解过程 Ｆｉｇ．１ Ｔｈｅ Ｄｅｃｏｍｐｏｓｉｔｉｏｎ ｏｆ ＭＮＳＤＦＢ ２．２ 低频系数融合规则 改 进 拉 普 拉 斯 能 量 和 （ ｓｕｍ⁃ｍｏｄｉｆｉｅｄ⁃Ｌａｐｌａ⁃ ｃｉａｎ， ＳＭＬ） 是 空 间 域 内 典 型 的 清 晰 度 评 价 指 标［ １２］ ，反映了图像的边缘特征，能恰当地表征图 像的聚焦特性和清晰度。 由于 ＳＭＬ 反映的是邻 域窗口内像素的能量和，ＳＭＬ 规则也能够体现局 部区域内多个像素的综合特征。 因此，本文在低 频系数中采用 ＳＭＬ 规则。 其中，拉普拉斯（ ｍｏｄｉ⁃ ｆｉｅｄ⁃Ｌａｐｌａｃｉａｎ，ＭＬ）定义为 ＭＬ ｌ （ｉ，ｊ） ＝ ２Ｃ ｌ (ｉ，ｊ) － Ｃ ｌ (ｉ － １，ｊ) － Ｃ ｌ (ｉ ＋ １，ｊ) ＋ ２Ｃ ｌ (ｉ，ｊ) － Ｃ ｌ (ｉ，ｊ ＋ １) － Ｃ ｌ (ｉ，ｊ ＋ １) （３） 式中： Ｃ ｌ （ｉ，ｊ） 表示经过 ＭＮＳＤＦＢ 变换后，在 ｌ 尺度 分解低频系数在（ ｉ， ｊ）处的像素值。 ＭＬ ｌ （ｉ，ｊ） 表示 其对应的 ＭＬ 值。 改进的拉普拉斯能量和（ＳＭＬ）定 义为 ＳＭＬ ｌ （ｉ，ｊ） ＝ ∑ Ｏ ｏ ＝ －Ｏ∑ Ｐ ｐ ＝ －Ｐ ［ＭＬ ｌ （ｉ ＋ ｏ，ｊ ＋ ｐ）］ ２ （４） 式中：ｏ 和 ｐ 表示大小为（２Ｏ＋１） ×（２Ｐ＋１）的窗口的 中心位置。 设 Ｃ ｌ Ａ（ｉ，ｊ） 和两幅源图像经过在 ｌ 尺度 ＭＮＳＤＦＢ 分解后得到的低频分解系数， Ｃ ｌ Ｆ（ｉ，ｊ） 为 融合后的系数。 由 ＳＭＬ ｌ （ｉ，ｊ） 得到的决策图如下： Ｍｌ（ｉ，ｊ） ＝ １， ＳＭＬ ｌ Ａ（ｉ，ｊ） ≥ ＳＭＬ ｌ Ｂ（ｉ，ｊ） ０， ＳＭＬ ｌ Ａ（ｉ，ｊ） ＜ ＳＭＬ ｌ { Ｂ（ｉ，ｊ） （５） 因此，可以根据决策图选择融合后的 ＭＮＳＤＦＢ 低频系数 Ｃ ｌ Ｆ（ｉ，ｊ） ： Ｃ ｌ Ｆ（ｉ，ｊ） ＝ Ｃ ｌ Ａ（ｉ，ｊ）， Ｍ ｌ （ｉ，ｊ） ＝ １ Ｃ ｌ Ｂ（ｉ，ｊ）， Ｍ ｌ { （ｉ，ｊ） ＝ ０ （６） ２．３ 高频系数融合规则 ＭＮＳＤＦＢ 变换中的高频细节信息中包含着丰富 的图像特征细节信息，高频信息中绝对值较大的系 数对应着一些突变，如图像的边缘、纹理等重要特征 信息。 由于图像的特征信息不是由单一的像素所表 征的，而 是 由 这 一 区 域 的 多 个 像 素 来 表 征 和 体 现［１３］ 。 空间频率就是利用了相邻像素的特征来表 示的行频率和列频率，来区分图像的清晰程度。 然 而，空间频率没有利用相邻像素间的方向信息，不能 有效地区分图像的清晰区域和模糊区域。 Ｄａｓ 提出 的 ＭＳＦ ［１４］综合利用了图像的方向信息，同时结合列 频率和行频率来作为图像的显著性特征准则。 相对 于空间频率方法，ＭＳＦ 能更加有效地区分图像的清 晰区域和模糊区域。 图像空间频率的计算方法为 ＳＦ（ｍ，ｎ） ＝ １ ＭＮ∑ Ｍ ｍ ＝ １∑ Ｎ ｎ ＝ １ Ｉｍ，ｎ － Ｉｍ，ｎ－１ ( ) ２ － Ｉｍ，ｎ － Ｉｍ－１，ｎ ( ) ２ [ ] （７） 式中： Ｉｍ，ｎ 是图像 Ｉ 中 ｍ 行 ｎ 列像素值。 图像的方 向信息（ｄｉｒｅｃｔｉｏｎ ｉｎｆｏｒｍａｔｉｏｎ，ＤＦ）可以按照式（８）进 行计算： ＤＦ（ｍ，ｎ） ＝ １ ＭＮ∑ Ｍ ｍ ＝ １∑ Ｎ ｎ ＝ １ Ｉｍ，ｎ － Ｉｍ－１，ｎ－１ ( ) ２ ＋ Ｉｍ，ｎ － Ｉｍ－１，ｎ－１ ( ) ２ [ ] （８） 因此，ＭＳＦ 可以表示为 ＭＳＦ（ｍ，ｎ） ＝ ＤＦ （ｍ，ｎ） ２ ＋ ＳＦ （ｍ，ｎ） ２ （９） 由于 ＤＴＲＦ 可以忽略空间域图像的纹理信息， 使得图像中具有相似特征的像素连接成为相似像素 值区域，可以更好分聚焦区域和模糊。 因此，在图像 的多尺度变换域可以利用变 ＤＴＲＦ 滤波器对 ＭＳＦ 第 ２ 期 任晓霞，等：多小波和 ＮＳＤＦＢ 组合域递归滤波多聚焦图像融合 ·２４３·

.244 智能系统学报 第11卷 进行滤波，利用DTRF滤波器输出去生成更为接近 本文融合方法如图2所示并描述如下： 理想状态的融合决策图。利用式(9)计算所提出的 l)将源图像A和B进行l尺度的Multiwavelet MNSDFB高频系数的MSF值，并将其作为DTRF滤 分解，得到包含多小波低频系数G(i,)、 波器的输人，可以得到： G(i》和多小波高频系数G(i》、G哈(i》。 R(i,j)=RF(MSF(i,j)) (10) 2)对G(i,)和G(i,》的3个高频系数进行 由R(i,)得到的决策图为 NSDB分解得到k方向子带系数C.(i,)和 1, R(i,j》≥R.(i,i) M(i.j)= (11) C.(i,),k为NSDFB分解的方向。 0 R(i,》<R路(i,) 3)对低频系数和高频系数分别按照相应的融 根据决策图选择融合后的系数C,(i,》: 合规则选择系数得到融合后的系数C(i,)和 D,(i,j), M(i,j)=1 C(ij》。 C(i,j》= (12) D(ij). ,(i,j)=0 4)对融合后系数C(i,)和C(i,)进行反 NSDFB变换得到融合后的多小波低频系数 式中：C,(i,》、C.(i,)和C(i,》分别是源多 聚焦图像A、B和融合后图像F在k方向的高频分 G(i,)和高频系数G(i,)。 解系数。 5)对小波系数G(i,)、G(i,)进行反多小 2.4 算法步骤 波变换(IGHM)得到融合后的图像F。 低频 图像 MNSDFB 系数 A 融合 高撅 低颜融 合规则 低频 系数 系数 反向 融合图像 MNSDFB F 低频 融合 高频融 图像 MNSDFB 系数 高频 B 高频 合规则 系数 系数 图2本文提出的融合方法 Fig.2 Proposed image method 3实验及结果分析 法相同NSCT分解，融合规则采用修改的拉普拉斯 能量和融合图像：3)Kumar提出的采用Cross bilater- 为验证本文算法的有效性与优越性，将本文提 al filter的方法1s),本实验中采用和文献[l5]完全 出方法与其他3种算法进行了多聚焦图像融合对 相同的参数设置。本文提出方法采用2个尺度 比。1)Das提出的采用MSF激发脉冲耦合神经 GHM的Multiwavelet分解，对每个高频多小波系数 网络(pulse couple neural network,PCNN)融合所有 采用4方向NSDFB分解，低频系数采用SML融合 NSCT子带系数的方法。其中，非子采样方向滤波 规则，其他所有高频子带采用本文提出DTRF结合 器设置为“cd”,金字塔滤波器采用“maxflat”滤波 MSF的融合规则选择系数。图3分别是进行融合 器，分解方向设置为[2344]：2)采用与第1种方 测试的Clock、Lab、Disk、Pepsi4对源多聚焦图像。 a)Clock左聚焦图像 (b)Lab左聚焦图像(c)Disk左聚焦图像(d)Pepsi左聚焦图像 (e)Clock右聚焦图像 (DLab右聚焦图像 (g)Disk右聚焦图像(h)Pepsiz左聚焦图像 图3第一组源多聚焦图像 Fig.3 The first group of source images

进行滤波，利用 ＤＴＲＦ 滤波器输出去生成更为接近 理想状态的融合决策图。 利用式（９）计算所提出的 ＭＮＳＤＦＢ 高频系数的 ＭＳＦ 值，并将其作为 ＤＴＲＦ 滤 波器的输入，可以得到： Ｒ ｈ，ｋ （ｉ，ｊ） ＝ ＲＦ（ＭＳＦ ｈ，ｋ（ｉ，ｊ）） （１０） 由 Ｒ ｈ，ｋ（ｉ，ｊ） 得到的决策图为 Ｍ ｈ，ｋ （ｉ，ｊ） ＝ １， Ｒ ｈ，ｋ Ａ （ｉ，ｊ） ≥ Ｒ ｈ，ｋ Ｂ （ｉ，ｊ） ０， Ｒ ｈ，ｋ Ａ （ｉ，ｊ） ＜ Ｒ ｈ，ｋ { Ｂ （ｉ，ｊ） （１１） 根据决策图选择融合后的系数 Ｃ ｈ，ｋ Ｆ （ｉ，ｊ） ： Ｃ ｈ，ｋ Ｆ （ｉ，ｊ） ＝ Ｄ ｈ，ｋ Ａ （ｉ，ｊ）， Ｍ ｈ，ｋ （ｉ，ｊ） ＝ １ Ｄ ｈ，ｋ Ｂ （ｉ，ｊ）， Ｍ ｈ，ｋ { （ｉ，ｊ） ＝ ０ （１２） 式中： Ｃ ｈ，ｋ Ａ （ｉ，ｊ） 、 Ｃ ｈ，ｋ Ｂ （ｉ，ｊ） 和 Ｃ ｈ，ｋ Ｆ （ｉ，ｊ） 分别是源多 聚焦图像 Ａ、Ｂ 和融合后图像 Ｆ 在 ｋ 方向的高频分 解系数。 ２．４ 算法步骤 本文融合方法如图 ２ 所示并描述如下： １）将源图像 Ａ 和 Ｂ 进行 ｌ 尺度的 Ｍｕｌｔｉｗａｖｅｌｅｔ 分解， 得 到 包 含 多 小 波 低 频 系 数 Ｇ ｌ Ａ (ｉ，ｊ) 、 Ｇ ｌ Ｂ (ｉ，ｊ) 和多小波高频系数 Ｇ ｈ Ａ (ｉ，ｊ) 、 Ｇ ｈ Ｂ (ｉ，ｊ) 。 ２）对 Ｇ ｈ Ａ (ｉ，ｊ) 和 Ｇ ｈ Ｂ (ｉ，ｊ) 的 ３ 个高频系数进行 ＮＳＤＦＢ 分 解 得 到 ｋ 方 向 子 带 系 数 Ｃ ｈ，ｋ Ａ (ｉ，ｊ) 和 Ｃ ｈ，ｋ Ｂ (ｉ，ｊ) ，ｋ 为 ＮＳＤＦＢ 分解的方向。 ３）对低频系数和高频系数分别按照相应的融 合规则 选 择 系 数 得 到 融 合 后 的 系 数 Ｃ ｌ Ｆ (ｉ，ｊ) 和 Ｃ ｈ，ｋ Ｆ (ｉ，ｊ) 。 ４）对融合后系数 Ｃ ｌ Ｆ (ｉ，ｊ) 和 Ｃ ｈ，ｋ Ｆ (ｉ，ｊ) 进行反 ＮＳＤＦＢ 变 换 得 到 融 合 后 的 多 小 波 低 频 系 数 Ｇ ｌ Ｆ (ｉ，ｊ) 和高频系数 Ｇ ｈ Ｆ (ｉ，ｊ) 。 ５）对小波系数 Ｇ ｌ Ｆ (ｉ，ｊ) 、 Ｇ ｈ Ｆ (ｉ，ｊ) 进行反多小 波变换（ＩＧＨＭ）得到融合后的图像 Ｆ。 图 ２ 本文提出的融合方法 Ｆｉｇ．２ Ｐｒｏｐｏｓｅｄ ｉｍａｇｅ ｍｅｔｈｏｄ ３ 实验及结果分析 为验证本文算法的有效性与优越性，将本文提 出方法与其他 ３ 种算法进行了多聚焦图像融合对 比。 １）Ｄａｓ ［１４］提出的采用 ＭＳＦ 激发脉冲耦合神经 网络（ ｐｕｌｓｅ ｃｏｕｐｌｅ ｎｅｕｒａｌ ｎｅｔｗｏｒｋ，ＰＣＮＮ） 融合所有 ＮＳＣＴ 子带系数的方法。 其中，非子采样方向滤波 器设置为“ ｃｄ”，金字塔滤波器采用“ ｍａｘｆｌａｔ” 滤波 器，分解方向设置为［２ ３ ４ ４］；２）采用与第 １ 种方 法相同 ＮＳＣＴ 分解，融合规则采用修改的拉普拉斯 能量和融合图像；３）Ｋｕｍａｒ 提出的采用 Ｃｒｏｓｓ ｂｉｌａｔｅｒ⁃ ａｌ ｆｉｌｔｅｒ 的方法［１５］ ，本实验中采用和文献［１５］ 完全 相同的参数设置。 本文提出方法采用 ２ 个尺度 ＧＨＭ 的 Ｍｕｌｔｉｗａｖｅｌｅｔ 分解，对每个高频多小波系数 采用 ４ 方向 ＮＳＤＦＢ 分解，低频系数采用 ＳＭＬ 融合 规则，其他所有高频子带采用本文提出 ＤＴＲＦ 结合 ＭＳＦ 的融合规则选择系数。 图 ３ 分别是进行融合 测试的 Ｃｌｏｃｋ、Ｌａｂ、Ｄｉｓｋ、Ｐｅｐｓｉ ４ 对源多聚焦图像。 图 ３ 第一组源多聚焦图像 Ｆｉｇ．３ Ｔｈｅ ｆｉｒｓｔ ｇｒｏｕｐ ｏｆ ｓｏｕｒｃｅ ｉｍａｇｅｓ ·２４４· 智 能 系 统 学 报 第 １１ 卷

第2期 任晓霞，等：多小波和NSDFB组合域递归滤波多聚焦图像融合 .245. Lab图像的融合结果如图4(a)~(d)所示，对 焦图像中的聚焦良好区域的图像信息，同时引入了 于人物头部区域来说，本文提出的方法具有更少的 更少的伪影信息。 伪影信息，而其他3个方法的图像在该区域具有更 除了主观的视觉评价外，本文采用互信息(M如 为明显的伪影。对于Disk图像的融合结果来说，图 tual information,MI)、Qv[6]对融合结果进行客观 5(a)和图5(b)中具有明显的伪影，而图5(d)中出 评价，M和QA这两个指标并不要求知道理想的 现了明显的边缘模糊现象，图5()这些现象最为不 融合图像（也称为参考图像），并且大量的应用在图 明显：为了更为明显地展示不同方法的融合效果，图 像融合客观评价中718)。M也称为相关嫡，用来评 4(e)~(h)和图5(e)~(h)给出了不同方法获得的 价源图像与融合图像之间的相似度，互信息值越大， 融合结果图像与源图像的残差图像。对于残差图像 说明融合的效果越好；Qr利用Sobel边缘检测来 来说，更少的残差信息说明融合图像的相应区域绝 衡量有多少源图像的边缘细节信息转移到了融合图 大多数来自于源图像中清晰的部分，很少或者没有 像，Q的值越大，说明融合的效果越好：表1给出 来自于另一幅源图像的模糊部分。图4(e)和图5 了4种融合方法得到融合图像的评价指标的计算结 (g)具有最为明显的残差信息，图4(f)、图4(g)、图 果。从表1看，本文提出的方法对于4对多聚焦图 5(e)和图5(f)的残差虽不太明显，但依然比图4 像都取得了最高的I和Q值，这表明所提出方 (h)和图5(h)要明显得多。图4(h)和图5(h)几乎 法的融合结果相对于其他3种方法来说，融合了更 没有残留信息。以上分析表明：本文提出的融合方 多源图像的信息。 法相对于其他3种方法，能够获得更多来自源多聚 表1不同融合方法的客观评价 Table 1 Objective evaluation of fused images with different methods 图像评价指标Das方法 NSCT SML Kumar方法 本文方法 MI 7.7478 7.6351 7.4774 8.1208 Lab OARVE 0.7183 0.7338 0.7321 0.7416 MI 7.2022 6.8510 6.6735 7.7612 Disk OARIF 0.6992 0.7139 0.6950 0.7184 MI 7.2205 7.4496 7.4139 7.7855 Clock OE 0.6477 0.6743 0.6873 0.6889 MI 7.7260 7.4017 7.2282 7.5135 Pepsi OARVE 0.7592 0.7768 0.7798 0.7825 (a)Das方法 (b)NSCT SMI (c)Kumar方法 (d)本文方法 (e)图(a)与图3(e)的差 ()图(b)与图3(e)的差(g)图(c)与图3(e)的差(h)图(d)与图3(e)的差 图4Lb图像融合实验 Fig.4 Image fusion experiment on Lab image 76 (a)Das方法 (b)NSCT SML (c)Kumar方法 (d)本文方法

Ｌａｂ 图像的融合结果如图 ４（ ａ） ～ （ ｄ） 所示，对 于人物头部区域来说，本文提出的方法具有更少的 伪影信息，而其他 ３ 个方法的图像在该区域具有更 为明显的伪影。 对于 Ｄｉｓｋ 图像的融合结果来说，图 ５（ａ）和图 ５（ｂ）中具有明显的伪影，而图 ５（ｄ）中出 现了明显的边缘模糊现象，图 ５（ｄ）这些现象最为不 明显；为了更为明显地展示不同方法的融合效果，图 ４（ｅ） ～ （ｈ）和图 ５（ ｅ） ～ （ ｈ）给出了不同方法获得的 融合结果图像与源图像的残差图像。 对于残差图像 来说，更少的残差信息说明融合图像的相应区域绝 大多数来自于源图像中清晰的部分，很少或者没有 来自于另一幅源图像的模糊部分。 图 ４（ ｅ） 和图 ５ （ｇ）具有最为明显的残差信息，图 ４（ｆ）、图 ４（ｇ）、图 ５（ｅ）和图 ５（ ｆ） 的残差虽不太明显，但依然比图 ４ （ｈ）和图 ５（ｈ）要明显得多。 图 ４（ｈ）和图 ５（ｈ）几乎 没有残留信息。 以上分析表明：本文提出的融合方 法相对于其他 ３ 种方法，能够获得更多来自源多聚 焦图像中的聚焦良好区域的图像信息，同时引入了 更少的伪影信息。 除了主观的视觉评价外，本文采用互信息（Ｍｕ⁃ ｔｕａｌ ｉｎｆｏｒｍａｔｉｏｎ， ＭＩ）、Ｑ ＡＢ／ Ｆ［１６］ 对融合结果进行客观 评价，ＭＩ 和 Ｑ ＡＢ／ Ｆ这两个指标并不要求知道理想的 融合图像（也称为参考图像），并且大量的应用在图 像融合客观评价中［１７⁃１８］ 。 ＭＩ 也称为相关熵，用来评 价源图像与融合图像之间的相似度，互信息值越大， 说明融合的效果越好；Ｑ ＡＢ／ Ｆ利用 Ｓｏｂｅｌ 边缘检测来 衡量有多少源图像的边缘细节信息转移到了融合图 像，Ｑ ＡＢ／ Ｆ的值越大，说明融合的效果越好；表 １ 给出 了 ４ 种融合方法得到融合图像的评价指标的计算结 果。 从表 １ 看，本文提出的方法对于 ４ 对多聚焦图 像都取得了最高的 ＭＩ 和 Ｑ ＡＢ／ Ｆ值，这表明所提出方 法的融合结果相对于其他 ３ 种方法来说，融合了更 多源图像的信息。 表 １ 不同融合方法的客观评价 Ｔａｂｌｅ １ Ｏｂｊｅｃｔｉｖｅ ｅｖａｌｕａｔｉｏｎ ｏｆ ｆｕｓｅｄ ｉｍａｇｅｓ ｗｉｔｈ ｄｉｆｆｅｒｅｎｔ ｍｅｔｈｏｄｓ 图像 评价指标 Ｄａｓ 方法 ＮＳＣＴ＿ＳＭＬ Ｋｕｍａｒ 方法 本文方法 Ｌａｂ ＭＩ ７．７４７ ８ ７．６３５ １ ７．４７７ ４ ８．１２０ ８ Ｑ ＡＢ／ Ｆ ０．７１８ ３ ０．７３３ ８ ０．７３２ １ ０．７４１ ６ Ｄｉｓｋ ＭＩ ７．２０２ ２ ６．８５１ ０ ６．６７３ ５ ７．７６１ ２ Ｑ ＡＢ／ Ｆ ０．６９９ ２ ０．７１３ ９ ０．６９５ ０ ０．７１８ ４ Ｃｌｏｃｋ ＭＩ ７．２２０ ５ ７．４４９ ６ ７．４１３ ９ ７．７８５ ５ Ｑ ＡＢ／ Ｆ ０．６４７ ７ ０．６７４ ３ ０．６８７ ３ ０．６８８ ９ Ｐｅｐｓｉ ＭＩ ７．７２６ ０ ７．４０１ ７ ７．２２８ ２ ７．５１３ ５ Ｑ ＡＢ／ Ｆ ０．７５９ ２ ０．７７６ ８ ０．７７９ ８ ０．７８２ ５ 图 ４ Ｌａｂ 图像融合实验 Ｆｉｇ．４ Ｉｍａｇｅ ｆｕｓｉｏｎ ｅｘｐｅｒｉｍｅｎｔ ｏｎ Ｌａｂ ｉｍａｇｅ 第 ２ 期 任晓霞，等：多小波和 ＮＳＤＦＢ 组合域递归滤波多聚焦图像融合 ·２４５·

.246 智能系统学报 第11卷 (e)图(a)与图3b)的差 ()图(b)与图3(b)的差(g)图(c)与图3(b)的差(h)图(d)与图3(b)的差 图5Disk图像融合实验 Fig.5 Image fusion experiment on Disk image 由于MI和Q这两种图像融合评价方法仅 表2不同融合方法的客观评价 仅衡量融合中的信息量和边缘保留量。为了对 Table 2 Objective evaluation of fused images with dif- 融合结果图像的灰度失真、结构畸变等方面进行 ferent methods 评估，本文另外采用了SSM1]和VIF20]两个指 图像 评价指标Das方法 Kumar方法 本文方法 标进行了客观评价。鉴于上述两个指标为有参 SSIM 0.9875 0.9879 0.9913 考图像评价方法，利用了的Lena、Barbara、Peppers VIF 0.9334 0.8949 0.9459 3幅测试图像，对这3幅图像的左右两个部分分 Lena MI 5.9953 6.2832 6.5310 别进行高斯模糊，结果如图6所示。对图6中3 OARF 0.6876 0.7158 0.7181 对源图像采用本文方法、Das方法和Kumar方法 SSIM 0.9957 0.9966 0.9987 分别进行融合，融合结果的客观评价结果如表2 VIF 0.9770 0.9588 0.9926 所示，由表2可知，本文方法相对于另外两种方 Barbara M 6.2859 6.6144 6.7473 法来说，不仅MI和QF这两个指标取得较大值， OARE 0.7582 0.7691 0.7734 SSIM和VIF两个指标也得到最大值。 SSIM 0.9940 0.9935 0.9941 VIF 0.9712 0.9431 0.9775 Peppers M 6.7859 7.0801 7.4432 DAWE 0.6739 0.6848 0.6876 最后，表3给出了本文提出的MNSDFB变换方 法和NSCT方法在融合不同大小图像的运行时间对 比，由于本文试验中所有源图像只有480×640、512× (a)左聚焦图像1 (b)右聚焦图像1 512两种尺寸。因此，本文只选择了Clock图像和 Lb图像进行试验。对于上述两种融合方法来说， 低频均采用了均值规则，高频均采用极大值规则。 进行实验的计算机的CPU为四核2.5GHz,内存为 4GB。实验结果表明，所提出MNSDFB变换方法相 对于NSCT方法具有更快的运算速度。 综上所述，本文的算法得到的融合图像能够有 (c)左聚焦图像2 (d)右聚焦图像2 效选择源图像中清晰区域的图像信息，对图像边缘 以及细节信息表述得更为清晰突出，图像灰度失真 和结构畸变较小，主观目视效果更好，并且也获得了 最好的客观评价指标。 表3 NSDFB与NSCT融合不同图像的运算时间对比 Table 3 Comparisons on computation time with NSDFB and NSCT method (e)左聚焦图像3 ()右聚焦图像3 图像 尺寸 NSCT方法 本文方法 图6第二组源多聚焦图像 Clock 512×512 34.73 24.74 Fig.6 The second group of source images Lab 480×640 36.20 25.54

图 ５ Ｄｉｓｋ 图像融合实验 Ｆｉｇ．５ Ｉｍａｇｅ ｆｕｓｉｏｎ ｅｘｐｅｒｉｍｅｎｔ ｏｎ Ｄｉｓｋ ｉｍａｇｅ 由于 ＭＩ 和 Ｑ ＡＢ ／ Ｆ 这两种图像融合评价方法仅 仅衡量融合中的信息量和边缘保留量。 为了对 融合结果图像的灰度失真、结构畸变等方面进行 评估，本文另外采用了 ＳＳＩＭ ［ １９］ 和 ＶＩＦ ［ ２０］ 两个指 标进行了客观评价。 鉴于上述两个指标为有参 考图像评价方法，利用了的 Ｌｅｎａ、Ｂａｒｂａｒａ、Ｐｅｐｐｅｒｓ ３ 幅测试图像，对这 ３ 幅图像的左右两个部分分 别进行高斯模糊，结果如图 ６ 所示。 对图 ６ 中 ３ 对源图像采用本文方法、Ｄａｓ 方法和 Ｋｕｍａｒ 方法 分别进行融合，融合结果的客观评价结果如表 ２ 所示，由表 ２ 可知，本文方法相对于另外两种方 法来说，不仅 ＭＩ 和 Ｑ ＡＢ ／ Ｆ 这两个指标取得较大值， ＳＳＩＭ 和 ＶＩＦ 两个指标也得到最大值。 图 ６ 第二组源多聚焦图像 Ｆｉｇ．６ Ｔｈｅ ｓｅｃｏｎｄ ｇｒｏｕｐ ｏｆ ｓｏｕｒｃｅ ｉｍａｇｅｓ 表 ２ 不同融合方法的客观评价 Ｔａｂｌｅ ２ Ｏｂｊｅｃｔｉｖｅ ｅｖａｌｕａｔｉｏｎ ｏｆ ｆｕｓｅｄ ｉｍａｇｅｓ ｗｉｔｈ ｄｉｆ⁃ ｆｅｒｅｎｔ ｍｅｔｈｏｄｓ 图像 评价指标 Ｄａｓ 方法 Ｋｕｍａｒ 方法 本文方法 Ｌｅｎａ ＳＳＩＭ ０．９８７ ５ ０．９８７ ９ ０．９９１ ３ ＶＩＦ ０．９３３ ４ ０．８９４ ９ ０．９４５ ９ ＭＩ ５．９９５ ３ ６．２８３ ２ ６．５３１ ０ Ｑ ＡＢ／ Ｆ ０．６８７ ６ ０．７１５ ８ ０．７１８ １ Ｂａｒｂａｒａ ＳＳＩＭ ０．９９５ ７ ０．９９６ ６ ０．９９８ ７ ＶＩＦ ０．９７７ ０ ０．９５８ ８ ０．９９２ ６ ＭＩ ６．２８５ ９ ６．６１４ ４ ６．７４７ ３ Ｑ ＡＢ／ Ｆ ０．７５８ ２ ０．７６９ １ ０．７７３ ４ Ｐｅｐｐｅｒｓ ＳＳＩＭ ０．９９４ ０ ０．９９３ ５ ０．９９４ １ ＶＩＦ ０．９７１ ２ ０．９４３ １ ０．９７７ ５ ＭＩ ６．７８５ ９ ７．０８０ １ ７．４４３ ２ Ｑ ＡＢ／ Ｆ ０．６７３ ９ ０．６８４ ８ ０．６８７ ６ 最后，表 ３ 给出了本文提出的 ＭＮＳＤＦＢ 变换方 法和 ＮＳＣＴ 方法在融合不同大小图像的运行时间对 比，由于本文试验中所有源图像只有 ４８０×６４０、５１２× ５１２ 两种尺寸。 因此，本文只选择了 Ｃｌｏｃｋ 图像和 Ｌａｂ 图像进行试验。 对于上述两种融合方法来说， 低频均采用了均值规则，高频均采用极大值规则。 进行实验的计算机的 ＣＰＵ 为四核 ２．５ ＧＨｚ，内存为 ４ ＧＢ。 实验结果表明，所提出 ＭＮＳＤＦＢ 变换方法相 对于 ＮＳＣＴ 方法具有更快的运算速度。 综上所述，本文的算法得到的融合图像能够有 效选择源图像中清晰区域的图像信息，对图像边缘 以及细节信息表述得更为清晰突出，图像灰度失真 和结构畸变较小，主观目视效果更好，并且也获得了 最好的客观评价指标。 表 ３ ＮＳＤＦＢ 与 ＮＳＣＴ 融合不同图像的运算时间对比 Ｔａｂｌｅ ３ Ｃｏｍｐａｒｉｓｏｎｓ ｏｎ ｃｏｍｐｕｔａｔｉｏｎ ｔｉｍｅ ｗｉｔｈ ＮＳＤＦＢ ａｎｄ ＮＳＣＴ ｍｅｔｈｏｄ ｓ 图像 尺寸 ＮＳＣＴ 方法 本文方法 Ｃｌｏｃｋ ５１２×５１２ ３４．７３ ２４．７４ Ｌａｂ ４８０×６４０ ３６．２０ ２５．５４ ·２４６· 智 能 系 统 学 报 第 １１ 卷

第2期 任晓霞，等：多小波和NSDFB组合域递归滤波多聚焦图像融合 ·247. image J.Acta optica sinica,2005,25(9):1176-1180. 4 结束语 [7]WANG Zhihui,ZHAO Baojun,SHEN Tingzhi.Image fusion 图像的多尺度分解和融合规则是决定多聚焦图 based on multiwavelet transformation and PCNN[J.Journal 像融合结果的关键因素，本文利用多小波的高频子 of information computational science,2008,5(3):1397- 带系数进行NSDFB分解构建了MNSDFB变换，低频 1406. 系数采用SML融合规则，高频系数采用域变换递归 [8]朱四荣，王迎春.基于多小波变换的多聚焦图像融合 滤波器融合规则。实验结果表明：与NSCT变换改 [J].计算机工程与应用，2010,46(6)：169-170,226. ZHU Sirong,WANG Yingchun.Multi-focus image fusion 进空间频率激发PCNN的融合规则、NSCT变换采用 method based on multi-wavelet transform[J].Computer en- SML的融合方法和Cross bilateral filter融合方法相 gineering and applications,2010,46(6):169-170,226. 比较，本文所提出的方法获取的融合图像保留丰富 [9]DA CUNHA A L,ZHOU Jianping,DO M N.The nonsub- 图像细节信息，获得了更好的视觉效果，并且采用互 sampled contourlet transform:Theory,design,and applica- 信息、QF、SSM、VIF等客观指标进行比较，取得 tions[J].IEEE transactions on image processing,2006,15 了较好的客观融合评价结果。多小波和MNSDFB (10):3089-3101. 变换依然存在冗余问题，因此，如何减少本文提出方 「10]米德伶，冯鹏，魏彪，等.一种基于清晰度计算的NSCT 法的冗余度是未来的研究方向。 域多聚焦图像融合算法[J].光学与光电技术，2010,8 (2):7-10. 参考文献： MI Deling,FENG Peng,WEI Biao,et al.A multi-focus image fusion algorithm based on the calculation of clarity in [1]BAI Xiangzhi,ZHANG Yu,ZHOU Fugen,et al.Quadtree- NSCT domain [J].Optics optoelectronic technology, based multi-focus image fusion using a weighted focus-meas- 2010,8(2)：7-10. ure[J].Information fusion,2015,22:105-118. [11]GASTAL E S L,OLIVEIRA MM.Domain transform for [2]许良凤，林辉，胡敏.基于差分进化算法的多模态医学 edge-aware image and video processing[J].ACM transac- 图像融合[J].电子测量与仪器学，2013,27(2)：110- tions on graphics,2011,30(4):69. 114 [12]GENG Peng,HUANG Min,LIU Shuaiqi,et al.Multifocus XU Liangfeng,LIN Hui,HU Min.Multi-modality medicine image fusion method of Ripplet transform based on cycle image fusion based on differential evolution algorithm [J]. spinning[J].Multimedia tools and applications,2014:1- Journal of electronic measurement and instrument,2013,27 11,doi:10.1007/s11042-014-1942.1. (2):110-114 [13]屈小波，闫敬文，杨贵德.改进拉普拉斯能量和的尖锐 [3]李华锋，柴毅，张晓阳.提升静态小波域内多聚焦图像 频率局部化Contourlet域多聚焦图像融合方法[J].光 融合算法[J].光电工程，2011,38(3)：131-137,144. 学精密工程，2009,17(5)：1203-1212 LI Huafeng,CHAI Yi,ZHANG Xiaoyang.Multi-focus im- QU Xiaobo,YAN Jingwen,YANG Guide.Multifocus im- age fusion algorithm in lifting stationary wavelet domain[J]. age fusion method of sharp frequency localized Contourlet Opto-electronic engineering,2011,38(3):131-137,144. transform domain based on sum-modified-Laplacian [J]. [4]LIU Xuan,ZHOU Yue,WANG Jiajun.Image fusion based Optics and precision engineering,2009,17 (5):1203- on shearlet transform and regional features[J].AEU-interna- 1212. tional journal of electronics and communications,2013,68 [14]DAS S,KUNDU M K.NSCT-based multimodal medical (6):471-477. image fusion using pulse-coupled neural network and modi- [5]苗启广，王宝树.基于小波变换与局部能量的多聚焦图 fied spatial frequency[J].Medical biological engineer- 像融合[J].计算机科学，2005,32(2)：229-232. ing&computing.,2012,50(10):1105-1114. MIAO Qiguang,WANG Baoshu.Multi-focus image fusion [15]KUMAR B K S.Image fusion based on pixel significance u- based on wavelet transform and local energy[J].Computer sing cross bilateral filter[]]Signal,image and video pro- science,.2005,32(2):229-232. cessing,2015,9(5):1193-1204. [6]宗思光，王江安.多尺度形态算子融合图像滤波技术及 [16]GANGAPURE V N,BANERJEE S,CHOWDHURY A S. 滤波质量评价[J].光学学报，2005,25(9)：1176-1180. Steerable local frequency based multispectral multifocus im- ZONG Siguang,WANG Jiang'an.Multi-scale mathematical age fusion[]].Information fusion,2015,23:99-115. morphology fusion filtering and the evaluation of the filtered [17]李俊峰，姜晓丽，戴文战.基于提升小波变换的医学图

４ 结束语 图像的多尺度分解和融合规则是决定多聚焦图 像融合结果的关键因素，本文利用多小波的高频子 带系数进行 ＮＳＤＦＢ 分解构建了 ＭＮＳＤＦＢ 变换，低频 系数采用 ＳＭＬ 融合规则，高频系数采用域变换递归 滤波器融合规则。 实验结果表明：与 ＮＳＣＴ 变换改 进空间频率激发 ＰＣＮＮ 的融合规则、ＮＳＣＴ 变换采用 ＳＭＬ 的融合方法和 Ｃｒｏｓｓ ｂｉｌａｔｅｒａｌ ｆｉｌｔｅｒ 融合方法相 比较，本文所提出的方法获取的融合图像保留丰富 图像细节信息，获得了更好的视觉效果，并且采用互 信息、 Ｑ ＡＢ／ Ｆ 、ＳＳＩＭ、ＶＩＦ 等客观指标进行比较，取得 了较好的客观融合评价结果。 多小波和 ＭＮＳＤＦＢ 变换依然存在冗余问题，因此，如何减少本文提出方 法的冗余度是未来的研究方向。 参考文献： ［１］ＢＡＩ Ｘｉａｎｇｚｈｉ， ＺＨＡＮＧ Ｙｕ， ＺＨＯＵ Ｆｕｇｅｎ， ｅｔ ａｌ． Ｑｕａｄｔｒｅｅ⁃ ｂａｓｅｄ ｍｕｌｔｉ⁃ｆｏｃｕｓ ｉｍａｇｅ ｆｕｓｉｏｎ ｕｓｉｎｇ ａ ｗｅｉｇｈｔｅｄ ｆｏｃｕｓ⁃ｍｅａｓ⁃ ｕｒｅ［Ｊ］． Ｉｎｆｏｒｍａｔｉｏｎ ｆｕｓｉｏｎ， ２０１５， ２２： １０５⁃１１８． ［２］许良凤， 林辉， 胡敏． 基于差分进化算法的多模态医学 图像融合［ Ｊ］． 电子测量与仪器学， ２０１３， ２７（ ２）： １１０⁃ １１４． ＸＵ Ｌｉａｎｇｆｅｎｇ， ＬＩＮ Ｈｕｉ， ＨＵ Ｍｉｎ． Ｍｕｌｔｉ⁃ｍｏｄａｌｉｔｙ ｍｅｄｉｃｉｎｅ ｉｍａｇｅ ｆｕｓｉｏｎ ｂａｓｅｄ ｏｎ ｄｉｆｆｅｒｅｎｔｉａｌ ｅｖｏｌｕｔｉｏｎ ａｌｇｏｒｉｔｈｍ ［ Ｊ］． Ｊｏｕｒｎａｌ ｏｆ ｅｌｅｃｔｒｏｎｉｃ ｍｅａｓｕｒｅｍｅｎｔ ａｎｄ ｉｎｓｔｒｕｍｅｎｔ， ２０１３， ２７ （２）： １１０⁃１１４． ［３］李华锋， 柴毅， 张晓阳． 提升静态小波域内多聚焦图像 融合算法［Ｊ］． 光电工程， ２０１１， ３８（３）： １３１⁃１３７， １４４． ＬＩ Ｈｕａｆｅｎｇ， ＣＨＡＩ Ｙｉ， ＺＨＡＮＧ Ｘｉａｏｙａｎｇ． Ｍｕｌｔｉ⁃ｆｏｃｕｓ ｉｍ⁃ ａｇｅ ｆｕｓｉｏｎ ａｌｇｏｒｉｔｈｍ ｉｎ ｌｉｆｔｉｎｇ ｓｔａｔｉｏｎａｒｙ ｗａｖｅｌｅｔ ｄｏｍａｉｎ［Ｊ］． Ｏｐｔｏ⁃ｅｌｅｃｔｒｏｎｉｃ ｅｎｇｉｎｅｅｒｉｎｇ， ２０１１， ３８（３）： １３１⁃１３７， １４４． ［４］ＬＩＵ Ｘｕａｎ， ＺＨＯＵ Ｙｕｅ， ＷＡＮＧ Ｊｉａｊｕｎ． Ｉｍａｇｅ ｆｕｓｉｏｎ ｂａｓｅｄ ｏｎ ｓｈｅａｒｌｅｔ ｔｒａｎｓｆｏｒｍ ａｎｄ ｒｅｇｉｏｎａｌ ｆｅａｔｕｒｅｓ［Ｊ］． ＡＥＵ⁃ｉｎｔｅｒｎａ⁃ ｔｉｏｎａｌ ｊｏｕｒｎａｌ ｏｆ ｅｌｅｃｔｒｏｎｉｃｓ ａｎｄ ｃｏｍｍｕｎｉｃａｔｉｏｎｓ， ２０１３， ６８ （６）： ４７１⁃４７７． ［５］苗启广， 王宝树． 基于小波变换与局部能量的多聚焦图 像融合［Ｊ］． 计算机科学， ２００５， ３２（２）： ２２９⁃２３２． ＭＩＡＯ Ｑｉｇｕａｎｇ， ＷＡＮＧ Ｂａｏｓｈｕ． Ｍｕｌｔｉ⁃ｆｏｃｕｓ ｉｍａｇｅ ｆｕｓｉｏｎ ｂａｓｅｄ ｏｎ ｗａｖｅｌｅｔ ｔｒａｎｓｆｏｒｍ ａｎｄ ｌｏｃａｌ ｅｎｅｒｇｙ［ Ｊ］． Ｃｏｍｐｕｔｅｒ ｓｃｉｅｎｃｅ， ２００５， ３２（２）： ２２９⁃２３２． ［６］宗思光， 王江安． 多尺度形态算子融合图像滤波技术及 滤波质量评价［Ｊ］． 光学学报， ２００５， ２５（９）： １１７６⁃１１８０． ＺＯＮＧ Ｓｉｇｕａｎｇ， ＷＡＮＧ Ｊｉａｎｇ＇ａｎ． Ｍｕｌｔｉ⁃ｓｃａｌｅ ｍａｔｈｅｍａｔｉｃａｌ ｍｏｒｐｈｏｌｏｇｙ ｆｕｓｉｏｎ ｆｉｌｔｅｒｉｎｇ ａｎｄ ｔｈｅ ｅｖａｌｕａｔｉｏｎ ｏｆ ｔｈｅ ｆｉｌｔｅｒｅｄ ｉｍａｇｅ［Ｊ］． Ａｃｔａ ｏｐｔｉｃａ ｓｉｎｉｃａ， ２００５， ２５（９）： １１７６⁃１１８０． ［７］ＷＡＮＧ Ｚｈｉｈｕｉ， ＺＨＡＯ Ｂａｏｊｕｎ， ＳＨＥＮ Ｔｉｎｇｚｈｉ． Ｉｍａｇｅ ｆｕｓｉｏｎ ｂａｓｅｄ ｏｎ ｍｕｌｔｉｗａｖｅｌｅｔ ｔｒａｎｓｆｏｒｍａｔｉｏｎ ａｎｄ ＰＣＮＮ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ ｉｎｆｏｒｍａｔｉｏｎ ＆ ｃｏｍｐｕｔａｔｉｏｎａｌ ｓｃｉｅｎｃｅ， ２００８， ５（３）： １３９７⁃ １４０６． ［８］朱四荣， 王迎春． 基于多小波变换的多聚焦图像融合 ［Ｊ］． 计算机工程与应用， ２０１０， ４６（６）： １６９⁃１７０， ２２６． ＺＨＵ Ｓｉｒｏｎｇ， ＷＡＮＧ Ｙｉｎｇｃｈｕｎ． Ｍｕｌｔｉ⁃ｆｏｃｕｓ ｉｍａｇｅ ｆｕｓｉｏｎ ｍｅｔｈｏｄ ｂａｓｅｄ ｏｎ ｍｕｌｔｉ⁃ｗａｖｅｌｅｔ ｔｒａｎｓｆｏｒｍ［ Ｊ］． Ｃｏｍｐｕｔｅｒ ｅｎ⁃ ｇｉｎｅｅｒｉｎｇ ａｎｄ ａｐｐｌｉｃａｔｉｏｎｓ， ２０１０， ４６（６）： １６９⁃１７０， ２２６． ［９］ＤＡ ＣＵＮＨＡ Ａ Ｌ， ＺＨＯＵ Ｊｉａｎｐｉｎｇ， ＤＯ Ｍ Ｎ． Ｔｈｅ ｎｏｎｓｕｂ⁃ ｓａｍｐｌｅｄ ｃｏｎｔｏｕｒｌｅｔ ｔｒａｎｓｆｏｒｍ： Ｔｈｅｏｒｙ， ｄｅｓｉｇｎ， ａｎｄ ａｐｐｌｉｃａ⁃ ｔｉｏｎｓ［Ｊ］． ＩＥＥＥ ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ ｉｍａｇｅ ｐｒｏｃｅｓｓｉｎｇ， ２００６， １５ （１０）： ３０８９⁃３１０１． ［１０］米德伶， 冯鹏， 魏彪， 等． 一种基于清晰度计算的 ＮＳＣＴ 域多聚焦图像融合算法［Ｊ］． 光学与光电技术， ２０１０， ８ （２）： ７⁃１０． ＭＩ Ｄｅｌｉｎｇ， ＦＥＮＧ Ｐｅｎｇ， ＷＥＩ Ｂｉａｏ， ｅｔ ａｌ． Ａ ｍｕｌｔｉ⁃ｆｏｃｕｓ ｉｍａｇｅ ｆｕｓｉｏｎ ａｌｇｏｒｉｔｈｍ ｂａｓｅｄ ｏｎ ｔｈｅ ｃａｌｃｕｌａｔｉｏｎ ｏｆ ｃｌａｒｉｔｙ ｉｎ ＮＳＣＴ ｄｏｍａｉｎ ［ Ｊ ］． Ｏｐｔｉｃｓ ＆ ｏｐｔｏｅｌｅｃｔｒｏｎｉｃ ｔｅｃｈｎｏｌｏｇｙ， ２０１０， ８（２）： ７⁃１０． ［１１］ＧＡＳＴＡＬ Ｅ Ｓ Ｌ， ＯＬＩＶＥＩＲＡ Ｍ Ｍ． Ｄｏｍａｉｎ ｔｒａｎｓｆｏｒｍ ｆｏｒ ｅｄｇｅ⁃ａｗａｒｅ ｉｍａｇｅ ａｎｄ ｖｉｄｅｏ ｐｒｏｃｅｓｓｉｎｇ［ Ｊ］． ＡＣＭ ｔｒａｎｓａｃ⁃ ｔｉｏｎｓ ｏｎ ｇｒａｐｈｉｃｓ， ２０１１， ３０（４）： ６９． ［１２］ＧＥＮＧ Ｐｅｎｇ， ＨＵＡＮＧ Ｍｉｎ， ＬＩＵ Ｓｈｕａｉｑｉ， ｅｔ ａｌ． Ｍｕｌｔｉｆｏｃｕｓ ｉｍａｇｅ ｆｕｓｉｏｎ ｍｅｔｈｏｄ ｏｆ Ｒｉｐｐｌｅｔ ｔｒａｎｓｆｏｒｍ ｂａｓｅｄ ｏｎ ｃｙｃｌｅ ｓｐｉｎｎｉｎｇ［Ｊ］． Ｍｕｌｔｉｍｅｄｉａ ｔｏｏｌｓ ａｎｄ ａｐｐｌｉｃａｔｉｏｎｓ， ２０１４： １⁃ １１， ｄｏｉ： １０．１００７ ／ ｓ１１０４２⁃０１４⁃１９４２⁃１． ［１３］屈小波， 闫敬文， 杨贵德． 改进拉普拉斯能量和的尖锐 频率局部化 Ｃｏｎｔｏｕｒｌｅｔ 域多聚焦图像融合方法［ Ｊ］． 光 学精密工程， ２００９， １７（５）： １２０３⁃１２１２． ＱＵ Ｘｉａｏｂｏ， ＹＡＮ Ｊｉｎｇｗｅｎ， ＹＡＮＧ Ｇｕｉｄｅ． Ｍｕｌｔｉｆｏｃｕｓ ｉｍ⁃ ａｇｅ ｆｕｓｉｏｎ ｍｅｔｈｏｄ ｏｆ ｓｈａｒｐ ｆｒｅｑｕｅｎｃｙ ｌｏｃａｌｉｚｅｄ Ｃｏｎｔｏｕｒｌｅｔ ｔｒａｎｓｆｏｒｍ ｄｏｍａｉｎ ｂａｓｅｄ ｏｎ ｓｕｍ⁃ｍｏｄｉｆｉｅｄ⁃Ｌａｐｌａｃｉａｎ ［ Ｊ ］． Ｏｐｔｉｃｓ ａｎｄ ｐｒｅｃｉｓｉｏｎ ｅｎｇｉｎｅｅｒｉｎｇ， ２００９， １７ （ ５）： １２０３⁃ １２１２． ［１４］ ＤＡＳ Ｓ， ＫＵＮＤＵ Ｍ Ｋ． ＮＳＣＴ⁃ｂａｓｅｄ ｍｕｌｔｉｍｏｄａｌ ｍｅｄｉｃａｌ ｉｍａｇｅ ｆｕｓｉｏｎ ｕｓｉｎｇ ｐｕｌｓｅ⁃ｃｏｕｐｌｅｄ ｎｅｕｒａｌ ｎｅｔｗｏｒｋ ａｎｄ ｍｏｄｉ⁃ ｆｉｅｄ ｓｐａｔｉａｌ ｆｒｅｑｕｅｎｃｙ［ Ｊ］． Ｍｅｄｉｃａｌ ＆ ｂｉｏｌｏｇｉｃａｌ ｅｎｇｉｎｅｅｒ⁃ ｉｎｇ ＆ ｃｏｍｐｕｔｉｎｇ， ２０１２， ５０（１０）： １１０５⁃１１１４． ［１５］ＫＵＭＡＲ Ｂ Ｋ Ｓ． Ｉｍａｇｅ ｆｕｓｉｏｎ ｂａｓｅｄ ｏｎ ｐｉｘｅｌ ｓｉｇｎｉｆｉｃａｎｃｅ ｕ⁃ ｓｉｎｇ ｃｒｏｓｓ ｂｉｌａｔｅｒａｌ ｆｉｌｔｅｒ［Ｊ］． Ｓｉｇｎａｌ， ｉｍａｇｅ ａｎｄ ｖｉｄｅｏ ｐｒｏ⁃ ｃｅｓｓｉｎｇ， ２０１５， ９（５）： １１９３⁃１２０４． ［１６］ＧＡＮＧＡＰＵＲＥ Ｖ Ｎ， ＢＡＮＥＲＪＥＥ Ｓ， ＣＨＯＷＤＨＵＲＹ Ａ Ｓ． Ｓｔｅｅｒａｂｌｅ ｌｏｃａｌ ｆｒｅｑｕｅｎｃｙ ｂａｓｅｄ ｍｕｌｔｉｓｐｅｃｔｒａｌ ｍｕｌｔｉｆｏｃｕｓ ｉｍ⁃ ａｇｅ ｆｕｓｉｏｎ［Ｊ］． Ｉｎｆｏｒｍａｔｉｏｎ ｆｕｓｉｏｎ， ２０１５， ２３： ９９⁃１１５． ［１７］李俊峰， 姜晓丽， 戴文战． 基于提升小波变换的医学图 第 ２ 期 任晓霞，等：多小波和 ＮＳＤＦＢ 组合域递归滤波多聚焦图像融合 ·２４７·

.248. 智能系统学报 第11卷 像融合[J].中国图象图形学报，2014,19(11)：1639- 作者简介： 1648. 任晓霞，女，1982年生，讲师，主要研 LI Junfeng,JIANG Xiaoli,DAI Wenzhan.Medical image 究方向为计算机技术、图形图像技术。 fusion based on lifting wavelet transform[J].Journal of im- age and graphics,2014,19(11):1639-1648. [l8]高继森，董亚楠，沈瑜，等.基于改进Tetrolet变换的图 像融合算法研究[J].计算机科学，2015,42(5)：320- 323. 孙秀明，男，1978年生，副教授，主 GAO Jisen,DONG Yanan,SHEN Yu,et al.Research of 要研究方向为图像去噪、图像融合以及 image fusion algorithm based on improved Tetrolet transform 为数字化资源开发。 []Computer science,2015,42(5):320-323. [19]WANG Zhou,BOVIK A C,SHEIKH H R,et al.Image quality assessment:from error visibility to structural simi- larity[J].IEEE transactions on image processing,2004, 耿鹏，男，1979年生，副教授，主要 13(4):600-612. 研究方向为图像融合、小波分析。主持 [20]SHEIKH H R,BOVIK A C.Image information and visual 河北省自然科学基金项目、河北省教育 quality[J].IEEE transactions on image processing,2006. 厅高等学校科学研究项目等多项，发表 15(2):430-444. 学术论文多篇，其中被SCI检索4篇，被 EI检索6篇。 EEE第20届智能工程系统国际会议 IEEE 20th Jubilee International Conference on Intelligent Engineering Systems 2016 Authors are welcome to submit original and unpublished papers and attend the IEEE 20th Jubilee International Con- ference on Intelligent Engineering Systems 2016(INES 2016)to be held on June 30-July 2 2016 in Hungary. Topics include but not limited to: Artificial Intelligence in Engineering Reasoning,Learning,Decision Making,Knowledge Based Systems,Expert Systems;CAD/CAM/CAE Systems:Product Modeling,Shape Modeling,Manufacturing Process Planning;Communica- tions Software and Systems in Engineering:Design Methodologies and Tools,Object-oriented,UML,Software Engineer- ing;Computational Intelligence in Engineering:Machine Learning,Genetic Algorithms,Neural Nets,Fuzzy Systems, Fuzzy and Neuro-fuzzy Control;Intelligent Manufacturing Systems:Production Planning and Scheduling,Rapid Prototy- ping,Flexible Manufacturing Systems,Collaborative Engineering,Concurrent Engineering;Intelligent Mechatronics and Robotics Systems:Control,Perception and Recognition,Sensing and Sensor Data Fusion,Intelligent Sensors,Intelligent Motion Control,Service Robots Intelligent Signal Processing;Intelligent Transportation Systems:Navigation Systems,On- board Systems,Real-time Traffic Control;Man-Machine Systems:Human Computer Interaction,Multimedia Communica- tions,Advanced Computer Graphics,Virtual Reality;Ontologies and Semantic Engineering:Ontology,Thesaurus,Disam- biguation,Semantic Inference,Natural Language Interaction;Systems Engineering:Systems Analysis,Systems Methodolo- gy,Self-Organizing Systems,Systems Integration,Large Scale Systems,Systems Simulation,Diagnosis and Performance Monitoring;Web Engineering:Intelligent Information Retrieval and Recommendation,Textual and Visual Interaction and Interfaces,User Behavior Analysis and Modeling,Data Modeling,Information Visualization,Knowledge Acquisition. Website:http://www.ines-conf.org/ines-conf/2016index.html

像融合［Ｊ］． 中国图象图形学报， ２０１４， １９（１１）： １６３９⁃ １６４８． ＬＩ Ｊｕｎｆｅｎｇ， ＪＩＡＮＧ Ｘｉａｏｌｉ， ＤＡＩ Ｗｅｎｚｈａｎ． Ｍｅｄｉｃａｌ ｉｍａｇｅ ｆｕｓｉｏｎ ｂａｓｅｄ ｏｎ ｌｉｆｔｉｎｇ ｗａｖｅｌｅｔ ｔｒａｎｓｆｏｒｍ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ ｉｍ⁃ ａｇｅ ａｎｄ ｇｒａｐｈｉｃｓ， ２０１４， １９（１１）： １６３９⁃１６４８． ［１８］高继森， 董亚楠， 沈瑜， 等． 基于改进 Ｔｅｔｒｏｌｅｔ 变换的图 像融合算法研究［ Ｊ］． 计算机科学， ２０１５， ４２（５）： ３２０⁃ ３２３． ＧＡＯ Ｊｉｓｅｎ， ＤＯＮＧ Ｙａｎａｎ， ＳＨＥＮ Ｙｕ， ｅｔ ａｌ． Ｒｅｓｅａｒｃｈ ｏｆ ｉｍａｇｅ ｆｕｓｉｏｎ ａｌｇｏｒｉｔｈｍ ｂａｓｅｄ ｏｎ ｉｍｐｒｏｖｅｄ Ｔｅｔｒｏｌｅｔ ｔｒａｎｓｆｏｒｍ ［Ｊ］． Ｃｏｍｐｕｔｅｒ ｓｃｉｅｎｃｅ， ２０１５， ４２（５）： ３２０⁃３２３． ［１９］ ＷＡＮＧ Ｚｈｏｕ， ＢＯＶＩＫ Ａ Ｃ， ＳＨＥＩＫＨ Ｈ Ｒ， ｅｔ ａｌ． Ｉｍａｇｅ ｑｕａｌｉｔｙ ａｓｓｅｓｓｍｅｎｔ： ｆｒｏｍ ｅｒｒｏｒ ｖｉｓｉｂｉｌｉｔｙ ｔｏ ｓｔｒｕｃｔｕｒａｌ ｓｉｍｉ⁃ ｌａｒｉｔｙ［ Ｊ］． ＩＥＥＥ ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ ｉｍａｇｅ ｐｒｏｃｅｓｓｉｎｇ， ２００４， １３（４）： ６００⁃６１２． ［２０］ＳＨＥＩＫＨ Ｈ Ｒ， ＢＯＶＩＫ Ａ Ｃ． Ｉｍａｇｅ ｉｎｆｏｒｍａｔｉｏｎ ａｎｄ ｖｉｓｕａｌ ｑｕａｌｉｔｙ［Ｊ］． ＩＥＥＥ ｔｒａｎｓａｃｔｉｏｎｓ ｏｎ ｉｍａｇｅ ｐｒｏｃｅｓｓｉｎｇ， ２００６， １５（２）： ４３０⁃４４４． 作者简介： 任晓霞，女，１９８２ 年生，讲师，主要研 究方向为计算机技术、图形图像技术。 孙秀明，男，１９７８ 年生，副教授， 主 要研究方向为图像去噪、图像融合以及 为数字化资源开发。 耿鹏，男，１９７９ 年生，副教授，主要 研究方向为图像融合、小波分析。 主持 河北省自然科学基金项目、河北省教育 厅高等学校科学研究项目等多项，发表 学术论文多篇，其中被 ＳＣＩ 检索 ４ 篇，被 ＥＩ 检索 ６ 篇。 ＩＥＥＥ 第 ２０ 届智能工程系统国际会议 ＩＥＥＥ ２０ｔｈ Ｊｕｂｉｌｅｅ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｉｎｔｅｌｌｉｇｅｎｔ Ｅｎｇｉｎｅｅｒｉｎｇ Ｓｙｓｔｅｍｓ ２０１６ Ａｕｔｈｏｒｓ ａｒｅ ｗｅｌｃｏｍｅ ｔｏ ｓｕｂｍｉｔ ｏｒｉｇｉｎａｌ ａｎｄ ｕｎｐｕｂｌｉｓｈｅｄ ｐａｐｅｒｓ ａｎｄ ａｔｔｅｎｄ ｔｈｅ ＩＥＥＥ ２０ｔｈ Ｊｕｂｉｌｅｅ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎ⁃ ｆｅｒｅｎｃｅ ｏｎ Ｉｎｔｅｌｌｉｇｅｎｔ Ｅｎｇｉｎｅｅｒｉｎｇ Ｓｙｓｔｅｍｓ ２０１６ （ＩＮＥＳ ２０１６） ｔｏ ｂｅ ｈｅｌｄ ｏｎ Ｊｕｎｅ ３０－Ｊｕｌｙ ２ ２０１６ ｉｎ Ｈｕｎｇａｒｙ． Ｔｏｐｉｃｓ ｉｎｃｌｕｄｅ ｂｕｔ ｎｏｔ ｌｉｍｉｔｅｄ ｔｏ： Ａｒｔｉｆｉｃｉａｌ Ｉｎｔｅｌｌｉｇｅｎｃｅ ｉｎ Ｅｎｇｉｎｅｅｒｉｎｇ ： Ｒｅａｓｏｎｉｎｇ， Ｌｅａｒｎｉｎｇ， Ｄｅｃｉｓｉｏｎ Ｍａｋｉｎｇ， Ｋｎｏｗｌｅｄｇｅ Ｂａｓｅｄ Ｓｙｓｔｅｍｓ， Ｅｘｐｅｒｔ Ｓｙｓｔｅｍｓ；ＣＡＤ／ ＣＡＭ／ ＣＡＥ Ｓｙｓｔｅｍｓ： Ｐｒｏｄｕｃｔ Ｍｏｄｅｌｉｎｇ， Ｓｈａｐｅ Ｍｏｄｅｌｉｎｇ， Ｍａｎｕｆａｃｔｕｒｉｎｇ Ｐｒｏｃｅｓｓ Ｐｌａｎｎｉｎｇ；Ｃｏｍｍｕｎｉｃａ⁃ ｔｉｏｎｓ Ｓｏｆｔｗａｒｅ ａｎｄ Ｓｙｓｔｅｍｓ ｉｎ Ｅｎｇｉｎｅｅｒｉｎｇ： Ｄｅｓｉｇｎ Ｍｅｔｈｏｄｏｌｏｇｉｅｓ ａｎｄ Ｔｏｏｌｓ， Ｏｂｊｅｃｔ⁃ｏｒｉｅｎｔｅｄ， ＵＭＬ， Ｓｏｆｔｗａｒｅ Ｅｎｇｉｎｅｅｒ⁃ ｉｎｇ；Ｃｏｍｐｕｔａｔｉｏｎａｌ Ｉｎｔｅｌｌｉｇｅｎｃｅ ｉｎ Ｅｎｇｉｎｅｅｒｉｎｇ： Ｍａｃｈｉｎｅ Ｌｅａｒｎｉｎｇ， Ｇｅｎｅｔｉｃ Ａｌｇｏｒｉｔｈｍｓ， Ｎｅｕｒａｌ Ｎｅｔｓ， Ｆｕｚｚｙ Ｓｙｓｔｅｍｓ， Ｆｕｚｚｙ ａｎｄ Ｎｅｕｒｏ⁃ｆｕｚｚｙ Ｃｏｎｔｒｏｌ；Ｉｎｔｅｌｌｉｇｅｎｔ Ｍａｎｕｆａｃｔｕｒｉｎｇ Ｓｙｓｔｅｍｓ： Ｐｒｏｄｕｃｔｉｏｎ Ｐｌａｎｎｉｎｇ ａｎｄ Ｓｃｈｅｄｕｌｉｎｇ， Ｒａｐｉｄ Ｐｒｏｔｏｔｙ⁃ ｐｉｎｇ， Ｆｌｅｘｉｂｌｅ Ｍａｎｕｆａｃｔｕｒｉｎｇ Ｓｙｓｔｅｍｓ， Ｃｏｌｌａｂｏｒａｔｉｖｅ Ｅｎｇｉｎｅｅｒｉｎｇ， Ｃｏｎｃｕｒｒｅｎｔ Ｅｎｇｉｎｅｅｒｉｎｇ；Ｉｎｔｅｌｌｉｇｅｎｔ Ｍｅｃｈａｔｒｏｎｉｃｓ ａｎｄ Ｒｏｂｏｔｉｃｓ Ｓｙｓｔｅｍｓ： Ｃｏｎｔｒｏｌ， Ｐｅｒｃｅｐｔｉｏｎ ａｎｄ Ｒｅｃｏｇｎｉｔｉｏｎ， Ｓｅｎｓｉｎｇ ａｎｄ Ｓｅｎｓｏｒ Ｄａｔａ Ｆｕｓｉｏｎ， Ｉｎｔｅｌｌｉｇｅｎｔ Ｓｅｎｓｏｒｓ， Ｉｎｔｅｌｌｉｇｅｎｔ Ｍｏｔｉｏｎ Ｃｏｎｔｒｏｌ， Ｓｅｒｖｉｃｅ Ｒｏｂｏｔｓ Ｉｎｔｅｌｌｉｇｅｎｔ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｉｎｇ；Ｉｎｔｅｌｌｉｇｅｎｔ Ｔｒａｎｓｐｏｒｔａｔｉｏｎ Ｓｙｓｔｅｍｓ： Ｎａｖｉｇａｔｉｏｎ Ｓｙｓｔｅｍｓ， Ｏｎ⁃ ｂｏａｒｄ Ｓｙｓｔｅｍｓ， Ｒｅａｌ⁃ｔｉｍｅ Ｔｒａｆｆｉｃ Ｃｏｎｔｒｏｌ；Ｍａｎ⁃Ｍａｃｈｉｎｅ Ｓｙｓｔｅｍｓ： Ｈｕｍａｎ Ｃｏｍｐｕｔｅｒ Ｉｎｔｅｒａｃｔｉｏｎ， Ｍｕｌｔｉｍｅｄｉａ Ｃｏｍｍｕｎｉｃａ⁃ ｔｉｏｎｓ， Ａｄｖａｎｃｅｄ Ｃｏｍｐｕｔｅｒ Ｇｒａｐｈｉｃｓ， Ｖｉｒｔｕａｌ Ｒｅａｌｉｔｙ；Ｏｎｔｏｌｏｇｉｅｓ ａｎｄ Ｓｅｍａｎｔｉｃ Ｅｎｇｉｎｅｅｒｉｎｇ： Ｏｎｔｏｌｏｇｙ， Ｔｈｅｓａｕｒｕｓ， Ｄｉｓａｍ⁃ ｂｉｇｕａｔｉｏｎ， Ｓｅｍａｎｔｉｃ Ｉｎｆｅｒｅｎｃｅ， Ｎａｔｕｒａｌ Ｌａｎｇｕａｇｅ Ｉｎｔｅｒａｃｔｉｏｎ；Ｓｙｓｔｅｍｓ Ｅｎｇｉｎｅｅｒｉｎｇ： Ｓｙｓｔｅｍｓ Ａｎａｌｙｓｉｓ， Ｓｙｓｔｅｍｓ Ｍｅｔｈｏｄｏｌｏ⁃ ｇｙ， Ｓｅｌｆ⁃Ｏｒｇａｎｉｚｉｎｇ Ｓｙｓｔｅｍｓ， Ｓｙｓｔｅｍｓ Ｉｎｔｅｇｒａｔｉｏｎ， Ｌａｒｇｅ Ｓｃａｌｅ Ｓｙｓｔｅｍｓ， Ｓｙｓｔｅｍｓ Ｓｉｍｕｌａｔｉｏｎ， Ｄｉａｇｎｏｓｉｓ ａｎｄ Ｐｅｒｆｏｒｍａｎｃｅ Ｍｏｎｉｔｏｒｉｎｇ；Ｗｅｂ Ｅｎｇｉｎｅｅｒｉｎｇ： Ｉｎｔｅｌｌｉｇｅｎｔ Ｉｎｆｏｒｍａｔｉｏｎ Ｒｅｔｒｉｅｖａｌ ａｎｄ Ｒｅｃｏｍｍｅｎｄａｔｉｏｎ， Ｔｅｘｔｕａｌ ａｎｄ Ｖｉｓｕａｌ Ｉｎｔｅｒａｃｔｉｏｎ ａｎｄ Ｉｎｔｅｒｆａｃｅｓ， Ｕｓｅｒ Ｂｅｈａｖｉｏｒ Ａｎａｌｙｓｉｓ ａｎｄ Ｍｏｄｅｌｉｎｇ， Ｄａｔａ Ｍｏｄｅｌｉｎｇ， Ｉｎｆｏｒｍａｔｉｏｎ Ｖｉｓｕａｌｉｚａｔｉｏｎ， Ｋｎｏｗｌｅｄｇｅ Ａｃｑｕｉｓｉｔｉｏｎ． Ｗｅｂｓｉｔｅ：ｈｔｔｐ： ／ ／ ｗｗｗ．ｉｎｅｓ⁃ｃｏｎｆ．ｏｒｇ ／ ｉｎｅｓ⁃ｃｏｎｆ ／ ２０１６ｉｎｄｅｘ．ｈｔｍｌ ·２４８· 智 能 系 统 学 报 第 １１ 卷