d2bldwt9_7.m

来自「matlab下的提升小波程序」· M 代码 · 共 28 行

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 
% Contact me : [e-mail],http://guohanwei.51.net,QQ:461414909

 [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:xx);
         det_sig=img(l,xx+1:2*xx);
         img(l,1:2*xx)= bldwt9_7(app_sig,det_sig);
     end
     xx=2*xx;
     yy=2*yy;
 end