利用函数imnoise2处理噪声污染和spfilt处理滤波器

f=imread('0505.tif');
subplot(3,2,1),imshow(f)

[M,N]=size(f);%M和N分别显示图像f的行数和列数
 R=imnoise2('salt & pepper',M,N,0.1,0);%概率为0.1的胡椒污染图像
 c=find(R==0);%寻找R=0，并赋值为255
 gp=f;
 gp(c)=255;
 subplot(3,2,2),imshow(gp)
 
 R=imnoise2('salt & pepper',M,N,0,0.1);%盐粒污染图像
 c=find(R==1);%寻找R=1，并赋值为255
 gs=f;
 gs(c)=255;
 subplot(3,2,3),imshow(gs)
 
 fp=spfilt(gp,'amean',3,3);%使用算术平均滤波器得到的图像
 subplot(3,2,4),imshow(fp)
 
 fpmax=spfilt(gp,'max',3,3);%使用最大滤波器得到的图像
 subplot(3,2,5),imshow(fpmax)
 
 fpmin=spfilt(gs,'min',3,3);%使用最小滤波器得到的图像
 subplot(3,2,6),imshow(fpmin)