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).