高斯噪声和椒盐噪声

clear all;clc;close all;
im=imread('1.jpg');
% o=im;
[w, h ,z]=size(im);
im=rgb2gray(im);
o=im;%原图备份
figure;seq=0;
seq=seq+1;subplot(2,2,seq);imshow(im);title('原图');pos=axis;xlabel('','position',[(pos(1)+pos(2))/2 pos(4)]);

%高斯噪声
gau=normrnd(0,0.1,w,h);
im1=im2double(im);%转成double
im=im1+gau;
seq=seq+1;subplot(2,2,seq);imshow(im);title('高斯噪声');pos=axis;xlabel('','position',[(pos(1)+pos(2))/2 pos(4)]);

%椒盐噪声
im=o;
k1=0.1;
k2=0.3;
a1=rand(w,h)<k1;
a2=rand(w,h)<k2;
im(a1&a2)=0;%通过同维矩阵可以选取元素!
im(a1&~a2)=255;
seq=seq+1;subplot(2,2,seq);imshow(im);title('椒盐噪声');pos=axis;xlabel('','position',[(pos(1)+pos(2))/2 pos(4)]);