快答辩了，求助基于离散余弦变换的数字水印matlab仿真源码

eagle217 · 发表于 2006-5-29 10:00

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

各位兄弟姐妹,小弟跪求如题所示源代码！若有高手也帮小弟调试一下程序，不胜感激！从水印提取以后就运行不好~~若有好的建议发到我邮箱：<a href="mailteagle217@tom.com target="_blank" >eagle217@tom.com 真诚感谢 </A> %读入原始图像，设置参数 trueImage=imread('hua.jpg'); trueImage=rgb2gray(trueImage); alfa=.1;LENGTH=2500; subplot(2,2,1); imshow(trueImage); title('original image');
 %对原图像进行DCT变换 dctF1=dct2(trueImage); subplot(2,2,2); imshow(log(abs(dctF1)),[]); title('DCT coefficient matrix'); [m,n]=size(dctF1);
 %产生水印序列并对其进行排序 randn('state',10); waterMark1=randn(LENGTH,1); subplot(2,2,3); plot(waterMark1); title('watermark sequence'); [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(x)=Y1(x)*(1+alfa*Y0(k)); k=k-1; else M(x)=Y1(x); end x=x-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('embedded image');

%水印提取逆过程是嵌入过程的逆过程 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%function watermark_detect(image,Y1,I0,waterMark1) %image:嵌入水印的图像 %Y1:原始图像的序列排序 %I0:原始水印的排列序号 %waterMark1):原始水印序列 %对嵌入水印图像进行DCT变换 dctW1=dct2(image); %找出幅值较大的系数 B=dctW1(1); [Y2,I2]=sort(B); [m1,n1]=size(dctW1); y=m1*n1; k=length(waterMark1); N0=zeros(k,1); %提取水印序列 while k>=1 N0(k)=(Y2(y)-Y1(y))/alfa/Y1(y); k=k-1; y=y-1; end k=length(waterMark1); waterMark2=zeros(k,1); for i=1:k waterMark2(I0(i))=N0(i); end %选取50个测试序列，其中第10个位提取出的水印 figure; for i=1:50 if i==10 waterMark=waterMark2; else waterMark=rand(k,1); end %计算各个序列与原水印序列的相关值 c=waterMark'*waterMark1/sqrt(waterMark'*watermark); stem(i,c); hold on; end %攻击实验 disp('input your choice according to the following image processing operation:'); disp('0--exit'); disp('1--smoothing patterns'); disp('2--adding uniform noise');%添加噪声 disp('3--adding fiter [10,10]');%滤波 disp('4--cutting part of the image');%剪切 disp('5--10 quality JPEG compression');%压缩 disp('6--rotate 45');%旋转 
d=input('please input your choice:') while d~=0 switch d case 1 watermark_detect(idctF1,Y1,I0,waterMark1); case 2 WImage2=idctF1; noise0=10*rand(size(WImage2)); WImage2=WImage2+noise0; figure; inshow(WImage2,[]); title('adding uniform noise'); watermark_detect(WImage2,Y1,I0,waterMark1); case 3 WImage3=idctF1; H=fspecial('gaussian',[10,10]'); WImage3=imfilter(WImage3,H); figure; imshow(WImage,[]); title('through filter [10,10]'); watermark_detect(WImage3,Y1,I0,waterMark1); case 4 WImage4=idctF1; WImage4(1:128,1:128)=256; figure; imshow(WImage4); title('cutting part of the image'); watermark_detect(WImage4,Y1,I0,waterMark1); case 5 Wimage5=idctF1; Wimage5=im2double(WImage5); cnum=10; dctm=dctmtx(8); P1=dctm; P2=dctm.'; imageDCT=blkproc(WImage5,[8,8],'P1*x*P2'dctm,dctm.'); DCTvar=im2col(imageDCT,[8,8],'distinct').'; n=size(DCTvar,1); DCTvar=(sum(DCTvar.*DCTvar)-(sum(DCTvar)/n).^2)/n; [dum,order]=sort(DCTvar); cnum=64-cnum; mask=ones(8,8); mask(order(1:cnum))=zeros(1,cnum); im88=zeros(9,9); im88(1:8,1:8)=mask; im128128=kron(im88(1:8,1:8),ones(16)); dctm=dctmtx(8); P1=dctm.'; P2=mask(1;8,1:8); P3=dctm; WImage5=blkproc(imageDCT,[8,8],'P1*(x.*P2)*P3',dctm.',mask(1:8,1:8),dctm); figure; imshow(WImage5); title('JPEG Image'); watermark_detect(WImage5,Y1,I0,waterMark1); case 6 WImage6=idctF1; WImage6=imrotate(WImage,45,'bilinear','crop'); figure; imshow(WImage6); title('rotate 45'); watermark_detect(WImage6,Y1,I0,waterMark1); case 0 break; otherwise error('you have input a valid value'); end d=input('please input your choice:'); end

happy · 发表于 2006-5-29 11:33

最好先把出现的问题说明一下

账号		自动登录	找回密码
密码			我要加入

快答辩了，求助基于离散余弦变换的数字水印matlab仿真源码

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（eagle217）快答辩了，求助基于离散余弦变换的...