MATLAB实现Anroid置乱与复原

Anroid.m

% 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;