Segmentation of Pulmonary Vascular Trees from Thoracic 3D CT Images

This paper describes an algorithm for extracting pulmonary vascular trees (arteries plus veins) from three-dimensional (3D) thoracic computed tomographic (CT) images. The algorithm integrates tube enhancement filter and traversal approaches which are based on eigenvalues and eigenvectors of a Hessian matrix to extract thin peripheral segments as well as thick vessels close to the lung hilum. The resultant algorithm was applied to a simulation data set and 44 scans from 22 human subjects imaged via multidetector-row CT (MDCT) during breath holds at 85% and 20% of their vital capacity. A quantitative validation was performed with more than 1000 manually identified points selected from inside the vessel segments to assess true positives (TPs) and 1000 points randomly placed outside of the vessels to evaluate false positives (FPs) in each case. On average, for both the high and low volume lung images, 99% of the points was properly marked as vessel and 1% of the points were assessed as FPs. Our hybrid segmentation algorithm provides a highly reliable method of segmenting the combined pulmonary venous and arterial trees which in turn will serve as a critical starting point for further quantitative analysis tasks and aid in our overall goal of establishing a normative atlas of the human lung

Hemorrhage Detection and Analysis in Traumatic Pelvic Injuries

Traumatic pelvic injuries associated with high-energy pelvic fractures are life-threatening injuries. Extensive bleeding is relatively common with pelvic fractures. However, bleeding is especially prevalent with high-energy fractures. Hemorrhage remains the major cause of death that occur within the first 24 hours after a traumatic pelvic injury. Emergent-life saving treatment is required for high-energy pelvic fractures associated with hemorrhage. A thorough understanding of potential sources of bleeding within a short period is essential for diagnosis and treatment planning. Computed Tomography (CT) images have been widely in use in identifying the potential sources of bleeding. A pelvic CT scan contains a large number of images. Analyzing each slice in a scan via simple visual inspection is very time consuming. Time is a crucial factor in emergency medicine. Therefore, a computer-assisted pelvic trauma decision-making system is advantageous for assisting physicians in fast and accurate decision making and treatment planning. The proposed project presents an automated system to detect and segment hemorrhage and combines it with the other extracted features from pelvic images and demographic data to provide recommendations to trauma caregivers for diagnosis and treatment. The first part of the project is to develop automated methods to detect arteries by incorporating bone information. This part of the project merges bone edges and segments bone using a seed growing technique. Later the segmented bone information is utilized along with the best template matching to locate arteries and extract gray level information of the located arteries in the pelvic region. The second part of the project focuses on locating the source of hemorrhage and its segmentation. The hemorrhage is segmented using a novel rule based hemorrhage segmentation approach. This approach segments hemorrhage through hemorrhage matching, rule optimization, and region growing. Later the position of hemorrhage in the image and the volume of the hemorrhage are determined to analyze hemorrhage severity. The third part of the project is to automatically classify the outcome using features extracted from the medical images and patient medical records and demographics. A multi-stage feature selection algorithm is used to select the predominant features among all the features. Finally, boosted logistic model tree is used to classify the outcome. The methods are tested on CT images of traumatic pelvic injury patients. The hemorrhage segmentation and classification results seem promising and demonstrate that the proposed method is not only capable of automatically segmenting hemorrhage and classifying outcome, but also has the potential to be used for clinical applications. Finally, the project is extended to abdominal trauma and a novel knowledge based heuristic technique is used to detect and segment spleen from the abdominal CT images. This technique is tested on a limited number of subjects and the results are promising

Lung nodule modeling and detection for computerized image analysis of low dose CT imaging of the chest.

