DCT数据隐藏实验报告(1)

1、中南大学Central South University信息隐藏实验报告学院：信息科学与工程学院班级：信息安全1201学号：0909121724姓名:吕秋言2015年6月时间：实验二：离散余弦变换(DCT)算法隐藏信息一：实验目的：掌握在频率域隐藏信息，在木实验屮采用离散余弦变换(DCT)算法隐藏信息，同时理解DCT在图像处理屮的应用原理。二：实验内容：(1)预备知识：掌握DCT的工作原理a二维离散余弦变换矩阵形式：正变换：F=DfD#反变换：f=D，FD产生 DCT 矩阵的 MATLAB 函数：D=dctmtx(N);b. 二维DCT变换。格式：B = dct2 (A)B = dct2 (A

2、, m, n)B = dct2 (A, m n)说明:B = dct2 (A)i+算A的DCT变换B, A与B的大小相同;B = dct2 (A, m, n)和B二dct2 (A, m n)通过对A补0或剪裁，使B的大小为mxrioc. DCT反变换。格式：B = idct2 (A)B = idct2 (A, m, n)B = idct2 (A, m n)d. 计算DCT变换矩阵。格式：D = dctmtx (n)说明:D = dctmtx (n)返回-个nxn的DCT变换矩阵，输出矩阵D为double类型。F=DfD*说明二维余弦正反变换在Matlab中的实现。 RGB=imread(aut

3、umn.tif)； l=rgb2gray(RGB); figure(l); imshow(l); figure(2); J=dct2(l); imshow(log(abs(J)4); colormap(jet(64),colorbar; figure(3); J(abs(J)I_dctlW(p,q)=l;elseW(p,q)=0;endendendfigure;subplot( 1,2,1);imshow(RJl);title(低通滤波攻击工subplot( 1,2,2);imshow(W,);title(1|像中提取的水印);%WJh 7K印相彳以比校J=double(J);W=double

4、(W);sumJ=0;sumW=0;for j=l:Nfor i=l :NsumJ=sumJ+J(j,i)*W(j,i);sumW=sumW+W(j,i)*W(j,i)；endendCH=(sumW-sumJ)/sumJ基于DTC数字水印算法的Matlab实现源代码:二256;%原图像长度232;%水印图像长度K二出I二zeros (M, M);II=zeros (K, K);B=zeros(M, M);Idct=zeros(K, K);D二zeros (M, M);%读取原图像I=imread ( 33. png);subplot (2, 2, 1);%显示原图像imshow(I);titl

5、eC原图像)；%产生水印序列randn(，seed*, 10);mark=randn(1024,1);subplot (2, 2, 2);plot (mark);titleC水印序列)；%嵌入水印T=l;for m=l:Nfor n二1:Nx=(mT)*K+l;y二(nT)*K+l;II=I (x:x+K-l, y:y+K-1);%将原图分成 8*8 的子块Idct二dct2(II)珊对子块进行DCT变换if x=l&y=lalfa=0. 002;elsealfa=0;endB=Idct*仃 +alfa*mark (T) ;%嵌入水印Bidct=idct2 (B);%进行 DCT 反变换I (

6、x:x+KT, y:y+K-l)=Bidct;T二T+l;endendsubplot (2, 2, 3);imshow(I) ；%显示嵌入水印后的图像title(，tu) ; imwrite(T,嵌入后的图像.bmp);%进行相关性比较figure;for i=l:50if i=10;mark2=markl，;elsemark2=randn(1024, 1);end%计算相关值c=(mark2, *mark)/sqrt(mark2, *mark2);stem(i, c);hold on;endRGB=imread(*psb.tif);I=rgb2gray(RGB);subplot( 1,2,1);imshow(I);title(ystx);J=dct2(I);sub