DCT域图像水印技术.doc

《DCT域图像水印技术.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、DCT域图像水印技术%1.实验目的：使用MATLAB软件，熟悉DCT域图像水印技术，并学会该技术将特定的内容嵌入图像中。%1.实验原理：离散余弦变化(DiscreteConsineTransform)简称DCT。任何连续的实对称函数的傅立叶变换屮貝含有余弦项，因此余弦变换与傅立叶变换一样有明确的物理意义，DCT变换避免了傅立叶变换屮的复数运算，它是基于实数的正交变换。DCT变换矩阵的基向量很近似于Toeplitz矩阵(系数矩阵对称且沿着与主对角线平行的任意一对角线上的元素都相等)的特征向量，而Toeplit

2、z矩阵乂体现了人类语言及图像信号的和关特性。故DCT常常被认为是对语音和图像编码的最佳变换，同吋DCT算法较易于在数字信号处理器屮快速实现，因此它H前在图像编码屮占有重要的地位。二维DCT变换是H前最常使用的有损数字图像床缩系统——JPEG系统的核心。在本次实验中我们用了matlab内嵌的函数dct2()和idct2(),来进行det变化和反变换。大大简化的试验的过程。与空域图像水印相比，DCT域图像水印对压缩，滤波和其他-•些数字处理算子具有更强的稳健性，同时乂与常用的图像压缩标准JPEG兼容，因而得到了

3、广泛的重视，基于DCT的数字水印技术是FI前水印技术屮研究的最多，垠深入,而且也是最成熟的。本次试验为了简化过程，所以没有进行通常的图像分块DCT的过程，而是直接寻找DCT变换值小于门限的位置，并在这些位置小嵌入水印。以图像Lena为载体图像，大小为256X256。shuiyin.txt为存放水印的文件(复旦大学数字水印2006年12月22口)。DCT嵌入水印及水印提取的过程如下：1•读原始图像和水印文件。2.图像进行DCT变换。3.根据水印的字符多少來确定DCT门限值的大小。4.根据DCT门限值来确定图像

4、屮要嵌入水印的位置。5.在变换后的图像屮嵌入水印。6.用之前得到的嵌入水印的位置提取出水印。%1.实验内容：以下是进行DCT域图像隐藏的程序。DCT_suiyinqnmclcfid=fopen(,shuiyin.txt,；rt,);%打开shuiyin.txt文件，对其进行读操作[A,count]=fread(fid);%将标识符为fid的文件以二进制的方式全部读入A中%count中存放的是A的大小fclose(fid);%关闭文件RGB=imread(fLenna.bmpf);figure(I),imsh

5、ow(RGB);I=rgb2gray(RGB);figure(2)4mshow(I);J=dct2(I);[mn]=size(J);装入图像%显示图像%将真彩图像转换为灰度图像%显示灰度图像%进行余弦变换计算矩阵J的大小%%%%%%%%%%%%%%%%%%%%%%%%%%%根据shuiyin文件中的字符数，计算DCT变换的门限值的大小，并放在num变量中ifcount<=20num=0.00I;elseifcount<=44num=0.002;elseifcount<=l16num=0.005;elseif

6、count<=224num=0.01;elsenum=0.1;endendendend%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%记录DCT变换值小于num的位置%并将DCT变换值小于num的元素设为0j=0；fori=l:m*nifabs(J(i))

7、%%%%%%%%%%%%%%%%%%%%%%将A中字符逐一嵌入J中fori=l:countJ(s(i))=A(i);end%%%%%%%%%%%%%%%%%%%%%figure(4),imshow(K);显示嵌入水印后的图勺imwrite(K,‘result.jpg；jpg')；fori=l:countJ(s(i))end提取20个嵌入的数字四.结论与分析:Lenna.bmptempl.jpgtemp2.jpgresult.jpgtempl.jpg为直接由Lenna.bmp转化而来的灰度图像。temp2.j

8、pg为将templ.jpg中DCT变换值小于num的元素设为0后的图像。result.jpg为嵌入水印后的结果图像。可以看出，后三幅图片用肉眼完全看不到差别。说明DCT值小于一个较小的数的数F1很小，基本对图片不起作用，大概要将当DCT值小于10的都置为0时，图像可以看到比较明显的模糊。