数字水印功能实现了两种嵌入水印的方法嵌入高斯序列和以可见.docVIP

数字水印 功能：实现了两种嵌入水印的方法：嵌入高斯序列和以可见方式嵌入固定文字。实现了两种攻击方法：剪切和对图像重置。 嵌入高斯序列方法中计算了相似度，并画出了相似度检测图；计算了峰值信噪比。 在命令窗口输入Menu，会出现如下窗口： （一） 首先嵌入高斯序列，点击之后在命令窗口出现提示：“input the key please(0-100):”和“input the a please”。输入一个符合要求的数字作为密钥，并以此为种子产生高斯序列。输入34和0.6之后： 命令窗口显示： ans = The lodge intensity is 0.6 ans = The similarity is 5696.54 ans = The PSNR is 25.0525 第一列是嵌入深度，第二列是相似度，第三列是峰值信噪比。 嵌入水印前后的图像及相似度检测图如下所示： 峰值信噪比随着潜入深度的变化： （二） 在总菜单上点击‘嵌入姓名’会嵌入我的名字： 图像如下所示： （三） 选择‘攻击水印’会出现一个子菜单，如下所示： 本菜单实现剪切方法攻击图像，剪切右半部分： 首先出现提示“请输入剪切的比例（百分数0—100） :”输入一个数字，如5 出现： 如果出入30： 剪切下半部分，同样输入5： 点击返回上级菜单之后返回到总菜单。 （四） 选择“置入高斯序列噪声” 选择抖动攻击： 选择滤波攻击： 选择旋转攻击，输入旋转度数15： 点击退出，退出程序。 代码部分 总菜单： %Menu.m k=menu( 请 选 择 水 印 输 入 方 法 , 嵌 入 高 斯 序 列 , 嵌 入 姓 名 , 攻 击 水 印 , 退出 ); switch k case 1, ren1 case 2, ren2 case 3, Attracks otherwise exit end; 嵌入水印，提取水印，作图 %ren1.m %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%嵌入水印 clear close all I=imread(lena.bmp); figure; subplot(1,3,1); imshow(I); title(the initial); key=input(input the key please: ); save keydata key; randn(seed,key);%%%%%%%%%%%%%产生高斯序列 W=randn(1,1024); W=mod(W,256); save gause W; a=input(input the a please: ); save adata a; k=1; for i=1:32 for j=1:32 F{i,j}=I( ((i-1)*8+1):(i*8) , ((j-1)*8+1):(j*8) ); D{i,j}=dct2(F{i,j}); D1{i,j}=D{i,j}; D1{i,j}(1,1)=D{i,j}(1,1)+a*W(k); k=k+1; end end M=zeros(256,256); for i=1:32 for j=1:32 D2{i,j}=idct2(D1{i,j}); M( ((i-1)*8+1):(i*8) , ((j-1)*8+1):(j*8) )=D2{i,j}; end end %%%%%%%%%%%%%%%%%%%%%两组for循环构成嵌入算法 M=round(M); M=M/255; subplot(1,3,2); imshow(M); title(With Print); imwrite(M,lena-Wprint.bmp); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%提取水印 P=imread(lena-Wprint.bmp); I=imread(lena.bmp); ww=zeros(1,1024); k=1; for i=1:32 for j=1:32 F1{i,j}=I( ((i-1)*8+1):(i*8) , ((j-1)*8+1):(j*8) );