Liver segmentation using 3D CT scans.

Master of Science in Computer Science. University of KwaZulu-Natal, Durban, 2018.Abstract available in PDF file

Computational methods for the analysis of functional 4D-CT chest images.

Medical imaging is an important emerging technology that has been intensively used in the last few decades for disease diagnosis and monitoring as well as for the assessment of treatment effectiveness. Medical images provide a very large amount of valuable information that is too huge to be exploited by radiologists and physicians. Therefore, the design of computer-aided diagnostic (CAD) system, which can be used as an assistive tool for the medical community, is of a great importance. This dissertation deals with the development of a complete CAD system for lung cancer patients, which remains the leading cause of cancer-related death in the USA. In 2014, there were approximately 224,210 new cases of lung cancer and 159,260 related deaths. The process begins with the detection of lung cancer which is detected through the diagnosis of lung nodules (a manifestation of lung cancer). These nodules are approximately spherical regions of primarily high density tissue that are visible in computed tomography (CT) images of the lung. The treatment of these lung cancer nodules is complex, nearly 70% of lung cancer patients require radiation therapy as part of their treatment. Radiation-induced lung injury is a limiting toxicity that may decrease cure rates and increase morbidity and mortality treatment. By finding ways to accurately detect, at early stage, and hence prevent lung injury, it will have significant positive consequences for lung cancer patients. The ultimate goal of this dissertation is to develop a clinically usable CAD system that can improve the sensitivity and specificity of early detection of radiation-induced lung injury based on the hypotheses that radiated lung tissues may get affected and suffer decrease of their functionality as a side effect of radiation therapy treatment. These hypotheses have been validated by demonstrating that automatic segmentation of the lung regions and registration of consecutive respiratory phases to estimate their elasticity, ventilation, and texture features to provide discriminatory descriptors that can be used for early detection of radiation-induced lung injury. The proposed methodologies will lead to novel indexes for distinguishing normal/healthy and injured lung tissues in clinical decision-making. To achieve this goal, a CAD system for accurate detection of radiation-induced lung injury that requires three basic components has been developed. These components are the lung fields segmentation, lung registration, and features extraction and tissue classification. This dissertation starts with an exploration of the available medical imaging modalities to present the importance of medical imaging in today’s clinical applications. Secondly, the methodologies, challenges, and limitations of recent CAD systems for lung cancer detection are covered. This is followed by introducing an accurate segmentation methodology of the lung parenchyma with the focus of pathological lungs to extract the volume of interest (VOI) to be analyzed for potential existence of lung injuries stemmed from the radiation therapy. After the segmentation of the VOI, a lung registration framework is introduced to perform a crucial and important step that ensures the co-alignment of the intra-patient scans. This step eliminates the effects of orientation differences, motion, breathing, heart beats, and differences in scanning parameters to be able to accurately extract the functionality features for the lung fields. The developed registration framework also helps in the evaluation and gated control of the radiotherapy through the motion estimation analysis before and after the therapy dose. Finally, the radiation-induced lung injury is introduced, which combines the previous two medical image processing and analysis steps with the features estimation and classification step. This framework estimates and combines both texture and functional features. The texture features are modeled using the novel 7th-order Markov Gibbs random field (MGRF) model that has the ability to accurately models the texture of healthy and injured lung tissues through simultaneously accounting for both vertical and horizontal relative dependencies between voxel-wise signals. While the functionality features calculations are based on the calculated deformation fields, obtained from the 4D-CT lung registration, that maps lung voxels between successive CT scans in the respiratory cycle. These functionality features describe the ventilation, the air flow rate, of the lung tissues using the Jacobian of the deformation field and the tissues’ elasticity using the strain components calculated from the gradient of the deformation field. Finally, these features are combined in the classification model to detect the injured parts of the lung at an early stage and enables an earlier intervention

Segmentation of Juxtapleural Pulmonary Nodules Using a Robust Surface Estimate

An algorithm was developed to segment solid pulmonary nodules attached to the chest wall in computed tomography scans. The pleural surface was estimated and used to segment the nodule from the chest wall. To estimate the surface, a robust approach was used to identify points that lie on the pleural surface but not on the nodule. A 3D surface was estimated from the identified surface points. The segmentation performance of the algorithm was evaluated on a database of 150 solid juxtapleural pulmonary nodules. Segmented images were rated on a scale of 1 to 4 based on visual inspection, with 3 and 4 considered acceptable. This algorithm offers a large improvement in the success rate of juxtapleural nodule segmentation, successfully segmenting 98.0% of nodules compared to 81.3% for a previously published plane-fitting algorithm, which will provide for the development of more robust automated nodule measurement methods

CAD system for lung nodule analysis.

