clear all
imgorig = imread("111-0.jpg");
img = rgb2gray(imgorig);
imgns = imnoise(img, "gaussian", 0, 0.08);
x = length(img(:,1));
y = length(img(1,:));
for i=2:x-1
 for j=2:y-1
 imgcl(i, j) =  imgns(i, j)/9+ imgns(i-1,j-1)/9+ imgns(i+1,j+1)/9+ imgns(i+1,j)/9+ imgns(i-1,j)/9+ imgns(i,j+1)/9+ imgns(i,j-1)/9+ imgns(i-1,j+1)/9+ imgns(i+1,j-1)/9;
 imgcl2(i, j) = median(
 [imgns(i, j),
 imgns(i-1,j-1),
 imgns(i+1,j+1),
 imgns(i+1,j),
 imgns(i-1,j),
 imgns(i,j+1),
 imgns(i,j-1),
 imgns(i-1,j+1),
 imgns(i+1,j-1)]);
 endfor
endfor
subplot(1, 3, 1), subimage(imgns);
subplot(1, 3, 2), subimage(imgcl);
subplot(1, 3, 3), subimage(imgcl2);