www.gusucode.com > 小波变换的数字水印matlab源码程序 > new_wavelet/gongji.m

    addwaterimage1=imread('Add_Water_Image.bmp');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%加椒盐噪声
addwaterimage2=IMNOISE(addwaterimage1,'salt & pepper',0.1);
% addwaterimage3=IMNOISE(addwaterimage1,'salt & pepper',0.02);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%中值滤波
addwaterimage4=medfilt2(addwaterimage1,[3,3]);
% addwaterimage5=medfilt2(addwaterimage1,[1,5]);
addwaterimage6=IMNOISE(addwaterimage1,'gaussian',0,0.1);
% addwaterimage7=IMNOISE(addwaterimage1,'gaussian',0,0.005);
imwrite(uint8(addwaterimage2),'salt0.1.bmp');
% imwrite(uint8(addwaterimage3),'salt0.02.bmp');
imwrite(uint8(addwaterimage4),'midle33.bmp');
% imwrite(uint8(addwaterimage5),'midle15.bmp');
imwrite(uint8(addwaterimage6),'g0.1.bmp');
% imwrite(uint8(addwaterimage7),'g0.005.bmp');

IM=imread('barbara_512.bmp');
PSNR1=psnr(IM,addwaterimage2)
% PSNR2=psnr(IM,addwaterimage3)
PSNR3=psnr(IM,addwaterimage4)
% PSNR4=psnr(IM,addwaterimage5)
PSNR5=psnr(IM,addwaterimage6)
% PSNR6=psnr(IM,addwaterimage7)