% 自适应中值滤波
clc,clear,close all % 清理命令区、清理工作区、关闭显示图形
warning off % 消除警告
feature jit off % 加速代码运行
[filename ,pathname]=…
uigetfile({‘.bmp’;'.jpg’;},‘选择图片’); % 选择图片路径
str = [pathname filename]; % 合成路径+文件名
im = imread(str); % 原图像
im = imnoise(im,‘gaussian’,0,1e-3); % 原图像 + 白噪声

im1 = adapmedian_filter( im,3 ); % 应用自适应中值滤波
figure(‘color’,[1,1,1])
subplot(121),imshow(im,[]),title(‘original image’)
colormap(jet) % 颜色
shading interp % 消隐
subplot(122),imshow(im1,[]),title(‘自适应中值滤波’)
colormap(jet) % 颜色
shading interp % 消隐
function Z = adapmedian_filter(X, Smax)
% 函数对输入图像进行邻域窗口大小可变的自适应中值滤波
% 函数输入
% X：输人二维图像矩阵
% smax：中值滤波邻域窗日的最大值，必须是大于1的奇数
% 函数输出
% Z：输出图像矩阵，数据类型与输人相同
% smax必须是大于1的奇数
if (Smax<=1) | (Smax/2 == round(Smax/2)) | (Smax ~= round(Smax))
error(‘smax必须是大于1的奇数’)
end
if size(X,3)~=1
error(‘图像应该为2维矩阵’)
end
if ~isa(X,‘double’)
X = double(X)/255; % 数据类型
end
[M, N] = size(X);
% 初始化
Z = X; Z(:)=0;
LevelA = false(size(X)); % 初始化，同X的 全0逻辑矩阵
% 开始滤波
for k=3:2:Smax
Zmin = ordfilt2(X,1,ones(k,k),‘symmetric’ ); % 排序滤波
Zmax = ordfilt2(X,k*k,ones(k,k), ‘symmetric’); % 排序滤波
Zmed = medfilt2(X,[k,k], ‘symmetric’); % 中值滤波
% 判断是否进入B层
LevelB = (Zmed>Zmin)&(Zmax>Zmed)&LevelA; % 判断A层
ZB = (X > Zmin) & (Zmax > X);
outputZxy = LevelB & ZB; % 交运算
outputZmed = LevelB & ~ZB; % 交运算
Z(outputZxy) = X(outputZxy); % 赋值
Z(outputZmed) = Zmed(outputZmed); % 赋值
LevelA = LevelA | LevelB; % 非运算
if all(LevelA(😃)
break;
end
end
Z(~LevelA) = Zmed(~LevelA); % 赋值
Z = im2uint8(Z); % 类型转化
end