| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 1260 人关注过本帖
标题:[求助]有哪位高手帮我看看这个程序
取消只看楼主 加入收藏
maidy2006
Rank: 1
等 级:新手上路
帖 子:3
专家分:0
注 册:2007-4-9
收藏
 问题点数:0 回复次数:1 
[求助]有哪位高手帮我看看这个程序
我搞了很久都没搞懂,程序如下面,特别是后半段,哪位高手帮我详细的解一下,谢谢了.clc
clear
%读入原图象
trueImage=imread('autumn.tif');
alfa=.1;
LENGTH=2500;
subplot(2,3,1);
imshow(trueImage);
title('原始图象');
%对原图象进行DCT变换
dctF1=rgb2gray(trueImage);
subplot(2,3,2);
imshow(dctF1);
title('灰度图象');
I=dct2(dctF1);
subplot(2,3,3)
imshow(log(abs(I)),[ ]);
title('DCT cofficient matrix');
[m,n]=size(I);
%产生水印序列并对其排序
z=1:LENGTH;
randn('state',10);
watermark1=randn(LENGTH,1);
subplot(2,3,4)
[Y0,I0]=sort(watermark1);
plot(I0,Y0,'k')
title('watermark seqence');
%找出水印嵌入位置(幅值较大的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,3,5);
imshow(idctF1,[ ]);
title('嵌入水印后的图象');
end
搜索更多相关主题的帖子: title 
2007-04-09 18:50
maidy2006
Rank: 1
等 级:新手上路
帖 子:3
专家分:0
注 册:2007-4-9
收藏
得分:0 

希望哪位高手帮帮忙!!!!!!!!!!!

2007-04-10 13:31
快速回复:[求助]有哪位高手帮我看看这个程序
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.032800 second(s), 8 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved