Local bone enhancement fuzzy clustering for segmentation of MR trabecular bone images

Abstract

Segmentation of trabecular bone from magnetic resonance (MR) images is a challenging task due to spatial resolution limitations, signal-to-noise ratio constraints, and signal intensity inhomogeneities. This article examines an alternative approach to trabecular bone segmentation using partial membership segmentation termed fuzzy C-means clustering incorporating local second order features for bone enhancement (BE-FCM) at multiple scales. This approach is meant to allow for a soft segmentation that accounts for partial volume effects while suppressing the influence of noise. A soft segmentation method was developed and evaluated on three different sets of data; interscan reproducibility was evaluated on six test-retest in vivo MR scans of the proximal femur, correlation between MR and HR-pQCT measurements was evaluated on 49 in vivo scans from the distal tibia, and the potential for fracture discrimination was evaluated using MR scans of calcaneus specimens from 15 participants with and 15 participants without vertebral fracture. The algorithm was compared to fuzzy clustering using the intensity as the only feature (I-FCM) and a dual thresholding algorithm. The metric evaluated was bone volume over total volume (BV/TV) within user-defined regions of interest. BE-FCM had a higher interscan reproducibility (rms CV: 2.0%) compared to I-FCM (5.6%) and thresholding (4.2%), and expressed higher correlation to HR-pQCT data (r = 0.79, p < 10(-11)) compared to I-FCM (r = 0.74, p < 10(-8)) and thresholding (r = 0.70, p < 10(-6)). BE-FCM was also the method that was best able to differentiate between a control and a vertebral fracture group at a 95% significance level. The results suggest that trabecular bone segmentation by BE-FCM can provide a precise BV/TV measurement that is sensitive to pathology. The segmentation method may become useful in MR imaging-based quantification of bone microarchitecture.