数字图像处理（自适应中值滤波）

自适应中值滤波

常规的中值滤波在噪声密度不是很大的情况下（噪声出现的概率小于0.2时），性能很好，但是当概率较大时，常规的中值滤波处理完之后仍然存在噪声并且丢失了细节。自适应中值滤波器在进行滤波处理时，会根据某些条件来改变S的尺寸。滤波器输出的是一个数值，该数值用来表示（x,y）的像素值。

算法原理

算法的原理如书中p209中所示：


其中要使用到一个重要函数ordfilt2():

Matlab代码

clc;
clear;
f=imread('Kobe.jpg');
image_gray=rgb2gray(f);%得到灰度图像
ff =image_gray;
alreadyProcessed = false(size(image_gray));%生成逻辑非的矩阵
% 迭代.
Smax=7;
for k = 3:2:Smax
   zmin = ordfilt2(image_gray, 1, ones(k, k), 'symmetric');
   zmax = ordfilt2(image_gray, k * k, ones(k, k), 'symmetric');
   zmed = medfilt2(image_gray, [k k], 'symmetric');
   
   processUsingLevelB = (zmed > zmin) & (zmax > zmed) & ...
       ~alreadyProcessed; %需要转到B步骤的像素
   zB = (image_gray > zmin) & (zmax > image_gray);
   outputZxy  = processUsingLevelB & zB;%满足步骤AB的输出原值对应像素的位置
   outputZmed = processUsingLevelB & ~zB;%满足A不满足B的输出中值对应的像素位置
   ff(outputZxy) = image_gray(outputZxy);
   ff(outputZmed) = zmed(outputZmed);
   
   alreadyProcessed = alreadyProcessed | processUsingLevelB;
   if all(alreadyProcessed(:))
      break;
   end
end
ff(~alreadyProcessed) = zmed(~alreadyProcessed);
f1=imnoise(image_gray,'salt & pepper',0.3);%添加椒盐噪声后的图像
f2=medfilt2(f1,[3,3]);%中值滤波后的图像
subplot(2,2,1);
imshow(image_gray);
title('原图');
subplot(2,2,2);
imshow(f1);
title('椒盐噪声污染后的图像');
subplot(2,2,3);
imshow(f2);
title('中值滤波');
subplot(2,2,4);
imshow(ff);
title('自适应中值滤波');

实验效果展示