wabsract.m

用matlab实现的DCT水印处理

function aw=wabsract(imw,delta)
% 分块 DCT 和 量化的盲水印提取
%------------------------NOTES---------------------------------------------
% Written by : zj  - 01/05/07
% Contact    : [e-mail]
% CpoyRight@2007-2010
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[imwrow,imwcol]=size(imw);
if imwrow<=1||imwcol<=1
    error('输入测试图像数据必须是一个矩阵');
end 
fun1=@dct2;
imw_dct=blkproc(imw,[8 8],fun1);
numblk=imwrow/8;
for k=1:numblk
    for p=1:numblk
         x=(k-1)*8+1; y=(p-1)*8+1;
         cur_dct00=imw_dct(x,y);
         akp=floor(cur_dct00/delta);   % akp 是量化值
         if mod(akp,2)==1
             aw(k,p)=1;
         else
             aw(k,p)=0;
         end
    end
end