图像处理之像素级环形扩散

刚入门学习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("还原图");

运行截图：图像处理之像素级环形扩散
后续还会继续学习，尽量将其完全实现。。。。。。