數字圖像處理實驗之中值濾波器

代碼如下：

主函數：

% 輸入圖像
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
           

實驗結果如下：