Question

谁能帮我解释这段MATLAB代码、

clear all; close all; clc;
M=512;%原图像长度
N=64; %水印长度
[filename1,pathname]=uigetfile('*.*','select the image');
image1=imread(num2str(filename1));
subplot(2,2,1);imshow(image1); title('original image'); % orginal image for watermarking
image1=double(image1);%，matlab将图象中的数据存储为double精度型
imagew=imread('水印.jpg');
subplot(2,2,2);imshow(imagew);title('original watermark'); %original watermark
%嵌入水印
[ca,ch,cv,cd]=dwt2(image1,'db1');%使用指定的小波基函数 'db1' 对二维信号 image1 进行二维
%离散小波变幻,ca，ch,cv,cd 分别为近似分量、水平细节分量、垂直细节分量和对角细节分量
[cas,chs,cvs,cds]=dwt2(ca,'db1');
for i=1:N
for j=1:N
if imagew(i,j)==0
a=-1;
else
a=1;
end
cas(i,j)=cas(i,j)*(1+a*0.03);
end
end
IM=idwt2(cas,chs,cvs,cds,'db1') ;%二维离散小波反变换，重构图像
markedimage=double(idwt2(IM,ch,cv,cd,'db1'));
%显示嵌入后水印图像
subplot(2,2,3);colormap(gray(256));image(markedimage);title('marked image');
imwrite(markedimage,gray(256),'watermarked.bmp','bmp');
%提取水印
image1=imread(num2str(filename1));image1=double(image1);
imaged=imread('watermarked.bmp');
[ca,ch,cv,cd] = dwt2(image1,'db1');
[cas,chs,cvs,cds]=dwt2(ca,'db1');
[caa,chh,cvv,cdd]=dwt2(imaged,'db1');
[caas,chhs,cvvs,cdds]=dwt2(caa,'db1');
for p=1:N
for q=1:N
a=caas(p,q)/cas(p,q)-1;
if a<0
W(p,q)=0;
else
W(p,q)=255;
end
end
end
%显示提取的水印
subplot(2,2,4);
colormap(gray(256));image(W);title('从含水印图像中提取的水印');
imwrite(W,gray(256),'watermark.bmp','bmp');

若以下回答无法解决问题，邀请你更新回答

Answer 1

一个简单的提取图像中水印的程序，你要解释什么？

追问

这段每句是什么意思

追答

意思是对于imagew上面每个像素元，如果该像素元值为0，则把cas矩阵对应位置的数乘以0.97，反之则乘以1.03

Answer 2

你在逗我们么。。。哪儿有代码啊

追问

你在逗我么、在补充问题里有代码.......