clear all; close all;

img = imread('PictureBW.bmp');

imshow(img);
title('ORYGINALNY');
figure


saltpepper = imnoise(img, 'salt & pepper', 0.02);
imshow(saltpepper);
title('SÓL I PIEPRZ');
figure

mediana1x1 = medfilt2(saltpepper, [1 1]);
imshow(mediana1x1);
title('MEDIANA 1x1');
figure

mediana3x3 = medfilt2(saltpepper, [3 3]);
imshow(mediana3x3);
title('MEDIANA 3x3');
figure

mediana5x5 = medfilt2(saltpepper, [5 5]);
imshow(mediana5x5);
title('MEDIANA 5x5');
figure
mediana7x7 = medfilt2(saltpepper, [7 7]);
imshow(mediana7x7);
title('MEDIANA 7x7');

medianowa = [2,1,2;1,5,1;2,1,2];










oryginalZmedniana3x3 = medfilt2(img, [3 3]);
oryginalZmedniana5x5 = medfilt2(img, [5 5]);
oryginalZmedniana7x7 = medfilt2(img, [7 7]);


Aerror3x3 = 0;
AAerror5x5 = 0;
AAAerror7x7 = 0;
ERROR1x1 = 0;
ERROR3x3 = 0;
ERROR5x5 = 0;
ERROR7x7 = 0;

for m = 1:512
    for n = 1:512
        ERROR1x1 = (1/(512*512)) * (sum( sum( ((img(m,n) - mediana1x1(m,n))^2) )));
        ERROR3x3 = ERROR3x3 + (1/(512*512)) * (sum( sum( ((img(m,n) - mediana3x3(m,n))^2) )));
        ERROR5x5 = ERROR5x5 + (1/(512*512)) * (sum( sum( ((img(m,n) - mediana5x5(m,n))^2) )));
        ERROR7x7 = ERROR7x7 + (1/(512*512)) * (sum( sum( ((img(m,n) - mediana7x7(m,n))^2) )));

        Aerror3x3 = Aerror3x3 + (1/(512*512)) * (sum( sum( ((img(m,n) - oryginalZmedniana3x3(m,n))^2) )));
        AAerror5x5 = AAerror5x5 + (1/(512*512)) * (sum( sum( ((img(m,n) - oryginalZmedniana5x5(m,n))^2) )));
        AAAerror7x7 = AAAerror7x7 + (1/(512*512)) * (sum( sum( ((img(m,n) - oryginalZmedniana7x7(m,n))^2) )));
    end
end