Podstawowe operacje morfologiczne

close all;
clear all;

IM_1 = imread('ertka.bmp');
se = strel('square',3);
erIM_1 = imerode(IM_1,se);
erIM_2 = imerode(erIM_1,se);
erIM_3 = imerode(erIM_2,se);

figure(1)
subplot(2,2,1)
imshow(IM_1)
title('Obraz originalny')
subplot(2,2,2)
imshow(erIM_1)
title('Obraz po pierwszej erozji')
subplot(2,2,3)
imshow(erIM_2)
title('Obraz po drugiej erozji')
subplot(2,2,4)
imshow(erIM_3)
title('Obraz po trzeciej erozji')

%---------------------------------------

IM_3 = imread('buska.bmp');
se2 = [0 1 0,
       0 0 0,
       0 1 0]
bezIM_3 = imerode(IM_3,se2);
figure(2)
subplot(1,2,1)
imshow(IM_3)
title('Obraz originalny')
subplot(1,2,2)
imshow(bezIM_3)
title('Obraz bez')

%---------------------------------

IM_2 = imread('ertka.bmp');
se1 = strel('square',3);
erIM_2 = imerode(IM_2,se1);
dyIM_2 = imdilate(IM_2,se1);
opIM_2 = imopen(IM_2,se1);
zaIM_2 = imclose(IM_2,se1);
samIM_2 = imerode(IM_2,se1);

temp = imdilate(samIM_2,se1)
temp1 = imdilate(temp,se1)
temp2 = imdilate(temp1,se1)

figure(3)
subplot(2,3,1)
imshow(IM_2)
title('Obraz originalny')
subplot(2,3,2)
imshow(erIM_2)
title('Obraz po erozji')
subplot(2,3,3)
imshow(dyIM_2)
title('Obraz po dylatacji')
subplot(2,3,4)
imshow(opIM_2)
title('Obraz po otwarciu')
subplot(2,3,5)
imshow(zaIM_2)
title('Obraz po zamknieciu')
subplot(2,3,6)
imshow(temp2)
title('SAMO RT')

%------------------------------------

IM_4 = imread('hom.bmp');

SE1=[0 1 0; 1 1 1; 0 1 0];
SE2=[1 0 1; 0 0 0; 1 0 1];

BW2 = bwhitmiss(IM_4,SE1,SE2)

figure(4)
subplot(1,2,1)
imshow(IM_4)
title('Obraz originalny')
subplot(1,2,2)
imshow(BW2)
title('Obraz po operacji transformacji trafi - nie trafi')