Lung cancer is the deadliest type of known cancer in the United States, claiming hundreds of thousands of lives each year. However, despite the high mortality rate, the 5-year survival rate after resection of Stage 1A non–small cell lung cancer is currently in the range of 62%– 82% and in recent studies even 90%. Patient survival is highly correlated with early detection. Computed Tomography (CT) technology services the early detection of lung cancer tremendously by offering a minimally invasive medical diagnostic tool. Some early types of lung cancer begin with a small mass of tissue within the lung, less than 3 cm in diameter, called a nodule. Most nodules found in a lung are benign, but a small population of them becomes malignant over time. Expert analysis of CT scans is the first step in determining whether a nodule presents a possibility for malignancy but, due to such low spatial support, many potentially harmful nodules go undetected until other symptoms motivate a more thorough search. Computer Vision and Pattern Recognition techniques can play a significant role in aiding the process of detecting and diagnosing lung nodules. This thesis outlines the development of a CAD system which, given an input CT scan, provides a functional and fast, second-opinion diagnosis to physicians. The entire process of lung nodule screening has been cast as a system, which can be enhanced by modern computing technology, with the hopes of providing a feasible diagnostic tool for clinical use. It should be noted that the proposed CAD system is presented as a tool for experts—not a replacement for them. The primary motivation of this thesis is the design of a system that could act as a catalyst for reducing the mortality rate associated with lung cancer

Characterisation and correction of respiratory-motion artefacts in cardiac PET-CT

Respiratory motion during cardiac Positron Emission Tomography (PET) Computed Tomography (CT) imaging results in blurring of the PET data and can induce mismatches between the PET and CT datasets, leading to attenuation-correction artefacts. The aim of this project was to develop a method of motion-correction to overcome both of these problems. The approach implemented was to transform a single CT to match the frames of a gated PET study, to facilitate respiratory-matched attenuation-correction, without the need for a gated CT. This is benecial for lowering the radiation dose to the patient and in reducing PETCT mismatches, which can arise even in gated studies. The heart and diaphragm were identied through phantom studies as the structures responsible for generating attenuation-correction artefacts in the heart and their motions therefore needed to be considered in transforming the CT. Estimating heart motion was straight-forward, due to its high contrast in PET, however the poor diaphragm contrast meant that additional information was required to track its position. Therefore a diaphragm shape model was constructed using segmented diaphragm surfaces, enabling complete diaphragm surfaces to be produced from incomplete and noisy initial estimates. These complete surfaces, in combination with the estimated heart motions were used to transform the CT. The PET frames were then attenuation-corrected with the transformed CT, reconstructed, aligned and summed, to produce motion-free images. It was found that motion-blurring was reduced through alignment, although benets were marginal in the presence of small respiratory motions. Quantitative accuracy was improved from use of the transformed CT for attenuation-correction (compared with no CT transformation), which was attributed to both the heart and the diaphragm transformations. In comparison to a gated CT, a substantial dose saving and a reduced dependence on gating techniques were achieved, indicating the potential value of the technique in routine clinical procedures

Automatic Approach for Lung Segmentation with Juxta-Pleural Nodules from Thoracic CT Based on Contour Tracing and Correction

This paper presents a fully automatic framework for lung segmentation, in which juxta-pleural nodule problem is brought into strong focus. The proposed scheme consists of three phases: skin boundary detection, rough segmentation of lung contour, and pulmonary parenchyma refinement. Firstly, chest skin boundary is extracted through image aligning, morphology operation, and connective region analysis. Secondly, diagonal-based border tracing is implemented for lung contour segmentation, with maximum cost path algorithm used for separating the left and right lungs. Finally, by arc-based border smoothing and concave-based border correction, the refined pulmonary parenchyma is obtained. The proposed scheme is evaluated on 45 volumes of chest scans, with volume difference (VD) 11.15±69.63 cm3, volume overlap error (VOE) 3.5057±1.3719%, average surface distance (ASD) 0.7917±0.2741 mm, root mean square distance (RMSD) 1.6957±0.6568 mm, maximum symmetric absolute surface distance (MSD) 21.3430±8.1743 mm, and average time-cost 2 seconds per image. The preliminary results on accuracy and complexity prove that our scheme is a promising tool for lung segmentation with juxta-pleural nodules

Quantifying Airway Dilatation in the Lungs from Computed Tomography

Non CF bronchiectasis and idiopathic pulmonary fibrosis (IPF) are pulmonary diseases characterised by the abnormal and permanent dilatation of the airways. Computed tomography (CT) is used in clinical practice to diagnose and monitor patients with the disease. Currently, analysis of the scans is performed by manual inspection and there is no established computerised method to quantify the enlargement of airways. I developed a pipeline to quantify the cross-sectional area for a given airway track. Using an airway segmentation, my proposed algorithm measures the area at contiguous intervals along the airway arclength from the Carina to the most distal point visible on CT. I showed the use of the data generated from the pipeline in two applications. First, I proposed a novel tapering measure as the gradient of a linear regression between a logarithmic area against the arclength. The measurement was applied to airways affected by bronchiectasis. Second, I used Bayesian Changepoint Detection (BCD) with the area measurements to locate the progression of IPF along the airway track. The proposed pipeline was applied to a set of clinically acquired scans. I show a statistical difference (p = 3.4×10−4 ) in the tapering measurement between bronchiectatic (n = 53) and controlled (n = 39) airways. In addition, I report a statistical difference (p = 7.2×10−3 ) in the change in measurement between airways remaining healthy (n = 14) and airways that have become bronchiectatic (n = 5). I show the tapering measurement is reproducible independent to voxel size, CT reconstruction, and radiation dose. Using BCD, I show on simulated data (n = 14) my proposed method can detect the progression of IPF within 2.5mm. Finally, using results from BCD, I present a novel measure of IPF progression as the percentage volume change in the diseased region of the airways