噪声模型

常见的噪声模型有(z是噪声值，m表示均值，s2表示方差)：

1、高斯噪声:m=a , s2=b^2

2、均匀分布噪声

3、伽马噪声 m=b/a , s2=b/a^2

4、指数分布噪声 m=1/a , s2=1/a^2

5、瑞利噪声 :

6、椒盐噪声
其中盐表示亮点，椒表示暗点。

几种常见的噪声的概率密度函数的matlab代码：

x=-5:0.2:5;
subplot(3,2,1);
y1=normpdf(x,0,1);
plot(x,y1);
title('高斯噪声图像'); 

subplot(3,2,2);   
y2=unifpdf(x,-2,4);
plot(x,y2);
title('均匀噪声图像');
    
subplot(3,2,3);
y3=zeros(size(x));
y3(1)=0.7;
y3(end)=0.7;
plot(x,y3);
title('椒盐噪声图像'); 

subplot(3,2,4);        
y4=raylpdf(x,2);       
plot(x,y4);
title('瑞利噪声图像');        
        
subplot(3,2,5);       
y5=exppdf(x,1);  
plot(x,y5);
title('指数噪声图像');          
        
subplot(3,2,6);  
y6=gampdf(x,2,4);       
plot(x,y6);
title('伽马噪声图像')

如何给一幅图像添加噪声？
例：将图像添加这六种常见的噪声并观察前后的效果。

I=imread('dragon.jpg');%读取图像
x=rgb2gray(I);%转换为灰度图像
subplot(2,4,1);%分区绘图
imshow(x);%显示图像
title('原图');%给图像加标题
%添加高斯噪声
J=imnoise(x,'gaussian',0.02);%在图像中加入高斯噪声
subplot(2,4,2);
imshow(J,[]);
title('高斯加噪后');
%添加椒盐噪声
K=imnoise(x,'salt & pepper',0.02);%在图像中加入椒盐噪声
subplot(2,4,3);
imshow(K,[]);
title('椒盐加噪后');
[M,N]=size(x)
%添加均匀分布噪声
R=unifrnd(-10,20,M,N);%得到均匀分布的随机数
H=double(x)+R;%添加噪声
H=uint8(round(R));
subplot(2,4,4);
imshow(H,[])
title('均匀分布加噪后');
%添加指数噪声
R1=exprnd(20,M,N);%得到指数分布的随机数
O=double(x)+R1;%添加指数噪声
O=uint8(round(R1));
subplot(2,4,5);
imshow(O,[])
title('指数加噪后');
%添加瑞利噪声
R2=raylrnd(20,M,N);%得到瑞利分布的随机数
P=double(x)+R2;%添加瑞利噪声
P=uint8(round(R2));
subplot(2,4,6);
imshow(P,[])
title('瑞利加噪后');
%添加伽马噪声
R3=gamrnd(-10,20,M,N);%得到伽马分布的随机数
Q=double(x)+R3;%添加伽马噪声
Q=uint8(round(R3));
subplot(2,4,7);
imshow(Q,[])
title('伽马加噪后');

例：在灰度图像中加入均值为0，方差为0.02的高斯噪声，通过100次相加求平均的方法去除所加的噪声，比较显示原图、加噪图与去噪图。

I=imread('dragon.jpg');%读取图像
x=rgb2gray(I);%转换为灰度图像
subplot(1,3,1);%分区绘图
imshow(x);%显示图像
title('原图');%给图像加标题
J=imnoise(x,'gaussian',0,0.02);%在图像中加入高斯噪声
subplot(1,3,2);
imshow(J);
title('加噪后');
%100次相加求平均去噪声
H1=zeros(size(x));
for i=1:100
    J=imnoise(x,'gaussian',0,0.02);
    H1=H1+double(J);
end
H=H1/100;
subplot(1,3,3);
imshow(uint8(H));
title('去噪后');