clear all;
close all;
I1=imread('D:\matlabtu\car1.bmp');
h=ones(4,4)/16;
I2=imfilter(I1,h); %ͨ¹ýËÄÁÚÓò·½·¨µÃµ½µÍ·Ö±æÂʵÄͼÏñ
figure(1);
imshow(I1),title('ÔʼͼÏñ');
figure(2);
imshow(I2),title('µÍ·Ö±æÂʵÄͼÏñ');
%ÓÃË«ÏßÐÔ²åÖµ·½·¨»ñµÃ²åֵͼÏñY1
[Y1,map]=imresize(I2,2,'bilinear');
%ÓÃ×î½üÁÚÓò²åÖµµÃµ½ÁÚÓò²åֵͼÏñY2
[Y2,map]=imresize(I2,2);
[c,s]=wavedec2(Y1,2,'haar');
sizey1=size(Y1);
%´ÓС²¨·Ö½âµÄ½á¹¹[c£¬s]ÖÐÌáÈ¡Y1µÚÒ»²ãµÄ³ß¶ÈϵÊýºÍС²¨ÏµÊý
Xa1=appcoef2(c,s,'haar',1);
Xh1=detcoef2('h',c,s,1);
Xv1=detcoef2('v',c,s,1);
Xd1=detcoef2('d',c,s,1);
ded1=[Xa1,Xh1,Xv1,Xd1];
nbcol=size(map,1);
[c,s]=wavedec2(Y1,2,'haar');
sizey2=size(Y2);
%´ÓС²¨·Ö½âµÄ½á¹¹[c£¬s]ÖÐÌáÈ¡Y2µÚÒ»²ãµÄ³ß¶ÈϵÊýºÍС²¨ÏµÊý
Xa2=appcoef2(c,s,'haar',1);
Xh2=detcoef2('h',c,s,1);
Xv2=detcoef2('v',c,s,1);
Xd2 =detcoef2('d',c,s,1);
ded1=[Xa2,Xh2,Xv2,Xd2];
nbcol=size(map,1);
Y=idwt2(Xa2,Xh1,Xv1,Xd1,'haar');
nbcol=size(map,1);
figure(3);
imshow(uint8(Y)),title('»ñµÃµÄ³¬·Ö±æÂÊͼÏñ');
% err = Y - I2;
% err = err( ;
% PSNRdb = 20 * log10(256/sqrt(mean(err .^2)));
ÓÐûÓп´µÄ¶®µÄ£¬µÚÈýÕÅͼƬ³ö²»À´
¸ø³öµÄ´íÎóÊÇ
![ÇóÖú »ùÓÚС²¨ºÍ²åÖµµÄ³¬·Ö±æÂÊͼÏñÖØ½¨Ëã·¨]()
b60b070828381f30b7282a1aaa014c086e06f070.jpg |