Weighted Median Filter

clear all; close all; clc;
a=imread('rice.png');
a=imnoise(a,'salt & pepper',0.2);
a=im2double(a);
[m n]=size(a);
si=2;
w=[1 1 1 1 1; 1 2 2 2 1;1 2 3 2 1 ;1 2 2 2 1;1 1 1 1 1];
for i=1:m
for j=1:n
con=0;
for k= -si:si
for p=-si:si
if (i+k>0 && j+p>0) &&(i+k < m && j+p < n)
for wi=1:w(3+k,3+p)
con=con+1;
s(con)=a(i+k,j+p);
end
end
end
end
b(i,j)=median(s);
end
end
figure,imshow(a),title('Original Image');
figure,imshow(b),title('Image After Weightedmedian');