yousun.rar_DCT watermark _compression lossy_水印

%有损压缩情形下水印的嵌入与提取 A=imread('lena.bmp'); subplot(2,2,1); imshow(A); A1=dct2(A); %对原始图像进行dct变换A1 [M,N]=size(A); [C,Q]=sort(A1); %对A调用函数进行变换 B=imread('dragon.jpg'); %B作为水印图像，进行嵌入 subplot(2,2,2); imshow(B); B1=dct2(B); %对水印图像进行dct变换B1 [m,n]=size(B); a1=0.01; %所设定的系数 for i=1:m %得到嵌入水印的图像的DCT变换 D(i,j) for j=1:n D1(i,j)=C(M-m+i,N-n+j)*(1+a1*B1(i,j)); end end for i=1:(M-m) for j=1:N D(i,j)=C(i,j); end end for i=(M-m+1):M for j=1:(N-n) D(i,j)=C(i,j); end for j=(N-n+1):N D(i,j)=D1(i-M+m,j-N+n); end end %对D(i,j)，利用记录矩阵Q进行逆排序，得到复原的DCT系数矩阵R for j=1:N for i=1:M R(Q(i,j),j)=D(i,j); end end E=idct2(R); %再进行取逆DCT，得到嵌入水印的图像E E1=E; E=uint8(E); subplot(2,2,3); imshow(E); R1=dct2(E1); %水印图像E，在图像压缩的情形下提取水印 for i=1:M for j=1:N if abs(R1(i,j))<1 R1(i,j)=0; end end end for j=1:N %水印提取,利用记录矩阵Q，首先恢复为SORT函数前的情形 for i=1:M F(i,j)=R1(Q(i,j),j); end end for i=1:m for j=1:n F1(i,j)=F(M-m+i,N-n+j); end end for i=1:m for j=1:n G1(i,j)=(F1(i,j)-C(M-m+i,N-n+j))/C(M-m+i,N-n+j)/a1; end end G=idct2(G1); G=uint8(G); %提取的水印为G subplot(2,2,4); imshow(G); k=1; %水印B与G的评价，将两个二维矩阵化为两个一维数组w1， w2 for i=1:m for j=1:n w1(k)=B(i,j); w2(k)=G(i,j); k=k+1; end end c=dot(w1,w2); %通常意义下的评价方法得到g1，g2，g3 c1=dot(w1,w1); c2=dot(w2,w2); g1=sqrt(c^2/(c1*c2)); g2=c/c1; g3=c/c2; for i=1:(k-1) %利用文章提到的改进方法进行评价得到p1 e(i)=1; end a=max(w1); b=min(w1); m1=(a+b)/2; for i=1:(m*n) s(i)=w1(i)-m1*e(i); t(i)=w2(i)-m1*e(i); end d=dot(s,t); d1=dot(s,s); d2=dot(t,t); p1=d/(sqrt(d1)*sqrt(d2));