embed.m

数字水印嵌入和提取算法 包含时空域和变换域的算法的matlab实现

close all;
clear all;

data=double(imread('E:\watermark\use picture\lena512.bmp','bmp'));%double(imread('cameraman.tif'));
% a=data(51:306,101:356,:);
% a1=a(:,:,1);
% a2=a(:,:,2);
% a3=a(:,:,3);
% k1=0.229.*a1+0.587.*a2+0.114.*a3;
% k2=0.596.*a1-0.275.*a2-0.321.*a3;
% k3=0.212.*a1-0.523.*a2+0.311.*a3;

watermark=imread('E:\watermark\use picture\picture\1234.bmp','bmp');
water=im2bw(watermark,0.4);
% water=zeros(32,32);
% for i=1:2
%     for j=1:32
%         water((i-1)*16+1:(i-1)*16+16,j)=[0,0,1,1,0,1,0,0,0,1,1,0,1,0,0,0];
%     end
% end

[a,h,v,d]=lwt2(data,'bior4.4');                   %9.7 lifting filter  dwt2(data,'db1');  %
figure(1);
showtu(a,h,v,d);

[m,n]=size(d);
[block,av]=fenkuai(d);
block1=block;
t0=min(min(av));
t1=max(max(av));

block1=embeded(block,av,t0,t1,water);
d1=makeup(block1);
k1=ilwt2(a,h,v,d1,'bior4.4');
% a1=1.0753.*k1+1.0276.*k2+0.6665.*k3;
% a2=1.0753.*k1-0.1996.*k2-0.6002.*k3;
% a3=1.0753.*k1-1.0362.*k2+1.7511.*k3;
% data22=zeros(256,256,3);
% data22(:,:,1)=a1;
% data22(:,:,2)=a2;
% data22(:,:,3)=a3;

figure(2);
showtu(data,uint8(k1),water,0);
imwrite(uint8(k1),'E:\watermark\use picture\256.bmp','bmp');