Bone image segmentation

A surrogate mechanostatistical microstructural model to inform whole hip cortical bone remodelling

Multiscale models of bone remodelling which incorporate clinically important features at a wide range of spatial scales are computationally intractable. In addition, there is a large gap in existing bone remodelling literature which deals with Haversian-level cortical bone remodelling behaviour and its effect on whole organ bone strength. In this study, we present an integration between continuum mechanics and surrogate modelling as an efficient circumvention to the problem of intractability, linked with an intricate 3D cortical bone remodelling algorithm which captures realistic microanatomical features at the Haversian level. The surrogate model predicted clear differentiation in long-term bone strength changes in clinically significant cases with extremely high computational efficiency.

Segmentation of histology slides of cortical bone using pulse coupled neural networks optimized by particle-swarm optimization

The aim of this study is to automatically discern the micro-features in histology slides of cortical bone using pulse coupled neural networks (PCNN). To the best knowledge of the authors, utilizing PCNN in such an application has not been reported in the literature and, as such, constitutes a novel application. The network parameters are optimized using particle swarm optimization (PSO) where the PSO fitness function was introduced as the entropy and energy of the bone micro-constituents extracted from a training image. Another novel contribution is combining the above with the method of adaptive threshold (T) where the PCNN algorithm is repeated until the best threshold T is found corresponding to the maximum variance between two segmented regions. To illustrate the quality of resulting segmentation according to this methodology, a comparison of the entropy/energy obtained of each pulse is reported. Suitable quality metrics (precision rate, sensitivity, specificity, accuracy, and dice) were used to benchmark the resulting segments against those found by a more traditional method namely K-means. The quality of the segments revealed by this methodology was found to be of much superior quality. Another testament to the quality of this methodology was that the images resulting from testing pulses were found to be of similarly good quality to those of the training images.

Automated segmentation of acetabulum and femoral head from 3-D CT images

This paper describes several new methods and software for automatic segmentation of the pelvis and the femur, based on clinically obtained multislice computed tomography (CT) data. The hip joint is composed of the acetabulum, cavity of the pelvic bone, and the femoral head. In vivo CT data sets of 60 actual patients were used in the study. The 120 (60 x 2) hip joints in the data sets were divided into four groups according to several key features for segmentation. Conventional techniques for classification of bony tissues were first employed to distinguish the pelvis and the femur from other CT tissue images in the hip joint. Automatic techniques were developed to extract the boundary between the acetabulum and the femoral head. An automatic method was built up to manage the segmentation task according to image intensity of bone tissues, size, center, shape of the femoral heads, and other characters. The processing scheme consisted of the following five steps: 1) preprocessing, including resampling 3-D CT data by a modified Sinc interpolation to create isotropic volume and to avoid Gibbs ringing, and smoothing the resulting images by a 3-D Gaussian filter; 2) detecting bone tissues from CT images by conventional techniques including histogram-based thresholding and binary morphological operations; 3) estimating initial boundary of the femoral head and the joint space between the acetabulum and the femoral head by a new approach utilizing the constraints of the greater trochanter and the shapes of the femoral head; 4) enhancing the joint space by a Hessian filter; and 5) refining the rough boundary obtained in step 3) by a moving disk technique and the filtered images obtained in step 4). The above method was implemented in a Microsoft Windows software package and the resulting software is freely available on the Internet. The feasibility of this method was tested on the data sets of 60 clinical cases (5000 CT images).

Segmentation of avascular necrosis of the femoral head using 3-D MR images

Avascular necrosis of the femoral head (ANFH) is a common clinical disorder in the orthopedic field. Traditional approaches to study the extent of ANFH rely primarily on manual segmentation of clinical magnetic resonance images (MRI). However, manual segmentation is insufficient for quantitative evaluation and staging of ANFH. This paper presents a new computerized approach for segmentation of necrotic lesions of the femoral head. The segmentation method consists of several steps including histogram based thresholding, 3-D morphological operations, oblique data reconstruction, and 2-D ellipse fitting. The proposed technique is rapid and efficient. In addition, it is available as a Microsoft Windows free software package on the Internet. Feasibility of the method is demonstrated on the data sets of 30 patients (1500 MR images).