%Read cameraman
I = imread('cameraman.tif');
se = strel('ball',5,5);

figure('Name', 'original');
imshow(I);

figure('Name', 'Basic morphological operations 5x5 ball');
Ierode = imerode(I,se);
Idilate = imdilate(I,se);
Iclose = imclose(I,se);
Iopen = imclose(I, se);
subplot(2,2,1), imshow(Ierode), title('Erode');
subplot(2,2,2), imshow(Idilate), title('Dilate');
subplot(2,2,3), imshow(Iopen), title('Open');
subplot(2,2,4), imshow(Iclose), title('Close');

%Look at different structuring elements
se1 = strel('square',11);      % 11-by-11 square
se2 = strel('line',10,45);     % line, length 10, angle 45 degrees
se3 = strel('disk',3);         % disk, radius 3
se4 = strel('ball',15,5);      % ball, radius 15, height 5

figure('Name', 'Erode');
Ierode1 = imerode(I,se1);
Ierode2 = imerode(I,se2);
Ierode3 = imerode(I,se3);
Ierode4 = imerode(I,se4);
subplot(2,2,1), imshow(Ierode1), title('11x11 square');
subplot(2,2,2), imshow(Ierode2), title('45 deg line len 10');
subplot(2,2,3), imshow(Ierode3), title('Disk, radius 3');
subplot(2,2,4), imshow(Ierode4), title('Ball, radius 15, height 5');


figure('Name', 'Dilate');
Idilate1 = imdilate(I,se1);
Idilate2 = imdilate(I,se2);
Idilate3 = imdilate(I,se3);
Idilate4 = imdilate(I,se4);
subplot(2,2,1), imshow(Idilate1), title('11x11 square');
subplot(2,2,2), imshow(Idilate2), title('45 deg line len 10');
subplot(2,2,3), imshow(Idilate3), title('Disk, radius 3');
subplot(2,2,4), imshow(Idilate4), title('Ball, radius 15, height 5');

figure('Name', 'Open');
Iopen1 = imopen(I,se1);
Iopen2 = imopen(I,se2);
Iopen3 = imopen(I,se3);
Iopen4 = imopen(I,se4);
subplot(2,2,1), imshow(Iopen1), title('11x11 square');
subplot(2,2,2), imshow(Iopen2), title('45 deg line len 10');
subplot(2,2,3), imshow(Iopen3), title('Disk, radius 3');
subplot(2,2,4), imshow(Iopen4), title('Ball, radius 15, height 5');


figure('Name', 'Close');
Iclose1 = imclose(I,se1);
Iclose2 = imclose(I,se2);
Iclose3 = imclose(I,se3);
Iclose4 = imclose(I,se4);
subplot(2,2,1), imshow(Iclose1), title('11x11 square');
subplot(2,2,2), imshow(Iclose2), title('45 deg line len 10');
subplot(2,2,3), imshow(Iclose3), title('Disk, radius 3');
subplot(2,2,4), imshow(Iclose4), title('Ball, radius 15, height 5');