A comparative study of automatic thresholding approaches for 3D x-ray microtomography of trabecular bone.

PURPOSE This paper presents a comparative study of automatic thresholding algorithms for segmenting trabecular bone volume in x-ray microtomography (μCT). METHODS First, a preprocessing stage was established, which considered noise reduction by applying anisotropic diffusion filtering and contrast enhancement by using morphological top-hats. Next, four automatic thresholding algorithms were implemented: clustering, maximum entropy, moment preservation, and concavity-based. These approaches analyze the preprocessed 3D μCT image histogram to optimize some parameters to find the best gray-level threshold. Thirty-eight vertebra bone samples were acquired from 19 normal Wistar rats, specifically the L3 and L4 vertebrae. The μCT images were acquired with a microfocus x-ray device at 100 slices/sample. Next, three human operators segmented the entire 3D μCT images manually to establish ground-truth segmentations so as to associate the segmentation problem with perceptual grouping. The normalized probabilistic Rand index (NPRI) was used to quantify the agreement between each computerized segmentation and the corresponding set of three ground-truth segmentations. Hence, the NPRI value should tend toward unity for an acceptable performance. Finally, a statistical analysis was done to determine which thresholding approach achieved the best performance. Besides, 3D morphometric indices were also measured. RESULTS The Games-Howell test (α = 0.05) was used to compare the equality of means from the NPRI results considering the four thresholding algorithms (multiple comparisons). This statistical analysis indicated that the clustering and moment preservation techniques performed similarly, with NPRI values of 0.594 ± 0.126 and 0.607 ± 0.127, respectively. CONCLUSIONS The main advantage of computerized segmentation is that it is fully automatic; that is, no interaction with the user is required. Thus, the method could be considered objective. Besides, the proposed preprocessing stage plays an important role in enhancing the μCT image quality to achieve better separation between the background volume and the trabecular bone volume.

[1]  Ralph Müller,et al.  Guidelines for assessment of bone microstructure in rodents using micro–computed tomography , 2010, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[2]  Andrew K. C. Wong,et al.  A new method for gray-level picture thresholding using the entropy of the histogram , 1985, Comput. Vis. Graph. Image Process..

[3]  Tiexiang Wen,et al.  Scale selection for morphological top-hat transformation based on mutual information , 2010, 2010 3rd International Congress on Image and Signal Processing.

[4]  Sharmila Majumdar,et al.  A Local Adaptive Threshold Strategy for High Resolution Peripheral Quantitative Computed Tomography of Trabecular Bone , 2007, Annals of Biomedical Engineering.

[5]  P. Filzmoser,et al.  Statistical Data Analysis Explained , 2008 .

[6]  P. Revell,et al.  Histomorphometry of bone. , 1983, Journal of clinical pathology.

[7]  Joaquim Teixeira de Assis,et al.  Bone quality analysis using X-ray microtomography and microfluorescence. , 2012, Applied radiation and isotopes : including data, instrumentation and methods for use in agriculture, industry and medicine.

[8]  Azriel Rosenfeld,et al.  Histogram concavity analysis as an aid in threshold selection , 1983, IEEE Transactions on Systems, Man, and Cybernetics.

[9]  G. Matsopoulos,et al.  Active segmentation of micro-CT trabecular bone images , 2011, 2011 10th International Workshop on Biomedical Engineering.

[10]  Theocharis Antoniadis,et al.  Bone labelling on micro-magnetic resonance images , 1999, Medical Image Anal..

[11]  Jitendra Malik,et al.  Scale-Space and Edge Detection Using Anisotropic Diffusion , 1990, IEEE Trans. Pattern Anal. Mach. Intell..

[12]  Nels A Hanson,et al.  Alternative approach to assessment of bone quality using micro-computed tomography. , 2004, Bone.

[13]  Guillaume Lavoué,et al.  A comparative study of existing metrics for 3D-mesh segmentation evaluation , 2010, The Visual Computer.

[14]  Guido Gerig,et al.  Nonlinear anisotropic filtering of MRI data , 1992, IEEE Trans. Medical Imaging.

[15]  Willi A Kalender,et al.  Ring artifact correction for high-resolution micro CT , 2009, Physics in medicine and biology.

