clc clear all orig_image=imread('D:\Aruna Research Papers\AMBTC_Interpolation papers\IP-1\Images 512\housee512.jpg'); sm_image=rgb2gray(orig_image); %orig_image1=[85,86,105,86; 84,95,75,82; 85,90,80,86; 86,86,84,81] % sm_image=[149,53,48,90; 113,103,44,66; 77,114,47,44; 73,132,98,48] orig_image1=sm_image; for i=2:1:511 for j=2:1:511 sm= sm_image(i,j)*0.6+ sm_image(i-1,j)*0.1+sm_image(i+1,j)*0.1+sm_image(i,j-1)*0.1+sm_image(i,j+1)*0.1; orig_image1(i,j)=sm; end end ip_image1=double(orig_image1); ip_image2=ip_image1; ip_image3=ip_image1; % AMBTC Image ip_image4=ip_image1; % Stego Image thr1=10; thr2=25; % % thr1=20; % thr2=35; % % thr1=30; % thr2=45; % % thr1=40; % thr2=55; % % thr1=50; % thr2=65; % Capacity=0; for i=1:4:512 for j=1:4:512 Average=round((ip_image1(i,j)+ip_image1(i,j+1)+ip_image1(i,j+2)+ip_image1(i,j+3)+ip_image1(i+1,j)+ip_image1(i+1,j+1)+ip_image1(i+1,j+2)+ip_image1(i+1,j+3)+ip_image1(i+2,j)+ip_image1(i+2,j+1)+ip_image1(i+2,j+2)+ip_image1(i+2,j+3)+ip_image1(i+3,j)+ip_image1(i+3,j+1)+ip_image1(i+3,j+2)+ip_image1(i+3,j+3))/16); x=0; y=0; count=0; if (ip_image1(i,j)thr1 && (High-Low)<=thr2) %(High-Low)>thr1 Capacity=Capacity+8; xx=0; yy=0; count=0; h1=8; h2=8; %%%%%%%%%%% ip_image2(i:i+3,j:j+3) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Less complex block%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%h1 a1=0; x1=255; b=round(a1+(x1-a1)*rand(1)); ccc=dec2bin(b,8); if (dec2bin(ip_image2(i,j),1)==ccc(1)) h1=h1-1; end if (dec2bin(ip_image2(i,j+2),1)==ccc(2)) h1=h1-1; end if (dec2bin(ip_image2(i+1,j),1)==ccc(3)) h1=h1-1; end if (dec2bin(ip_image2(i+1,j+2),1)==ccc(4)) h1=h1-1; end if (dec2bin(ip_image2(i+2,j),1)==ccc(5)) h1=h1-1; end if (dec2bin(ip_image2(i+2,j+2),1)==ccc(6)) h1=h1-1; end if (dec2bin(ip_image2(i+3,j),1)==ccc(7)) h1=h1-1; end if (dec2bin(ip_image2(i+3,j+2),1)==ccc(8)) h1=h1-1; end %%%%%%%%%%%%%%%%%%%%%%%%%%%%H2 if (dec2bin(ip_image2(i,j+1),1)==ccc(1)) h2=h2-1; end if (dec2bin(ip_image2(i,j+3),1)==ccc(2)) h2=h2-1; end if (dec2bin(ip_image2(i+1,j+2),1)==ccc(3)) h2=h2-1; end if (dec2bin(ip_image2(i+1,j+3),1)==ccc(4)) h2=h2-1; end if (dec2bin(ip_image2(i+2,j+1),1)==ccc(5)) h2=h2-1; end if (dec2bin(ip_image2(i+2,j+3),1)==ccc(6)) h2=h2-1; end if (dec2bin(ip_image2(i+3,j+1),1)==ccc(7)) h2=h2-1; end if (dec2bin(ip_image2(i+3,j+3),1)==ccc(8)) h2=h2-1; end if h1

