% Arnold置乱变换
function imgn=Arnold(img,a,b,t)
%{
img 载体图像
a,b 置乱矩阵参数
t 置乱次数
%}
%读取图像尺寸,Arnold变换针对方阵,所以一般m=n
% m和n为图像的行和列
[m,n,p]=size(img);
%方阵边长
N=m;
%置乱
imgn=zeros(m,n);
for i=1:t %置乱次数
for j=1:m
for k=1:n
%由于matlab索引不能为0,所以为防止取余后索引错误,将坐标系变换0~N-1,并+1
% x,y为变换后的坐标
x=mod((j-1)+a*(k-1),N)+1;
y=mod(b*(j-1)+(a*b+1)*(k-1),N)+1;
imgn(j,k)=img(x,y);
end
end
img=imgn;
end
imgn=uint8(imgn);
return;
reAnroid.m
% Arnold置乱还原
function imgn=reArnold(img,a,b,t)
%{
img 载体图像
a,b 置乱矩阵参数
t 置乱次数
%}
%读取图像尺寸,Arnold变换针对方阵,所以一般m=n
% m和n为图像的行和列
[m,n,p]=size(img);
%方阵边长
N=m;
%还原
imgn=zeros(m,n);
for i=1:t %置乱次数
for j=1:m
for k=1:n
%由于matlab索引不能为0,所以为防止取余后索引错误,将坐标系变换0~N-1,并+1
% x,y为变换后的坐标
x=mod((a*b+1)*(j-1)+(-a)*(k-1),N)+1;
y=mod(-b*(j-1)+(k-1),N)+1;
imgn(j,k)=img(x,y);
end
end
img=imgn;
end
return;