wavelet decomposition process code in matlabSolution WAVELET

wavelet decomposition process code in matlab???

Solution

WAVELET DECOMPOSITION PROCESS :-

function [coef,scf]=mywavedec2(x,N, wname)
coef={cA_N;cV_N;cH_N;cD_N;cV_N-1;cH_N-1;cD_N-1;… …;cV_1;cH_1;cD_1}
[Lo_D,Hi_D] = wfilters(wname,\'d\');
lf=length(Lo_D); imshow(x);title(\'Original Image\');
[r,c]=size(x);
xlabel([\'Size : \',num2str(r),\'*\',num2str(c)]);
xd=double(x);              
[rx,cx]=size(x);
[o2sa,f1sa,rsx]=sizcoef([rx,cx],lf,N); a=[o2sa,f1sa,rsx]
coef=[]; scf=[rx,cx];
for i=1:N   
          [cA,cV,cH,cD]=mydwt2(xd,Lo_D,Hi_D);
           xd=cA;
           outmp={cV;cH;cD};     
           scf=[size(cV);scf];    
           coef=[outmp;coef];       
end
coef=[cA;coef];
scf=[size(cA);scf];
plotcoef(N,wname,coef);
subplot(N+1,3,2);
yt=uint8(coef{1});
[yrow,ycol]=size(yt);
imshow(yt);title( [\'Approximation A\',num2str(N)]);
xlabel([\'Size : \',num2str(yrow),\'*\',num2str(ycol)]);
titllist={[\'Vertical Detail V\'];[\'Horizontal Detail H\'];
[\'Diagonal Detail D\']};
pn=2;
for pr=1:N   
     for pc=1:3      
             subplot(N+1,3,pn+2);     
              yt=[ ];       
               yt=uint8(coef{pn});       
              [yrow,ycol]=size(yt);     
              imshow(yt);title([ titllist{pc},num2str(N-pr+1)]);
            xlabel([\'Size : \',num2str(yrow),\'*\',num2str(ycol)]);
              if mod(pn+2,3)==1          
            ylabel([\'Level \',num2str(N-pr+1)]);       
                end      
              pn=pn+1;   
                 end
                end

function [o2sa,f1sa,rsx]=sizcoef(sx,lf,N)^4
o2sa=sx;
f1sa=[0 0];
osx=sx;
for i=1:N  
ot=osx+lf-1;   
sa=floor(ot/2);   
o2sa=[o2sa;ot];   
f1sa=[f1sa;sa];  
osx=sa;
end
rsx=zeros(N+1,2);
rsx(1,:)=sx;
addln=0;
for j=2:N+1                       
aln=floor((addln+lf-1)/2);    
tsx=f1sa(j,:)-aln;   
rsx(j,:)=tsx;   
addln=aln;
end
function y=downspl(x)
N=length(x);      
M=floor(N/2);     
i=1:M;
y(i)=x(2*i)