function imw=wembed(im,w,delta,numblk)
%  分塊 DCT 和 量化的盲水印嵌入
%------------------------NOTES---------------------------------------------
% Written by : zj  - 01/05/07
% Contact    : zxxyj0506@sina.com
% CpoyRight@2007-2010
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fun1=@dct2;
im_dct=blkproc(im,[8 8],fun1);
for k=1:numblk
    for p=1:numblk
         x=(k-1)*8+1; y=(p-1)*8+1;
         cur_dct00=im_dct(x,y);
         akp=round(cur_dct00/delta);   % akp 是量化值
         if mod(akp+w(k,p),2)==1
             im_dct(x,y)=(akp-1/2)*delta;
         else
             im_dct(x,y)=(akp+1/2)*delta;
         end
    end
end
fun2=@idct2;
imw=blkproc(im_dct,[8 8],fun2);