小波实验报告dwt水印

《小波实验报告dwt水印》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、一、题目：dwt水印二、目的：编程实现水印的嵌入与还原三、算法及其实现：水印的嵌入与还原小波变换具有空间-频率的多尺度性，对阁像的分解可以连续地从低分辨率过渡到高分辨率，也可以对阁像整体进行变换，DWT还具柯多分辨率特性，可以很好地旬人类视觉系统(HVS)特性相匹配。本算法将载体图像进行小波分解，利用公式实现水印信息的嵌入。将2值图像中0和1f-lPi參k的值变成-1和1。11一。将水印入宿主图形£耍的小波系数上。水印的提取完全是水印嵌入的逆过程。水印提取是为了从待检测图像中恢复出己绎存在的水印信息.其过程是水印嵌入过程的逆过程*贫先对嵌入水印的图像做三展高敵小波变換.得到低播重要系数

2、值LL>O.則W=1:否贝根据密钥产生二值伪随机序列.由达个二值伪随m序歹>J对W做伪随杉L处理后.得到水印图像W。四、实现工具：Matlab五、程序代码：(1)dwtl_cmbcd.m%嵌入源码clcclearall;%保存开始吋间start_time=cputime;k:20;blocksize=8;%设置块的人小%读入原图像file_name=,zimu.bmp,;cover_object=double(imread(file_name));%原阁像的行数与列数Mc=size(cover

3、_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%将随机数发牛器的状态置为1

4、100randn(’state’，l100);%产生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,cH

5、1,cV1,cD1J=dwt2(cover_object(x:x+blocksize-l，y:y+blocksize-1),'haar');a=cD1;%如果message_pad(kk)==0if(message一vector(kk)==0)%且(2,2)<(2,3)，交换它们if(cD1(2,2)(2,3),交换它们if(cD1(2,2)>=cD1(

6、2,3))temp=cD1(2,3);cDl(2,3)=cD1(2,2);cDl(2,2)=temp;endend%检查(2,2),(2,3)的差是否〉二kifcD1(2,2)>cD1(2,3)ifcDl(2,2)-cDl(2,3)

7、size-1)=idwt2(cA1，cHl，cVl，cD1，liaar'，[Mc，Nc]);%移动到下一块if(x+blocksize)>=Ncx=l;y=y+blocksize;elsex=x+block$izc;endwaitbar(kk/length(message_vector)，h);endclose(h);%转换为uint8并写入dwt_watermarked.bmpwatermarked_image_uint=uint8(wate