dct.m

扩频水印技术的MATLAB程序

clear
%
%读入原图像，设置参数
trueImage=imread('G:\dct\lena.bmp');
alfa=.1;
LENGTH=2500;
subplot(2,2,1);
imshow(trueImage);
title('原始图象');
%
%对原图象进行DCT变换
dctF1=dct2(trueImage);
subplot(2,2,2);
imshow(log(abs(dctF1)),[ ]);
title('DCT cofficient matrix');
[m,n]=size(dctF1);
%
%产生水印序列并对其排序
randn('state',10);
watermark1=randn(LENGTH,1);
subplot(2,2,3);
plot(watermark1);
title('watermark seqence')
[Y0,I0]=sort(watermark1);
%
%找出水印嵌入位置（幅值较大的n个频域成分）
A=dctF1(:);
[Y1,I1]=sort(A);
x=m*n;
k=LENGTH;
M=zeros(x,1);
%
%修改幅值较大的n个频域成分的幅值，嵌入水印（因为两个问题不同，所以有两个注释符）
for i=1:x
   if k>=1
    M(i)=Y1(i)*(1+alfa*Y0(k));
    k=k-1;
   else
    M(i)=Y1(i);
   end
   i=i+1;
end
N=zeros(x,1);
x=m*n;
for i=1:x
    N(I1(i))=M(i);
end
a=1;
for j=1:n
  for i=1:m
    dctF2(i,j)=N(a);
    a=a+1;
  end
end
%
%DCT反变换，得到嵌入水印的图象
idctF1=idct2(dctF2);
subplot(2,2,4);
imshow(idctF1,[]);
title('嵌入水印后的图像');