how to use the arithenco in this code for each block in image?
Info
Diese Frage ist geschlossen. Öffnen Sie sie erneut, um sie zu bearbeiten oder zu beantworten.
Ältere Kommentare anzeigen
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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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)))< 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');
% -----
4 Kommentare
Image Analyst
am 4 Mär. 2013
What is "arithenco"?
Walter Roberson
am 5 Mär. 2013
Bearbeitet: Walter Roberson
am 5 Mär. 2013
a routine that implements arithmetic encoding
Walter Roberson
am 6 Mär. 2013
This was closed as being a duplicate. What modifications did you make to it such that it is no longer a duplicate ?
Antworten (0)
Diese Frage ist geschlossen.
Produkte
am 4 Mär. 2013
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
