数字图像处理实验之中值滤波器

代码如下：
主函数：

% 输入图像
img_1 = imread('Fig0335(a)(ckt_board_saltpep_prob_pt05).tif');
[r, c] = size(img_1);
% 滤波器大小
m = 3;

% 作用参考我的上一篇博客
img_1s = m_10_Smooth_create_img(img_1, r, c, m - 1);
img_2 = m_10_Smooth_process_img(img_1s, r, c, m);

% 中值滤波处理
img_3 = m_12_Median_Filter_process_img(img_1, r, c, m);

subplot(131), imshow(img_1);
subplot(132), imshow(img_2);
subplot(133), imshow(img_3);

m_12_Median_Filter_process_img.m 代码：

function img_2 = m_12_Median_Filter_process_img(img_1, r, c, m)

% 创建一个填充边图像
img_1s = m_10_Smooth_create_img(img_1, r, c, m - 1);

% 创建新图像
img_2 = zeros(r, c);
% 逐个像素处理
for x = 1 : r
    for y = 1 : c
        % 将滤波器中的所有像素保存起来，用于求中位数
        arr = zeros(1, m*m);
        for i = 1:m
            for j = 1:m
                arr((i - 1) * m + j) = img_1s(x + i - 1, y + j - 1);
            end
        end
        % 求中位数
        img_2(x, y) = median(arr);
    end
end

img_2 = im2uint8(mat2gray(img_2));
end

实验结果如下：