Program for removing salt and pepper noise using median filter

Program for removing salt and pepper noise using median filter

Program:

clc;

a=imread('D:\Flower.jpg');   

b=rgb2gray(a);

ab = imnoise(b,'salt & pepper',0.02);   %adding noise

d=double(ab);

e=d;

[row,col]=size(d);

     for x=2:1:row-1

        for y=2:1:col-1

            a1=[a(x-1,y-1) a(x-1,y) a(x-1,y+1)...

                a(x,y-1) a(x,y) a(x,y+1)...

                a(x+1,y-1) a(x+1,y) a(x+1,y+1)];

            a2=sort(a1);

            med= a2(5);

            f(x,y)=med;

     end

end

figure(1);

imshow(ab);

title('salt and pepper noise ');

figure(2);

imshow(f);

title('salt and pepper noise removed ');

Output: