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

PURPOSE 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. METHODS 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. RESULTS 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. CONCLUSIONS 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.

[1]  J. C. Dunn,et al.  A Fuzzy Relative of the ISODATA Process and Its Use in Detecting Compact Well-Separated Clusters , 1973 .

[2]  S. Cummings,et al.  Bone density at various sites for prediction of hip fractures , 1993, The Lancet.

[3]  S. Majumdar,et al.  Correlation of Trabecular Bone Structure with Age, Bone Mineral Density, and Osteoporotic Status: In Vivo Studies in the Distal Radius Using High Resolution Magnetic Resonance Imaging , 1997, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[4]  Carl-Fredrik Westin,et al.  Using local 3D structure for segmentation of bone from computer tomography images , 1997, Proceedings of IEEE Computer Society Conference on Computer Vision and Pattern Recognition.

[5]  Alejandro F. Frangi,et al.  Muliscale Vessel Enhancement Filtering , 1998, MICCAI.

[6]  S. Majumdar,et al.  In Vivo High Resolution MRI of the Calcaneus: Differences in Trabecular Structure in Osteoporosis Patients , 1998, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[7]  Alan C. Evans,et al.  A nonparametric method for automatic correction of intensity nonuniformity in MRI data , 1998, IEEE Transactions on Medical Imaging.

[8]  Scott N. Hwang,et al.  Estimating voxel volume fractions of trabecular bone on the basis of magnetic resonance images acquired in vivo , 1999 .

[9]  S. Majumdar,et al.  Processing and Analysis of In Vivo High-Resolution MR Images of Trabecular Bone for Longitudinal Studies: Reproducibility of Structural Measures and Micro-Finite Element Analysis Derived Mechanical Properties , 2002, Osteoporosis International.

[10]  Aly A. Farag,et al.  A modified fuzzy c-means algorithm for bias field estimation and segmentation of MRI data , 2002, IEEE Transactions on Medical Imaging.

[11]  Ego Seeman,et al.  Changes in bone mineral density explain little of the reduction in vertebral or nonvertebral fracture risk with anti-resorptive therapy. , 2004, Bone.

[12]  J. Koenderink.  The structure of images , 2004, Biological Cybernetics.

[13]  Claus Christiansen,et al.  Consensus development conference: Prophylaxis and treatment of osteoporosis , 2005, Osteoporosis International.

[14]  Maxime Descoteaux,et al.  Bone Enhancement Filtering: Application to Sinus Bone Segmentation and Simulation of Pituitary Surgery , 2005, MICCAI.

[15]  Roland Krug,et al.  Feasibility of in vivo structural analysis of high-resolution magnetic resonance images of the proximal femur , 2005, Osteoporosis International.

[16]  M. Jergas,et al.  Accurate assessment of precision errors: How to measure the reproducibility of bone densitometry techniques , 2005, Osteoporosis International.

[17]  Felix W. Wehrli,et al.  A novel local thresholding algorithm for trabecular bone volume fraction mapping in the limited spatial resolution regime of in vivo MRI , 2005, IEEE Transactions on Medical Imaging.

[18]  Sharmila Majumdar,et al.  Characterization of trabecular bone structure from high-resolution magnetic resonance images using fuzzy logic. , 2006, Magnetic resonance imaging.

[19]  S. Majumdar,et al.  Autocalibrating parallel imaging of in vivo trabecular bone microarchitecture at 3 Tesla , 2006, Magnetic resonance in medicine.

[20]  F. Wehrli.  Structural and functional assessment of trabecular and cortical bone by micro magnetic resonance imaging , 2007, Journal of magnetic resonance imaging : JMRI.

[21]  S. Majumdar,et al.  Assessment of trabecular bone structure comparing magnetic resonance imaging at 3 Tesla with high-resolution peripheral quantitative computed tomography ex vivo and in vivo , 2008, Osteoporosis International.

[22]  S. Majumdar,et al.  In Vivo Determination of Bone Structure in Postmenopausal Women: A Comparison of HR‐pQCT and High‐Field MR Imaging , 2007, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[23]  Sharmila Majumdar,et al.  Evaluation of correction methods for coil-induced intensity inhomogeneities and their influence on trabecular bone structure parameters from MR images. , 2009, Medical physics.

[24]  G. H. van Lenthe,et al.  Non-invasive bone competence analysis by high-resolution pQCT: an in vitro reproducibility study on structural and mechanical properties at the human radius. , 2009, Bone.

[25]  Baowei Fei,et al.  A modified fuzzy C-means classification method using a multiscale diffusion filtering scheme , 2009, Medical Image Anal..