小波实验报告wt水印.docVIP

一、 题目：dwt水印 二、 目的：编程实现水印的嵌入与还原 三、 算法及其实现：水印的嵌入与还原 小波变换具有空间-频率的多尺度性，对阁像的分解可以连续地从低分辨率过渡到高分 辨率，也可以对阁像整体进行变换，DWT还具柯多分辨率特性，可以很好地旬人类视觉系 统(HVS)特性相匹配。 本算法将载体图像进行小波分解，利用公式实现水印信息的嵌入。将2值图像中0和1 f-l Pi 參 k 的值变成-1和1。1 1 一 。将水印入宿主图形￡耍的小波系数上。 水印的提取完全是水印嵌入的逆过程。 水印提取是为了从待检测图像中恢复出己绎存在的水 印信息.其过程是水印嵌入过程的逆过程*贫先对嵌入水印 的图像做三展高敵小波变換.得到低播重要系数值 LL\x.y) 根据标志矩阵Aw获得原始图像的低频業要系 数值 若 LU,(x,y)—LL,(^y O .則 W = 1 :否 贝根据密钥产生二值伪随机序列. 贝根据密钥产生二值伪随机序列. 由达个二值伪 随m序歹J对W做伪随杉L处理后.得到水印图像W。 四、 实现工具：Matlab 五、 程序代码： (1) dwtl_cmbcd.m %嵌入源码 clc clear all; %保存开始吋间 start_time=cputime; k:20; blocksize=8; %设置块的人小 %读入原图像 file_name=,zimu.bmp,; cover_object=double(imread(file_name)); %原阁像的行数与列数 Mc=size(cover_object, 1); %原图像行数 Nc=size(cover_object,2); %原图像列数 %最人嵌入信息撒 max_message=Mc*Nc/(blocksizeA2); %读入水印图像 file_name=,c.bmp; mcssagc=doublc(imrcad(filc_namc)) %%水印图像的行数与列数 Mm=size(message, 1); % 水印行数 Nm=size(message,2); %水印列数 mcssagc_vcctor=rc$hapc(mcssagc,l %检查水印信总是否过人 if (Mm*Nm max_message) erroif水印太大 end %将随机数发牛器的状态置为1100 randn(’state’，l 100); %产生watermarked_image，并写入原图信思 watermarked」mage=cover一object; %将图像分块嵌入 % 当(2,2) (2,3)且 message_pad(kk)=O % 当(2,2)(2,3)且 message_pad(kk)=l %%经过分析可以得出结论：在提取水印吋，如果cDl(2,2)〉cDl(2,3)便足嵌入了水印的黑色 像素， %%反之则是嵌入了內色像素 x=l; y=i; h=waitbar(0/嵌入水印，请等待); for (kk = 1 :length(message_vector)) %对块进行DWT变换 IcA1,cH1,cV1,cD1J = dwt2(cover_object(x:x+blocksize-l，y :y+blocksize-1 ),haar); a=cD 1; % 如果 message_pad(kk)== 0 if (message一vector(kk) == 0) %且(2,2) (2,3)，交换它们 if (cD 1 (2,2) cD 1(2,3)) tcmp=cD 1(2,3)； cDl(2,3)=cDl(2,2); cDl(2,2)=temp; end % 如果 message_pad(kk) == 1, elseif (message_vector(kk) == 1) %且(2,2) (2,3),交换它们 if (cD 1(2,2) = cD 1(2,3)) temp=cD 1(2,3); cDl(2,3)=cD 1(2,2); cDl(2,2)=temp; end end %检查(2,2) , (2,3)的差是否〉二k if cD 1 (2,2) cD 1(2,3) ifcDl(2,2)-cDl(2,3)k cDl(2,2)=cDl(2,2)+(k/2); cDl(2,3)=cDl(2,3)-(k/2); end else ifcDl(2,3)-cDl(2,2)k cDl(2,3)=cDl(2,3)+(k/2); cDl(2,2)=cDl(2,2)-(k/2); end end %IDWT watermarked_image(x: x+blocksize-1，y:y+blocksize-1 )= idwt2(cA 1，cHl，cV l，cD 1，liaar，[Mc，Nc