f11_8.m

机械测试信号处理 附书的matlab源码。

%装入原图像1
load woman;
X1=X;
map1=map;
subplot(221);
%画出woman图像
image(X1);
colormap(map1);
title('原图像1');
%axis square

%装入原图像2
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);
%画出wbarb图像
image(X2);
colormap(map2);
title('原图像2');
%axis square 

[c1,s1]=wavedec2(X1,2,'sym4');
sizec1=size(c1);
for I=1:sizec1(2)
     c1( I )=1.2*c1( I );
end

[c2,s2]=wavedec2(X2,2,'sym4');
c=c1+c2;
c=0.5*c;
s=s1+s2;
s=0.5*s;
xx=waverec2(c,s,'sym4');

subplot(223);
%画出融合后的图像
image(xx);
title('融合图像');
%axis square