close all
clc
%%% Original
original = rgb2gray(imread('fabolous.png'));
original = im2bw(255-original,0.5);
profe=imread('fabolous.png');
imshow(original)
%"Disk"
disk3 = imerode(original,strel('disk',3));
disk5 = imerode(original,strel('disk',5));
disk7 = imerode(original,strel('disk',7));
% "Diamond"
diamond3 = imerode(original,strel('diamond',3));
diamond5 = imerode(original,strel('diamond',5));
diamond7 = imerode(original,strel('diamond',7));
% "Square"
square3 = imerode(original,strel('square',3));
square5 = imerode(original,strel('square',5));
square7 = imerode(original,strel('square',7));
% Dilatación
% "Disk"
disk3di = imdilate(original,strel('disk',3));
disk5di = imdilate(original,strel('disk',5));
disk7di = imdilate(original,strel('disk',7));
% "Diamond"
diamond3di = imdilate(original,strel('diamond',3));
diamond5di = imdilate(original,strel('diamond',5));
diamond7di = imdilate(original,strel('diamond',7));
% "Square"
square3di = imdilate(original,strel('square',3));
square5di = imdilate(original,strel('square',5));
square7di = imdilate(original,strel('square',7));
% Apertura
% "Disk"
disk3open = imopen(original,strel('disk',3));
disk5open = imopen(original,strel('disk',5));
disk7open = imopen(original,strel('disk',7));
% "Diamond"
diamond3open = imopen(original,strel('diamond',3));
diamond5open = imopen(original,strel('diamond',5));
diamond7open = imopen(original,strel('diamond',7));
% "Square"
square3open = imopen(original,strel('square',3));
square5open = imopen(original,strel('square',5));
square7open = imopen(original,strel('square',7));
% Cierre
% "Disk"
disk3close = imclose(original,strel('disk',3));
disk5close = imclose(original,strel('disk',5));
disk7close = imclose(original,strel('disk',7));
% "Diamond"
diamond3close = imclose(original,strel('diamond',3));
diamond5close = imclose(original,strel('diamond',5));
diamond7close = imclose(original,strel('diamond',7));
% "Square"
square3close = imclose(original,strel('square',3));
square5close = imclose(original,strel('square',5));
square7close = imclose(original,strel('square',7));
%SKEL
BW5 = bwmorph(disk5,'skel',Inf);
BW6 = bwmorph(diamond5,'skel',Inf);
BW7 = bwmorph(square5,'skel',Inf);
BW8 = bwmorph(disk5di,'skel',Inf);
BW9 = bwmorph(diamond5di,'skel',Inf);
BW10 = bwmorph(square5di,'skel',Inf);
BW11 = bwmorph(disk5open,'skel',Inf);
BW12= bwmorph(diamond5open,'skel',Inf);
BW13 = bwmorph(square5open,'skel',Inf);
BW14 = bwmorph(disk5close,'skel',Inf);
BW15 = bwmorph(diamond5close,'skel',Inf);
BW16 = bwmorph(square5close,'skel',Inf);
figure
subplot(4,3,1),imshow(BW5),title('Original');
subplot(4,3,2),imshow(BW6),title('Erosionada');
subplot(4,3,3),imshow(BW7),title('Dilatada');
subplot(4,3,4),imshow(BW8),title('Abierta');
subplot(4,3,5),imshow(BW9),title('Cerrada');
subplot(4,3,6),imshow(BW10),title('Equeleto');
subplot(4,3,7),imshow(BW11),title('Original');
subplot(4,3,8),imshow(BW12),title('Erosionada');
subplot(4,3,9),imshow(BW13),title('Dilatada');
subplot(4,3,10),imshow(BW14),title('Abierta');
subplot(4,3,11),imshow(BW15),title('Cerrada');
subplot(4,3,12),imshow(BW16),title('Equeleto');

figure
subplot(3,2,1),imshow(original),title('Original');
subplot(3,2,2),imshow(erosionada),title('Erosionada');
subplot(3,2,3),imshow(dilatada),title('Dilatada');
subplot(3,2,4),imshow(abierta),title('Abierta');
subplot(3,2,5),imshow(cerrada),title('Cerrada');
subplot(3,2,6),imshow(BW2),title('Equeleto');
figure
subplot(3,1,1),imshow(disk3),title('"Erosionada-Disk3"');
subplot(3,1,2),imshow(disk5),title('"Erosionada-Disk5"');
subplot(3,1,3),imshow(disk7),title('"Erosionada-Disk7"');
figure
subplot(3,1,1),imshow(diamond3),title('"Erosionada-Diamond3"');
subplot(3,1,2),imshow(diamond5),title('"Erosionada-Diamond5"');
subplot(3,1,3),imshow(diamond7),title('"Erosionada-Diamond7"');
figure
subplot(3,1,1),imshow(square3),title('"Erosionada-Square3"');
subplot(3,1,2),imshow(square5),title('"Erosionada-Square5"');
subplot(3,1,3),imshow(square7),title('"Erosionada-Square7"');
%%%
figure
subplot(3,1,1),imshow(disk3di),title('"Dilatada-Disk3"');
subplot(3,1,2),imshow(disk5di),title('"Dilatada-Disk5"');
subplot(3,1,3),imshow(disk7di),title('"Dilatada-Disk7"');
 figure
subplot(3,1,1),imshow(diamond3di),title('"Dilatada-Diamond3"');
subplot(3,1,2),imshow(diamond5di),title('"Dilatada-Diamond5"');
subplot(3,1,3),imshow(diamond7di),title('"Dilatada-Diamond7"');
 figure
subplot(3,1,1),imshow(square3di),title('"Dilatada-Square3"');
subplot(3,1,2),imshow(square5di),title('"Dilatada-Square5"');
subplot(3,1,3),imshow(square7di),title('"Dilatada-Square7"');

 figure
subplot(3,1,1),imshow(disk3open),title('"Abierta-Disk3"');
subplot(3,1,2),imshow(disk5open),title('"Abierta-Disk5"');
subplot(3,1,3),imshow(disk7open),title('"Abierta-Disk7"');
 figure
subplot(3,1,1),imshow(diamond3open),title('"Abierta-Diamond3"');
subplot(3,1,2),imshow(diamond5open),title('"Abierta-Diamond5"');
subplot(3,1,3),imshow(diamond7open),title('"Abierta-Diamond7"');
 figure
subplot(3,1,1),imshow(square3open),title('"Abierta-Square3"');
subplot(3,1,2),imshow(square5open),title('"Abierta-Square5"');
subplot(3,1,3),imshow(square7open),title('"Abierta-Square7"');
figure
subplot(3,1,1),imshow(disk3close),title('"Cerrada-Disk3"');
subplot(3,1,2),imshow(disk5close),title('"Cerrada-Disk5"');
subplot(3,1,3),imshow(disk7close),title('"Cerrada-Disk7"');
figure
subplot(3,1,1),imshow(diamond3close),title('"Cerrada-Diamond3"');
subplot(3,1,2),imshow(diamond5close),title('"Cerrada-Diamond5"');
subplot(3,1,3),imshow(diamond7close),title('"Cerrada-Diamond7"');
figure
subplot(3,1,1),imshow(square3close),title('"Cerrada-Square3"');
subplot(3,1,2),imshow(square5close),title('"Cerrada-Square5"');
subplot(3,1,3),imshow(square7close),title('"Cerrada-Square7"');
C5=imfill(BW16,'holes');
plot(C5)
figure
plot(C5)
imshow(C5)