- [X,map]=imread('lina.bmp');
- X1=X;
- map1=map;
- subplot(221);
- image(X1);
- colormap(map1);
- title('lina');
- axis square
- [X,map]=imread('couple.bmp');
- X2=X;
- map2=map;
- subplot(222);
- image(X2);
- colormap(map2);
- title('couple');
- 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;
- xx=waverec2(c,s1,'sym4');
- subplot(223);
- image(xx);
- title('融合图象');
- axis square