ex2.m

for image processing

% Exercise 2
clear
clc

f = [20 20 20 30 40 40 40 40 ;
     20 20 20 30 40 40 40 40 ;
     20 20 20 30 40 40 40 40 ;
     30 30 30 30 30 30 30 30 ;
     40 40 40 30 20 20 20 20 ;
     40 40 40 30 20 50 20 20 ;
     40 40 40 30 20 40 20 20 ;
     40 40 40 30 20 20 20 20];

figure, imshow(uint8((f-20)*(255/30)));
title('Scaled Original Image');

mask_dim = 3;
W = [0 1 0 ; 1 1 1 ; 0 1 0]/5;

% Median Filter
dim = size(f);
k = 1;
for i = 1 : dim(1)-mask_dim+1
    for j = 1 : dim(2)-mask_dim+1
        for m = 1+i-1 : mask_dim+i-1
            for n = 1+j-1 : mask_dim+j-1
                a(k) = f(m,n);
                k = k + 1;
            end
        end
        a = sort(a,'ascend');
        if rem(length(a),2) == 0
            med = (a(length(a)/2) + a(length(a)/2 + 1)) / 2;
        else
            med = a((length(a)-1)/2 + 1);
        end
        g1(i,j) = med;
        k = 1;
    end
end
figure, imshow(uint8((g1-20)*(255/20)));
title('Scaled Median Filtered Image');

% Mask W Filter
g2 = imagefilter(f,W);
figure, imshow(uint8((f-12)*(255/28)));
title('Scaled Mask W Filtered Image');