利用椒盐噪声只有 2 个值的事实，中值滤波器在消除它方面做得很好。

• 使用中值过滤器过滤图像
• 查找中值滤波器完成大部分工作的像素
• 绘制这些像素的直方图

• 计算 0 和 1 的分数

  threshold = 0.2; % adjust to perform best on your images
  im = imread('../lena.bmp');
  imd = double(im);
  imd = imd./max(max(imd));
  imdg = imnoise(imd,'gaussian');
  imdsp = imnoise(imd,'salt & pepper');
  subplot(2,2,1), imshow(imdg)
  subplot(2,2,2), imshow(imdsp)
  imdgf = medfilt2(imdg);
  imdspf = medfilt2(imdsp);
  g = imdg.-imdgf;
  sp = imdsp.-imdspf;
  gm = abs(g)>threshold;
  spm = abs(sp)>threshold;
  subplot(2,2,3), hist(imdg(gm))
  subplot(2,2,4), hist(imdsp(spm))
  waitforbuttonpress