数字水印学习之混沌学习

xiaoer_1982

%%混沌置乱代码
  clear;
  clc;
  %os原始图像的路径
  [filename pathname]=uigetfile('*.bmp;*.tiff;*.tif', '读入图像');
  os=[pathname filename];
  o=imread(os);
  % figure,imshow(o);
  %计算原始图像的大小
  [m n]=size(o);
  %l=zeros(m,n);
  l=linspace(0,0,m*n);
  l(1)=0.98;
  for i=2:m*n
  l(i)=1-2*l(i-1)*l(i-1);
  end
  % %对生成的混沌序列进行排序
  % [lsort lindex]=sort(l);
  % t=linspace(0,0,m*n);
  % t(1)=0.98;
  % for i=1:m*n
  % t(i+1)=1-2*t(i)*t(i);
  % end
  % for i=1:m*n
  % if (l(i)>=0)
  % l(i)=1;
  % else
  % l(i)=0;
  % end
  % end
  % l=~l;
  % subplot(2,2,2),imshow(l),title('混沌图像');
  [lsort,lindex]=sort(l);
  %ozl =original zhi luan 表示置乱后的原始图像
  ozl=zeros(m,n);
  for i=1:m*n
  ozl(i)=o(lindex(i));
  end
  % ozl=reshape(ozl,m,n);
  %%oo =original fu yuan
  ofy=zeros(m,n);
  for i=1:m*n
  ofy(lindex(i))=ozl(i);
  end
  h= waitbar(0,'程序处理中，请耐心等待。。。');
  for i=1:100, % computation here %
  waitbar(i/100)
  end
  close(h) ;
  subplot(1,3,1),imshow(o),title('原始图像');
  subplot(1,3,2),imshow(ozl,[]),title('加密图像');
  subplot(1,3,3),imshow(ofy,[]),title('复原图像');
  %%%%
  clear;
  clc;
  %需要注意本程序实现对二进制文件的加密，对于非二进制文件不能采用此程序进行加密！！！！
  [filename pathname]=uigetfile('*.bmp;*.tiff;*.tif', '读入图像');
  os=[pathname filename];
  %os原始图像的路径
  o=imread(os);
  % figure,imshow(o);
  %计算原始图像的大小
  [m n]=size(o);
  l=zeros(m,n);
  l(1)=0.98;
  for i=2:m*n
  l(i)=1-2*l(i-1)*l(i-1);
  end
  % %对生成的混沌序列进行排序
  % [lsort lindex]=sort(l);
  % t=linspace(0,0,m*n);
  % t(1)=0.98;
  % for i=1:m*n
  % t(i+1)=1-2*t(i)*t(i);
  % end
  for i=1:m*n
  if (l(i)>=0)
  l(i)=1;
  else
  l(i)=0;
  end
  end
  % l=~l;
  ei=xor(o,l);
  di=xor(ei,l);
  h= waitbar(0,'程序处理中，请耐心等待。。。');
  for i=1:100, % computation here %
  waitbar(i/100)
  end
  close(h)
  subplot(2,2,1),imshow(o),title('原始图像');
  subplot(2,2,2),imshow(l),title('混沌图像');
  % figure,imshow(l);
  %ei表示加密图像encryption image
  subplot(2,2,3),imshow(ei),title('加密图像');
  %di表示解密图像dencryption image
  subplot(2,2,4),imshow(di,[]),title('解密图像');
  % figure,imshow(di);
  clear;
  clc;
  % o=imread('image/lena.bmp');
  %%%%%%%%%%%%读入混沌初始值
  prompt={'混沌初始值 '};
  title='初值';
  lines=1;
  def={'98'};
  mya=inputdlg(prompt,title,lines,def,'on');
  key11=['0.',mya{1}];
  key=str2num(key11);
  %%%%%%%%%%%%
  [filename pathname]=uiputfile({'*.bmp';'*.tiff';'*.tif'}, '保存图像');
  os=[pathname filename];
  m=512;
  n=512;
  l=zeros(m,n);
  l(1)=key;
  for i=2:m*n
  l(i)=1-2*l(i-1)*l(i-1);
  end
  %二值化的时候需要注意，此处的临界值是0，而不是0.5
  key=0;
  for i=1:m*n
  if (l(i)>=key)
  l(i)=0;
  else
  l(i)=1;
  end
  % waitbar((i/(m*n)),h);
  end
  l=~l;
  % imwrite(l,'image/wm.bmp');
  %进度条
  h= waitbar(0,'程序处理中，请耐心等待。。。');
  for i=1:100, % computation here %
  waitbar(i/100);
  end
  close(h) ;
  imwrite(l,os);
  % handles.figure;
  subplot(1,1,1),imshow(l),title('混沌图像');
  % test=imread(os);
  % figure,imshow(test,[]);
  clear;
  clc;
  % key=0.98
  %key为解密密钥
  %os原始图像的路径
  [filename pathname]=uigetfile('*.bmp;*.tiff;*.tif', '读入需要解密的图像');
  os=[pathname filename];
  ozl=imread(os);
  % figure,imshow(ozl,[]);
  %%%%%%%%%%%%%设置混沌初始值
  prompt={'混沌初始值 '};
  title='初值';
  lines=1;
  def={'98'};
  mya=inputdlg(prompt,title,lines,def,'on');
  key11=['0.',mya{1}];
  key=str2num(key11);
  %%%%%%%%%%%%设置混沌初始值
  %计算原始图像的大小
  [m n]=size(ozl);
  %l=zeros(m,n);
  l=linspace(0,0,m*n);
  l(1)=key;
  for i=2:m*n
  l(i)=1-2*l(i-1)*l(i-1);
  end
  % %对生成的混沌序列进行排序
  % [lsort lindex]=sort(l);
  % t=linspace(0,0,m*n);
  % t(1)=0.98;
  % for i=1:m*n
  % t(i+1)=1-2*t(i)*t(i);
  % end
  % for i=1:m*n
  % if (l(i)>=0)
  % l(i)=1;
  % else
  % l(i)=0;
  % end
  % end
  % l=~l;
  % subplot(2,2,2),imshow(l),title('混沌图像');
  [lsort,lindex]=sort(l);
  %ozl =original zhi luan 表示置乱后的原始图像
  % ozl=zeros(m,n);
  % for i=1:m*n
  % ozl(i)=o(lindex(i));
  % end
  % ozl=reshape(ozl,m,n);
  %%oo =original fu yuan
  ofy=zeros(m,n);
  for i=1:m*n
  ofy(lindex(i))=ozl(i);
  end
  h= waitbar(0,'程序处理中，请耐心等待。。。');
  for i=1:100, % computation here %
  waitbar(i/100)
  end
  close(h) ;
  subplot(1,2,1),imshow(ozl,[]);
  %,title('需要解密的图像');
  set(gca,'Title',text('String','需要解密的图像'))
  subplot(1,2,2),imshow(ofy,[]);
  %,title('解密后图像');
  set(gca,'Title',text('String','解密后图像'))
  % [filename pathname]=uiputfile({'*.bmp';'*.tiff';'*.tif'}, '保存加密后图像');
  % os=[pathname filename];
  % imwrite(ozl,os);
  % subplot(1,3,3),imshow(ofy,[]),title('复原图像');
  clear;
  clc;
  %key为加密密钥
  %os原始图像的路径
  [filename pathname]=uigetfile('*.bmp;*.tiff;*.tif', '读入需要加密的图像');
  os=[pathname filename];
  [filename pathname]=uiputfile({'*.bmp';'*.tiff';'*.tif'}, '保存加密后的图像');
  ozls=[pathname filename];
  % key=0.98;
  %%%%%%%%%%%%读入混沌初始值
  prompt={'混沌初始值 '};
  title='初值';
  lines=1;
  def={'98'};
  mya=inputdlg(prompt,title,lines,def,'on');
  key11=['0.',mya{1}];
  key=str2num(key11);
  %%%%%%%%%%%%
  o=imread(os);
  % figure,imshow(o);
  %计算原始图像的大小
  [m n]=size(o);
  %l=zeros(m,n);
  l=linspace(0,0,m*n);
  l(1)=key;
  for i=2:m*n
  l(i)=1-2*l(i-1)*l(i-1);
  end
  % %对生成的混沌序列进行排序
  % [lsort lindex]=sort(l);
  % t=linspace(0,0,m*n);
  % t(1)=0.98;
  % for i=1:m*n
  % t(i+1)=1-2*t(i)*t(i);
  % end
  % for i=1:m*n
  % if (l(i)>=0)
  % l(i)=1;
  % else
  % l(i)=0;
  % end
  % end
  % l=~l;
  % subplot(2,2,2),imshow(l),title('混沌图像');
  [lsort,lindex]=sort(l);
  %ozl =original zhi luan 表示置乱后的原始图像
  ozl=zeros(m,n);
  for i=1:m*n
  ozl(i)=o(lindex(i));
  end
  % ozl=reshape(ozl,m,n);
  %%oo =original fu yuan
  % ofy=zeros(m,n);
  % for i=1:m*n
  % ofy(lindex(i))=ozl(i);
  % end
  % h= waitbar(0,'程序处理中，请耐心等待。。。');
  % for i=1:100, % computation here %
  % waitbar(i/100)
  % end
  % close(h) ;
  %首先保存需要加密的图像
  ozl=uint8(ozl);
  ozl=double(ozl);
  imwrite(ozl,ozls);
  subplot(1,2,1)
  imshow(o)
  % title('原始图像')
  set(gca,'Title',text('String','原始图像'))
  subplot(1,2,2),imshow(ozl,[]);
  set(gca,'Title',text('String','加密后图像'))
  % title('加密后图像');
  % [filename pathname]=uiputfile({'*.bmp';'*.tiff';'*.tif'}, '保存加密后图像');
  % os=[pathname filename];
  % % [filename pathname]=uiputfile({'*.bmp';'*.tiff';'*.tif'}, '保存图像');
  % % os=[pathname filename];
  % imwrite(ozl,os);
  % subplot(1,3,3),imshow(ofy,[]),title('复原图像');
  % [filename pathname]=uiputfile({'*.bmp';'*.tiff';'*.tif'}, '保存加密后的图像');
  % os=[pathname filename];
  % %进度条
  % % h= waitbar(0,'程序处理中，请耐心等待。。。');
  % % for i=1:100, % computation here %
  % % waitbar(i/100);
  % % end
  % % close(h) ;
  % ozl=uint8(ozl);
  % imwrite(ozl,os);
  % figure,imshow(ozl,[]);
  % test=imread(os);
  % figure,imshow(test,[]);