1039049-

基于LSB的数字水印技术

%-------------------- 添加椒盐噪声----------------
%读入伪装图像
M=imread('伪装图像.bmp');           
figure(1)
subplot(121) 
imshow(M)                        
title('第1位平面伪装图像')

P2=imnoise(M,'salt & pepper',0.03); 
subplot(122)
imshow(P2)                          
title('添加椒盐噪声的伪装图像');

%---------椒盐算术均值滤波--------
l=[1 1 1                            
   1 1 1 
   1 1 1]; 
l=l/9;                           
k=conv2(P2,l);                        
figure(2)
subplot(221)
imshow(k,[])
title('算数均值滤波后的椒盐伪装图像')

%imwrite(uint8(k),'算数均值滤波图像.bmp')
%Img_cover = double(imread('算数均值滤波图像.bmp'));
%[nrow,ncol] = size(Img_cover);

%-------------------保存并读入伪装图像---------------
imwrite(uint8(k),'椒盐算数均值滤波图像.bmp');
k=imread('椒盐算数均值滤波图像.bmp');
I_water=imread('sy.bmp');%读入水印图像

%-----------------------提取算法---------------------
M=8;
water_detect=zeros(nrow,ncol);
for i=1:nrow
    for j=1:ncol
    Img_bin=dec2bin(k(i,j),8);%把一个十进制数转换成二进制数
        if Img_bin(N)=='0'
           water_detect(i,j)=0;
        else
            water_detect(i,j)=1;   
        end                        
    end
end
%结果显示
subplot(222);
imshow(I_water);%显示原水印图像
title('原水印图像256*256')
subplot(223);
imshow(water_detect);
title('均值提取出的水印图像256*256')