中值滤波器处理椒盐噪声

clc;
close all;
clear all;
%读取原始图像
im = imread(‘Fig0318(b)(ckt-board-slt-pep-both-0pt2).tif’);
im_noise_salt = imnoise(im,‘salt & pepper’,0.05); % 加入椒盐噪声
% im_noise_gaussian = imnoise(im,‘gaussian’);
%定义邻域尺寸
n1 = 2; m1 = 2n1+1;
n2 = 2; m2 = 2n2+1;
k = floor(m1m2/2)+1;
[h,l,c] = size(im_noise_salt); % [h,l,c] = size(im_noise_gaussian);
H = zeros(h,l);
t = zeros(n1,n2);
for i=n1+1:h-n1
for j=n2+1:l-n2
Neiborhood = im_noise_salt(i-n1:i+n1,j-n2:j+n2); 得到33邻域
t = Neiborhood(😃;
s = sort(t);
H(i,j) = s(k);
end
end
im_filtered2 = uint8(H);
subplot(2,2,1),imshow(im);title(‘原图’);
subplot(2，2，2),imshow(im_noise_salt);title(‘加入椒盐噪声污染0.1后的图像’);
subplot(2，2，3),imshow(im_filtered2);title(‘椒盐噪声图像进行中值滤波后的图像’);