剛入門學習matlab,自己根據《位級同步置亂擴散和像素級環形擴散圖像加密算法》這個文章實作了一下像素級環形擴散這一小節的代碼,和大家分享一下,能力有限,有不足之處還望見諒。代碼如下:
clear all
clc
X=imread('lena.jpg');
X=double(X);
figure()
subplot(221)
imshow(X,[]);
title("原圖");
S=size(X,2);
%密鑰生成
% % % % % % % %
x2=0.4;
L2=0.3;
x2 = mod(x2+sum(sum(X))/256,1);
for i=1:300
x2=F(x2,L2);
end
y3=zeros(1,S*S);
y4=zeros(1,S*S);
y3(1)=x2;
for i=1:S*S-1
y3(i+1)=F(y3(i),L2);
end
y4(1)=y3(end);
for i=1:S*S-1
y4(i+1)=F(y4(i),L2);
end
%混沌序列标準化
z3=floor(mod(y3.*10^14,256));
z4=floor(mod(y4.*10^14,256));
% % % % % % % %
H=reshape(flipud(rot90(X)),1,S^2);
H1=[];
H1(1)=bitxor(bitxor(sum(H),H(1)),z3(1));
for i=1:S^2-1
H1(i+1)=bitxor(bitxor(H(i+1),H1(i)),z3(i+1));
end
H1=reshape(H1,S,S);
H1=flipud(fliplr(H1));
V=reshape(H1,1,S^2);
E=[];
E(1)=bitxor(bitxor(V(1),sum(V)),z4(1));
for i=1:S^2-1
E(i+1)=bitxor(bitxor(E(i),V(i+1)),z4(i+1));
end
e=reshape(E,S,S);
subplot(222)
imshow(e,[]);
title("像素級環形擴散後");
% % % % % % % % %
% 還原
% % % % % % % % %
V1=[];
V1(1)=bitxor(bitxor(E(1),z4(1)),sum(V));
for i=1:S^2-1
V1(i+1)=bitxor(bitxor(E(i+1),z4(i+1)),E(i));
end
h22=reshape(V1,S,S);
h11=reshape(fliplr(flipud(h22)),1,S^2);
H2=[];
H2(1)=bitxor(bitxor(z3(1),h11(1)),sum(H));
for i=1:S^2-1
H2(i+1)=bitxor(bitxor(z3(i+1),h11(i+1)),h11(i));
end
H2=reshape(H2,S,S);
X1=rot90(flipud(H2),3);
subplot(223)
imshow(X1,[]);
title("還原圖");
運作截圖:
後續還會繼續學習,盡量将其完全實作。。。。。。