=0 && diff1<=7) Capacity=Capacity+3; a1=0; x1=7; b1=round(a1+(x1-a1)*rand(1)); diff2=b1; diff=abs(diff1-diff2); if mod(diff,2)==1 Low=Low + ceil(diff/2); High=High- floor(diff/2); else Low=Low+ floor(diff/2); High=High- ceil(diff/2); end elseif (diff1>=8 && diff1<=15) Capacity=Capacity+3; a1=0; x1=7; b1=round(a1+(x1-a1)*rand(1)); diff2=7+b1; diff=abs(diff1-diff2); if mod(diff,2)==1 Low=Low + ceil(diff/2); High=High- floor(diff/2); else Low=Low+ floor(diff/2); High=High- ceil(diff/2); end elseif (diff1>=16 && diff1<=31) Capacity=Capacity+4; a1=0; x1=15; b1=round(a1+(x1-a1)*rand(1)); diff2=16+b1; diff=abs(diff1-diff2); if mod(diff,2)==1 Low=Low + ceil(diff/2); High=High- floor(diff/2); else Low=Low+ floor(diff/2); High=High- ceil(diff/2); end elseif (diff1>=32 && diff1<=63) Capacity=Capacity+5; a1=0; x1=31; b1=round(a1+(x1-a1)*rand(1)); diff2=32+b1; diff=abs(diff1-diff2); if mod(diff,2)==1 Low=Low + ceil(diff/2); High=High- floor(diff/2); else Low=Low+ floor(diff/2); High=High- ceil(diff/2); end elseif (diff1>=64 && diff1<=127) Capacity=Capacity+6; a1=0; x1=63; b1=round(a1+(x1-a1)*rand(1)); diff2=64+b1; diff=abs(diff1-diff2); if mod(diff,2)==1 Low=Low + ceil(diff/2); High=High- floor(diff/2); else Low=Low+ floor(diff/2); High=High- ceil(diff/2); end else a1=0; x1=127; Capacity=Capacity+7; b1=round(a1+(x1-a1)*rand(1)); diff2=128+b1; diff=abs(diff1-diff2); if mod(diff,2)==1 Low=Low + ceil(diff/2); High=High- floor(diff/2); else Low=Low+ floor(diff/2); High=High- ceil(diff/2); end end %%%%%%%%%%%%%%%%%%%%%%%% if (ip_image2(i,j)==0) ip_image4(i,j)=Low; else ip_image4(i,j)=High; end if (ip_image2(i,j+1)==0) ip_image4(i,j+1)=Low; else ip_image4(i,j+1)=High; end if (ip_image2(i,j+2)==0) ip_image4(i,j+2)=Low; else ip_image4(i,j+2)=High; end if (ip_image2(i,j+3)==0) ip_image4(i,j+3)=Low; else ip_image4(i,j+3)=High; end %%%%%%%%%%%%%%%%%%%%% if (ip_image2(i+1,j)==0) ip_image4(i+1,j)=Low; else ip_image4(i+1,j)=High; end if (ip_image2(i+1,j+1)==0) ip_image4(i+1,j+1)=Low; else ip_image4(i+1,j+1)=High; end if (ip_image2(i+1,j+2)==0) ip_image4(i+1,j+2)=Low; else ip_image4(i+1,j+2)=High; end if (ip_image2(i+1,j+3)==0) ip_image4(i+1,j+3)=Low; else ip_image4(i+1,j+3)=High; end %%%%%%%%%%%%%%%%%%%%% if (ip_image2(i+2,j)==0) ip_image4(i+2,j)=Low; else ip_image4(i+2,j)=High; end if (ip_image2(i+2,j+1)==0) ip_image4(i+2,j+1)=Low; else ip_image4(i+2,j+1)=High; end if (ip_image2(i+2,j+2)==0) ip_image4(i+2,j+2)=Low; else ip_image4(i+2,j+2)=High; end if (ip_image2(i+2,j+3)==0) ip_image4(i+2,j+3)=Low; else ip_image4(i+2,j+3)=High; end %%%%%%%%%%%%%%%%%%%%% if (ip_image2(i+3,j)==0) ip_image4(i+3,j)=Low; else ip_image4(i+3,j)=High; end if (ip_image2(i+3,j+1)==0) ip_image4(i+3,j+1)=Low; else ip_image4(i+3,j+1)=High; end if (ip_image2(i+3,j+2)==0) ip_image4(i+3,j+2)=Low; else ip_image4(i+3,j+2)=High; end if (ip_image2(i+3,j+3)==0) ip_image4(i+3,j+3)=Low; else ip_image4(i+3,j+3)=High; end end end end figure,imshow(ip_image1); figure,imshow(ip_image3); figure,imshow(ip_image4); origImg = double(orig_image1); distImg = double(ip_image3); stegoImg = double(ip_image4); [M N] = size(origImg); error = origImg - distImg; MSE = sum(sum(error .* error)) / (M * N); AMBTC_PSNR = 10*log(255*255/MSE) / log(10) error = origImg - stegoImg; MSE = sum(sum(error .* error)) / (M * N); OrgStg_PSNR = 10*log(255*255/MSE) / log(10) error = distImg - stegoImg; MSE = sum(sum(error .* error)) / (M * N); AmbtcStg_PSNR = 10*log(255*255/MSE) / log(10) Capacity