Volume 9 Supplement 1
Unified wavelet and gaussian filtering for segmentation of CT images; application in segmentation of bone in pelvic CT images
© Vasilache et al; licensee BioMed Central Ltd. 2009
Published: 3 November 2009
The analysis of pelvic CT scans is a crucial step for detecting and assessing the severity of Traumatic Pelvic Injuries. Automating the processing of pelvic CT scans could impact decision accuracy, decrease the time for decision making, and reduce health care cost. This paper discusses a method to automate the segmentation of bone from pelvic CT images. Accurate segmentation of bone is very important for developing an automated assisted-decision support system for Traumatic Pelvic Injury diagnosis and treatment.
The automated method for pelvic CT bone segmentation is a hierarchical approach that combines filtering and histogram equalization, for image enhancement, wavelet analysis and automated seeded region growing. Initial results of segmentation are used to identify the region where bone is present and to target histogram equalization towards the specific area. Speckle Reducing Anisotropic Didffusion (SRAD) filter is applied to accentuate the desired features in the region. Automated seeded region growing is performed to refine the initial bone segmentation results.
The proposed method automatically processes pelvic CT images and produces accurate segmentation. Bone connectivity is achieved and the contours and sizes of bones are true to the actual contour and size displayed in the original image. Results are promising and show great potential for fracture detection and assessing hemorrhage presence and severity.
Preliminary experimental results of the automated method show accurate bone segmentation. The novelty of the method lies in the unique hierarchical combination of image enhancement and segmentation methods that aims at maximizing the advantages of the combined algorithms. The proposed method has the following advantages: it produces accurate bone segmentation with maintaining bone contour and size true to the original image and is suitable for automated bone segmentation from pelvic CT images.
Traumatic injuries are leading cause of death for patients with ages up to 45 years. Every year, on U.S territory alone, a staggering total of four million years of potential life are cut short due to traumatic injuries . Forty percent of the patients with fatal traumatic injuries die before even reaching the emergency room . Motor vehicle crashes account for 48% to 68% of traumatic injuries. Among traumatic injuries, Traumatic Pelvic Injury and associated complications, such as internal hemorrhage, infected hematomas, multi-organ failure and blood clots traveling to the brain, result in a mortality rate ranging from 8.6% to 50% . Even when the injury is not fatal it is usually the cause of life long disabilities . Computer-aided systems can impact trauma decision making by increasing decision accuracy and reducing time for decision making. Such effects can, in turn, lead to improved health care standards, better resource allocation and lower health care costs.
Traumatic injury decision making is extremely time sensitive therefore a decrease in decision making time is extremely significant for patient treatment. Because traumatic injuries are associated with specific causes and treatment methods, the chance of occurrence for fatality or long-term disability can be avoided or reduced by making more accurate decisions in the trauma unit . Although several computer-assisted trauma decision making systems already exist, the majority of the systems rely solely on patient demographics to find similar cases in trauma databases and provide a recommendation based on these cases. As a result the recommendations might not be accurate or specific enough for the purpose of practical implementation .
Detecting the presence and extent of the fracture is an important step in assessing the severity of a pelvic injury. Subtle fractures, as fractures of the acetabulum, hip displacement or presence of hemorrhage, are assessed based on the analysis of Computed Tomography (CT) images, as CT is more detailed than X-Ray. Since the number of CT slices can be quite large, it is important to develop a computer aided decision making system to analyze the data, because in general only a small portion of the dataset becomes important in establishing a diagnostic .
Segmentation of CT images poses a series of typical challenges. The density of cortical bone is significantly different than that of cancellous bone, as a result the appearance and overall grey level of the two types of bone tissue are very different. Cortical bone is bright and smooth while cortical bone has a spongy texture and an intensity similar to that of the surrounding tissue. It is for this reason that the bone density cannot be uniformly characterized . Additionally, CT scans typically includes tens of slices in which bones assume different shapes and positions depending on the patient vertical and horizontal alignment during the scanning procedure. A number of slices include joint regions where the normal distance between bones is very small. Such small distances can result in separate bones erroneously being merged with each other during segmentation.
In this paper, an automated method to segment bone from pelvic CT images is presented. Another portion of our work addresses hemorrhage detection in pelvic CT images using advanced image processing techniques. The following paragraphs will briefly discuss several existing segmentation techniques used for medical image segmentation.
X-ray computed tomography (CT) images are typically affected by partial volume effects which cause the transitions from bone to air or bone to fat to be less noticeable . Such faded transitions result in similar grey level values of bone and surrounding soft tissue. A popular approach in dealing with low contrast images and noisy edges is using methods that are based on deformable models (DMM) . DMMs allow knowledge about adjacent structures to be incorporated in the model and therefore are capable of producing more accurate segmentation . Deformable Model Methods (DMM) [9–13] are parametric segmentation methods that use closed curves or surfaces that can be deformed under internal and external forces. The intrinsic advantage of such methods is directly generating closed curves while some of the disadvantages include poor convergence towards concave boundaries and sensitivity to initialization. Another category of popular segmentation techniques are based on watershed transform. Watershed segmentation techniques [6, 14–17] are gradient based methods that the image based on its topology. Watershed methods can often result in over-segmentation of the image. Approaches to minimize over-segmentation include: marker-based watershed and low pass filtering of the image prior to watershed segmentation.
Region Growing methods [18–20] are very popular for the segmentation of medical images. In their classic form such methods segment an image by iteratively adding neighboring pixels, that satisfy certain similarity constraints, to the initial seeds. The major advantages of Region Growing are the simple concept, implementation and possibility of using multiple similarity criteria. Region Growing methods share the common disadvantage of being dependent on initialization - quality of segmentation is dependent of initialization quality.
Region Competition segmentation [21, 22] combines statistical properties of region growing methods with geometrical features of deformable models. Region Competition methods are also dependent on initialization and the precision of boundary detection is dependent on the size of the sampling window as well as the signal to noise ratio.
Level Set Methods (LSM) [23–26] are based on defining a moving contour as the zero-level set of a time-evolving scalar function defined over a regular grid. The curve is deformed according to the solution given by a set of partial differential equations (PDEs). As to how the initial curve can evolve, there are several different approaches: some apply a global update of the function while others are based on local update techniques.
Create a binary version of the raw image. All the pixels with grey level higher than zero become one. Pixels with grey value zero remain zero in the new image;
Apply morphological operations to the binary image so that the different objects in the image are separated;
Use blob analysis to select the object that has the largest area;
Detect the edge of the object. Create an alternate edge with better symmetry;
Create a mask that follows the edge created at step 4.
Initial mask formation
Apply a Gaussian filter to the original image;
Apply Wavelet analysis to the resulting image. Reconstruct the image using only the approximation matrix;
Calculate the means and deviation of the grey scale values that differ from zero, i.e. the pixels that are not background;
Calculate threshold value as the sum of means and standard deviation of the grey scale values calculated at Step 3;
- 5.Create an initial binary mask for bone by thresholding the low pass image;
Determine the bounding box for the bone region. Constrain following segmentation to the identified bounding box;
Perform histogram equalization on the image.
Forming seeds close to contour
Determine the optimal size of the filter window based on σ - the standard deviation of the grey level values of the image pixels;
Apply an optimal two dimensional Gaussian filter to a copy of the original image;
- 3.Obtain the intermediary bone mask through binary multiplication of the initial bone mask and the Gaussian filtered image:(6)
where I O (x, y; t) is the original image, I(x, y; t) is the output image, Ω is the image support, ∂Ω is the border of Ω, is the outer normal to ∂Ω and c(q) is the instantaneous coefficient of variation.
For each seed in the seeded bone mask the neighbors in an n × n neighborhood are identified (in this study n = 3);
For each seed in the seeded bone mask the values of the gradient along the eight possible directions are calculated in the n × n neighborhood;
For each of the neighbors in the nxn neighborhood their mxm neighborhood is identified (in this study m = 9);
Average grey level value of the neighbors in the mxm neighborhood is calculated;
- 5.A decision of adding a neighbor of a seed to the seeded bone mask is made based on two conditions:
he neighbor grey level value is greater than the average grey level value of the neighbors in the m × m neighborhood;
the gradient value corresponding to the neighbor is greater or equal to -1. This is a conservative criteria is meant to keep the seeds from growing outside of the edges of the bone and minimizes the risk of separate bones being merged.
The results provided by the proposed method will be discussed in what follows.
The testing CT dataset is courtesy of the Carolinas Healthcare System (CHS). CT scans pertain to Traumatic Pelvic Injury patients. The testing dataset consists of approximately 200 axial pelvic CT images from different segments of the CT scan.
Results of the proposed method
The algorithm has an approximate processing time of 30 seconds for each CT image, from initial input of the image to providing segmentation results. Such a performance suggests a considerable reduction in processing time compared to semi-automated segmentation methods. Semi-automated segmentation methods often require manual selection of seeds and many times such selection is performed in a repetitive process to ensure adequacy of seeds.
Conclusion and future work
The paper is describing a method for automatic bone segmentation from pelvic CT images. Accurately segmenting CT pelvic images in an automatic manner is an important initial step for assessing the presence of fracture in the context of an assisted decision-making system that provides recommendations for trauma care. The system will integrate demographic information, physiological data, fracture detection, assesment of hemorrhage severity to provide caregivers with recommendations regarding patient diagnosis and treatment. Image analysis for hemorhage detection already provided promising results. A larger database has recently been made available for testing future work includes even more extensive testing of the method.
List of abbreviations used
Speckle Reducing Anisotropic Difusion
Deformable Model Methods
Level Set Methods
partial diferential equations
This material is based upon work supported by the National Science Foundation under Grant No. IIS0758410. The database used in this paper is courtesy of Carolinas Healthcare System. The authors would like to thank Ms. Rebecca Smith for her feedback.
This article has been published as part of BMC Medical Informatics and Decision Making Volume 9, Supplement 1, 2009: 2008 International Workshop on Biomedical and Health Informatics. The full contents of the supplement are available online at http://www.biomedcentral.com/1472-6947/9?issue=S1.
- Broderick JD, McKenna DP: Injury Control. Trauma Management: An Emergency Medicine Approach. Edited by: Ferrera PC, Colucciello SA, Marx J, Verdile V. 2001, St. Louis Missouri: C.V. MosbyGoogle Scholar
- CCarrico J, Holcomb JB, Chaundry IH: Scientific priorities and strategic planning for resuscitation research and Life saving therapy following traumatic injury: Report of the PULSE Trauma Work Group. Academic Emergency Medicine. 2002, 9 (6): 621-626.View ArticleGoogle Scholar
- of Maryland National Study Center for Trauma/EMS U: Lower extremity injuries among restrained vehicle occupants.
- Ben-Menachem Y, Coldwell D, Young JW, Burgess AR: Hemorrhage associated with pelvic fractures: causes, diagnosis and emergent management. American Journal of Roentgenology. 1991, 157 (5): 1005-1014.View ArticlePubMedGoogle Scholar
- Ji SY, Smith R, Huynh T, Najarian K: A comparative analysis of multi-level computer-assisted decision making systems for traumatic injuries. BMC Medical Informatics and Decision Making. 2009, 9 (Suppl 1): S6-PubMed CentralView ArticlePubMedGoogle Scholar
- Straka M, Cruz AL, Kochl A, Sramek M, Groller E, Fleischmann D: 3D watershed transform combined with a probabilistic atlas for medical image segmentation. Proceedings of MIT 2003. 2003, 1-8.Google Scholar
- Sebastian T, Tek H, Crisco J, Kimia B: Segmentation of Carpal Bones from CT Images Using Skeletally Coupled Deformable Models. Medical Image Analysis. 2003, 7: 21-45.View ArticlePubMedGoogle Scholar
- Souza A, Udupa JK, Saha PK: Volume Rendering in the Presence of Partial Volume Effects. IEEE TRANSACTIONS ON MEDICAL IMAGING. 2005, 24 (2): 223-233.View ArticlePubMedGoogle Scholar
- Kainmueller D, Lamecker H, Zachow S, Hege HC: Coupling Deformable Models for Multi-object Segmentation. Lecture Notes in Computer Science. Edited by: Bello F, Edwards E. 2008, Springer-Verlag Berlin Heidelberg, 69-78.Google Scholar
- Nyul LG, Kanyo J, Mate E, Geza Makay EB, Fidrich M, Kuba A: Method for Automatically Segmenting the Spinal Cord and Canal from 3D CT Images. Lecture Notes in Computer Science. Edited by: Gagalowicz A, Philips W. 2005, Springer-Verlag Berlin Heidelberg, 456-463.Google Scholar
- Burnett SSC, Starkschall G, Stevens CW, Liao Z: A deformable-model approach to semi-automatic segmentation of CT images demonstrated by application to the spinal canal. Medical Physics. 2004, 31 (2): 251-263.View ArticlePubMedGoogle Scholar
- Sylvia R, Antonio GJ, Raphael P, Mariano A: Automatic segmentation of jaw tissues in CT using active appearance models and semi-automatic landmarking. Proceedings of MICCAI - International Conference on Medical Image Computing and Computer-Assisted Intervention 2006. 2006, 167-174.Google Scholar
- Pham D, Xu C, Prince J: A Survey of Current Methods in Medical Image Segmentation. Annual Review of Biomedical Engineering. 2000, 2: 315-338.View ArticlePubMedGoogle Scholar
- Shojaii R, Alirezaie J, Babyn P: Automatic lung segmentation in ct images using watershed transform. Proceedings of IEEE International Conference on Image Processing 2005. 2005Google Scholar
- Sramek M, Dimitrov L: Segmentation of tomographic data by hierarchical watershed transform. Medical Informatics and Technologies. 2002, 3: MI161-MI169.Google Scholar
- Mancas M, Gosselin B: Towards an automatic tumor segmentation using iterative watersheds. Proceedings of SPIE. 2004, 5370: 1598-1608.View ArticleGoogle Scholar
- Peter Z, Boussone V, Bergote C, Peyrin F: Aconstrained region growing approach based onwatershed for the segmentation of lowcontrast structures in bone micro-CT images. Pattern Recognition. 2008, 41: 2358-2368.View ArticleGoogle Scholar
- Mancas M, Gosselin B: Segmentation using a region-growing thresholding. Proceedings of SPIE - Image Processing. 2005, 5672:Google Scholar
- Mehnert A, Jackway P: An improved seeded region growing algorithm. Pattern Recognition Letters. 1997, 18 (10): 1065-1071.View ArticleGoogle Scholar
- Adams R, Bischof L: Seeded Region Growing. IEEE Transactions on Pattern Analysis and Machine Intelligence. 1994, 16 (6): 641-647.View ArticleGoogle Scholar
- Zhu SC, Yuille A: Region competition: Unifying snakes, region growing and Bayes/MDL for multiband image segmentation. IEEE Transactions on Pattern Analysis and Machine Intelligence. 1996, 18 (9): 884-900.View ArticleGoogle Scholar
- Tek H, Akova F, Ayvaci A: Region competition via local watershed operators. IEEE CVPR Computer Vision and Pattern Recognition. 2005, 2: 361-368.Google Scholar
- Street E, Hadjiiski L, Sahiner B, Gujar S, Ibrahim M, Mukherji SK, Chan HP: Automated volume analysis of head and neck lesions on CT scans using 3D level set segmentation. Medical Physics. 2007, 34 (11): 4399-4408.PubMed CentralView ArticlePubMedGoogle Scholar
- Xie Q, Chen X, Ma L, Zhou Z: Segmentation for CT Image Based on Improved Level-Set Approach. CISP Congress on Image and Signal Processing 2008. 2008, 3: 725-728.View ArticleGoogle Scholar
- Hong-Kai Z, Osher S, Fedkiw R: Fast surface reconstruction using the level set method. Proceedings of IEEE Workshop on Variational and Level Set Methods in Computer Vision. 2001, 194-201.View ArticleGoogle Scholar
- Yonggang S, Karl W: A Real-Time Algorithm for the Approximation of Level-Set-Based Curve Evolution. IEEE Transactions on Image Processing. 2008, 17: 645-656.View ArticleGoogle Scholar
- Yongjian Y, Acton S: Speckle Reducing Anisotropic Diffusion. IEEE Transactions on Image Processing. 2002, 11 (11): 1260-1270.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.