![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](https://image.slidesharecdn.com/waveletdecompositionprocesscodeinmatlabsolutionwave-230708181101-766c9ca1/85/wavelet-decomposition-process-code-in-matlabSolution-WAVE-pdf-1-320.jpg)
![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);](https://image.slidesharecdn.com/waveletdecompositionprocesscodeinmatlabsolutionwave-230708181101-766c9ca1/85/wavelet-decomposition-process-code-in-matlabSolution-WAVE-pdf-2-320.jpg)
wavelet decomposition process code in matlabSolution WAVE.pdf
- 1. 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
- 2. 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);