MATLAN图像处理之图像复原（自适应中值滤波）

%自适应空间滤波器
%自适应中值滤波器算法如下
%Sxy表示一个将被处理的，中心在（x,y）处的子图像
%Zmax 表示 Sxy中最大的亮度值
%Zmin 表示 Sxy中最小的亮度值
%Zmed 表示 Sxy中的亮度均值
%zxy  表示 坐标（x,y）处的亮度值

%该自适应中值滤波算法工作在两个层面 表示为  Level A  Level B
% Level A:若Zmin<Zmed<Zmax  则转向Level B
%         否则增加窗口尺寸
%         若窗口尺寸<=Smax,重复Level
%         否侧输出 Zmed
% Level B:若Zmin<Zxy<Zmax 则输出Zxy
%         否则输出 Zmed
% 其中  Smax表示允许的最大自适应滤波窗口的大小
% LevelA最后一步的另一种选择是输出Zxy来代替中值
% 这将产生一个稍微清楚的结果，但却可能探测不到与椒盐盐粒噪声值相同的
% 内涵与常熟背景中的椒盐噪声
clc
clear
f=imread('D:\亮度调节与空间滤波\5.jpg');
f=f(:,:,1);
[M,N]=size(f);
R=imnoise2('salt & pepper',M,N,0.1,0);%概率为0.1的椒盐噪声
c=find(R==0);
gp=f;
gp(c)=0;
R1=imnoise2('salt & pepper',M,N,0,0.1);%盐粒噪声
c=find(R1==1);
gs=f;
gs(c)=255;

f1=adpmedian(gp,7);
f2=adpmedian(gs,7);

subplot(3,2,1)
imshow(f)
subplot(3,2,3)
imshow(gp)
subplot(3,2,4)
imshow(gs)
subplot(3,2,5)
imshow(f1)
subplot(3,2,6)
imshow(f2)