dete.m

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

date=imread('E:\watermark\use picture\256.bmp','bmp');
% a1=date(:,:,1);
% a2=date(:,:,2);
% a3=date(:,:,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);
watermark2=imread('E:\watermark\use picture\sign32','bmp');
waterqa=im2bw(watermark2,0.4);
% dat=imnoise(date,'salt & pepper',0.02);   % add noise
% data(130:160,100:180)=0;           % cut
% % data=midfilt2
% data2=imrotate(date,-1,'bilinear','crop');
% dat=imrotate(data2,1,'bilinear','crop');
% % % % data3=ones(256,256)*255;
% % % % imshow(data3)
% % % for q=0.5:0.1:1.5
% %    data3=imrotate(data2,1,'bilinear','crop');
% %    dat=ones(256,256)*255;
%     dat(5:244,5:244)=dat(15:254,15:254);
 figure(4)
 imshow(date)
% dat=ones(256,256)*255;

%  data=date;

% dat=ones(256,256)*255;
% dat(1:226,1:216)=data3;
[a,h,v,d]=lwt2(double(date),'bior4.4');
[m,n]=size(d);
[block,av]=fenkuai(d);
t0=min(min(av));
t1=max(max(av));
water2=ones(32,32);
for i=1:2:64
    for j=1:2:64
        if abs(av(i,j)-av(i+1,j)-t0)<=abs(av(i,j)-av(i+1,j)-t1)
            water2(fix(i/2)+1,fix(j/2)+1)=0;
        else
            water2(fix(i/2)+1,fix(j/2)+1)=1;
        end
    end
end
figure(3);
showtu(date,date,water,water2);
% k1(1,q*10)=err(water,water2);
% k2(1,q*10)=err(waterqa,water2);
%   end