Image_processing/Task6_code.m at master · murtuzaal2/Image_processing · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
clear;

u0 = imread('test2.png');
u0 = im2double(u0);
u0 = imnoise(u0,'salt & pepper',0.1);
[m,n] = size(u0);
u = u0;

for i = 3:m-1
    for j = 3:n-2

        mean = 0;

        for x1 = -1:1:1
            for x2 = -1:1:1
                mean = mean+ u(i+x1,j+x2);
            end
        end

       for x1 = -1:1:1
            for x2 = -1:1:1

                B = ((u(i+x1,j+x2) == 1)||(u(i+x1,j+x2) == 0));

                if (B)
                    u(i+x1,j+x2) = (mean - u(i+x1,j+x2))/8;
                end

            end
        end


    end
end

figure
subplot(1,3,1),
imshow(u0),title('Noised(S&P)');
subplot(1,3,2),
imshow(u),title('by Average method');
subplot(1,3,3),
imshow(medfilt2(u0)),title('Using "medfilt2"');