Problems with Matlab Projects? You may face many Problems, but do not worry we are ready to solve your Problems. All you need to do is just leave your Comments. We will assure you that you will find a solution to your project along with future tips. On Request we will Mail you Matlab Codes for Registered Members of this site only, at free service...Follow Me.

image compression code


PERFORMANCE EVALUATION OF WAVELET BASED IMAGE COMPRESSION USING EZW & SPIHT ALGORITHMS


1 MAIN MODULE CODE (main.m)

GParm.figure = figure('Color',[1 1 1], ...% white color
'MenuBar','none', ...
'Colormap',gray(256), ...
'Name','EZW & SPIHT WAVELET INTERFACE', ... % figure window name
'Visible','on', 'Resize','off', ... % initially not drawn
'NumberTitle','off', ...
'ShareColor','off', ...
'RendererMode','manual','Renderer','painters', ...
'Units','pixels', ...
'Position',[130 115 800 580], ... % window position
'WindowStyle','normal', ...
'Pointer','arrow');
%%% CLOSE BUTTON
posbox = [130 300 550 75] ;
GParm.close = uicontrol( ...
'Parent',GParm.figure, ...
'Style','Pushbutton', ...
'Background',[.5 .6 .8], ...
'BusyAction','queue','Interruptible','off', ...
'Units','pixel', ...
'String','WAVELET BASED IMAGE COMPRESSION USING EZW & SPIHT ALGORITHM', ...
'Position',posbox,....
'Callback','cont' ...
);
posbox = [150 110 150 30] ;
GParm.close = uicontrol( ...
'Parent',GParm.figure, ...
'Style','Pushbutton', ...
'BusyAction','queue','Interruptible','off', ...
'Units','pixel', ...
'String','CONTINUE', ...
'Position',posbox, ...
'Callback','menu' ...
);
posbox = [500 110 150 30] ;
GParm.close = uicontrol( ...
'Parent',GParm.figure, ...
'Style','Pushbutton', ...
'BusyAction','queue','Interruptible','off', ...
'Units','pixel', ...
'String','CLOSE', ...
'Position',posbox, ...
'Callback','close1' ...
);
2 SUB MODULE CODE (cont.m)
function cont
figure;imshow('info.bmp');
SUB MODULE CODE (menu.m)
function menu
GParm.figure = figure('Color',[1 1 1], ...% white color
'MenuBar','none', ...
'Colormap',gray(256), ...
'Name','EMBEDDED ZERO-TREE WAVELET INTERFACE', ...% figure window name
'Visible','on', 'Resize','off', ... % initially not drawn
'NumberTitle','off', ...
'ShareColor','off', ...
'RendererMode','manual','Renderer','painters', ...
'Units','pixels', ...
'Position',[130 115 800 580], ... % window position
'WindowStyle','normal', ...
'Pointer','arrow');
%%% CLOSE BUTTON
posbox = [630 100 150 30] ;
GParm.close = uicontrol( ...
'Parent',GParm.figure, ...
'Style','Pushbutton', ...
'BusyAction','queue','Interruptible','off', ...
'Units','pixel', ...
'String','CLOSE', ...
'Position',posbox, ...
'Callback','close1' ...
);
posbox = [30 100 150 30] ;
GParm.close = uicontrol( ...
'Parent',GParm.figure, ...
'Style','Pushbutton', ...
'BusyAction','queue','Interruptible','off', ...
'Units','pixel', ...
'String','READ INPUT', ...
'Position',posbox, ...
'Callback','read' ...
) ;
posbox = [230 100 150 30] ;
GParm.close = uicontrol( ...
'Parent',GParm.figure, ...
'Style','Pushbutton', ...
'BusyAction','queue','Interruptible','off', ...
'Units','pixel', ...
'String','EZW IMPLEMENTATION', ...
'Position',posbox, ...
'Callback','menu1(f,map);' ...
);
posbox = [430 100 150 30] ;
GParm.close = uicontrol( ...
'Parent',GParm.figure, ...
'Style','Pushbutton', ...
'BusyAction','queue','Interruptible','off', ...
'Units','pixel', ...
'String','SPIHT IMPLEMENTATION', ...
'Position',posbox, ...
'Callback','menu2(f,map);' ...
);
SUB MODULE CODE (close.m)
function close1
close all;
SUB MODULE CODE (read.m)
function read
GParm.figure=figure('color',[1 1 1], ...%black color
'MenuBar','none', ...
'Colormap',gray(256), ...
'Name','READ INPUT IMAGE OF ANY FORMAT', ... %figure window name
'Visible','on','Resize','off', ... %initially not drawn
'NumberTitle','off', ...
'Units','pixels', ...
'Position',[130 115 795 580], ... % window position
'WindowStyle','normal', ...
'Pointer','arrow');
%FORMAT BUTTONS
posbox=[290 137 200 40];
GParm.close=uicontrol( ...
'parent',GParm.figure, ...
'Style','Pushbutton', ...
'Background',[.9 .7 .7], ...
'BusyAction','queue', ...%'Interruptible','0ff',
'Units','pixel', ...
'String','Figure of bmp format', ...
'Position',posbox, ....
'callback','close;[f,map]=bmp;' ...
);
posbox=[40 175 200 40];
GParm.close=uicontrol( ...
'parent',GParm.figure, ...
'Style','Pushbutton', ...
'Background',[.9 .7 .7], ...
'BusyAction','queue', ...%'Interruptible','0ff',
'Units','pixel', ...
'String','Figure of jpg format', ...
'Position',posbox, ....
'callback','close;[f,map]=jpg;' ...
);
posbox=[540 175 200 40];
GParm.close=uicontrol( ...
'parent',GParm.figure, ...
'Style','Pushbutton', ...
'Background',[.9 .7 .7], ...
'BusyAction','queue', ...%'Interruptible','0ff',
'Units','pixel', ...
'String','Figure of tif format', ...
'Position',posbox, ....
'callback','close;[f,map]=tif;' ...
);
posbox=[40 100 200 40];
GParm.close=uicontrol( ...
'parent',GParm.figure, ...
'Style','Pushbutton', ...
'Background',[.9 .7 .7], ...
'BusyAction','queue', ...%'Interruptible','0ff',
'Units','pixel', ...
'String','Figure of gif format', ...
'Position',posbox, ....
'callback','close;[f,map]=gif;' ...
);
posbox=[540 100 200 40];
GParm.close=uicontrol( ...
'parent',GParm.figure, ...
'Style','Pushbutton', ...
'Background',[.9 .7 .7], ...
'BusyAction','queue', ...%'Interruptible','0ff',
'Units','pixel', ...
'String','Figure of png format', ...
'Position',posbox, ....
'callback','close;[f,map]=png;' ...
);
SUB MODULE CODE (read.m)
function read
GParm.figure=figure('color',[1 1 1], ...%black color
'MenuBar','none', ...
'Colormap',gray(256), ...
'Name','READ INPUT IMAGE OF ANY FORMAT', ... %figure window name
'Visible','on','Resize','off', ... %initially not drawn
'NumberTitle','off', ...
'Units','pixels', ...
'Position',[130 115 795 580], ... % window position
'WindowStyle','normal', ...
'Pointer','arrow');
%FORMAT BUTTONS
posbox=[290 137 200 40];
GParm.close=uicontrol( ...
'parent',GParm.figure, ...
'Style','Pushbutton', ...
'Background',[.9 .7 .7], ...
'BusyAction','queue', ...%'Interruptible','0ff',
'Units','pixel', ...
'String','Figure of bmp format', ...
'Position',posbox, ....
'callback','close;[f,map]=bmp;' ...
);
posbox=[40 175 200 40];
GParm.close=uicontrol( ...
'parent',GParm.figure, ...
'Style','Pushbutton', ...
'Background',[.9 .7 .7], ...
'BusyAction','queue', ...%'Interruptible','0ff',
'Units','pixel', ...
'String','Figure of jpg format', ...
'Position',posbox, ....
'callback','close;[f,map]=jpg;' ...
);
posbox=[540 175 200 40];
GParm.close=uicontrol( ...
'parent',GParm.figure, ...
'Style','Pushbutton', ...
'Background',[.9 .7 .7], ...
'BusyAction','queue', ...%'Interruptible','0ff',
'Units','pixel', ...
'String','Figure of tif format', ...
'Position',posbox, ....
'callback','close;[f,map]=tif;' ...
);
posbox=[40 100 200 40];
GParm.close=uicontrol( ...
'parent',GParm.figure, ...
'Style','Pushbutton', ...
'Background',[.9 .7 .7], ...
'BusyAction','queue', ...%'Interruptible','0ff',
'Units','pixel', ...
'String','Figure of gif format', ...
'Position',posbox, ....
'callback','close;[f,map]=gif;' ...
);
posbox=[540 100 200 40];
GParm.close=uicontrol( ...
'parent',GParm.figure, ...
'Style','Pushbutton', ...
'Background',[.9 .7 .7], ...
'BusyAction','queue', ...%'Interruptible','0ff',
'Units','pixel', ...
'String','Figure of png format', ...
'Position',posbox, ....
'callback','close;[f,map]=png;' ...
);
SUB MODULE CODE (menu1m)
function menu1(s,map);
K=[];
GParm.figure = figure('Color',[1 1 1], ...% black color
'MenuBar','none', ...
'Colormap',gray(256), ...
'Name','EZW WAVELET INTERFACE', ... % figure window name
'Visible','on', 'Resize','off', ... % initially not drawn
'NumberTitle','off', ...
'ShareColor','off', ...
'RendererMode','manual','Renderer','painters', ...
'Units','pixels', ...
'Position',[130 115 800 580], ... % window position
'WindowStyle','normal', ...
'Pointer','arrow');
%%% CLOSE BUTTON
posbox = [630 100 150 30] ;
GParm.close = uicontrol( ...
'Parent',GParm.figure, ...
'Style','Pushbutton', ...
'BusyAction','queue','Interruptible','off', ...
'Units','pixel', ...
'String','CLOSE', ...
'Position',posbox, ...
'Callback','close1' ...
);
posbox = [30 100 150 30] ;
GParm.close = uicontrol( ...
'Parent',GParm.figure, ...
'Style','Pushbutton', ...
'BusyAction','queue','Interruptible','off', ...
'Units','pixel', ...
'String','BLACK & WHITE PROCESS', ...
'Position',posbox, ...
'Callback','mat=ezw(f,map);' ...
);
posbox = [315 100 180 30] ;
GParm.close = uicontrol( ...
'Parent',GParm.figure, ...
'Style','Pushbutton', ...
'BusyAction','queue','Interruptible','off', ...
'Units','pixel', ...
'String','PSNR vs COMPRESSION RATIO', ...
'Position',posbox, ...
'Callback','analysis(mat)' ...
);
7 SUB MODULE CODE (menu2.m)
function menu2(s,map);
K=[];
GParm.figure = figure('Color',[1 1 1], ...% black color
'MenuBar','none', ...
'Colormap',gray(256), ...
'Name','SPIHT WAVELET INTERFACE', ... % figure window name
'Visible','on', 'Resize','off', ... % initially not drawn
'NumberTitle','off', ...
'ShareColor','off', ...
'RendererMode','manual','Renderer','painters', ...
'Units','pixels', ...
'Position',[130 115 800 580], ... % window position
'WindowStyle','normal', ...
'Pointer','arrow');
%%% CLOSE BUTTON
posbox = [630 100 150 30] ;
GParm.close = uicontrol( ...
'Parent',GParm.figure, ...
'Style','Pushbutton', ...
'BusyAction','queue','Interruptible','off', ...
'Units','pixel', ...
'String','CLOSE', ...
'Position',posbox, ...
'Callback','close1' ...
);
posbox = [30 100 150 30] ;
GParm.close = uicontrol( ...
'Parent',GParm.figure, ...
'Style','Pushbutton', ...
'BusyAction','queue','Interruptible','off', ...
'Units','pixel', ...
'String','BLACK & WHITE PROCESS', ...
'Position',posbox, ...
'Callback','mat=spiht(f,map);' ...
);
posbox = [315 100 180 30] ;
GParm.close = uicontrol( ...
'Parent',GParm.figure, ...
'Style','Pushbutton', ...
'BusyAction','queue','Interruptible','off', ...
'Units','pixel', ...
'String','PSNR vs COMPRESSION RATIO', ...
'Position',posbox, ...
'Callback','analysis(mat)' ...
);
SUB MODULE CODE (jpg.m)
function [f,map]=jpg;
[fname,path,findex] = uigetfile('*.jpg','Open a image file :',100,100) ;
[f,map]=imread(fname);
figure('position',[350 370 280 100]);
imshow(f,map);
title('ORIGNAL IMAGE');
return
SUB MODULE CODE (gif.m)
function [f,map]=gif;
[fname,path,findex] = uigetfile('*.gif','Open a image file :',100,100) ;
[f,map]=imread(fname);
figure('position',[350 370 280 100]);
imshow(f,map);
title('ORIGNAL IMAGE');
return
SUB MODULE CODE (gif.m)
function [f,map]=gif;
[fname,path,findex] = uigetfile('*.gif','Open a image file :',100,100) ;
[f,map]=imread(fname);
figure('position',[350 370 280 100]);
imshow(f,map);
title('ORIGNAL IMAGE');
return
SUB MODULE CODE (bmp.m)
function [f,map]=bmp;
[fname,path,findex] = uigetfile('*.bmp','Open a image file :',100,100) ;
[f,map]=imread(fname);
figure('position',[350 370 280 100]);
imshow(f,map);
title('ORIGNAL IMAGE');
return
SUB MODULE CODE (tif.m)
function [f,map]=tif;
[fname,path,findex] = uigetfile('*.tif','Open a image file :',100,100) ;
[f,map]=imread(fname);
figure('position',[350 370 280 100]);
imshow(f,map);
title('ORIGNAL IMAGE');
return
SUB MODULE CODE (png.m)
function [f,map]=png;
[fname,path,findex] = uigetfile('*.png','Open a image file :',100,100) ;
[f,map]=imread(fname);
figure('position',[350 370 280 100]);
imshow(f,map);
title('ORIGNAL IMAGE');
return
SUB MODULE CODE (ezw.m)
function mat=ezw(f,map);
%Conversion RGB type of image to gray scale image
if ndims(f)==3
s=rgb2gray(f);
f=s;
else
s=f;
end
[m n]=size(s); %Resolution of intesity image
%Conversion of any Resolution of image to 128 X 128 resolution image
while m>225 | n>225
if m>225 & n>225
s=dyaddown(s,1,'m'); %Down sampling of an image
elseif m>225
s=dyaddown(s,1,'r');
else
s=dyaddown(s,1,'c');
end
[m n]=size(s);
end
[m n]=size(s);
if m>128 & n>128
if m>153 & n>153
s=s(25:25+127,25:25+127); %To extract required part of image to get required resolution
elseif m>153
s=s(25:25+127,1:128);
elseif n>153
s=s(1:128,25:25+127);
else
s=s(1:128,1:128);
end
%Zero padding to the image matrix to get required resolution
elseif m<128>
m1=128-m;
n1=128-n;
for i=1:m1
for j=1:n
s(m+i,j)=0;
end
end
for i=1:128
for j=1:n1
s(i,n+j)=0;
end
end
else
if m<128
m1=128-m;
for i=1:m1
for j=1:n
s(m+i,j)=0;
end
end
else
n1=128-n;
for i=1:m
for j=1:n1
s(i,n+j)=0;
end
end
end
s=s(1:128,1:128);
end
f1=double(s);
%Conversion of any data class of image matrix to double data class
%Showing Original image of 128 X 128 resolution
figure('position',[100 500 200 70]);
imshow(mat2gray(f1)); %To convert image to Gray scale image
title('ORIGINAL IMAGE');
rr=0.05;
rrs=waitbar(0,'PROCESSING EZW CODING...........');
waitbar(.05);
%Applying Discrete Wavelet Transform upto 7 levels using mother wavelet Haar
[a1,h1,v1,d1]=dwt2(s,'haar');
n=256;
ch1=wcodemat(h1,n); %To quantize the image matrix to n(means 128) levels
cv1=wcodemat(v1,n);
cd1=wcodemat(d1,n);
[a2,h2,v2,d2]=dwt2(a1,'haar');
ch2=wcodemat(h2,n);
cv2=wcodemat(v2,n);
cd2=wcodemat(d2,n);
[a3,h3,v3,d3]=dwt2(a2,'haar');
ch3=wcodemat(h3,n);
cv3=wcodemat(v3,n);
cd3=wcodemat(d3,n);
ca3=wcodemat(a3,n);
s1=[ca3,ch3;cv3,cd3];
s2=[s1,ch2;cv2,cd2];
s3=[s2,ch1;cv1,cd1];
%Showing 3-Level decomposition image
figure('position',[400 500 200 70]);
imshow(uint8(s3));
title('3-LEVEL DECOMPOSITION IMAGE');
[a4,h4,v4,d4]=dwt2(a3,'haar');
[a5,h5,v5,d5]=dwt2(a4,'haar');
[a6,h6,v6,d6]=dwt2(a5,'haar');
[a7,h7,v7,d7]=dwt2(a6,'haar');
s1=[a7,h7;v7,d7];
s2=[s1,h6;v6,d6];
s3=[s2,h5;v5,d5];
s4=[s3,h4;v4,d4];
s5=[s4,h3;v3,d3];
s6=[s5,h2;v2,d2];
X=[s6,h1;v1,d1]; %7-Level Decomposition matrix
figure('position',[400 500 200 70]);
imshow(X);
title('7-LEVEL DECOMPOSITION IMAGE');
mag=abs(X);
mag=uint16(mag);
maximum=max(max(mag));
maximum=double(maximum);
U=log2(maximum);
V=ceil(U);
R=V;
sublist=[];
[xx,yy]=size(X);
global N
A=mapping(xx);
[m,n]=size(A);
k=1;
N=[];
for k=1:m*n;
flag=0;
for i=1:m;
if flag==1;
break;
end
for j=1:n;
if A(i,j)==k;
N=[N;i,j];
flag=1;
break;
end
end
end
end
order=1;
bb=[];
zz=0;
global N;
[m,n]=size(N);
XX=zeros(128);
sublists=[];
for I=1:V
threshold=(2^(V-I));
%Dominant Pass
[D,X,sublist] = dominantpass(X,threshold,sublist);
DD{order}=D;
%Subordinate pass
S= subordinatepass(sublist,threshold);
SS{order}=S;
order=order+1;
if threshold ==1,
break;
end
%%decoding
RR=zeros(128); % reference matrix RR;
[a,b]=size(DD{I});
% dominant pass
i=1; j=1;
while i<=m;
if j>b; % b is the length of dominant pass vector
break;
end
if RR(N(i,1),N(i,2))==1;
i=i+1;
else
if DD{I}(j)=='p';
XX(N(i,1),N(i,2))=threshold;
RR(N(i,1),N(i,2))=1; %reference matrix =1;
sublists=[sublists;N(i,1),N(i,2)];
end
if DD{I}(j)=='n';
XX(N(i,1),N(i,2))=-threshold;
RR(N(i,1),N(i,2))=1; %reference matrix =1;
sublists=[sublists;N(i,1),N(i,2)];
end
if DD{I}(j)=='z';
RR(N(i,1),N(i,2))=1; %reference matrix =1;
end
if DD{I}(j)=='t';
RR(N(i,1),N(i,2))=1; %reference matrix =1;
RR=checkchildren(i,RR);% all zerotree's descendants are set to 1.
end
i=i+1;j=j+1;
end
end
[xx,yy]=size(sublists);
threshold=threshold/2;
for i=1:xx;
if I==V-1 |threshold==1;
break;
end
if SS{I}(i)==1;
if XX(sublists(i,1),sublists(i,2))>0;
XX(sublists(i,1),sublists(i,2))= XX(sublists(i,1),sublists(i,2))+ threshold;
else
XX(sublists(i,1),sublists(i,2))= XX(sublists(i,1),sublists(i,2))-threshold;
end
end
end
rec=[];
rec=XX;
ra7=rec(1,1);
rh7=rec(1,2);
rv7=rec(2,1);
rd7=rec(2,2);
rh6=rec(1:2,3:4);
rv6=rec(3:4,1:2);
rd6=rec(3:4,3:4);
rh5=rec(1:4,5:8);
rv5=rec(5:8,1:4);
rd5=rec(5:8,5:8);
rh4=rec(1:8,9:16);
rv4=rec(9:16,1:8);
rd4=rec(9:16,9:16);
rh3=rec(1:16,17:32);
rv3=rec(17:32,1:16);
rd3=rec(17:32,17:32);
rh2=rec(1:32,33:64);
rv2=rec(33:64,1:32);
rd2=rec(33:64,33:64);
rh1=rec(1:64,65:128);
rv1=rec(65:128,1:64);
rd1=rec(65:128,65:128);
rs7=idwt2(ra7,rh7,rv7,rd7,'haar');
rs6=idwt2(rs7,rh6,rv6,rd6,'haar');
rs5=idwt2(rs6,rh5,rv5,rd5,'haar');
rs4=idwt2(rs5,rh4,rv4,rd4,'haar');
rs3=idwt2(rs4,rh3,rv3,rd3,'haar');
rs2=idwt2(rs3,rh2,rv2,rd2,'haar');
rs1=idwt2(rs2,rh1,rv1,rd1,'haar');
rs1=fix(rs1);
figure('position',[700 500 200 70]);
imshow(mat2gray(rs1))
title('COMPRESSED IMAGE');
OrigSize=131072;
qq=(size(S,2));
aa=((qq/8)+(size(D,2)));
bb=[bb aa];
cc=size(bb,2);
if cc==1
zz=bb(1);
else
for W=1:cc-1
zz=zz+bb(W);
end
end
RecSize=zz;
RecSize=fix(RecSize);
CompRatio=OrigSize/RecSize;
err=0;
S=logical(S);
for i=1:128
for j=1:128
diff=(rs1(i,j)-f1(i,j))^2;
err=err+diff;
end
end
mse=err/(128*128);
psnr=10*log10(255*255/mse);
mat=[];
mat=[mat;psnr,CompRatio];
sprintf('\t%d\t\t%f:%d',psnr,CompRatio,1)
rr=rr+0.05;
waitbar(rr);
end
waitbar(1);
close(rrs);
return
SUB MODULE CODE (mapping.m)
function A = mapping(n)
% To create a Morton Scan order matrix
if n == 2
A = [1 2; 3 4];
else
B = mapping(n/2);
A = [B B+(n/2)^2; B+(n/2)^2*2 B+(n/2)^2*3];
end
SUB MODULE CODE (dominantpass.m)
function [D,X,sublist] = dominantpass(X,threshold,sublist)
% Dominant pass function
D=[];
global N
[m,n]=size(X);
% X is the coefficients matrix
R=zeros(m); % matrix R is a reference matrix, same size as X; '0' means
%this coefficient is not a descendant from zerotree root;
[a,b]=size(N);
if abs(X(1,1))>=threshold; % X(1,1) is DC coefficient
sublist=[sublist, abs(X(1,1))];
if X(1,1)>0;
D=[D,'p'];
else D=[D,'n'];
end
X(1,1)=0;
else D=[D,'z'];
end
for k=2:4
if abs(X(N(k,1),N(k,2)))>=threshold,
sublist=[sublist, abs(X(N(k,1),N(k,2)))];
% append this significant coefficient to the subordinate list;
if X(N(k,1),N(k,2))>0, % determine the sign
D=[D,'p']; % >0,assign a "p"
else
D=[D,'n'];% <0,assign>
end
X(N(k,1),N(k,2))=0;
% the significent coefficients is replaced by a '0' in the coefficients matrix
else
% 2,3,4 has no parents,just check its descendants.
result = checkdescendants1( k,X,threshold,0);
if result==1,
D=[D,'z'];
else
D=[D,'t'];
R(N(k,1),N(k,2))=1;
end
end
end
for k=5:a;
if abs(X(N(k,1),N(k,2)))>=threshold,
sublist=[sublist, abs(X(N(k,1),N(k,2)))];
if X(N(k,1),N(k,2))>0, % determine the sign
D=[D,'p']; % >0,assign a "p"
else D=[D,'n'];% <0,assign>
end
X(N(k,1),N(k,2))=0;
% else zerotree = checkancestors( k,R,0);
else zerotree = checkancestors1( k,R,0);
% check coefficient is descended from zerotree root?
if zerotree==0, % Not descends from zerotree root
result = checkdescendants1( k,X,threshold,0);
% Check its has significant descendants?
if result==1,
D=[D,'z']; % isolated zero
else D=[D,'t'];% zerotree
R(N(k,1),N(k,2))=1;
% if zerotree, reference matrix coefficient=1
end
end
end
end
SUB MODULE CODE (checkancestors1.m)
function zerotree = checkancestors1(j,R,zerotree)
% zerotree=1 means coefficient is descended from zerotree root
% initial set zerotree=0
global N
i = floor((j+3)/4 );
while i>1,
if R(N(i,1),N(i,2))==1;
zerotree=1;
break;
end;
i = floor((i+3)/4 );
end
SUB MODULE CODE (checkdescendants.m)
function result = checkdescendants1(j,X,threshold,result)
% initial set result=0
% if the result=1, means that a coefficient has at least
% 1 significant descendant.
global N
[m,n]=size(N);
for i=(4*j-3):4*j;
if result==1 | i>m,
break;
end
if abs(X(N(i,1),N(i,2)))>=threshold
result=1;
break;
else
result=checkdescendants1(i,X,threshold,result);
end;
end;
SUB MODULE CODE (checkchildren.m)
function [RR]=checkchildren(j,RR)
% if a symbol 't' is encounted, then make all its descendants in reference
% matrix RR's components equal 1---ZEROTREES
global N
[m,n]=size(N);
for i=(4*j-3):4*j;
if i<=m,
RR(N(i,1),N(i,2))=1;
[RR]=checkchildren(i,RR);
end
end;
SUB MODULE CODE (subordinatepass.m)
function S = subordinatepass(sublist,threshold)
S=[];
[m,n]=size(sublist);
for i=1:n;
if sublist(1,i)>=threshold
%if sublist(1,i)>=threshold,
S=[S,1];
sublist(1,i)=sublist(1,i)-threshold;
else
S=[S,0];
end
end
SUB MODULE CODE (spiht.m)
function mat=spiht(f,map);
%Conversion RGB type of image to gray scale image
if ndims(f)==3
s=rgb2gray(f);
f=s;
else
s=f;
end
[m n]=size(s); %Resolution of intesity image
%Conversion of any Resolution of image to 128 X 128 resolution image
while m>225 | n>225
if m>225 & n>225
s=dyaddown(s,1,'m'); %Down sampling of an image
elseif m>225
s=dyaddown(s,1,'r');
else
s=dyaddown(s,1,'c');
end
[m n]=size(s);
end
[m n]=size(s);
if m>128 & n>128
if m>153 & n>153
s=s(25:25+127,25:25+127); %To extract required part of image to get required resolution
elseif m>153
s=s(25:25+127,1:128);
elseif n>153
s=s(1:128,25:25+127);
else
s=s(1:128,1:128);
end
%Zero padding to the image matrix to get required resolution
elseif m<128>
m1=128-m;
n1=128-n;
for i=1:m1
for j=1:n
s(m+i,j)=0;
end
end
for i=1:128
for j=1:n1
s(i,n+j)=0;
end
end
else
if m<128
m1=128-m;
for i=1:m1
for j=1:n
s(m+i,j)=0;
end
end
else
n1=128-n;
for i=1:m
for j=1:n1
s(i,n+j)=0;
end
end
end
s=s(1:128,1:128);
end
f1=double(s);
%Conversion of any data class of image matrix to double data class
%Showing Original image of 128 X 128 resolution
figure('position',[100 500 200 70]);
imshow(mat2gray(f1)); %To convert image to Gray scale image
title('ORIGINAL IMAGE');
wm=0.2;
wc=waitbar(0,'PROCESSING SPIHT CODING...........');
waitbar(.2);
%Applying Discrete Wavelet Transform upto 7 levels using mother wavelet Haar
[a1,h1,v1,d1]=dwt2(s,'haar');
n=256;
ch1=wcodemat(h1,n); %To quantize the image matrix to n(means 128) levels
cv1=wcodemat(v1,n);
cd1=wcodemat(d1,n);
[a2,h2,v2,d2]=dwt2(a1,'haar');
ch2=wcodemat(h2,n);
cv2=wcodemat(v2,n);
cd2=wcodemat(d2,n);
[a3,h3,v3,d3]=dwt2(a2,'haar');
ch3=wcodemat(h3,n);
cv3=wcodemat(v3,n);
cd3=wcodemat(d3,n);
ca3=wcodemat(a3,n);
s1=[ca3,ch3;cv3,cd3];
s2=[s1,ch2;cv2,cd2];
s3=[s2,ch1;cv1,cd1];
%Showing 3-Level decomposition image
figure('position',[400 500 200 70]);
imshow(uint8(s3));
title('3-LEVEL DECOMPOSITION IMAGE');
[a4,h4,v4,d4]=dwt2(a3,'haar');
[a5,h5,v5,d5]=dwt2(a4,'haar');
[a6,h6,v6,d6]=dwt2(a5,'haar');
[a7,h7,v7,d7]=dwt2(a6,'haar');
s1=[a7,h7;v7,d7];
s2=[s1,h6;v6,d6];
s3=[s2,h5;v5,d5];
s4=[s3,h4;v4,d4];
s5=[s4,h3;v3,d3];
s6=[s5,h2;v2,d2];
dec=[s6,h1;v1,d1]; %7-Level Decomposition matrix
figure('position',[400 500 200 70]);
imshow(dec);
%SPIHT algorithm implementation
%To know Threshold value
mag=abs(dec);
mag=uint16(mag);
maximum=max(max(mag));
maximum=double(maximum);
n=log2(maximum);
n=floor(n);
n1=n;
LSP=[]; %List of Significant Pixels
LIP=[]; %List of Insignificant Pixels
LIS=[]; %List of Insignificant sets
LIP1=[];
rec=zeros(128,128); %Reconstruction matrix
a=[]; %Transmission Bit Stream
mat=[]; %To store PSNR & Compression ratio
sprintf('\t%s\t\t%s','PSNR','Compression Ratio')
for i=1:n1
threshold=2^(n-i+1);
%Refinement pass
if i~=1
x=size(LSP,1);
for j=1:x
m=LSP(j,:);
p=m(1);
q=m(2);
v=abs(dec(p,q));
b1=dec2bin(v); %To convert decimal number to binary number
n2=size(b1,2);
n3=n-i+1;
n4=n2-n3;
b2=b1(n4);
b2=uint8(b2);
%To find Nth MSB bit
if b2==48
b3=0;
else
b3=1;
end
a=[a b3];
if rec(p,q)>=0
rec(p,q)=rec(p,q)+b3*2^(n-i+1);
else
rec(p,q)=rec(p,q)-b3*2^(n-i+1);
end
end
end
%Sorting pass
%Initialization of LIP & LIS
if i==1
LIP=[1,1;1,2;2,1;2,2];
LIS=[1,2,1;2,1,1;2,2,1;0,0,0];
end
x=size(LIP,1);
LIP1=LIP;
k=0;
for j=1:x
k=k+1;
m=LIP(j,:);
p=m(1);
q=m(2);
v=dec(p,q);
ab=abs(v);
if ab
a=[a 0];
else
a=[a 1];
x1=size(LIP1,1);
LIP1=[LIP1(1:k-1,:);LIP1(k+1:x1,:)];
k=k-1;
LSP=[LSP;p,q];
add=2^(n-i+1);
if v<0
a=[a 1];
add=-add;
else
a=[a 0];
end
rec(p,q)=add;
end
end
LIP=LIP1;
j=1;
while max(LIS(j,:))~=0
m=LIS(j,:);
p=m(1);
q=m(2);
t=m(3);
D=[2*p-1,2*q-1;2*p,2*q-1;2*p-1,2*q;2*p,2*q]; %Descedents or Children
D1=[];
for k=1:4
m1=D(k,:);
p1=m1(1);
q1=m1(2);
v=dec(p1,q1);
D1=[D1;v];
end
if t==2
L=[];
L1=[];
for k=1:4
m1=D(k,:);
p1=m1(1);
q1=m1(2);
L=[L;2*p1-1,2*q1-1;2*p1,2*q1-1;2*p1-1,2*q1;2*p1,2*q1]; %Grand children
end
for k=1:16
m1=L(k,:);
p1=m1(1);
q1=m1(2);
v=dec(p1,q1);
L1=[L1;v];
end
if max(abs(L1))
a=[a 0];
else
a=[a 1];
x=size(LIS,1);
LIS=[LIS(1:j-1,:);LIS(j+1:x,:)];
j=j-1;
LIS=LIS(1:x-2,:);
for k=1:4
m1=D(k,:);
p1=m1(1);
q1=m1(2);
LIS=[LIS;p1,q1,1];
end
x=size(LIS,1);
LIS=[LIS;0,0,0];
end
else
if max(abs(D1))
a=[a 0];
else
a=[a 1];
x=size(LIS,1);
LIS=[LIS(1:j-1,:);LIS(j+1:x,:)];
j=j-1;
if p<33>
LIS(x-1,:)=[p,q,2];
LIS(x,:)=[0,0,0];
end
for k=1:4
v=D1(k);
m1=D(k,:);
p1=m1(1);
q1=m1(2);
if abs(v)
a=[a 0];
LIP=[LIP;p1,q1];
else
a=[a 1];
LSP=[LSP;p1,q1];
add=2^(n-i+1);
if v<0
a=[a 1];
add=-add;
else
a=[a 0];
end
rec(p1,q1)=add;
end
end
end
end
j=j+1;
end
ra7=rec(1,1);
rh7=rec(1,2);
rv7=rec(2,1);
rd7=rec(2,2);
rh6=rec(1:2,3:4);
rv6=rec(3:4,1:2);
rd6=rec(3:4,3:4);
rh5=rec(1:4,5:8);
rv5=rec(5:8,1:4);
rd5=rec(5:8,5:8);
rh4=rec(1:8,9:16);
rv4=rec(9:16,1:8);
rd4=rec(9:16,9:16);
rh3=rec(1:16,17:32);
rv3=rec(17:32,1:16);
rd3=rec(17:32,17:32);
rh2=rec(1:32,33:64);
rv2=rec(33:64,1:32);
rd2=rec(33:64,33:64);
rh1=rec(1:64,65:128);
rv1=rec(65:128,1:64);
rd1=rec(65:128,65:128);
rs7=idwt2(ra7,rh7,rv7,rd7,'haar');
rs6=idwt2(rs7,rh6,rv6,rd6,'haar');
rs5=idwt2(rs6,rh5,rv5,rd5,'haar');
rs4=idwt2(rs5,rh4,rv4,rd4,'haar');
rs3=idwt2(rs4,rh3,rv3,rd3,'haar');
rs2=idwt2(rs3,rh2,rv2,rd2,'haar');
rs1=idwt2(rs2,rh1,rv1,rd1,'haar');
rs1=fix(rs1);
figure('position',[700 500 200 70]);
imshow(mat2gray(rs1));
title('COMPRESSED IMAGE');
OrigSize=131072;
RecSize=size(a,2);
RecSize=RecSize/8;
RecSize=fix(RecSize);
CompRatio=OrigSize/RecSize;
CompRatio=ceil(CompRatio);
err=0;
a=logical(a);
for i=1:128
for j=1:128
diff=(rs1(i,j)-f1(i,j))^2;
err=err+diff;
end
end
mse=err/(128*128);
psnr=10*log10(255*255/mse);
mat=[mat;psnr,CompRatio];
sprintf('\t%d\t\t%d:%d',psnr,CompRatio,1)
wm=wm+0.05;
waitbar(wm);
end
waitbar(1);
close(wc);
return
4.21 SUB MODULE CODE (analysis.m)
function analysis(mat)
psnr=mat(:,1);
cnr=mat(:,2);
x=size(psnr,1);
psnr=reshape(psnr,1,x);
cnr=reshape(cnr,1,x);
figure;plot(cnr,psnr);
xlabel('Compression Ratio');
ylabel('PSNR(in dB)');
Posted by bharadwaj
Labels:

