cwq5.m

matlab程序源码,这个是我和同学一起做的,大家可以参考看看怎么样的,谢谢大家,可以提议供大家学!

%图象融合
load woman; %装入原图像
X1=X;map1=map;
subplot(221);
image(X1);
colormap(map1);
title('woman');
axis square %画出woman 图像
load wbarb; %装入原图像
X2=X;map2=map;
for I =1:256
for j=1:256
if(X2(I, j)>100)
X2(I, j)=1.2*X2(I, j);
else
X2(I, j)=0.5*X2(I, j);
end
end
end
subplot(222);
image(X2);
colormap(map2);
title('wbarb');
axis square %画出wbarb 图像
%对上述二图像进行分解
[c1,s1]=wavedec2(X1,2,'sym4');
[c2,s2]=wavedec2(X2,2,'sym4');
%对分解的图像进行融合
c=c1+c2;
%应用融合系数进行图像重构并显示
xx=waverec2(c,s1,'sym4');
subplot(223);image(xx);
title('融合图象1');
sizec1=size(c1);
for I=1:sizec1(2)
 c1( I )=1.2*c1( I );
end
sizec2=size(c2);
for j=1:sizec2(2)
 c2( j)=0.8*c2( j );
end
%通过减小融合系数以减小图象宽度
c=0.5*(c1+c2);
%对融合系数进行图象重构
xxx=waverec2(c,s2,'sym4');
%显示重构结果
subplot(224);image(xxx);
title('融合图象2');