图像处理实验三 椒盐噪声

下载ftp上共享的图像文件完成如下的任务：

1、下载ftp上共享的图像文件Fig0318(a)(ckt-board-orig).tif，使用函数imnoise添

加椒盐噪声，噪声密度设为0.08 给出处理前后的图像结果对比

%%椒盐噪声

clc

clear

figure(1)

f= imread('I:\数字图像处理\IR\Fig0318(a)(ckt-board-orig).TIF');

subplot 211

imshow(f)

title('原始图像')

R=imnoise(f,'salt & pepper',0.08)

subplot 212

imshow(R);

title('噪声为0.08的椒盐噪声图像')

2、下载ftp上共享的图像文件Fig4.11(a).jpg，使用共享的自定义函数imnoise2（见

示例代码2），添加概率为0.15的胡椒噪声污染，再使用共享的自定义函数spfilt，使用反调和滤波器进行图像恢复，给出恢复前后的图像对比；类似的，当

Fig4.11(a).jpg被概率为0.15的盐粒噪声污染时，使用spfilt进行图像恢复，并

给出恢复前后的图像对比。

%%

clc

clear

figure(1)

f= imread('I:\数字图像处理\FD\I数字图像处理FD.JPG')

subplot 321

imshow(f);

title('原始图像')

PQ=size(f);

R=imnoise2( 'salt & pepper',PQ(1), PQ(2),0.15,0 ); %%胡椒噪声

c=find(R==0)

gp=f;

gp(c)=0;

subplot 323

imshow(gp);

title('胡椒图像')

R=imnoise2( 'salt & pepper',PQ(1), PQ(2),0,0.15 ); %%盐粒噪声

c=find(R==1)

gs=f;

gs(c)=255;

subplot 324

imshow(gs);

title('盐粒图像')

fp=spfilt(gp,'chmean',3,3,1.5)

subplot 325

imshow(fp);

title('胡椒图像恢复')

fs=spfilt(gs,'chmean',3,3,-1.5)

subplot 326

imshow(fs);

title('椒盐图像恢复')

将以上结果保存，写一个简单的实验报告，所有文件打包成压缩文件提交到ftp中的homework文件夹，命名规则：实验三+姓名+学号