基于混沌遗传算法和contourlet变换的医学图像融合

1、    【摘要】  目的： 改善传统医学图象融合方法对细节信息的丢失. 方法： 利用contourlet的多尺度、方向性和各向异性等优点，提出了一种基于局域特性匹配度融合的改进算法. 首先对原图进行contourlet分解，对分解各子带构造局域特征的匹配度，设定匹配度的阈值，对匹配度在不同阈值范围内的系数进行不同的加权融合；最后对融合系数进行contourlet逆变换得到融合图像. 在阈值的选取上，引入混沌遗传算法求解阈值的全局最优解. 结果： 运用传统小波变换和本文提出对方法对两组医学图象进行融合处理，客观评价参数表明本文提出的算法

2、效果更优. 结论： 该算法能够在保有原图信息的同时，有效的增强细节信息. 【关键词】  contourlet 医学图像 图像融合 遗传算法【Abstract】AIM: To improve the traditional image fusion algorithm in order to avoid the loss of  the detailed information in the processe of image fusion. METHODS: Utilizing the contourlets advantages of multisca

3、le,directionality and anisotropy,an advanced image fusion algorithm based on characteristic matching of region statistics in contourlet domain was proposed in this paper. First，the source images were transformed into contourlet domain. Then,we constructed the characteristic matching of each subband

4、and set a threshold for it. The coefficients whose characteristic matching was in different ranges applied different fusion rules. Finally, the fused image was obtained by taking the inverse contourlet transformation. Chaos genetic algorithm was adopted to optimize fusion image by searching optimal

5、solution of the threshold of characteristic matching. RESULTS: The comparison between the result of the traditional method and that of the proposed method showed that the later one was better. CONCLUSION: The proposed method can effectively preserve the information of the source images and enhance t

6、he details of the fused image.【Keywords】 contourlet; medical image; image fusion; genetic algorithms0   引言    传统的医学图像融合算法多采用基于小波变换，小波系数加权平均，或者阈值法融合，但二维可分小波仅能捕捉有限的方向信息，不能稀疏的表示含线或者面奇异的二维图像1-3. Do等4提出了一种稀疏的图像二维表示方法： contourlet变换. 该变换的“各向异性(anisotropy)”使得它能比小波变换更好的表现边缘特征，更加适合用来进行图像的处

7、理5. 并且小波变换下的传统方法有很大的局限，比如加权平均法中加权因子和基于特征匹配度融合方法中的阈值大小都不易确定6. 本文提出了一种新的基于匹配度融合算法：在contourlet变换多方向的优势下，阈值的确定采用混沌遗传算法. 混沌遗传算法是近年提出的一种优化算法，能够有效的防止遗传算法的“早熟”问题，使得遗传算法能够以更快的速度收敛，得到阈值的全局最优解.1    材料和方法1.1   材料   选用两组融合图象，一组为实验标准图象，另一组为腰椎冠状位的MRI-T1和MRI-T2成像，编码格式为DICOM，来源于西安红十字会医院

8、磁共振室.1.2   方法1.2.1   contourlet变换原理   contourlet变换也称塔型方向滤波器组(pyramidal directional filter bank,PDFB). contourlet 变换将多尺度分析和方向分析分拆进行：用拉普拉斯塔形滤波器结构(laplacian pyramid,LP)对图像进行多分辨率分解. 首先产生原始信号的一个低通采样逼近及原始图像与低通预测图像之间的一个差值图像，对得到的低通图像继续分解得到下一层的低通逼近和差值图像. 如此逐步滤波得到图像的多分辨率分解. 而后应用二

9、维方向滤波器组(directional filter bank,DFB)对分解得到的每一级高频分量在任意尺度上再分解得到 个方向子带，将分布在同方向上的奇异点合成为一个系数. 图1A给出了离散Contourlet 变换的滤波器组结构. 原始图像经PDFB结构分解得到一个低通图像和分布于多尺度多方向上的高频分量. 图1B为Contourlet 频域分解图实际应用中方向数随着尺度增大而增多.   A： 滤波器组结构图； B： 频率分解图.图1   contourlet变换（略）1.2.2   图像融合框架   本文的融

10、合方法的改进之一就是将contourlet 引入图像融合，可以利用其优良特性更好地提取原始图像中的几何特征，为融合图像提供更多的信息.1.2.2.1   contourlet分解   将精确配准的原图像进行Contourlet变换,得到相应的contourlet系数集合. 在融合中一般采用3层或4层分解. 设原图像为A和B，使用J级LP分解，每个尺度上的方向数分别为i(j)，则原图像的Contourlet分解过程为:    A(aJA,bJA,bJ-1A,b2A,b1A)(1)    B(aJB,bJB,bJ-1B

11、,b2B,b1B)(2)    bJid2Ji,d2JI,d1i,   i=A,B(3)    其中aJi是低频子带; bJi为尺度j上的方向子带集合; dki为尺度j上第k个方向的高频子带,k=1,2,2J.1.2.2.2   图像融合   对于分解后的低频子带和所有高频子带，使用基于区域特性的融合规则进行判别和融合处理,得到各尺度上融合后的contourlet系数.1.2.2.3   Contourlet重构   重构是分解的逆过程,对融合后的cont

12、ourlet系数进行contourlet逆变换，得到重构的融合图像，该图像包含原有多幅图像中的信息. 设 为融合后的结果，此过程可以表示为    (aJF,bJF,bJ-1F,b2F,b1F)F(4)1.2.3   融合规则   不同的融合算法都有其相应的应用价值. 目前还没有一种融合算法适合所有的图像融合. contourlet对线条有高度敏感性，这对细节质量的提高很有好处. 由于contourlet变换系数与其邻域有极强的依赖性，在子带系数处理中，我们选用了基于区域特征的方法，所以对若干个方向的高频分量，本文采用如下方法进行融合：

13、    （1） 选择一个矩形窗口L×K，如，计算contourlet分解各子带系数的均值. 子带中以点(n,m)为中心的区域均值为：    M-i(n,m)=1/L×K/nL/mKfi(n+n,m+m)(5)    i代表原图A或B，fi表示图A或图B的相应子带系数.     （2） 计算以点(n,m)为中心的区域方差i2(n,m)     2i(n,m)=1/L×K/nL/mK(fi(n+n,m+m)-M-I(n,m)2(6)  &#

14、160; （3） 计算以点(n,m)为中心图像A与图像B子带对应区域的协方差：    2(n,m)=1/L×K/nL/mK(fA(n+n,m+m)-M-A(n,m)×fB(n+n,m+m)-M-B(n,m)(7)    （4） 构造匹配度(n,m)：    (n,m)=2×2(n,m)/2A(n,m)+2B(n,m)(8)    （5） 选定阈值，其值大小一般在0.50.8之间. 最优值取决于待融合的图像. 本文用混沌遗传算法来确定阈值. 具体方法见后. 如果(n,m)，则融合结果中对应点值dKF(n,m)为（d2j为高频子带）：  &