d2bldwt9_7.m

SPIHT的MATLAB源代码,希望对有需要的朋友有所帮助.

function img=d2bldwt9_7(dwt_img,level)
%
% img=d2bldwt9_7(dwt_img,level)
%
% The forward lifting discrete 2-D wavelet transform (D2BLDWT) using the 9-7 Antonini filters 
 [x,y]=size(dwt_img);
 xx=x/(2^level);
 yy=y/(2^level);
 img=dwt_img;
 for i=1:level
     
     for  k=1:2*yy
         app_sig=img(1:xx,k);
         det_sig=img(xx+1:2*xx,k);
         img(1:2*xx,k)= bldwt9_7(app_sig,det_sig);
     end
     
     for l=1:2*xx
         
         app_sig=img(l,1:yy);
         det_sig=img(l,yy+1:2*yy);
         img(l,1:2*yy)= bldwt9_7(app_sig,det_sig);
     end
     xx=2*xx;
     yy=2*yy;
 end