-
Notifications
You must be signed in to change notification settings - Fork 1
/
main.m
82 lines (61 loc) · 2.25 KB
/
main.m
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
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
image = imread('images/kwiaty.png');
%se:
% 1
% 111
% 1
se = [1 0; -1 0; 0 1; 0 -1];
%podstawowe operacje
eroded_image = erode(image, se);
dilated_image = dilate(image, se);
opened_image = openn(image, se);
closed_image = closee(image, se);
figure('Name', 'erozja');
imshow(eroded_image);
imwrite(eroded_image, 'wynik_erozja.png');
figure('Name', 'dylacja');
imshow(dilated_image);
imwrite(dilated_image, 'wynik_dylacja.png');
figure('Name', 'otwarcie');
imshow(opened_image);
imwrite(opened_image, 'wynik_otwarcie.png');
figure('Name', 'zamknięcie');
imshow(closed_image);
imwrite(closed_image, 'wynik_zamkniecie.png');
%gradienty
gradient_erode = erode_gradient(image, se);
gradient_dilate = dilate_gradient(image, se);
gradient_morphology = morphology_gradient(image, se);
figure('Name', 'gradient przez erozje');
imshow(gradient_erode);
imwrite(gradient_erode, 'wynik_gradient_erozja.png');
figure('Name', 'gradient przez dylację');
imshow(gradient_dilate);
imwrite(gradient_dilate, 'wynik_gradient_dylacja.png');
figure('Name', 'gradient morfologiczny');
imshow(gradient_morphology);
imwrite(gradient_morphology, 'wynik_gradient_morfologiczny.png');
%top-hat
top_hat_image = top_hat(image, se);
bottom_hat_image = bottom_hat(image, se);
figure('Name', 'top hat');
imshow(top_hat_image);
imwrite(top_hat_image, 'wynik_top_hat.png');
figure('Name', 'bottom hat');
imshow(bottom_hat_image);
imwrite(bottom_hat_image, 'wynik_bottom_hat.png');
%rekonstrukcja przez dylację -- operacja trwa około 15 min !!!
mask = rgb2hsv(imread('images/paski.png'));
marker = rgb2hsv(imread('images/znacznik.png'));
reconstructed = dilation_reconstruction(marker, mask, se, 30);
imwrite(reconstructed, 'rekonstrukcja_dylacja_wynik.png');
%rekonstrukcja przez erozję -- operacja trwa około 40 min !!!
mask = imread('images/shapes.jpg');
marker = imread('images/znacznik_shapes.png');
reconstructed = erosion_reconstruction(marker, mask, se, 35);
imwrite(reconstructed, 'rekonstrukcja_erozja_wynik.png');
% otwarcie przez rekonstrukcję
opened_image = open_reconstruction(image, se);
imwrite(opened_image, 'otwarcie_rekonstrukcja_wynik.png');
% zamknięcie przez rekonstrukcję
closed_image = close_reconstruction(image, se);
imwrite(closed_image, 'zamkniecie_rekonstrukcja_wynik.png');