1039049-

基于LSB的数字水印技术

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

%添加高斯噪声
P1=imnoise(M,'gaussian',0.01);     
subplot(122) 
imshow(P1)                       
title('加高斯噪声的伪装图像');

%------对高斯躁声中值滤波-------
g=medfilt2(P1,[3,3]);                       
figure(2)
subplot(221)
imshow(g)
title('中值滤波后的高斯伪装图像');

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

%-------------------保存并读入伪装图像---------------
imwrite(uint8(g),'高斯中值滤波图像.bmp');
g=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(g(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')