离散傅里叶、离散余弦、离散小波水印实验报告.doc

实验报告
专业: 班级: 指导老师:
姓名: tcy 学号: 实验室:
实验名称: 变换域数字图像水印算法时间:

实验目的及要求
(一)实验目的
结合本课程图像处理的基本知识,在离散傅里叶变换域(DFT)、离散余弦变换域(DCT)、离散小波变换域(DWT)等进行数字图像水印算法的研究。
(二)实验要求
(DFT)、离散余弦变换域(DCT)、离散小波变换域(DWT)、轮廓波变换域(contourlet transfrom)、曲波变换域(curvelet transfrom)等进行数字图像水印的嵌入及提取;
,提取出水印,并同原始水印进行必要的比较。
二、实验设备(环境)及要求
1、支持Intel Pentium Ⅲ及其以上CPU,内存256MB以上、硬盘1GB以上容量的微机;软件配有Windows98/2000/XP操作系统及MATLAB软件;
2、实验过程中,务必分析实验结果,按要求写出实验报告。(建议同时网上提交电子版实验报告:yw6895@)
三、实验内容与步骤
1、程序实现: x=(p-1)*K+1; y=(q-1)*K+1;
I_dct=I(x:x+K-1,y:y+K-1);
I_dct1=dct2(I_dct);
if J(p,q)==0
alfa=-1;
else
alfa=1;
end
I_dct2=I_dct1+alfa*E;
I_dct=idct2(I_dct2);
D(x:x+K-1,y:y+K-1)=I_dct;
end
end
clear all
size=256; N=32;K=8; D=zeros(size);E=;
I=imread('');
subplot(2,2,1);imshow(I);title('原始公开图像');
I=double(I)/512;
J=imread(''); %读入水印
subplot(2,2,2);imshow(J);title('水印图像');
J=double(imresize(J,[64,64])); %对水印缩放
subplot(2,2,3);imshow(J);title('缩放水印图像');

%嵌入水印
for p=1:size/K
for q=1:size/K
for q=1:size/K
x=(p-1)*K+1;
y=(q-1)*K+1;
I1=I(x:x+K-1,y:y+K-1);
I3=RR(x:x+K-1,y:y+K-1);
I_dct1=dct2(I1);
I_dct2=dct2(I3);
if I_dct2>I_dct1
W(p,q)=1;
else
W(p,q)=0;
end
end
end
figure
subplot(1,2,1);imshow(RR,[]);
title('旋转攻击');
subplot(1,2,2);imshow(W,[]);
title('图像中提取的水印');

% 剪切攻击提取水印
for p=1:size/K
for q=1:size