examp16_17.m

《小波分析理论与Matlab7实现例程》一书里面各章节的源码

%装载原图像1
load bust
X1=X;
map1=map;
%画出woman图像
subplot(221);image(X1);
colormap(map1);title('原图像1');
axis square
%装载原图像2
load mask
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
%画出wbarb图像
subplot(222);
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