混沌系统在图像列的方向上加密

function chaoticDemoForColumnChange
%应用混沌序列对数字图像的每一列置乱
%程序设计：李立宗
%2012年3月3日
% http://blog.youkuaiyun.com/superdont
% lilizong【at】Gmail
clear;
clc;
[filename pathname]=uigetfile('*.bmp;*.tiff;*.tif', '读入图像');
os=[pathname filename];
%os原始图像的路径
o=imread(os);
% figure,imshow(o);

%计算原始图像的大小
[m n]=size(o);
key1=n;
l=zeros(1,key1);
bl=zeros(1,key1);
l(1)=0.539;
for i=2:key1
    l(i)=1-2*l(i-1)*l(i-1);
    bl(i)=rem(round(l(i)*100000),key1);
    
end
%加密图像
os=o;
for i=1:n
    %     i
    os(:,i)=circshift(os(:,i),bl(i)); 
 end
%开始解密图像
or=os;

for i=1:n
    or(:,i)=circshift(or(:,i),m-bl(i));
end
figure,
subplot(1,3,1),imshow(o),title('原始图像');
subplot(1,3,2),imshow(os,[]),title('加密图像');
subplot(1,3,3),imshow(or,[]),title('解密图像');
% 因为放在一个窗口内图像比较小，不便于测试，如下代码用于测试。
figure,imshow(o);
figure,imshow(os,[]);
figure,imshow(or,[]);