shamir叠像术 分类： 图像处理 2015-07-08 16:50 17人阅读 评论(1) 收藏

K=imread('SHNU.bmp');

height=info.Height;

width=info.Width;

A=zeros(height,width);

B=zeros(height,width);

C=zeros(height,width);

% 随机生成A,B,C矩阵

for i=1: height

for j=1: width

random=rand()*2;

if(random>1)

A(i,j)=1;B(i,j)=1;C(i,j)=1;

else

A(i,j)=0;B(i,j)=0;C(i,j)=0;

end

end

end

%%%%%生成子秘钥A,B,C图片

for i=1: height

for j=1: width

if(K(i,j)==0)

random=rand()*4;

if(random<1)

A(i,j)=0;A(i,j+1)=0;A(i+1,j)=1;A(i+1,j+1)=1;

B(i,j)=0;B(i+1,j)=0;B(i,j+1)=1;B(i+1,j+1)=1;

C(i,j)=0;C(i+1,j+1)=0;C(i,j+1)=1;C(i+1,j)=1;

else if(random<2)

A(i,j)=0;A(i,j+1)=0;A(i+1,j)=1;A(i+1,j+1)=1;

C(i,j+1)=0;C(i+1,j)=0;C(i,j)=1;C(i+1,j+1)=1;

B(i,j+1)=0;B(i+1,j+1)=0;B(i,j)=1;B(i+1,j)=1;

else if(random<3)

B(i,j)=0;B(i+1,j)=0;B(i,j+1)=1;B(i+1,j+1)=1;

A(i+1,j)=0;A(i+1,j+1)=0;A(i,j)=1;A(i,j+1)=1;

C(i,j+1)=0;C(i+1,j)=0;C(i,j)=1;C(i+1,j+1)=1;

else if(random<4)

B(i,j+1)=0;B(i+1,j+1)=0;B(i,j)=1;B(i+1,j)=1;

C(i,j)=0;C(i+1,j+1)=0;C(i,j+1)=1;C(i+1,j)=1;

A(i+1,j)=0;A(i+1,j+1)=0;A(i,j)=1;A(i,j+1)=1;

end

end

end

end

end

end

end

imwrite(A,'a.bmp','bmp');

imwrite(B,'b.bmp','bmp');

imwrite(C,'c.bmp','bmp');

figure

subplot(2,2,1),imshow(K);

subplot(2,2,2),imshow(A);

subplot(2,2,3),imshow(B);

subplot(2,2,4),imshow(C);

axes_handle = get(gcf, 'children');

axes(axes_handle(4)); title('original image');

axes(axes_handle(3)); title('image a');

axes(axes_handle(2)); title('image b');

axes(axes_handle(1)); title('image c');

%%%%隐藏信息的恢复

A1=imread('a.bmp');

B1=imread('b.bmp');

C1=imread('c.bmp');

random=rand()*3;

if(random<1)

K=A1;K12=A1;K13=A1;K23=A1;

else if(random<2)

K=B1;K12=B1;K13=B1;K23=B1;

else

K=C1;K12=C1;K13=C1;K23=C1;

end

end

for i=1: height-1

for j=1: width-1

if((xor(xor(A1(i,j),B1(i,j)),C1(i,j))==0)&&(xor(xor(A1(i,j+1),B1(i,j+1)),C1(i,j+1))==0)&&(xor(xor(A1(i+1,j),B1(i+1,j)),C1(i+1,j))==0)&&(xor(xor(A1(i+1,j+1),B1(i+1,j+1)),C1(i+1,j+1))==0))

K(i,j)=0;K(i,j+1)=0;K(i+1,j)=0;K(i+1,j+1)=0;

end

end

end

for i=1: height-1,

for j=1: width-1,

if((A1(i,j)||B1(i,j))==0&&(A1(i,j+1)||B1(i,j+1))==0&&(A1(i+1,j)||B1(i+1,j))==0&&(A1(i+1,j+1)||B1(i+1,j+1))==0)

K12(i,j)=0;K12(i,j+1)=0;K12(i+1,j)=0;K12(i+1,j+1)=0;

end

end

end

for i=1: height-1,

for j=1: width-1,

if((A1(i,j)||C1(i,j))==0&&(A1(i,j+1)||C1(i,j+1))==0&&(A1(i+1,j)||C1(i+1,j))==0&&(A1(i+1,j+1)||C1(i+1,j+1))==0)

K13(i,j)=0;K13(i,j+1)=0;K13(i+1,j)=0;K13(i+1,j+1)=0;

end

end

end

for i=1: height-1,

for j=1: width-1,

if((C1(i,j)||B1(i,j))==0&&(C1(i,j+1)||B1(i,j+1))==0&&(C1(i+1,j)||B1(i+1,j))==0&&(C1(i+1,j+1)||B1(i+1,j+1))==0)

K23(i,j)=0;K23(i,j+1)=0;K23(i+1,j)=0;K23(i+1,j+1)=0;

end

end

end

imwrite(K,'K.bmp','bmp')

imwrite(K12,'image_a+b.bmp','bmp')

imwrite(K13,'image_a+c.bmp','bmp')

imwrite(K23,'image_b+c.bmp','bmp')

figure

imshow(K);

title('recovery image');

figure

subplot(2,2,1),imshow(K12);

subplot(2,2,2),imshow(K13);

subplot(2,2,3),imshow(K23);

axes_handle = get(gcf, 'children');

axes(axes_handle(3)); title('image a+b');

axes(axes_handle(2)); title('image a+c');

axes(axes_handle(1)); title('image b+c');