Anatomically Anchored Template-Based Level Set Segmentation: Application to Quadriceps Muscles in MR Images from the Osteoarthritis Initiative

In this paper, we present a semi-automated segmentation method for magnetic resonance images of the quadriceps muscles. Our method uses an anatomically anchored, template-based initialization of the level set-based segmentation approach. The method only requires the input of a single point from the user inside the rectus femoris. The templates are quantitatively selected from a set of images based on modes in the patient population, namely, sex and body type. For a given image to be segmented, a template is selected based on the smallest Kullback-Leibler divergence between the histograms of that image and the set of templates. The chosen template is then employed as an initialization for a level set segmentation, which captures individual anatomical variations in the image to be segmented. Images from 103 subjects were analyzed using the developed method. The algorithm was trained on a randomly selected subset of 50 subjects (25 men and 25 women) and tested on the remaining 53 subjects. The performance of the algorithm on the test set was compared against the ground truth using the Zijdenbos similarity index (ZSI). The average ZSI means and standard deviations against two different manual readers were as follows: rectus femoris, 0.78 +/- 0.12; vastus intermedius, 0.79 +/- 0.10; vastus lateralis, 0.82 +/- 0.08; and vastus medialis, 0.69 +/- 0.16.