a = imread('PSWebImageBitMap0001.bmp',  'bmp');
b = a(:,:,2);

%J = imnoise(I,'salt & pepper',0.02);
c = medfilt2(b);

figure(3)
clf
imagesc(c)
colormap(gray)
title('Median filter')
