【图像融合】基于梯度能量、局部能量、 PCA三种融合规则实现MRI-CT图像融合附matlab代码

1 简介

​大多数图像融合算法只就图像的某一个特征进行融合,容易造成其他特征信息损失.针对此问题,提出了一种利用哈尔小波变换的特性,考虑梯度能量、局部能量、 PCA三种特征，采用区域方差与变换系数相结合的方法进行融合,最后经哈尔小波逆变换得到融合结果.

2 部分代码

clearx1=imread('ct.jpg');subplot(221);imshow(x1);title('ct图像');x2=imread('mri.jpg');subplot(222);imshow(x2);title('mri图像');X1=double(x1);X2=double(x2);%将两幅图像进行一级db4小波变换[ca1,ch1,cv1,cd1]=dwt2(X1,'db4');subplot(223);b1=uint8([ca1 ch1;cv1 cd1]);imshow(b1);title('ct图像一级db4小波分解');[ca2,ch2,cv2,cd2]=dwt2(X2,'db4');subplot(224);b2=uint8([ca2 ch2;cv2 cd2]);imshow(b2);title('mri图像一级db4小波分解');%按照一定融合规则进行融合ca=dipinchuli(ca1,ca2);ch=PCAfusion(ch1,ch2);cv=PCAfusion(cv1,cv2);cd=PCAfusion(cd1,cd2);figure;subplot(121);imshow(uint8([ca ch;cv cd]));title('融合后小波系数');X=idwt2(ca,ch,cv,cd,'db4');subplot(122);imshow(uint8(X));title('融合后重构图像');

3 仿真结果

4 参考文献

[1]张宏怡, 杨文, 林永平. 基于梯度和能量的图像融合算法[J]. 厦门理工学院学报, 2010, 18(1):5.

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

部分理论引用网络文献，若有侵权联系博主删除。