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)