% Median Filter function
% usage:
% [output_img] = imgMedianFilter ([input_img], mask_size)
% mask_size - size of matrix that goes through image finding average

function [img_out] = imgMedianFilter (img_in, mask_size)

i_max = size(img_in,1);         % X-size of input image
j_max = size(img_in,2);         % Y-size of input image

border = mask_size - 2;
center = round(mask_size^2/2);
img_in = imgResize(img_in,border);

for i=1:1:i_max
    for j=1:1:j_max
        A = img_in(i:i+mask_size-1,j:j+mask_size-1);
        B = A(:);
        C = sort(B);
        img_out(i,j) = C(center);
    end
end