[16]  P. Croucher,et al.  A computerised technique for the quantitative assessment of resorption cavities in trabecular bone. , 1990, Bone.

[17]  Tony F. Chan,et al.  Active contours without edges , 2001, IEEE Trans. Image Process..

[18]  F Peyrin,et al.  A method for the automatic characterization of bone architecture in 3D mice microtomographic images. , 2003, Computerized Medical Imaging and Graphics.

[19]  M Ding,et al.  Accuracy of cancellous bone volume fraction measured by micro-CT scanning. , 1999, Journal of biomechanics.

[20]  Chen Li,et al.  Automatic Thresholding of Micro-CT Trabecular Bone Images , 2008, 2008 International Conference on BioMedical Engineering and Informatics.

[21]  Godfried T. Toussaint,et al.  On the application of the convex hull to histogram analysis in threshold selection , 1983, Pattern Recognit. Lett..

[22]  Martial Hebert,et al.  Toward Objective Evaluation of Image Segmentation Algorithms , 2007, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[23]  Fuhui Long,et al.  Feature selection based on mutual information criteria of max-dependency, max-relevance, and min-redundancy , 2003, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[24]  Martial Hebert,et al.  Measures of Similarity , 2005, 2005 Seventh IEEE Workshops on Applications of Computer Vision (WACV/MOTION'05) - Volume 1.

[25]  Joaquim Teixeira de Assis,et al.  Non-Destructive Imaging Materials Investigation by Microfocus 3D X-Ray Computed Tomography , 2009, IEEE Transactions on Nuclear Science.

[26]  Scott T. Acton,et al.  Speckle reducing anisotropic diffusion , 2002, IEEE Trans. Image Process..

[27]  N. Otsu.  A threshold selection method from gray level histograms , 1979 .

[28]  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.

[29]  Julia F. Barrett,et al.  Artifacts in CT: recognition and avoidance. , 2004, Radiographics : a review publication of the Radiological Society of North America, Inc.

[30]  P. Rüegsegger,et al.  Direct Three‐Dimensional Morphometric Analysis of Human Cancellous Bone: Microstructural Data from Spine, Femur, Iliac Crest, and Calcaneus , 1999, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[31]  William E. Lorensen,et al.  Marching cubes: A high resolution 3D surface construction algorithm , 1987, SIGGRAPH.

[32]  I Lima,et al.  Ethanol bone evaluation using 3D microtomography. , 2008, Micron.

[33]  Michael J. Flynn,et al.  Microfocus X-ray sources for 3D microtomography , 1994 .

[34]  S Tetradis,et al.  MicroCT morphometry analysis of mouse cancellous bone: intra- and inter-system reproducibility. , 2011, Bone.

[35]  John F. Canny,et al.  A Computational Approach to Edge Detection , 1986, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[36]  T N Hangartner,et al.  Thresholding technique for accurate analysis of density and geometry in QCT, pQCT and microCT images. , 2007, Journal of musculoskeletal & neuronal interactions.

[37]  M. Bouhlel,et al.  Trabecular Bone Image Segmentation Using Iterative Watershed and Multi Resolution Analysis , 2011 .

[38]  F Guilak,et al.  Volume and surface area measurement of viable chondrocytes in situ using geometric modelling of serial confocal sections , 1994, Journal of microscopy.

[39]  Aly A. Farag,et al.  Segmentation of trabecular bones from Vertebral bodies in volumetric CT spine images , 2009, 2009 16th IEEE International Conference on Image Processing (ICIP).

[40]  D M Spengler,et al.  Multiplanar variations in the structural characteristics of cancellous bone. , 1994, Bone.

[41]  G. Tysarczyk-Niemeyer,et al.  Application of automatic image segmentation to tibiae and vertebrae from ovariectomized rats. , 1997, Bone.

[42]  Wen-Hsiang Tsai,et al.  Moment-preserving thresholding: a new approach , 1995 .

[43]  Françoise Peyrin,et al.  A constrained region growing approach based on watershed for the segmentation of low contrast structures in bone micro-CT images , 2008, Pattern Recognit..