how to use the arithenco command for each block in embedding process

Question

Revathi sri on 2 Mar 2013

0
Link

Direct link to this question

https://www.mathworks.com/matlabcentral/answers/65559-how-to-use-the-arithenco-command-for-each-block-in-embedding-process

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');
% -------