zhiluan.rar_水印预处理

clc close all clear all % 保存开始时间 % iTimes=4; %置乱次数 iTimes = input('Please input Arnold scramble times:'); % 读入水印图像 file_name='c:\6.jpg'; I=imread(file_name); I=rgb2gray(I); message=im2double(I); %%水印图像的行数与列数 Mm=size(message,1); %水印行数 Nm=size(message,2); %水印列数 %对水印图像进行arnold置乱 tempImg=message; for n=1:iTimes % 次数 for u=1:Mm for v=1:Nm temp=tempImg(u,v); ax=mod((u-1)+(v-1),Mm)+1; ay=mod((u-1)+2*(v-1),Nm)+1; outImg(ax,ay)=temp; end end tempImg=outImg; end zhiluantu=tempImg; message_vector1=reshape(message,1,Mm*Nm); %%%%%%%%%%%%%置乱复原 iTimes=60; tempImg=reshape(message_vector1(1:Mm*Nm),Mm,Nm); %tempImg=message_vector; message_arnold=tempImg; for n=1:iTimes % 次数 for u=1:Mm for v=1:Nm temp=tempImg(u,v); ax=mod((u-1)+(v-1),Mm)+1; ay=mod((u-1)+2*(v-1),Nm)+1; outImg(ax,ay)=temp; end end tempImg=outImg; end message=outImg; % % 显示提取水印与原始水印 figure subplot(2,2,1); imshow(message,[]); title('复原水印'); subplot(2,2,2); imshow(zhiluantu,[]); title('置乱水印'); subplot(2,2,3); imshow(message,[]); title('复原水印'); subplot(2,2,4); imshow(I); title('原始水印');