小波实验报告dwt水印

一、题目：dwt水印
二、目的：编程实现水印的嵌入与还原
三、算法及其实现：水印的嵌入与还原
小波变换具有空间-频率的多尺度性，对图像的分解可以连续地从低分辨率过渡到高分辨率，也可以对图像整体进行变换，DWT还具有多分辨率特性，可以很好地与人类视觉系统(HVS)特性相匹配。

本算法将载体图像进行小波分解，利用公式实现水印信息的嵌入。

将2值图像中0和1
的值变成-1和1。

将水印入宿主图形重要的小波系数上。

水印的提取完全是水印嵌入的逆过程。

四、实现工具：Matlab
五、程序代码：
（1）dwt1_embed.m
%嵌入源码
clc
clear all;
% 保存开始时间
start_time=cputime;
k=20;
blocksize=8; % 设置块的大小
% 读入原图像
file_name='zimu.bmp';
cover_object=double(imread(file_name));
% 原图像的行数与列数
Mc=size(cover_object,1); %原图像行数
Nc=size(cover_object,2); %原图像列数
% 最大嵌入信息量
max_message=Mc*Nc/(blocksize^2);
% 读入水印图像
file_name='c.bmp';
message=double(imread(file_name));
%%水印图像的行数与列数
Mm=size(message,1); %水印行数
Nm=size(message,2); %水印列数
message_vector=reshape(message,1,Mm*Nm);
% 检查水印信息是否过大
if (Mm*Nm> max_message)
error('水印太大')
end
%将随机数发生器的状态置为1100
randn('state',1100);
% 产生watermarked_image,并写入原图信息
watermarked_image=cover_object;
% 将图像分块嵌入
% 当(2,2) > (2,3) 且message_pad(kk)=0
% 当(2,2) < (2,3) 且message_pad(kk)=1
%%经过分析可以得出结论：在提取水印时，如果cD1(2,2)>cD1(2,3)便是嵌入了水印的黑色像素，
%%反之则是嵌入了白色像素
x=1;
y=1;
h=waitbar(0,'嵌入水印，请等待');
for (kk = 1:length(message_vector))
% 对块进行DWT变换
[cA1,cH1,cV1,cD1] = dwt2(cover_object(x:x+blocksize-1,y:y+blocksize-1),'haar');
a=cD1;
% 如果message_pad(kk)== 0
if (message_vector(kk) == 0)
% 且(2,2) < (2,3) ，交换它们
if (cD1(2,2) < cD1(2,3))
temp=cD1(2,3);
cD1(2,3)=cD1(2,2);
cD1(2,2)=temp;
end
% 如果message_pad(kk) == 1,
elseif (message_vector(kk) == 1)
% 且(2,2) > (2,3) ，交换它们
if (cD1(2,2) >= cD1(2,3))
temp=cD1(2,3);
cD1(2,3)=cD1(2,2);
cD1(2,2)=temp;
end
end
% 检查(2,2) ，(2,3)的差是否>= k
if cD1(2,2) > cD1(2,3)
if cD1(2,2) - cD1(2,3) < k
cD1(2,2)=cD1(2,2)+(k/2);
cD1(2,3)=cD1(2,3)-(k/2);
end
else
if cD1(2,3) - cD1(2,2) < k
cD1(2,3)=cD1(2,3)+(k/2);
cD1(2,2)=cD1(2,2)-(k/2);
end
end
%IDWT
watermarked_image(x:x+blocksize-1,y:y+blocksize-1)= idwt2(cA1,cH1,cV1,cD1,'haar',[Mc,Nc]);
% 移动到下一块
if (x+blocksize) >= Nc
x=1;
y=y+blocksize;
else
x=x+blocksize;
end
waitbar(kk/length(message_vector),h);
end
close(h);
% 转换为uint8并写入dwt_watermarked.bmp watermarked_image_uint=uint8(watermarked_image);
imwrite(watermarked_image_uint,'dwt_watermarked.bmp','bmp');
% 显示运行时间
elapsed_time=cputime-start_time,
% 计算psnr值
psnr=psnr(cover_object,watermarked_image),
% 显示图像
figure(1)
imshow(message);
title('水印');
figure(2)
subplot(1,2,1);
imshow(watermarked_image,[])
title('嵌入水印图像')
subplot(1,2,2);
imshow('zimu.bmp');
title('原图');
（2）dwt1_recover.m
%DWT源码
%提取源码
clear all;
% 保存开始时间
start_time=cputime;
blocksize=8; % 设置块的大小
% 读入嵌入水印图像
file_name='dwt_watermarked.bmp';
watermarked_image=double(imread(file_name));
% 嵌入水印图像的行数与列数
Mw=size(watermarked_image,1); %嵌入水印图像的行数Nw=size(watermarked_image,2); %嵌入水印图像的列数
% 最大嵌入信息量
max_message=Mw*Nw/(blocksize^2);
% 读入原始水印
file_name='c.bmp';
orig_watermark=double(imread(file_name));
% 原始水印的行数与列数
Mo=size(orig_watermark,1); %原始水印的行数
No=size(orig_watermark,2); %原始水印的列数
% 将图像分块提取水印
x=1;
y=1;
h=waitbar(0,'水印提取中，请等待');
for (kk = 1:max_message)
% 对块进行dwt变换
[cA1,cH1,cV1,cD1] = dwt2(watermarked_image(x:x+blocksize-1,y:y+blocksize-1),'haar');
% 如果cD1(2,2) > cD1(3,3) 那么message_vector(kk)=0
% 否则message_vector(kk)=1
if cD1(2,2)>cD1(3,3)
message_vector(kk)=0;
else
message_vector(kk)=1;
end
% 移动到下一块
if (x+blocksize) >= Mw
x=1;
y=y+blocksize;
else
x=x+blocksize;
end
waitbar(kk/max_message,h);
end
close(h);
% 将message重新排列
message_vector=message_vector(1,1:Mo*No); %%只需前Mo*No个
message=reshape(message_vector,Mo,No);
% 显示运行时间
elapsed_time=cputime-start_time,
% 显示提取水印与原始水印
figure(3)
subplot(1,2,1);
imshow(message,[]);
title('提取水印')
subplot(1,2,2);
imshow(orig_watermark,[])
title('原始水印');
六、运行结果：
1、把C嵌入到图形中
把c从水印图形中提取出来并对比
七：结果分析：。