please can you help , how to solve and apply the image compression for embedded image using arithmetic encoding in matlab?

clc;
clear all;
close all;
I1=imread('u1.jpg');
I=imresize(I1,[256 256]);
[r c]=size(I);
bs=32; % Block Size (64x64)

nob=(r/bs)*(c/bs); % Total number of 16 Blocks(4*4)

% Dividing the image into 64*64 Blocks
kk=0;
for i=1:(r/bs)
for j=1:(c/bs)
Block(:,:,kk+j)=I((bs*(i-1)+1:bs*(i-1)+bs),(bs*(j-1)+1:bs*(j-1)+bs));
end
kk=kk+(r/bs);
end
% -------------------------------------------------------------------------
% -------------------------------------------------------------------------
% Display all the blocks
subplot(8,8,1);
imshow(Block(:,:,1));
subplot(8,8,2);
imshow(Block(:,:,2));
subplot(8,8,3);
imshow(Block(:,:,3));
subplot(8,8,4);
imshow(Block(:,:,4));
subplot(8,8,5);
imshow(Block(:,:,5));
subplot(8,8,6);
imshow(Block(:,:,6));
subplot(8,8,7);
imshow(Block(:,:,7));
subplot(8,8,8);
imshow(Block(:,:,8));
subplot(8,8,9);
imshow(Block(:,:,9));
subplot(8,8,10);
imshow(Block(:,:,10));
subplot(8,8,11);
imshow(Block(:,:,11));
subplot(8,8,12);
imshow(Block(:,:,12));
subplot(8,8,13);
imshow(Block(:,:,13));
subplot(8,8,14);
imshow(Block(:,:,14));
subplot(8,8,15);
imshow(Block(:,:,15));
subplot(8,8,16);
imshow(Block(:,:,16));
subplot(8,8,17);
imshow(Block(:,:,17));
subplot(8,8,18);
imshow(Block(:,:,18));
subplot(8,8,19);
imshow(Block(:,:,19));
subplot(8,8,20);
imshow(Block(:,:,20));
subplot(8,8,21);
imshow(Block(:,:,21));
subplot(8,8,22);
imshow(Block(:,:,22));
subplot(8,8,23);
imshow(Block(:,:,23));
subplot(8,8,24);
imshow(Block(:,:,24));
subplot(8,8,25);
imshow(Block(:,:,25));
subplot(8,8,26);
imshow(Block(:,:,26));
subplot(8,8,27);
imshow(Block(:,:,27));
subplot(8,8,28);
imshow(Block(:,:,28));
subplot(8,8,29);
imshow(Block(:,:,29));
subplot(8,8,30);
imshow(Block(:,:,30));
subplot(8,8,31);
imshow(Block(:,:,31));
subplot(8,8,32);
imshow(Block(:,:,32));
subplot(8,8,33);
imshow(Block(:,:,33));
subplot(8,8,34);
imshow(Block(:,:,34));
subplot(8,8,35);
imshow(Block(:,:,35));
subplot(8,8,36);
imshow(Block(:,:,36));
subplot(8,8,37);
imshow(Block(:,:,37));
subplot(8,8,38);
imshow(Block(:,:,38));
subplot(8,8,39);
imshow(Block(:,:,39));
subplot(8,8,40);
imshow(Block(:,:,40));
subplot(8,8,41);
imshow(Block(:,:,41));
subplot(8,8,42);
imshow(Block(:,:,42));
subplot(8,8,43);
imshow(Block(:,:,43));
subplot(8,8,44);
imshow(Block(:,:,44));
subplot(8,8,45);
imshow(Block(:,:,45));
subplot(8,8,46);
imshow(Block(:,:,46));
subplot(8,8,47);
imshow(Block(:,:,47));
subplot(8,8,48);
imshow(Block(:,:,48));
subplot(8,8,49);
imshow(Block(:,:,49));
subplot(8,8,50);
imshow(Block(:,:,50));
subplot(8,8,51);
imshow(Block(:,:,51));
subplot(8,8,52);
imshow(Block(:,:,52));
subplot(8,8,53);
imshow(Block(:,:,53));
subplot(8,8,54);
imshow(Block(:,:,54));
subplot(8,8,55);
imshow(Block(:,:,55));
subplot(8,8,56);
imshow(Block(:,:,56));
subplot(8,8,57);
imshow(Block(:,:,57));
subplot(8,8,58);
imshow(Block(:,:,58));
subplot(8,8,59);
imshow(Block(:,:,59));
subplot(8,8,60);
imshow(Block(:,:,60));
subplot(8,8,61);
imshow(Block(:,:,61));
subplot(8,8,62);
imshow(Block(:,:,62));
subplot(8,8,63);
imshow(Block(:,:,63));
subplot(8,8,64);
imshow(Block(:,:,64));
<pre class="language-matlab">% ----------------------------------------------------------------------------------------------------------------------
% ----------------------------------------------------------------------------------------------------------------------
% Accessing individual Blocks
% figure;imshow(Block(:,:,1));% This shows u the first 64*64 Block in a figure window
% figure;imshow(Block(:,:,2));
% figure;imshow(Block(:,:,3));
% figure;imshow(Block(:,:,4));
% figure;imshow(Block(:,:,5));
% figure;imshow(Block(:,:,6));
% figure;imshow(Block(:,:,7));
%figure;imshow(Block(:,:,8));
% figure;imshow(Block(:,:,9));
% figure;imshow(Block(:,:,10));
% figure;imshow(Block(:,:,11));
% figure;imshow(Block(:,:,12));
% figure;imshow(Block(:,:,13));
% figure;imshow(Block(:,:,14));
% figure;imshow(Block(:,:,15));
% figure;imshow(Block(:,:,16));
% figure;imshow(Block(:,:,17));
% figure;imshow(Block(:,:,18));
% figure;imshow(Block(:,:,19));
% figure;imshow(Block(:,:,20));
% figure;imshow(Block(:,:,21));
% figure;imshow(Block(:,:,22));
% figure;imshow(Block(:,:,23));
% figure;imshow(Block(:,:,24));
% figure;imshow(Block(:,:,25));
% figure;imshow(Block(:,:,26));
% figure;imshow(Block(:,:,27));
% figure;imshow(Block(:,:,28));
% figure;imshow(Block(:,:,29));
% figure;imshow(Block(:,:,30));
% figure;imshow(Block(:,:,31));
% figure;imshow(Block(:,:,32));
% figure;imshow(Block(:,:,33));
% figure;imshow(Block(:,:,34));
% figure;imshow(Block(:,:,35));
% figure;imshow(Block(:,:,36));
% figure;imshow(Block(:,:,37));
% figure;imshow(Block(:,:,38));
% figure;imshow(Block(:,:,39));
% figure;imshow(Block(:,:,40));
% figure;imshow(Block(:,:,41));
% figure;imshow(Block(:,:,42));
% figure;imshow(Block(:,:,43));
% figure;imshow(Block(:,:,44));
% figure;imshow(Block(:,:,45));
% figure;imshow(Block(:,:,46));
% figure;imshow(Block(:,:,47));
% figure;imshow(Block(:,:,48));
% figure;imshow(Block(:,:,49));
% figure;imshow(Block(:,:,50));
% figure;imshow(Block(:,:,51));
% figure;imshow(Block(:,:,52));
% figure;imshow(Block(:,:,53));
% figure;imshow(Block(:,:,54));
% figure;imshow(Block(:,:,55));
% figure;imshow(Block(:,:,56));
% figure;imshow(Block(:,:,57));
% figure;imshow(Block(:,:,58));
% figure;imshow(Block(:,:,59));
% figure;imshow(Block(:,:,60));
% figure;imshow(Block(:,:,61));
% figure;imshow(Block(:,:,62));
% figure;imshow(Block(:,:,63));
% figure;imshow(Block(:,:,64));
% --------------------------------------------------------------------------------------------------------------------
% -------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA1,CH1,CV1,CD1] = lwt2(double(Block(:,:,1)),LS);
imlwt1=[abs(CA1) abs(CH1); abs(CV1) abs(CD1)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt1, []);
% title('2-D WAVELET');
CA1= uint8(abs(CA1));
% figure,imshow(CA1);
% title('RESULTANT IMAGE-1');
%----------------------------------------------------------------------
%----------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA2,CH2,CV2,CD2] = lwt2(double(Block(:,:,2)),LS);
imlwt2=[abs(CA2) abs(CH2); abs(CV2) abs(CD2)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt2, []);
% title('2-D WAVELET');
CA2= uint8(abs(CA2));
% figure,imshow(CA2);
% title('RESULTANT IMAGE-2');
%----------------------------------------------------------------------------
%------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA3,CH3,CV3,CD3] = lwt2(double(Block(:,:,3)),LS);
imlwt3=[abs(CA3) abs(CH3); abs(CV3) abs(CD3)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt3, []);
% title('2-D WAVELET');
CA3= uint8(abs(CA3));
% figure,imshow(CA3);
% title('RESULTANT IMAGE-3');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA4,CH4,CV4,CD4] = lwt2(double(Block(:,:,4)),LS);
imlwt4=[abs(CA4) abs(CH4); abs(CV4) abs(CD4)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt4, []);
% title('2-D WAVELET');
CA4= uint8(abs(CA4));
% figure,imshow(CA4);
% title('RESULTANT IMAGE-4');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA5,CH5,CV5,CD5] = lwt2(double(Block(:,:,5)),LS);
imlwt5=[abs(CA5) abs(CH5); abs(CV5) abs(CD5)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt5, []);
% title('2-D WAVELET');
CA5= uint8(abs(CA5));
% figure,imshow(CA5);
% title('RESULTANT IMAGE-5');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA6,CH6,CV6,CD6] = lwt2(double(Block(:,:,6)),LS);
imlwt6=[abs(CA6) abs(CH6); abs(CV6) abs(CD6)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt6, []);
% title('2-D WAVELET');
CA6= uint8(abs(CA6));
% figure,imshow(CA6);
% title('RESULTANT IMAGE-6');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA7,CH7,CV7,CD7] = lwt2(double(Block(:,:,7)),LS);
imlwt7=[abs(CA7) abs(CH7); abs(CV7) abs(CD7)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt7, []);
% title('2-D WAVELET');
CA7= uint8(abs(CA7));
% figure,imshow(CA7);
% title('RESULTANT IMAGE-7');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA8,CH8,CV8,CD8] = lwt2(double(Block(:,:,8)),LS);
imlwt8=[abs(CA8) abs(CH8); abs(CV8) abs(CD8)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt8, []);
% title('2-D WAVELET');
CA8= uint8(abs(CA8));
% figure,imshow(CA8);
% title('RESULTANT IMAGE-8');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA9,CH9,CV9,CD9] = lwt2(double(Block(:,:,9)),LS);
imlwt9=[abs(CA9) abs(CH9); abs(CV9) abs(CD9)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt9, []);
% title('2-D WAVELET');
CA9= abs(CA9);
% figure,imshow(CA9);
% title('RESULTANT IMAGE-9');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA10,CH10,CV10,CD10] = lwt2(double(Block(:,:,10)),LS);
imlwt10=[abs(CA10) abs(CH10); abs(CV10) abs(CD10)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt10, []);
% title('2-D WAVELET');
CA10= uint8(abs(CA10));
% figure,imshow(CA10);
% title('RESULTANT IMAGE-10');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA11,CH11,CV11,CD11] = lwt2(double(Block(:,:,11)),LS);
imlwt11=[abs(CA11) abs(CH11); abs(CV11) abs(CD11)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt11, []);
% title('2-D WAVELET');
CA11= uint8(abs(CA11));
% figure,imshow(CA11);
% title('RESULTANT IMAGE-11');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA12,CH12,CV12,CD12] = lwt2(double(Block(:,:,12)),LS);
imlwt12=[abs(CA12) abs(CH12); abs(CV12) abs(CD12)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt12, []);
% title('2-D WAVELET');
CA12= uint8(abs(CA12));
% figure,imshow(CA12);
% title('RESULTANT IMAGE-12');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA13,CH13,CV13,CD13] = lwt2(double(Block(:,:,13)),LS);
imlwt13=[abs(CA13) abs(CH13); abs(CV13) abs(CD13)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt13, []);
% title('2-D WAVELET');
CA13= uint8(abs(CA13));
% figure,imshow(CA13);
% title('RESULTANT IMAGE-13');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA14,CH14,CV14,CD14] = lwt2(double(Block(:,:,14)),LS);
imlwt14=[abs(CA14) abs(CH14); abs(CV14) abs(CD14)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt14, []);
% title('2-D WAVELET');
CA14= uint8(abs(CA14));
% figure,imshow(CA14);
% title('RESULTANT IMAGE-14');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA15,CH15,CV15,CD15] = lwt2(double(Block(:,:,15)),LS);
imlwt15=[abs(CA15) abs(CH15); abs(CV15) abs(CD15)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt15, []);
% title('2-D WAVELET');
CA15= uint8(abs(CA15));
% figure,imshow(CA15);
% title('RESULTANT IMAGE-15');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA16,CH16,CV16,CD16] = lwt2(double(Block(:,:,16)),LS);
imlwt16=[abs(CA16) abs(CH16); abs(CV16) abs(CD16)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt16, []);
% title('2-D WAVELET');
CA16= uint8(abs(CA16));
% figure,imshow(CA16);
% title('RESULTANT IMAGE-16');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA17,CH17,CV17,CD17] = lwt2(double(Block(:,:,17)),LS);
imlwt17=[abs(CA17) abs(CH17); abs(CV17) abs(CD17)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt17, []);
% title('2-D WAVELET');
CA17= uint8(abs(CA17));
% figure,imshow(CA17);
% title('RESULTANT IMAGE-17');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA18,CH18,CV18,CD18] = lwt2(double(Block(:,:,18)),LS);
imlwt18=[abs(CA18) abs(CH18); abs(CV18) abs(CD18)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt18, []);
% title('2-D WAVELET');
CA18= uint8(abs(CA18));
% figure,imshow(CA18);
% title('RESULTANT IMAGE-18');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA19,CH19,CV19,CD19] = lwt2(double(Block(:,:,19)),LS);
imlwt19=[abs(CA19) abs(CH19); abs(CV19) abs(CD19)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt19, []);
% title('2-D WAVELET');
CA19= uint8(abs(CA19));
% figure,imshow(CA19);
% title('RESULTANT IMAGE-19');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA20,CH20,CV20,CD20] = lwt2(double(Block(:,:,20)),LS);
imlwt20=[abs(CA20) abs(CH20); abs(CV20) abs(CD20)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt20, []);
% title('2-D WAVELET');
CA20= uint8(abs(CA20));
% figure,imshow(CA20);
% title('RESULTANT IMAGE-20');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA21,CH21,CV21,CD21] = lwt2(double(Block(:,:,21)),LS);
imlwt21=[abs(CA21) abs(CH21); abs(CV21) abs(CD21)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt21, []);
% title('2-D WAVELET');
CA21= uint8(abs(CA21));
% figure,imshow(CA21);
% title('RESULTANT IMAGE-21');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA22,CH22,CV22,CD22] = lwt2(double(Block(:,:,22)),LS);
imlwt22=[abs(CA22) abs(CH22); abs(CV22) abs(CD22)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt22, []);
% title('2-D WAVELET');
CA22= uint8(abs(CA22));
% figure,imshow(CA22);
% title('RESULTANT IMAGE-22');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA23,CH23,CV23,CD23] = lwt2(double(Block(:,:,23)),LS);
imlwt23=[abs(CA23) abs(CH23); abs(CV23) abs(CD23)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt23, []);
% title('2-D WAVELET');
CA23= uint8(abs(CA23));
% figure,imshow(CA23);
% title('RESULTANT IMAGE-23');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA24,CH24,CV24,CD24] = lwt2(double(Block(:,:,24)),LS);
imlwt24=[abs(CA24) abs(CH24); abs(CV24) abs(CD24)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt24, []);
% title('2-D WAVELET');
CA24= uint8(abs(CA24));
% figure,imshow(CA24);
% title('RESULTANT IMAGE-24');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA25,CH25,CV25,CD25] = lwt2(double(Block(:,:,25)),LS);
imlwt25=[abs(CA25) abs(CH25); abs(CV25) abs(CD25)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt25, []);
% title('2-D WAVELET');
CA25= uint8(abs(CA25));
% figure,imshow(CA25);
% title('RESULTANT IMAGE-25');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA26,CH26,CV26,CD26] = lwt2(double(Block(:,:,26)),LS);
imlwt26=[abs(CA26) abs(CH26); abs(CV26) abs(CD26)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt26, []);
% title('2-D WAVELET');
CA26= uint8(abs(CA26));
% figure,imshow(CA26);
% title('RESULTANT IMAGE-26');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA27,CH27,CV27,CD27] = lwt2(double(Block(:,:,27)),LS);
imlwt27=[abs(CA27) abs(CH27); abs(CV27) abs(CD27)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt27, []);
% title('2-D WAVELET');
CA27= uint8(abs(CA27));
% figure,imshow(CA27);
% title('RESULTANT IMAGE-27');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA28,CH28,CV28,CD28] = lwt2(double(Block(:,:,28)),LS);
imlwt28=[abs(CA28) abs(CH28); abs(CV28) abs(CD28)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt28, []);
% title('2-D WAVELET');
CA28= uint8(abs(CA28));
% figure,imshow(CA28);
% title('RESULTANT IMAGE-28');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA29,CH29,CV29,CD29] = lwt2(double(Block(:,:,29)),LS);
imlwt29=[abs(CA29) abs(CH29); abs(CV29) abs(CD29)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt29, []);
% title('2-D WAVELET');
CA29= uint8(abs(CA29));
% figure,imshow(CA29);
% title('RESULTANT IMAGE-29');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA30,CH30,CV30,CD30] = lwt2(double(Block(:,:,30)),LS);
imlwt30=[abs(CA30) abs(CH30); abs(CV30) abs(CD30)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt30, []);
% title('2-D WAVELET');
CA30= uint8(abs(CA30));
% figure,imshow(CA30);
% title('RESULTANT IMAGE-30');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA31,CH31,CV31,CD31] = lwt2(double(Block(:,:,31)),LS);
imlwt31=[abs(CA31) abs(CH31); abs(CV31) abs(CD31)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt31, []);
% title('2-D WAVELET');
CA31= uint8(abs(CA31));
% figure,imshow(CA31);
% title('RESULTANT IMAGE-31');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA32,CH32,CV32,CD32] = lwt2(double(Block(:,:,32)),LS);
imlwt32=[abs(CA32) abs(CH32); abs(CV32) abs(CD32)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt32, []);
% title('2-D WAVELET');
CA32= uint8(abs(CA32));
% figure,imshow(CA32);
% title('RESULTANT IMAGE-32');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA33,CH33,CV33,CD33] = lwt2(double(Block(:,:,33)),LS);
imlwt33=[abs(CA33) abs(CH33); abs(CV33) abs(CD33)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt33, []);
% title('2-D WAVELET');
CA33= uint8(abs(CA33));
% figure,imshow(CA33);
% title('RESULTANT IMAGE-33');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA34,CH34,CV34,CD34] = lwt2(double(Block(:,:,34)),LS);
imlwt34=[abs(CA34) abs(CH34); abs(CV34) abs(CD34)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt34, []);
% title('2-D WAVELET');
CA34= uint8(abs(CA34));
% figure,imshow(CA34);
% title('RESULTANT IMAGE-34');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA35,CH35,CV35,CD35] = lwt2(double(Block(:,:,35)),LS);
imlwt35=[abs(CA35) abs(CH35); abs(CV35) abs(CD35)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt35, []);
% title('2-D WAVELET');
CA35= uint8(abs(CA35));
% figure,imshow(CA35);
% title('RESULTANT IMAGE-35');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA36,CH36,CV36,CD36] = lwt2(double(Block(:,:,36)),LS);
imlwt36=[abs(CA36) abs(CH36); abs(CV36) abs(CD36)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt36, []);
% title('2-D WAVELET');
CA36= uint8(abs(CA36));
% figure,imshow(CA36);
% title('RESULTANT IMAGE-36');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA37,CH37,CV37,CD37] = lwt2(double(Block(:,:,37)),LS);
imlwt37=[abs(CA37) abs(CH37); abs(CV37) abs(CD37)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt37, []);
% title('2-D WAVELET');
CA37= uint8(abs(CA37));
% figure,imshow(CA37);
% title('RESULTANT IMAGE-37');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA38,CH38,CV38,CD38] = lwt2(double(Block(:,:,38)),LS);
imlwt38=[abs(CA38) abs(CH38); abs(CV38) abs(CD38)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt38, []);
% title('2-D WAVELET');
CA38= uint8(abs(CA38));
% figure,imshow(CA38);
% title('RESULTANT IMAGE-38');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA39,CH39,CV39,CD39] = lwt2(double(Block(:,:,39)),LS);
imlwt39=[abs(CA39) abs(CH39); abs(CV39) abs(CD39)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt39, []);
% title('2-D WAVELET');
CA39= uint8(abs(CA39));
% figure,imshow(CA39);
% title('RESULTANT IMAGE-39');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA40,CH40,CV40,CD40] = lwt2(double(Block(:,:,40)),LS);
imlwt40=[abs(CA40) abs(CH40); abs(CV40) abs(CD40)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt40, []);
% title('2-D WAVELET');
CA40= uint8(abs(CA40));
% figure,imshow(CA40);
% title('RESULTANT IMAGE-40');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA41,CH41,CV41,CD41] = lwt2(double(Block(:,:,41)),LS);
imlwt41=[abs(CA41) abs(CH41); abs(CV41) abs(CD41)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt41, []);
% title('2-D WAVELET');
CA41= uint8(abs(CA41));
% figure,imshow(CA41);
% title('RESULTANT IMAGE-41');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA42,CH42,CV42,CD42] = lwt2(double(Block(:,:,42)),LS);
imlwt42=[abs(CA42) abs(CH42); abs(CV42) abs(CD42)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt42, []);
% title('2-D WAVELET');
CA42= uint8(abs(CA42));
% figure,imshow(CA42);
% title('RESULTANT IMAGE-42');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA43,CH43,CV43,CD43] = lwt2(double(Block(:,:,43)),LS);
imlwt43=[abs(CA43) abs(CH43); abs(CV43) abs(CD43)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt43, []);
% title('2-D WAVELET');
CA43= uint8(abs(CA43));
% figure,imshow(CA43);
% title('RESULTANT IMAGE-43');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA44,CH44,CV44,CD44] = lwt2(double(Block(:,:,44)),LS);
imlwt44=[abs(CA44) abs(CH44); abs(CV44) abs(CD44)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt44, []);
% title('2-D WAVELET');
CA44= uint8(abs(CA44));
% figure,imshow(CA44);
% title('RESULTANT IMAGE-44');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA45,CH45,CV45,CD45] = lwt2(double(Block(:,:,45)),LS);
imlwt45=[abs(CA45) abs(CH45); abs(CV45) abs(CD45)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt45, []);
% title('2-D WAVELET');
CA45= uint8(abs(CA45));
% figure,imshow(CA45);
% title('RESULTANT IMAGE-45');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA46,CH46,CV46,CD46] = lwt2(double(Block(:,:,46)),LS);
imlwt46=[abs(CA46) abs(CH46); abs(CV46) abs(CD46)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt46, []);
% title('2-D WAVELET');
CA46= uint8(abs(CA46));
% figure,imshow(CA46);
% title('RESULTANT IMAGE-46');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA47,CH47,CV47,CD47] = lwt2(double(Block(:,:,47)),LS);
imlwt47=[abs(CA47) abs(CH47); abs(CV47) abs(CD47)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt47, []);
% title('2-D WAVELET');
CA47= uint8(abs(CA47));
% figure,imshow(CA47);
% title('RESULTANT IMAGE-47');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA48,CH48,CV48,CD48] = lwt2(double(Block(:,:,48)),LS);
imlwt48=[abs(CA48) abs(CH48); abs(CV48) abs(CD48)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt48, []);
% title('2-D WAVELET');
CA48= uint8(abs(CA48));
% figure,imshow(CA48);
% title('RESULTANT IMAGE-48');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA49,CH49,CV49,CD49] = lwt2(double(Block(:,:,49)),LS);
imlwt49=[abs(CA49) abs(CH49); abs(CV49) abs(CD49)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt49, []);
% title('2-D WAVELET');
CA49= uint8(abs(CA49));
% figure,imshow(CA49);
% title('RESULTANT IMAGE-49');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA50,CH50,CV50,CD50] = lwt2(double(Block(:,:,50)),LS);
imlwt50=[abs(CA50) abs(CH50); abs(CV50) abs(CD50)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt50, []);
% title('2-D WAVELET');
CA50= uint8(abs(CA50));
% figure,imshow(CA50);
% title('RESULTANT IMAGE-50');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA51,CH51,CV51,CD51] = lwt2(double(Block(:,:,51)),LS);
imlwt51=[abs(CA51) abs(CH51); abs(CV51) abs(CD51)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt51, []);
% title('2-D WAVELET');
CA51= uint8(abs(CA51));
% figure,imshow(CA51);
% title('RESULTANT IMAGE-51');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA52,CH52,CV52,CD52] = lwt2(double(Block(:,:,52)),LS);
imlwt52=[abs(CA52) abs(CH52); abs(CV52) abs(CD52)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt52, []);
% title('2-D WAVELET');
CA52= uint8(abs(CA52));
% figure,imshow(CA52);
% title('RESULTANT IMAGE-52');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA53,CH53,CV53,CD53] = lwt2(double(Block(:,:,53)),LS);
imlwt53=[abs(CA53) abs(CH53); abs(CV53) abs(CD53)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt53, []);
% title('2-D WAVELET');
CA53= uint8(abs(CA53));
% figure,imshow(CA53);
% title('RESULTANT IMAGE-53');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA54,CH54,CV54,CD54] = lwt2(double(Block(:,:,54)),LS);
imlwt54=[abs(CA54) abs(CH54); abs(CV54) abs(CD54)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt54, []);
% title('2-D WAVELET');
CA54= uint8(abs(CA54));
% figure,imshow(CA54);
% title('RESULTANT IMAGE-54');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA55,CH55,CV55,CD55] = lwt2(double(Block(:,:,55)),LS);
imlwt55=[abs(CA55) abs(CH55); abs(CV55) abs(CD55)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt55, []);
% title('2-D WAVELET');
CA55= uint8(abs(CA55));
% figure,imshow(CA55);
% title('RESULTANT IMAGE-55');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA56,CH56,CV56,CD56] = lwt2(double(Block(:,:,56)),LS);
imlwt56=[abs(CA56) abs(CH56); abs(CV56) abs(CD56)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt56, []);
% title('2-D WAVELET');
CA56= uint8(abs(CA56));
% figure,imshow(CA56);
% title('RESULTANT IMAGE-56');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA57,CH57,CV57,CD57] = lwt2(double(Block(:,:,57)),LS);
imlwt57=[abs(CA57) abs(CH57); abs(CV57) abs(CD57)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt57, []);
% title('2-D WAVELET');
CA57= uint8(abs(CA57));
% figure,imshow(CA57);
% title('RESULTANT IMAGE-57');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA58,CH58,CV58,CD58] = lwt2(double(Block(:,:,58)),LS);
imlwt58=[abs(CA58) abs(CH58); abs(CV58) abs(CD58)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt58, []);
% title('2-D WAVELET');
CA58= uint8(abs(CA58));
% figure,imshow(CA58);
% title('RESULTANT IMAGE-58');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA59,CH59,CV59,CD59] = lwt2(double(Block(:,:,59)),LS);
imlwt59=[abs(CA59) abs(CH59); abs(CV59) abs(CD59)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt59, []);
% title('2-D WAVELET');
CA59= uint8(abs(CA59));
% figure,imshow(CA59);
% title('RESULTANT IMAGE-59');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA60,CH60,CV60,CD60] = lwt2(double(Block(:,:,60)),LS);
imlwt60=[abs(CA60) abs(CH60); abs(CV60) abs(CD60)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt60, []);
% title('2-D WAVELET');
CA60= uint8(abs(CA60));
% figure,imshow(CA60);
% title('RESULTANT IMAGE-60');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA61,CH61,CV61,CD61] = lwt2(double(Block(:,:,61)),LS);
imlwt61=[abs(CA61) abs(CH61); abs(CV61) abs(CD61)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt61, []);
% title('2-D WAVELET');
CA61= uint8(abs(CA61));
% figure,imshow(CA61);
% title('RESULTANT IMAGE-61');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA62,CH62,CV62,CD62] = lwt2(double(Block(:,:,62)),LS);
imlwt62=[abs(CA62) abs(CH62); abs(CV62) abs(CD62)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt62, []);
% title('2-D WAVELET');
CA62= uint8(abs(CA62));
% figure,imshow(CA62);
% title('RESULTANT IMAGE-62');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA63,CH63,CV63,CD63] = lwt2(double(Block(:,:,63)),LS);
imlwt63=[abs(CA63) abs(CH63); abs(CV63) abs(CD63)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt63, []);
% title('2-D WAVELET');
CA63= uint8(abs(CA63));
% figure,imshow(CA63);
% title('RESULTANT IMAGE-63');
%----------------------------------------------------------------------------
LS = liftwave('cdf2.2','Int2Int');
[CA64,CH64,CV64,CD64] = lwt2(double(Block(:,:,64)),LS);
imlwt64=[abs(CA64) abs(CH64); abs(CV64) abs(CD64)]; % display the absolute value of 1D-IWT coeff
% figure,imshow(imlwt64, []);
% title('2-D WAVELET');
CA64= uint8(abs(CA64));
% figure,imshow(CA64);
% title('RESULTANT IMAGE-64');
%--------------------------------------------------------------------------
T=5;
[m n]=size(CA1);
for i=1:m
    for j=1:n
        if ((abs(CA1(i,j)))&lt; T)
            CA1(i,j)=CA1(i,j);
        elseif ((abs(CA1(i,j)))>= T)
            CA1(i,j)=(sign(CA1(i,j))*(floor((abs(CA1(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CAc1);
% title('COMPRESSED IMAGE-1');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA2);
for i=1:m
    for j=1:n
        if ((abs(CA2(i,j)))&lt; T)
            CA2(i,j)=CA2(i,j);
        elseif ((abs(CA2(i,j)))>= T)
            CA2(i,j)=(sign(CA2(i,j))*(floor((abs(CA2(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA2);
% title('COMPRESSED IMAGE-2');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=4;
[m n]=size(CA3);
for i=1:m
    for j=1:n
        if ((abs(CA3(i,j)))&lt; T)
            CA3(i,j)=CA3(i,j);
        elseif ((abs(CA3(i,j)))>= T)
            CA3(i,j)=(sign(CA3(i,j))*(floor((abs(CA3(i,j))-T)/2))+T);
        end
    end
end
figure,imshow(CA3);
title('COMPRESSED IMAGE-3');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=4;
[m n]=size(CA4);
for i=1:m
    for j=1:n
        if ((abs(CA4(i,j)))&lt; T)
            CA4(i,j)=CA4(i,j);
        elseif ((abs(CA4(i,j)))>= T)
            CA4(i,j)=(sign(CA4(i,j))*(floor((abs(CA4(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA4);
% title('COMPRESSED IMAGE-4');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA5);
for i=1:m
    for j=1:n
        if ((abs(CA5(i,j)))&lt; T)
            CA5(i,j)=CA5(i,j);
        elseif ((abs(CA5(i,j)))>= T)
            CA5(i,j)=(sign(CA5(i,j))*(floor((abs(CA5(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA5);
% title('COMPRESSED IMAGE-5');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA6);
for i=1:m
    for j=1:n
        if ((abs(CA6(i,j)))&lt; T)
            CA6(i,j)=CA6(i,j);
        elseif ((abs(CA6(i,j)))>= T)
            CA6(i,j)=(sign(CA6(i,j))*(floor((abs(CA6(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA6);
% title('COMPRESSED IMAGE-6');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA7);
for i=1:m
    for j=1:n
        if ((abs(CA7(i,j)))&lt; T)
            CA7(i,j)=CA7(i,j);
        elseif ((abs(CA7(i,j)))>= T)
            CA7(i,j)=(sign(CA7(i,j))*(floor((abs(CA7(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA7);
% title('COMPRESSED IMAGE-7');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA8);
for i=1:m
    for j=1:n
        if ((abs(CA8(i,j)))&lt; T)
            CA8(i,j)=CA8(i,j);
        elseif ((abs(CA8(i,j)))>= T)
            CA8(i,j)=(sign(CA8(i,j))*(floor((abs(CA8(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA8);
% title('COMPRESSED IMAGE-8');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA9);
for i=1:m
    for j=1:n
        if ((abs(CA9(i,j)))&lt; T)
            CA9(i,j)=CA9(i,j);
        elseif ((abs(CA9(i,j)))>= T)
            CA9(i,j)=(sign(CA9(i,j))*(floor((abs(CA9(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA9);
% title('COMPRESSED IMAGE-9');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA10);
for i=1:m
    for j=1:n
        if ((abs(CA10(i,j)))&lt; T)
            CA10(i,j)=CA10(i,j);
        elseif ((abs(CA10(i,j)))>= T)
            CA10(i,j)=(sign(CA10(i,j))*(floor((abs(CA10(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA10);
% title('COMPRESSED IMAGE-10');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA11);
for i=1:m
    for j=1:n
        if ((abs(CA11(i,j)))&lt; T)
            CA11(i,j)=CA11(i,j);
        elseif ((abs(CA11(i,j)))>= T)
            CA11(i,j)=(sign(CA11(i,j))*(floor((abs(CA11(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA11);
% title('COMPRESSED IMAGE-11');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA12);
for i=1:m
    for j=1:n
        if ((abs(CA12(i,j)))&lt; T)
            CA12(i,j)=CA12(i,j);
        elseif ((abs(CA12(i,j)))>= T)
            CA12(i,j)=(sign(CA12(i,j))*(floor((abs(CA12(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA12);
% title('COMPRESSED IMAGE-12');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA13);
for i=1:m
    for j=1:n
        if ((abs(CA13(i,j)))&lt; T)
            CA13(i,j)=CA13(i,j);
        elseif ((abs(CA13(i,j)))>= T)
            CA13(i,j)=(sign(CA13(i,j))*(floor((abs(CA13(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA13);
% title('COMPRESSED IMAGE-13');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA14);
for i=1:m
    for j=1:n
        if ((abs(CA14(i,j)))&lt; T)
            CA14(i,j)=CA14(i,j);
        elseif ((abs(CA14(i,j)))>= T)
            CA14(i,j)=(sign(CA14(i,j))*(floor((abs(CA14(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA14);
% title('COMPRESSED IMAGE-14');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA15);
for i=1:m
    for j=1:n
        if ((abs(CA1(i,j)))&lt; T)
            CA15(i,j)=CA15(i,j);
        elseif ((abs(CA15(i,j)))>= T)
            CA15(i,j)=(sign(CA15(i,j))*(floor((abs(CA15(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA15);
% title('COMPRESSED IMAGE-15');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA16);
for i=1:m
    for j=1:n
        if ((abs(CA16(i,j)))&lt; T)
            CA16(i,j)=CA16(i,j);
        elseif ((abs(CA16(i,j)))>= T)
            CA16(i,j)=(sign(CA16(i,j))*(floor((abs(CA16(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA16);
% title('COMPRESSED IMAGE-16');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA17);
for i=1:m
    for j=1:n
        if ((abs(CA17(i,j)))&lt; T)
            CA17(i,j)=CA17(i,j);
        elseif ((abs(CA17(i,j)))>= T)
            CA17(i,j)=(sign(CA17(i,j))*(floor((abs(CA17(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA17);
% title('COMPRESSED IMAGE-17');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA18);
for i=1:m
    for j=1:n
        if ((abs(CA18(i,j)))&lt; T)
            CA18(i,j)=CA18(i,j);
        elseif ((abs(CA18(i,j)))>= T)
            CA18(i,j)=(sign(CA18(i,j))*(floor((abs(CA18(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA18);
% title('COMPRESSED IMAGE-18');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA19);
for i=1:m
    for j=1:n
        if ((abs(CA19(i,j)))&lt; T)
            CA19(i,j)=CA19(i,j);
        elseif ((abs(CA19(i,j)))>= T)
            CA19(i,j)=(sign(CA19(i,j))*(floor((abs(CA19(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA19);
% title('COMPRESSED IMAGE-19');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=3;
[m n]=size(CA20);
for i=1:m
    for j=1:n
        if ((abs(CA20(i,j)))&lt; T)
            CA20(i,j)=CA20(i,j);
        elseif ((abs(CA20(i,j)))>= T)
            CA20(i,j)=(sign(CA20(i,j))*(floor((abs(CA20(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA20);
% title('COMPRESSED IMAGE-20');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=4;
[m n]=size(CA21);
for i=1:m
    for j=1:n
        if ((abs(CA21(i,j)))&lt; T)
            CA21(i,j)=CA21(i,j);
        elseif ((abs(CA21(i,j)))>= T)
            CA21(i,j)=(sign(CA21(i,j))*(floor((abs(CA21(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA21);
% title('COMPRESSED IMAGE-21');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=4;
[m n]=size(CA22);
for i=1:m
    for j=1:n
        if ((abs(CA22(i,j)))&lt; T)
            CA22(i,j)=CA22(i,j);
        elseif ((abs(CA22(i,j)))>= T)
            CA22(i,j)=(sign(CA22(i,j))*(floor((abs(CA22(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA22);
% title('COMPRESSED IMAGE-22');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA23);
for i=1:m
    for j=1:n
        if ((abs(CA23(i,j)))&lt; T)
            CA23(i,j)=CA23(i,j);
        elseif ((abs(CA23(i,j)))>= T)
            CA23(i,j)=(sign(CA23(i,j))*(floor((abs(CA23(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA23);
% title('COMPRESSED IMAGE-23');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA24);
for i=1:m
    for j=1:n
        if ((abs(CA24(i,j)))&lt; T)
            CA24(i,j)=CA24(i,j);
        elseif ((abs(CA24(i,j)))>= T)
            CA24(i,j)=(sign(CA24(i,j))*(floor((abs(CA24(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA24);
% title('COMPRESSED IMAGE-24');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA25);
for i=1:m
    for j=1:n
        if ((abs(CA25(i,j)))&lt; T)
            CA25(i,j)=CA25(i,j);
        elseif ((abs(CA25(i,j)))>= T)
            CA25(i,j)=(sign(CA25(i,j))*(floor((abs(CA25(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA25);
% title('COMPRESSED IMAGE-25');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=3;
[m n]=size(CA26);
for i=1:m
    for j=1:n
        if ((abs(CA26(i,j)))&lt; T)
            CA26(i,j)=CA26(i,j);
        elseif ((abs(CA26(i,j)))>= T)
            CA26(i,j)=(sign(CA26(i,j))*(floor((abs(CA26(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA26);
% title('COMPRESSED IMAGE-26');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=3;
[m n]=size(CA27);
for i=1:m
    for j=1:n
        if ((abs(CA27(i,j)))&lt; T)
            CA27(i,j)=CA27(i,j);
        elseif ((abs(CA27(i,j)))>= T)
            CA27(i,j)=(sign(CA27(i,j))*(floor((abs(CA27(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA27);
% title('COMPRESSED IMAGE-27');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=3;
[m n]=size(CA28);
for i=1:m
    for j=1:n
        if ((abs(CA28(i,j)))&lt; T)
            CA28(i,j)=CA28(i,j);
        elseif ((abs(CA28(i,j)))>= T)
            CA28(i,j)=(sign(CA28(i,j))*(floor((abs(CA28(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA28);
% title('COMPRESSED IMAGE-28');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA29);
for i=1:m
    for j=1:n
        if ((abs(CA29(i,j)))&lt; T)
            CA29(i,j)=CA29(i,j);
        elseif ((abs(CA29(i,j)))>= T)
            CA29(i,j)=(sign(CA29(i,j))*(floor((abs(CA29(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA29);
% title('COMPRESSED IMAGE-29');
% -----------------------------------------------------------------------------------------------------------------------
% -------------------------------------------------------------------------
T=5;
[m n]=size(CA30);
for i=1:m
    for j=1:n
        if ((abs(CA30(i,j)))&lt; T)
            CA30(i,j)=CA30(i,j);
        elseif ((abs(CA30(i,j)))>= T)
            CA30(i,j)=(sign(CA30(i,j))*(floor((abs(CA30(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA30);
% title('COMPRESSED IMAGE-30');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA31);
for i=1:m
    for j=1:n
        if ((abs(CA31(i,j)))&lt; T)
            CA31(i,j)=CA31(i,j);
        elseif ((abs(CA31(i,j)))>= T)
            CA31(i,j)=(sign(CA31(i,j))*(floor((abs(CA31(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA31);
% title('COMPRESSED IMAGE-31');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA32);
for i=1:m
    for j=1:n
        if ((abs(CA32(i,j)))&lt; T)
            CA32(i,j)=CA32(i,j);
        elseif ((abs(CA32(i,j)))>= T)
            CA32(i,j)=(sign(CA32(i,j))*(floor((abs(CA32(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA32);
% title('COMPRESSED IMAGE-32');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA33);
for i=1:m
    for j=1:n
        if ((abs(CA1(i,j)))&lt; T)
            CA33(i,j)=CA33(i,j);
        elseif CA33((abs(CA33(i,j)))>= T)
            CA33(i,j)=(sign(CA33(i,j))*(floor((abs(CA33(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA33);
% title('COMPRESSED IMAGE-33');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=4;
[m n]=size(CA34);
for i=1:m
    for j=1:n
        if ((abs(CA34(i,j)))&lt; T)
            CA34(i,j)=CA34(i,j);
        elseif ((abs(CA34(i,j)))>= T)
            CA34(i,j)=(sign(CA34(i,j))*(floor((abs(CA34(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA34);
% title('COMPRESSED IMAGE-34');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=3;
[m n]=size(CA35);
for i=1:m
    for j=1:n
        if ((abs(CA35(i,j)))&lt; T)
            CA35(i,j)=CA35(i,j);
        elseif ((abs(CA35(i,j)))>= T)
            CA35(i,j)=(sign(CA35(i,j))*(floor((abs(CA35(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA35);
% title('COMPRESSED IMAGE-35');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA36);
for i=1:m
    for j=1:n
        if ((abs(CA36(i,j)))&lt; T)
            CA36(i,j)=CA36(i,j);
        elseif ((abs(CA36(i,j)))>= T)
            CA36(i,j)=(sign(CA36(i,j))*(floor((abs(CA36(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA36);
% title('COMPRESSED IMAGE-36');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=4;
[m n]=size(CA37);
for i=1:m
    for j=1:n
        if ((abs(CA37(i,j)))&lt; T)
            CA37(i,j)=CA37(i,j);
        elseif ((abs(CA37(i,j)))>= T)
            CA37(i,j)=(sign(CA37(i,j))*(floor((abs(CA37(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA37);
% title('COMPRESSED IMAGE-37');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=4;
[m n]=size(CA38);
for i=1:m
    for j=1:n
        if ((abs(CA38(i,j)))&lt; T)
            CA38(i,j)=CA38(i,j);
        elseif ((abs(CA38(i,j)))>= T)
            CA38(i,j)=(sign(CA38(i,j))*(floor((abs(CA38(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA38);
% title('COMPRESSED IMAGE-38');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA39);
for i=1:m
    for j=1:n
        if ((abs(CA39(i,j)))&lt; T)
            CA39(i,j)=CA39(i,j);
        elseif ((abs(CA39(i,j)))>= T)
            CA39(i,j)=(sign(CA39(i,j))*(floor((abs(CA39(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA39);
% title('COMPRESSED IMAGE-39');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA40);
for i=1:m
    for j=1:n
        if ((abs(CA40(i,j)))&lt; T)
            CA40(i,j)=CA40(i,j);
        elseif ((abs(CA40(i,j)))>= T)
            CA40(i,j)=(sign(CA40(i,j))*(floor((abs(CA40(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA40);
% title('COMPRESSED IMAGE-40');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=4;
[m n]=size(CA41);
for i=1:m
    for j=1:n
        if ((abs(CA41(i,j)))&lt; T)
            CA41(i,j)=CA41(i,j);
        elseif ((abs(CA41(i,j)))>= T)
            CA41(i,j)=(sign(CA41(i,j))*(floor((abs(CA41(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA41);
% title('COMPRESSED IMAGE-41');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=4;
[m n]=size(CA42);
for i=1:m
    for j=1:n
        if ((abs(CA42(i,j)))&lt; T)
            CA42(i,j)=CA42(i,j);
        elseif ((abs(CA42(i,j)))>= T)
            CA42(i,j)=(sign(CA42(i,j))*(floor((abs(CA42(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA42);
% title('COMPRESSED IMAGE-42');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=4;
[m n]=size(CA43);
for i=1:m
    for j=1:n
        if ((abs(CA43(i,j)))&lt; T)
            CA43(i,j)=CA43(i,j);
        elseif ((abs(CA43(i,j)))>= T)
            CA43(i,j)=(sign(CA43(i,j))*(floor((abs(CA43(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA43);
% title('COMPRESSED IMAGE-43');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA44);
for i=1:m
    for j=1:n
        if ((abs(CA44(i,j)))&lt; T)
            CA44(i,j)=CA44(i,j);
        elseif ((abs(CA44(i,j)))>= T)
            CA44(i,j)=(sign(CA44(i,j))*(floor((abs(CA44(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA44);
% title('COMPRESSED IMAGE-44');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=4;
[m n]=size(CA45);
for i=1:m
    for j=1:n
        if ((abs(CA45(i,j)))&lt; T)
            CA45(i,j)=CA45(i,j);
        elseif ((abs(CA45(i,j)))>= T)
            CA45(i,j)=(sign(CA45(i,j))*(floor((abs(CA45(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA45);
% title('COMPRESSED IMAGE-45');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=4;
[m n]=size(CA46);
for i=1:m
    for j=1:n
        if ((abs(CA46(i,j)))&lt; T)
            CA46(i,j)=CA46(i,j);
        elseif ((abs(CA46(i,j)))>= T)
            CA46(i,j)=(sign(CA46(i,j))*(floor((abs(CA46(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA46);
% title('COMPRESSED IMAGE-46');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=3;
[m n]=size(CA47);
for i=1:m
    for j=1:n
        if ((abs(CA47(i,j)))&lt; T)
            CA47(i,j)=CA47(i,j);
        elseif ((abs(CA47(i,j)))>= T)
            CA47(i,j)=(sign(CA47(i,j))*(floor((abs(CA47(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA47);
% title('COMPRESSED IMAGE-47');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA48);
for i=1:m
    for j=1:n
        if ((abs(CA48(i,j)))&lt; T)
            CA48(i,j)=CA48(i,j);
        elseif ((abs(CA48(i,j)))>= T)
            CA48(i,j)=(sign(CA48(i,j))*(floor((abs(CA48(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA48);
% title('COMPRESSED IMAGE-48');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=4;
[m n]=size(CA49);
for i=1:m
    for j=1:n
        if ((abs(CA49(i,j)))&lt; T)
            CA49(i,j)=CA49(i,j);
        elseif ((abs(CA49(i,j)))>= T)
            CA49(i,j)=(sign(CA49(i,j))*(floor((abs(CA49(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA49);
% title('COMPRESSED IMAGE-49');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=4;
[m n]=size(CA50);
for i=1:m
    for j=1:n
        if ((abs(CA50(i,j)))&lt; T)
            CA50(i,j)=CA50(i,j);
        elseif ((abs(CA50(i,j)))>= T)
            CA50(i,j)=(sign(CA50(i,j))*(floor((abs(CA50(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA50);
% title('COMPRESSED IMAGE-50');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA51);
for i=1:m
    for j=1:n
        if ((abs(CA51(i,j)))&lt; T)
            CA51(i,j)=CA51(i,j);
        elseif ((abs(CA51(i,j)))>= T)
            CA51(i,j)=(sign(CA51(i,j))*(floor((abs(CA51(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA51);
% title('COMPRESSED IMAGE-51');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA52);
for i=1:m
    for j=1:n
        if ((abs(CA52(i,j)))&lt; T)
            CA52(i,j)=CA52(i,j);
        elseif ((abs(CA52(i,j)))>= T)
            CA52(i,j)=(sign(CA52(i,j))*(floor((abs(CA52(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA52);
% title('COMPRESSED IMAGE-52');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA53);
for i=1:m
    for j=1:n
        if ((abs(CA53(i,j)))&lt; T)
            CA53(i,j)=CA53(i,j);
        elseif ((abs(CA53(i,j)))>= T)
            CA53(i,j)=(sign(CA53(i,j))*(floor((abs(CA53(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA53);
% title('COMPRESSED IMAGE-53');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=4;
[m n]=size(CA54);
for i=1:m
    for j=1:n
        if ((abs(CA54(i,j)))&lt; T)
            CA54(i,j)=CA54(i,j);
        elseif ((abs(CA54(i,j)))>= T)
            CA54(i,j)=(sign(CA54(i,j))*(floor((abs(CA54(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA54);
% title('COMPRESSED IMAGE-54');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA55);
for i=1:m
    for j=1:n
        if ((abs(CA55(i,j)))&lt; T)
            CA55(i,j)=CA55(i,j);
        elseif ((abs(CA55(i,j)))>= T)
            CA55(i,j)=(sign(CA55(i,j))*(floor((abs(CA55(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA55);
% title('COMPRESSED IMAGE-55');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA56);
for i=1:m
    for j=1:n
        if ((abs(CA56(i,j)))&lt; T)
            CA56(i,j)=CA56(i,j);
        elseif ((abs(CA56(i,j)))>= T)
            CA56(i,j)=(sign(CA56(i,j))*(floor((abs(CA56(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA56);
% title('COMPRESSED IMAGE-56');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=5;
[m n]=size(CA57);
for i=1:m
    for j=1:n
        if ((abs(CA57(i,j)))&lt; T)
            CA57(i,j)=CA57(i,j);
        elseif ((abs(CA57(i,j)))>= T)
            CA57(i,j)=(sign(CA57(i,j))*(floor((abs(CA57(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA57);
% title('COMPRESSED IMAGE-57');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=3;
[m n]=size(CA58);
for i=1:m
    for j=1:n
        if ((abs(CA58(i,j)))&lt; T)
            CA58(i,j)=CA58(i,j);
        elseif ((abs(CA58(i,j)))>= T)
            CA58(i,j)=(sign(CA58(i,j))*(floor((abs(CA58(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA58);
% title('COMPRESSED IMAGE-58');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=3;
[m n]=size(CA59);
for i=1:m
    for j=1:n
        if ((abs(CA59(i,j)))&lt; T)
            CA59(i,j)=CA59(i,j);
        elseif ((abs(CA59(i,j)))>= T)
            CA59(i,j)=(sign(CA59(i,j))*(floor((abs(CA59(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA59);
% title('COMPRESSED IMAGE-59');
% -----------------------------------------------------------------------------------------------------------------------
% -----------------------------------------------------------------------------------------------------------------------
T=3;
[m n]=size(CA60);
for i=1:m
    for j=1:n
        if ((abs(CA60(i,j)))&lt; T)
            CA60(i,j)=CA60(i,j);
        elseif ((abs(CA60(i,j)))>= T)
            CA60(i,j)=(sign(CA60(i,j))*(floor((abs(CA60(i,j))-T)/2))+T);
        end
    end
end
% figure,imshow(CA60);
% title('COMPRESSED IMAGE-60');
% -----

[Walter Roberson]

duplicates http://www.mathworks.co.uk/matlabcentral/answers/65559-how-to-use-the-arithenco-command-for-each-block-in-embedding-process

[Image Analyst]

revathi, when you post two questions with huge blocks of code, apparently completely or mostly the same, and you edited both of them, how are we supposed to know which is the latest and corrected version?

[Walter Roberson]

This post had been closed as being a duplicate of another post. If it is to be re-opened as not being duplicate, you need to indicate clearly the differences between this post and the other one.