mmdBrain
Extract the lateral ventricle from an MRI image of the brain.

The MRI image of a brain slice is read.

Image a = mmReadgray("mribrain.tif");
mmShow(a);

a

The ventricle is enhanced using an opening with a disk of radius 10 followed by a reconstruction.

Image b = mmOpen(a, mmSedisk(10));
Image c = mmInfrec(b, a);
mmShow(b);
mmShow(c);

b	c

The result of the open by reconstruction is subtracted from the original image. Note that the three operations: open, reconstruction and the subtraction could be done at once using the mmOpenrecth (open by reconstruction top-hat) function. On the right, the enhanced ventricle is thresholded.

Image d = mmSubm(a, c);
mmShow(d);
Image e = mmThreshad(d, 50);
mmShow(e);

d	e

Finally, the ventricle is selected as the connected object with area larger than 70 pixels. For visualization purposes, the result of the segmentation is overlayed on the original brain image.

Image f = mmAreaopen(e, 70);
mmShow(f);
mmShow(a, f);

f	a, f