ant = imread("ant.jpg")
burnedAnt = createAntSegmentation(ant)
fig = figure();
ax(1) = axes('Units','normalized','Position', [ .1 .1 .4 .8]);
ax(2) = axes('Units','normalized','Position', [ .5 .1 .4 .8]);
imshow(ant,'Parent', ax(1))
imshow(burnedAnt, 'Parent', ax(2))
linkaxes(ax)
function burnedAnt = createAntSegmentation(ant)
greyant = rgb2gray(ant);
adaptedAnt = adaptthresh(greyant, 0.38,"ForegroundPolarity","dark");
BW = imbinarize(greyant,adaptedAnt);
BWopen = bwareaopen(~BW,800);
se = strel("disk",32);
BWdilate = imdilate(BWopen,se);
BWfilled = imfill(BWdilate,"holes");
burnedAnt = imoverlay(ant,~BWfilled,"k");
end