matlab数字水印 ！！！大神求助！

Question

clear;
close all;
M=imread('yin.jpg'); %读取水印图像
[rm,cm]=size(M);
I=imread('meizi.jpg'); %将图像矩阵转为double型
I=double(I)/255; %对矩阵计算处理要将其转化为double型，用unit8会溢出
[r,c]=size(I);
alpha=0.2;
k1=randn(1,8); %randn是均值为0方差为1的正态分布
k2=randn(1,8);
T=dctmtx(8); %产生二维DCT变换矩阵
%对图像分块
cda0=blkproc(I,[8 8],'P1*x*P2',T,T');
cda1=cda0;
for i=1:rm
for j=1:cm
x=(i-1)*8;y=(j-1)*8;
if M(i,j)==1
k=k1;
else
k=k2;
end
cda1(x+1,y+8)=cda0(x+1,y+8)+alpha*k(1);
cda1(x+2,y+7)=cda0(x+2,y+7)+alpha*k(2);
cda1(x+3,y+6)=cda0(x+3,y+6)+alpha*k(3);
cda1(x+4,y+5)=cda0(x+4,y+5)+alpha*k(4);
cda1(x+5,y+4)=cda0(x+5,y+4)+alpha*k(5);
cda1(x+6,y+3)=cda0(x+6,y+3)+alpha*k(6);
cda1(x+7,y+2)=cda0(x+7,y+2)+alpha*k(7);
cda1(x+8,y+1)=cda0(x+8,y+1)+alpha*k(8);
end
end
%嵌入水印后的图像
im_cover=blkproc(cda1,[8 8],'P1*x*P2',T',T);
%提取水印图像
dca1=blkproc(im_cover,[8 8],'P1*x*P2',T,T');
a=zeros(1,8);
for i=1:rm
for j=1:cm
x=(i-1)*8;y=(j-1)*8;
a(1)=dca1(x+1,y+8);
a(2)=dca1(x+2,y+7);
a(3)=dca1(x+3,y+6);
a(4)=dca1(x+4,y+5);
a(5)=dca1(x+5,y+4);
a(6)=dca1(x+6,y+3);
a(7)=dca1(x+7,y+2);
a(8)=dca1(x+8,y+1);
if corr2(a,k1)>corr2(a,k2)
mark(i,j)=1;
else
mark(i,j)=0;
end
end
end
figure;
subplot(221);imshow(I);
title('载体图像');
subplot(222);imshow(M);
title('原水印图像');
subplot(223);imshow(im_cover);
title('嵌入水印后图像');
subplot(224);imshow(mark);
title('提取的水印图像');

错误显示：
??? Subscripted assignment dimension mismatch.
Error in ==> blkproc at 89
aa(border(1)+(1:ma),border(2)+(1:na)) = a;
Error in ==> Untitled0 at 12
cda0=blkproc(I,[8 8],'P1*x*P2',T,T');

Answer 1

您好，当前的错误是由于您使用的图像为彩色，利用
I=rgb2gray(I);
转为灰度图像即可解决。
PS: 运行了您的程序，改灰度以后并不能出正确结果，因为后面的语句还有其他问题，望您检查一下程序算法。

追问

您好，谢谢回复。
程序我漏打了转成灰度的图像，实际上我有写，是的，错误还是如上的显示。
您能具体指出后面的错误吗？

追答

如果您使用了灰度图像，上面的错误是不会显示的。
取而代之是下面的错误：
??? Attempted to access cda0(1,344); index out of bounds because size(cda0)=[223,337].
Error in ==> Digital_Watermark at 29
cda1(x+1,y+8)=cda0(x+1,y+8)+alpha*k(1);

这个错误是说您的第一个循环会超出矩阵cda1 cda0的范围

追问

不是很明白为什么会出现这样的错误？
能给出修改的办法吗？

Answer 2

蒋门神是孟州的一个恶霸，强占了金眼彪施恩开的快活林酒店，施恩武艺不精，于是请武松帮忙，武松一路上一共喝了三四十碗酒，醉打蒋门神，帮施恩重新夺回快活林酒店。