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

2 次查看（过去 30 天）
Revathi sri
Revathi sri 2013-3-2
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');
% -------

请先登录，再回答此问题。

回答（1 个）

Walter Roberson
Walter Roberson 2013-3-2

类别

在 Help Center 和 File Exchange 中查找有关 Computer Vision with Simulink 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by