[求助]大家帮我看看这段程序错在哪？

zunny · 发表于 2006-4-19 13:58

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

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

x

程序如下： %读入二值水印图像，变换化为一个一维序列 Water=imread('D:\My Essay\草稿\H','bmp'); Encry=reshape(Water,1,1024);
 %读入载体图像Lena Lena=imread('D:\My Essay\草稿\Lena','bmp'); I=im2double(Lena);
%水印的嵌入，DCT变换 k=0; for i=1:32 for j=1:32 x=(i-1)*8+1; y=(j-1)*8+1; BLOCK=I(x:x+8-1,y:y+8-1); BLOCK_dct=dct2(BLOCK); %图像8、8分块然后DCT BLOCK_zig=zigzag(BLOCK_dct); %这个zigzag函数为自己编写，使之之字形排列，经验证无误的。同样下面的izigzag是把之字形后的恢复，也无误。 BLOCK_wtrz=BLOCK_zig; k=k+1; if Encry(k)==0 e=-1; else e=1; end BLOCK_wtrz(2)=BLOCK_zig(2)*(1+0.01*e); BLOCK_wtr=izigzag(BLOCK_wtrz); BLOCK_wtrd=idct2(BLOCK_wtr); Iw(x:x+8-1,y:y+8-1)=BLOCK_wtrd; end end
 %水印析出 count=0; for i=1:32 for j=1:32 count=count+1; x=(i-1)*8+1; y=(j-1)*8+1; iBLOCK=I(x:x+8-1,y:y+8-1); iBLOCKw=Iw(x:x+8-1,y:y+8-1); iBLOCK_dct=dct2(iBLOCK); iBLOCKw_dct=dct2(iBLOCKw); iBLOCK_zig=zigzag(iBLOCK_dct); iBLOCKw_zig=zigzag(iBLOCKw_dct); if iBLOCKw_zig(2)<iBLOCK_zig(2) Decry(count)=0; else Decry(count)=1; end end end 理论上，应该是Encry（水印数据）和Decry（从图中析出后的水印数据）相等的，可实际执行之后，二者变化很大。 有一些现象，如Encry是double array；而Decry却为logical array. 请问，怎么改动才可使两者一样呢。

zunny · 发表于 2006-4-19 16:38

怎么没有人回答啊，随便说点看法也好啊，起码会有点灵感出来

zunny · 发表于 2006-4-21 11:39

我再把它顶起来，看看是否有人可以给出解决方案阿

tracy_dan · 发表于 2006-4-21 13:01

太难了

happy · 发表于 2006-4-21 22:11

能否把你的图片传上来

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

[求助]大家帮我看看这段程序错在哪？

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

回复：（zunny）[求助]大家帮我看看这段程序错在哪？