数字水印复习.doc

数字水印整理的作业题和书本知识归纳总结，希望对29号参加考试的同学会 作业一：读入一幅图像，截取100*200的矩形区域，并旋转。 I=imread(xinxing.jpg); //读入图像 imshow(I) //显示图像 title(原始图像) //显示图像的名字，可以不写 G=rgb2gray(I); //将彩色图像转换成灰度图像 imshow(G) //显示图像 title(灰度图像) //显示图像的名字，可以不写 L=imcrop(G,[20 0 100 200]); //截取100*200的矩形区域，20 和 0是截取的位置，可以自己设定，但要保证截取的图片足够大 title(100*200矩形区域) //显示图像的名字，可以不写 M=imrotate(L,30); //逆时针旋转30°，如果顺时针需在30前加一个负号 imshow(M) //显示图像 title(旋转30°) //显示图像的名字，可以不写 作业二：对一幅图像求傅里叶变换和逆变换 I=imread(xiang.jpg); imshow(I) title(原始图像) G=fft2(I); imshow(G) title(傅里叶变换) H=fftshift(G); imshow(H) title(傅里叶逆变换) I=imread(xiang.jpg); //读入彩图 K=rgb2gray(I); //将读入的彩色图像阮欢成灰度图像 imshow(K) title(灰度图像) [cA1,cH1,cV1,cD1]=dwt2(K,db1); //固定语句，只有K可以改，跟灰度图像的名称一样就行 //一下四句是显示图像语句，后面的title语句都可以不写，221表示4个图像分两行两列显示，1、2、3、4分别代表第一幅图，第二幅图。。。。。。第四幅图 subplot(221);imshow(cA1,[]);title(近似水平系数) subplot(222);imshow(cH1,[]);title(水平细节系数) subplot(223);imshow(cV1,[]);title(垂直细节系数) subplot(224);imshow(cD1,[]);title(对角细节系数) K=idwt2(cA1,cH1,cV1,cD1,wname,[ ]) //我也不懂 作业三：将水印图像嵌入到原始载体中，并提取出水印信息 a=imread(xi.jpg); figure,imshow(a) b=imread(xii.png); [M,N]=size(b) for i=1:M for j=1:N if double(b(i,j))==0 w(i,j)=0; else w(i,j)=1; end end end //for 语句是将水印信息二值化 s=a; for i=1:M for j=1:N s(i,j)=bitset(a(i,j),4,w(i,j)); //将水印图像嵌入到载体的第四位，注意a(i,j)是指原始载体，w(i,j)是指水印图像，一定要分清二者的位置 end end imwrite(s,4wei.bmp); figure,imshow(s) for i=1:M for j=1:N mm(i,j)=bitget(uint8(s(i,j)),4);//提取水印信息，注意上面的水印嵌入到第几位，提取时一定要是相应的位数 if double(mm(i,j))==1 mm(i,j)=255; else mm(i,j)=0; end end end imwrite(mm,A4.bmp); figure,imshow(mm) 作业四：比较嵌入水印后的图像和原始图像的各种指标： //本次作业的图像是在作业三的基础上完成的，而且需要是灰度图像，所以先把读入的图像转换成灰度图像再计算 b1=imread(b1.jpg); a1=rgb2gray(b1); b2=imread(b2.jpg); an=rgb2gray(b2); a1 = double(a1); an=double(an); [M,N]=size(a1); a1n=(sum(sum(a1.^2))); a2n=(a1-an).^2; n=(sum(sum(a2n))); SNR=10*l