www.pudn.com > enhancement+fusing.rar > Wavelet_fusing.m

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图像 
 
[cl,sl]=wavedec2(X1,2,'sym4'); 
sizec1=size(cl); 
for I=1:sizec1(2) 
    cl( I )=1.2*cl( I ); 
end 
 
[c2,s2]=wavedec2(X2,2,'sym4'); 
c=cl+c2; 
c=0.5*c; 
xx=waverec2(c,sl,'sym4'); 
subplot(223); 
 
image(xx); 
title('融合图象'); 
axis square                                   % 画出融合后的图像