MATLAB排序滤波

% 排序滤波
clc,clear,close all % 清理命令区、清理工作区、关闭显示图形
warning off % 消除警告
feature jit off % 加速代码运行
im = imread(‘brain.bmp’); % 原图像
im = imnoise(im,‘gaussian’,0,1e-3); % 原图像 + 白噪声

im1 = ordfilt2(im, 1, true(5)); % 应用排序滤波
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 A = ord_filt2(im,order,domain)

if ~isa(im,‘double’) % 是否为double类型
A = double(im)/255;
end

domainSize = size(domain); % 维数
center = floor((domainSize + 1) / 2);
[r,c] = find(domain);
r = r - center(1); % 以中心分开
c = c - center(2); % 以中心分开
padSize = max(max(abs®), max(abs©)); % 求距离中心点最大长度
A = padarray(A, padSize * [1 1], 0, ‘both’); % 周围全部用0填充，向左右增2列0，上下增两行0
% A =
% 1 3 4
% 2 3 4
% 3 4 5
% B = padarray(A, 2 * [1 1], 0, ‘both’)
% 0 0 0 0 0 0 0
% 0 0 0 0 0 0 0
% 0 0 1 3 4 0 0
% 0 0 2 3 4 0 0
% 0 0 3 4 5 0 0
% 0 0 0 0 0 0 0
% 0 0 0 0 0 0 0
Ma = size(A,1); % 行
offsets = c*Ma + r;
% 确保offsets有效
if ~isreal(offsets) || any(floor(offsets) ~= offsets) || any(~isfinite(offsets))
error(message(‘offsets无效’))
end

% 排序
B = ord_filt2(A, order, offsets);
end