DCT域图像数字水印算法.docVIP

DCT域图像数字水印算法 设计目的 利用所学的计算机语言进行编程，实现图像的DCT域变换，并利用相应的数字水印算法实现水印的嵌入提取。通过对该实验的完成，来提高对语言的掌握程度，提高自身的编程水平。在了解该算法的同时，可以提出一些自己的观点，来加深对该问题的了解。 设计任务 1.了解数字图像水印的概念、特点、分类、及其应用； 2.掌握数字图像水印的原理； 3.理解数字图像水印嵌入和检测过程； 4.理解数字图像水印的攻击； 利用C语言或MATLAB编程，在DCT域的低频、中频、高频段分别实现数字图像水印的嵌入和提取。 三．设计原理 在DCT域,不同的DCT 系数作为水印载体对水印的鲁棒性有不同的影响。为使水印具有较好的稳健性,用来嵌入水印的DCT 系数应满足如下条件: (1) 在经过常见信号处理和噪声干扰后仍能很好地保留,即这些DCT系数不应过多地为信号处理和噪声干扰所改变; (2) 具有较大的感觉容量以便嵌入水印后不会引起视觉质量的明显改变。 第一个要求是为了保证水印在嵌入图像后有较大的稳健性。当加入水印的DCT 系数被改变较小时,水印便可能被保留,这是显然的。第二个要求是同时针对不可见性和稳健性而提出的,较大的感觉容量意味着在主观视觉效果不变的前提下有较大的改变裕度,这也意味着可以嵌入较强的水印信号。 水印嵌入位置在哪才能保证其稳健性？Cox提出了水印信息应当嵌入到在视觉系统感觉上最重要的分量上，即低频系数上。也有学者提出在保证水印不可见性的前提下，把水印嵌入到分量上，理由是系数比低频系数大得多，对一个较大数作较小的改变，并不会引起太大的变化。为了兼顾水印的稳健性和不可见性，也有文献提出将水印嵌入到DCT中频系数中。 本实验主要验证将水印嵌入到载体图像的系数、低频系数、中频系数和高频系数中的性能。 直流系数嵌入水印实验原理 选择二值图像作为水印信息,根据水印图像的二值性选择不同的嵌入系数,并将载体图像进行的分块,将数字水印的值直接植入到载体灰度图像的DCT变换域中直流系数（指的是DCT系数最左上角系数）,实现水印的嵌入。具体方法如下: 1) 读取原始公开图像和二值水印图像； 2) 将原始公开图像分割为互不覆盖的图像块，,,对 进行以DCT变换,得到； 3) 取二值水印图像中的一个元素嵌入原始公开图像块的直流系数中；嵌入公式如下：。其中是嵌入强度，是嵌入位置，直流系数。 4) 对嵌入水印信息后的图像块进行反DCT 变换； 5) 合并图像块，得到嵌入二值水印后的图像。 水印提取过程可以描述如下: 1) 读取原始公开图像，和嵌入水印后的图像； 2) 将原始公开图像和嵌入水印后的图像分别完全分割为互不覆盖的图像块和,, ,并对他们进行DCT变换； 3) 将相对应的图像块和，进行如下的运算: 4) 用得到的值进行如下的判定: 其中,p和q为水印图像的大小 ,即提取的水印图像。 低频系数嵌入水印 低频系数嵌入水印的原理与直流系数相同，只是嵌入到DCT系数的低频位置。 中频系数嵌入水印 中频系数嵌入水印的原理与直流系数相同，只是嵌入到DCT系数的低频位置。 高频系数嵌入水印 高频系数嵌入水印的原理与直流系数相同，只是嵌入到DCT系数的高频位置。 实验中，对于加入水印后的图像的检测采用峰值信噪比(PSNR),对于水印的相似性检测用归一化互相关系数(NC)，计算公式如下： (7) 其中，为原始图像，为含水印图像，M,N为图像的行，列 （8） 其中，为原始水印，为提取出的水印。 四．实验内容 查阅相关数字水印综述文章（至少3篇），熟悉数字图像水印的概念、特点、分类、及其应用，掌握一个完整的数字水印系统的基本组成。 MATLAB编程，根据上述原理，在DCT域的直流，低频、中频、高频段分别实现数字图像水印的嵌入和提取。 编写峰值信噪比PSNR及归一化相关系数NC的子函数，对算法的性能进行客观评价。 编写含水印图像的常见攻击算法（压缩，剪切，加噪，滤波，旋转等），计算攻击后图像的峰值信噪比及提取的水印的归一化相关系数，观察其变化。 改变水印的嵌入域，分别将水印嵌入于图像DCT系数的直流，低频，中频和高频部分，比较水印嵌入后图像的效果及抗攻击性能。 实验子程序简介 B=dct2(A,M,N) 表示对二维图像矩阵A进行DCT变换 其中，A表示要变换的图像，M和N是可选参数，表示填充后的图像矩阵大小。B表示变换后得到的系数矩阵。 2.imcrop 图像剪切函数 3.imrotate 图像旋转函数 4.Imnoise 图像加噪函数 思考题 在保证含水印图像透明性的前提下，计算出的PSNR值为多少