From a computerized image analysis prospective, early diagnosis of lung cancer involves detection of doubtful nodules and classification into different pathologies. The detection stage involves a detection approach, usually by template matching, and an authentication step to reduce false positives, usually conducted by a classifier of one form or another; statistical, fuzzy logic, support vector machines approaches have been tried. The classification stage matches, according to a particular approach, the characteristics (e.g., shape, texture and spatial distribution) of the detected nodules to common characteristics (again, shape, texture and spatial distribution) of nodules with known pathologies (confirmed by biopsies). This thesis focuses on the first step; i.e., nodule detection. Specifically, the thesis addresses three issues: a) understanding the CT data of typical low dose CT (LDCT) scanning of the chest, and devising an image processing approach to reduce the inherent artifacts in the scans; b) devising an image segmentation approach to isolate the lung tissues from the rest of the chest and thoracic regions in the CT scans; and c) devising a nodule modeling methodology to enhance the detection rate and lend benefits for the ultimate step in computerized image analysis of LDCT of the lungs, namely associating a pathology to the detected nodule. The methodology for reducing the noise artifacts is based on noise analysis and examination of typical LDCT scans that may be gathered on a repetitive fashion; since, a reduction in the resolution is inevitable to avoid excessive radiation. Two optimal filtering methods are tested on samples of the ELCAP screening data; the Weiner and the Anisotropic Diffusion Filters. Preference is given to the Anisotropic Diffusion Filter, which can be implemented on 7x7 blocks/windows of the CT data. The methodology for lung segmentation is based on the inherent characteristics of the LDCT scans, shown as distinct bi-modal gray scale histogram. A linear model is used to describe the histogram (the joint probability density function of the lungs and non-lungs tissues) by a linear combination of weighted kernels. The Gaussian kernels were chosen, and the classic Expectation-Maximization (EM) algorithm was employed to estimate the marginal probability densities of the lungs and non-lungs tissues, and select an optimal segmentation threshold. The segmentation is further enhanced using standard shape analysis based on mathematical morphology, which improves the continuity of the outer and inner borders of the lung tissues. This approach (a preliminary version of it appeared in [14]) is found to be adequate for lung segmentation as compared to more sophisticated approaches developed at the CVIP Lab (e.g., [15][16]) and elsewhere. The methodology developed for nodule modeling is based on understanding the physical characteristics of the nodules in LDCT scans, as identified by human experts. An empirical model is introduced for the probability density of the image intensity (or Hounsfield units) versus the radial distance measured from the centroid – center of mass - of typical nodules. This probability density showed that the nodule spatial support is within a circle/square of size 10 pixels; i.e., limited to 5 mm in length; which is within the range that the radiologist specify to be of concern. This probability density is used to fill in the intensity (or Hounsfield units) of parametric nodule models. For these models (e.g., circles or semi-circles), given a certain radius, we calculate the intensity (or Hounsfield units) using an exponential expression for the radial distance with parameters specified from the histogram of an ensemble of typical nodules. This work is similar in spirit to the earlier work of Farag et al., 2004 and 2005 [18][19], except that the empirical density of the radial distance and the histogram of typical nodules provide a data-driven guide for estimating the intensity (or Hounsfield units) of the nodule models. We examined the sensitivity and specificity of parametric nodules in a template-matching framework for nodule detection. We show that false positives are inevitable problems with typical machine learning methods of automatic lung nodule detection, which invites further efforts and perhaps fresh thinking into automatic nodule detection. A new approach for nodule modeling is introduced in Chapter 5 of this thesis, which brings high promise in both the detection, and the classification of nodules. Using the ELCAP study, we created an ensemble of four types of nodules and generated a nodule model for each type based on optimal data reduction methods. The resulting nodule model, for each type, has lead to drastic improvements in the sensitivity and specificity of nodule detection. This approach may be used as well for classification. In conclusion, the methodologies in this thesis are based on understanding the LDCT scans and what is to be expected in terms of image quality. Noise reduction and image segmentation are standard. The thesis illustrates that proper nodule models are possible and indeed a computerized approach for image analysis to detect and classify lung nodules is feasible. Extensions to the results in this thesis are immediate and the CVIP Lab has devised plans to pursue subsequent steps using clinical data