2 comments:

Unknown said...

I run this code but i keep getting an error message, as in :: (??? Error using ==> figure
Invalid property found.
Object Name : figure
Property Name : 'ShareColor'.

Error in ==> main at 1
GParm.figure = figure('Color',[1 1 1], ...% white color)

Please help me out. Thanks

July 1, 2014 at 7:44 AM  
Unknown said...

I run this code but i keep getting an error response, as in:: (??? Error using ==> figure
Invalid property found.
Object Name : figure
Property Name : 'ShareColor'.

Error in ==> main at 1
GParm.figure = figure('Color',[1 1 1], ...% white color)

Please help me out with the errors. You can send it to my E-mail (necs.mobolaji@gmail.com) Thanks

July 1, 2014 at 7:46 AM  

Post a Comment

Subscribe to: Post Comments (Atom)

Recent Comments

Popular Matlab Topics

  • Voice Conversion in matlab
    Introduction to Methods for Voice Conversion Speech processing is currently a key focus for many researchers in the area of DSP. In this pr...
  • Xray Imaging using Pseudo coloring Matlab
    ABSTRACT             Samples with the same color or appearance in the visible region of the electromagnetic spectrum, which are therefore i...
  • Radar Jammer in Matlab
    Project Objective The basic objective was to design a program in Matlab that simulates how police radar guns work and various methods of ho...
  • Programming in LabVIEW MathScript-For Loops
    The For Loop The for loop is one way to make LABVIEW MATHSCRIPT repeat a series of computations using different values. The for loop has t...
  • Hand Gesture Recognition Using Neural Network
    Chapter 1: Introduction 1.1 Introduction For the past few years, the common input computer devices did not change a lot. This means, the...
  • Iris Detection matlab
    Motivation behind Iris Detection Passwords are Bad Pin numbers, email passwords, credit card numbers, and protected premises access numbers...
  • 600 MATLAB Source Codes
    MATLAB Source Codes Click on below links args , a program which shows how to count and report command line arguments; arpack , a library...
  • Permanent Magnet Synchronous Motor Control code
    Introduction The objective of this module is to use the Embedded Target for the C2000, and the Digital Motor Control Library (DMClib) to co...
  • ANALYSIS OF FINITE WORD-LENGTH EFFECTS
    ANALYSIS OF FINITE WORD-LENGTH EFFECTS The MATLAB function a2dR given below can be used to generate the decimal equivalent of the binary r...
  • Filtering Periodic Signals using matlab
    Introduction In this lab, we will look at the effect of filtering signals using a frequency domain implementation of an LTI system, i.e., m...

Share your knowledge - help others

Submit your Papers, Projects & Presentations, Distribute to all over the world through this site for free on your name.Send your data to bharadwaj@matlabcodes.com .Kindly Request your code only through comments, not through mail.Other requests are ignored.

Crazy over Matlab Projects ? - Join Now - Follow Me

Sites U Missed to Visit ?

Related Posts Plugin for WordPress, Blogger...

Latest Articles

Loading...

Special Search For Matlab Projects

MATLAB PROJECTS

counter

Bharadwaj. Powered by Blogger.