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

:..实验报告专业:班级:指导老师:姓名:tcy学号:实验室:实验名称:变换域数字图像水印算法时间:一、实验目的及要求(一)实验目的结合本课程图像处理的基本知识,在离散傅里叶变换域(DFT)、离散余弦变换域(DCT)、离散小波变换域(DWT)等进行数字图像水印算法的研究。(二)(DFT)、离散余弦变换域(DCT)、离散小波变换域(DWT)、轮廓波变换域(contourlettransfrom)、曲波变换域(curvelettransfrom)等进行数字图像水印的嵌入及提取;,提取出水印,并同原始水印进行必要的比较。二、实验设备(环境)及要求1、支持IntelPentiumⅢ及其以上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);ifJ(p,q)==0alfa=-1;elsealfa=1;endI_dct2=I_dct1+alfa*E;I_dct=idct2(I_dct2);D(x:x+K-1,y:y+K-1)=I_dct;endendclearallsize=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('缩放水印图像');%嵌入水印forp=1:size/Kforq=1:size/Kforq=1:size/Kx=(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);ifI_dct2>I_dct1W(p,q)=1;elseW(p,q)=0;endendendfiguresubplot(1,2,1);imshow(RR,[]);title('旋转攻击');subplot(1,2,2);imshow(W,[]);title('图像中提取的水印');%剪切攻击提取水印forp=1:size/Kforq=1:size/Kx=(p-1)*K+1;y=(q-1)*K+1;I1=I(x:x+K-1,y:y+K-1);I4=PP(x:x+K-1,y:y+K-1);I_dct1=dct2(I1);I_dct2=dct2(I4);ifI_dct2>I_dct1W(p,q)=1;elseW(p,q