自适应中值滤波器matlab实现

将下面代码直接贴入matlab中，并将读入图像修改成自己机子上的，就可以运行了。可以按照“%%”顺序分步来运行

%% function 自适应中值滤波器

%%%%%%%%%%%%%%%

%实现两个功能：

%1.对高密度的椒盐噪声有好的滤除效果；

%2.滤波时减少对图像的模糊；

%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%

%原理：

%1.椒盐噪声概率越大，滤波器窗口需越大。故若滤波器窗口随噪声概率自适应变化，才能有好的滤除效果

%2.为减少对图像的模糊，需在得出原图像值并非椒盐噪声点时，保留原图像值不变；

%3.椒盐噪声点的特点：该点的值为该点领域上的最大或最小；%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%

%步骤（得到图像中某点(x,y)（即窗口中心点）的值的步骤）：

%1.设定一个起始窗口，以及窗口的最大尺寸；

%2.（此步用于确定窗口大小）对窗口内像素排序，判断中值是否是噪声点，若不是，继续第3步，若是，转到第5步；

%3.判断中心点是否是噪声点，若不是，则输出该点的值（即图像中该点的原值不变）；若是，则输出中值；

%4.窗口尺寸增大，若新窗口尺寸小于设定好的最大值，重复第2步，若大于，则滤波器输出前一个窗口的中值；

%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%

%参数说明：

%被噪声污染的图像（即退化图像也即待处理图像）:Inoise

%滤波器输出图像：Imf

%起始窗口尺寸:nmin*nmin(只取奇数)，窗口尺寸最大值：nmax*nmax %图像大小:Im*In

%窗口内图像的最大值Smax,中值Smed,最小值Smin

%%%%%%%%%%%%%%%%%%%%

%%

clear

clf

%% 读入图像I

I=imread('e:/photo/cat.jpg');

%转化为灰度图Ig

Ig=rgb2gray(I);

%被密度为0.2的椒盐噪声污染的图像Inoise

Inoise=imnoise(Ig,'salt & pepper',0.2);

%或者是被方差为0.2的高斯噪声污染的图像Inoise

%Inoise=imnoise(Ig,'gaussian',0.2);

%显示原图的灰度图Ig和噪声图像Inoise

subplot(2,2,1),imshow(Ig);xlabel('a.原始灰度图像');

subplot(2,2,2),imshow(Inoise);xlabel('b.被噪声污染的图像');

%% 定义参数

%获取图像尺寸:Im,In

[Im,In]=size(Inoise);

%起始窗口尺寸:nmin*nmin(窗口尺寸始终取奇数)

nmin=3;

%最大窗口尺寸:nmax*nmax

nmax=9;

%定义复原后的图像Imf

Imf=Inoise;

%为了处理到图像的边界点，需将图像扩充

%因为窗口尺寸是弹性的，所以将Inoise固定扩充到最大:I_ex[(Im+(nmax-1))*(In+(nmax-1))]

I_ex=[zeros((nmax-1)/2,In+(nmax-1));zeros(Im,(nmax-

1)/2),Inoise,zeros(Im,(nmax-1)/2);zeros((nmax-1)/2,In+(nmax-1))];

%% 自适应滤波过程

%遍历图像Inoise中的每一点

for x=1:Im

for y=1:In

for n=nmin:2:nmax

%图像Inoise中的某点(x,y)的领域Sxy，对应在I_ex中为(x+[(nmax-1)/2-(n-1)/2]:x+[(nmax-1)/2-(n-1)/2]+(n-1),y+(nmax-1)/2-(n-

1)/2:y+[(nmax-1)/2-(n-1)/2]+(n-1))

Sxy=I_ex(x+(nmax-1)/2-(n-1)/2:x+(nmax-1)/2+(n-1)/2,y+(nmax-1)/2-(n-1)/2:y+(nmax-1)/2+(n-1)/2);

Smax=max(max(Sxy));%求出窗口内像素的最大值

Smin=min(min(Sxy));%求出窗口内像素的最小值

Smed=median(median(Sxy));%求出窗口内像素的中值

%判断中值是否是噪声点

if Smed>Smin && Smed

%若中值既大于最小值又小于最大值，则不是

%是，则退出该if语句，增大窗口尺寸，再次判断

%不是，则判断该点的原值是不是噪声点

if Imf(x,y)<=Smin || Imf(x,y)>=Smax

%若该点的原值既大于最小值又小于最大值，则不是

%不是，则输出原值，即不作处理

%是，则输出中值

Imf(x,y)=Smed;

end

break

%有输出则不再进行循环判断

end

end

%当n=max时，输出中值

Imf(x,y)=Smed;

end

end

subplot(2,2,3),imshow(Imf);xlabel('d.最大尺寸为9的自适应中值滤波器的滤波效果');

%% 与普通中值滤波器的对比

Imf1=medfilt2(Inoise,[3,3]);

Imf2=medfilt2(Imf1,[3,3]);

subplot(2,2,4),imshow(Imf2);xlabel('e.尺寸为3的普通中值滤波器两次滤波效果'); 下面是运行结果：