mmdfila
Detect Filarial Worms.

A microscopic gray-scale image, with two filarial worms, is read.

a = mmreadgray('danaus.tif');

mmshow(a);

a

The Close by Reconstruction Top-Hat operator is applied to regularize the image background.

b = mmcloserecth(a,mmsebox(5));

mmshow(b);

b

The gray-scale opening by the elementary cross is applied to remove narrow objects.

c = mmopen(b);

mmshow(c);

c

The gray-scale area open operator is applied to remove small objects.

d = mmareaopen(c,200);

mmshow(d);

d

The threshold operator is applied to extract a reduced set of structures that include the two worms present in the image.

e = mmthreshad(d,50);

Warning: converting image from scalar to uint8

mmshow(e);

e

The objective of the sequence of transformations, that begin with the homotopic skeleton, is to eliminateg the structures that are not worms. The information used for the filtering is that the worms are longer than any other structure found.

f = mmthin(e);

mmshow(f);

f

The first 12 points of the skeleton branches, counting from their extremities, are eliminated. The structures that were not eliminated will be the markers for extracting the two worms.

g = mmthin(f,mmendpoints, 12);

mmshow(g);

g

The binary reconstruction operator is applied to reconstruct the binary image produced by the threshold from the marker image.

h = mminfrec(g,e);

mmshow(h);

h

The structures extracted are overlaid to the input gray-scale image.

mmshow(a,h);

a,h