12-14.m

matlab图像处理源码

load tartan;
X1=X;map1=map;
load sinsin;
X2=X;map2=map;
% 打开图像
subplot(121)
image(X1),colormap(map1);
title('图像tartan')
subplot(122)
image(X2),colormap(map2);
title('图像sinsin')

%%%    图像直接融合
figure,subplot(131)
image((X1+X2)/2),colormap(map2);
% 在空域内图像直接融合
title('两图像直接相加融合')

%%%    图像直接融合
figure,subplot(131)
image((X1+X2)/2),colormap(map2);
% 在空域内图像直接融合
title('两图像直接相加融合')

%%%    图像小波变换融合
[C1,L1]=wavedec2(X1,2,'sym4');
[C2,L2]=wavedec2(X2,2,'sym4');
% 分别对上述两图像进行分解
C=C1+C2;
% 对分解系数进行融合
X=waverec2(C,L1,'sym4');
% 融合后的信号进行图像重构
subplot(133)
image(X/2),colormap(map2);
% 显示小波变换的图像融合
title('两图像小波变换融合')