Computer-aided segmentation and estimation of indices in brain CT scans

The importance of neuro-imaging as one of the biomarkers for diagnosis and prognosis of pathologies and traumatic cases is well established. Doctors routinely perform linear measurements on neuro-images to ascertain severity and extent of the pathology or trauma from significant anatomical changes. However, it is a tedious and time consuming process and manually assessing and reporting on large volume of data is fraught with errors and variation. In this paper we present a novel technique for segmentation of significant anatomical landmarks using artificial neural networks and estimation of various ratios and indices performed on brain CT scans. The proposed method is efficient and robust in detecting and measuring sizes of anatomical structures on non-contrast CT scans and has been evaluated on images from subjects with ages between 5 to 85 years. Results show that our method has average ICC of ≥0.97 and, hence, can be used in processing data for further use in research and clinical environment

CT diagnosis of early stroke : the initial approach to the new CAD tool based on multiscale estimation of ischemia

Background: Computer aided diagnosis (CAD) becomes one of the most important diagnostic tools for urgent states in cerebral stroke and other life-threatening conditions where time plays a crucial role. Routine CT is still diagnostically insufficient in hyperacute stage of stroke that is in the therapeutic window for thrombolytic therapy. Authors present computer assistant of early ischemic stroke diagnosis that supports the radiologic interpretations. A new semantic-visualization system of ischemic symptoms applied to noncontrast, routine CT examination was based on multiscale image processing and diagnostic content estimation. Material/Methods: Evaluation of 95 sets of examinations in patients admitted to a hospital with symptoms suggesting stroke was undertaken by four radiologists from two medical centers unaware of the final clinical findings. All of the consecutive cases were considered as having no CT direct signs of hyperacute ischemia. At the first test stage only the CTs performed at the admission were evaluated independently by radiologists. Next, the same early scans were evaluated again with additional use of multiscale computer-assistant of stroke (MulCAS). Computerized suggestion with increased sensitivity to the subtle image manifestations of cerebral ischemia was constructed as additional view representing estimated diagnostic content with enhanced stroke symptoms synchronized to routine CT data preview. Follow-up CT examinations and clinical features confirmed or excluded the diagnosis of stroke constituting 'gold standard' to verify stroke detection performance. Results: Higher AUC (area under curve) values were found for MulCAS -aided radiological diagnosis for all readers and the differences were statistically significant for random readers-random cases parametric and non-parametric DBM MRMC analysis. Sensitivity and specificity of acute stroke detection for the readers was increased by 30% and 4%, respectively. Conclusions: Routine CT completed with proposed method of computer assisted diagnosis provided noticeable better diagnosis efficiency of acute stroke according to the rates and opinions of all test readers. Further research includes fully automatic detection of hypodense regions to complete assisted indications and formulate the suggestions of stroke cases more objectively. Planned prospective studies will let evaluate more accurately the impact of this CAD tool on diagnosis and further treatment in patients suffered from stroke. It is necessary to determine whether this method is possible to be applied widely

Computer aided assessment of CT scans of traumatic brain injury patients

A thesis submitted in partial fulfilment for the degree of Doctor of PhilosophyOne of the serious public health problems is the Traumatic Brain Injury, also known as silent epidemic, affecting millions every year. Management of these patients essentially involves neuroimaging and noncontrast CT scans are the first choice amongst doctors. Significant anatomical changes identified on the neuroimages and volumetric assessment of haemorrhages and haematomas are of critical importance for assessing the patients’ condition for targeted therapeutic and/or surgical interventions. Manual demarcation and annotation by experts is still considered gold standard, however, the interpretation of neuroimages is fraught with inter-observer variability and is considered ’Achilles heel’ amongst radiologists. Errors and variability can be attributed to factors such as poor perception, inaccurate deduction, incomplete knowledge or the quality of the image and only a third of doctors confidently report the findings. The applicability of computer aided dianosis in segmenting the apposite regions and giving ’second opinion’ has been positively appraised to assist the radiologists, however, results of the approaches vary due to parameters of algorithms and manual intervention required from doctors and this presents a gap for automated segmentation and estimation of measurements of noncontrast brain CT scans. The Pattern Driven, Content Aware Active Contours (PDCAAC) Framework developed in this thesis provides robust and efficient segmentation of significant anatomical landmarks, estimations of their sizes and correlation to CT rating to assist the radiologists in establishing the diagnosis and prognosis more confidently. The integration of clinical profile of the patient into image segmentation algorithms has significantly improved their performance by highlighting characteristics of the region of interest. The modified active contour method in the PDCAAC framework achieves Jaccard Similarity Index (JI) of 0.87, which is a significant improvement over the existing methods of active contours achieving JI of 0.807 with Simple Linear Iterative Clustering and Distance Regularized Level Set Evolution. The Intraclass Correlation Coefficient of intracranial measurements is >0.97 compared with radiologists. Automatic seeding of the initial seed curve within the region of interest is incorporated into the method which is a novel approach and alleviates limitation of existing methods. The proposed PDCAAC framework can be construed as a contribution towards research to formulate correlations between image features and clinical variables encompassing normal development, ageing, pathological and traumatic cases propitious to improve management of such patients. Establishing prognosis usually entails survival but the focus can also be extended to functional outcomes, residual disability and quality of life issues

A Rapid Segmentation-Insensitive "Digital Biopsy" Method for Radiomic Feature Extraction: Method and Pilot Study Using CT Images of Non-Small Cell Lung Cancer.

Quantitative imaging approaches compute features within images' regions of interest. Segmentation is rarely completely automatic, requiring time-consuming editing by experts. We propose a new paradigm, called "digital biopsy," that allows for the collection of intensity- and texture-based features from these regions at least 1 order of magnitude faster than the current manual or semiautomated methods. A radiologist reviewed automated segmentations of lung nodules from 100 preoperative volume computed tomography scans of patients with non-small cell lung cancer, and manually adjusted the nodule boundaries in each section, to be used as a reference standard, requiring up to 45 minutes per nodule. We also asked a different expert to generate a digital biopsy for each patient using a paintbrush tool to paint a contiguous region of each tumor over multiple cross-sections, a procedure that required an average of <3 minutes per nodule. We simulated additional digital biopsies using morphological procedures. Finally, we compared the features extracted from these digital biopsies with our reference standard using intraclass correlation coefficient (ICC) to characterize robustness. Comparing the reference standard segmentations to our digital biopsies, we found that 84/94 features had an ICC >0.7; comparing erosions and dilations, using a sphere of 1.5-mm radius, of our digital biopsies to the reference standard segmentations resulted in 41/94 and 53/94 features, respectively, with ICCs >0.7. We conclude that many intensity- and texture-based features remain consistent between the reference standard and our method while substantially reducing the amount of operator time required

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

Automatic detection and classification of nasopharyngeal carcinoma on PET/CT with support vector machine

Purpose: Positron emission tomography/computed tomography (PET/CT) has established values for imaging of head and neck cancers, including the nasopharyngeal carcinoma (NPC), utilizing both morphologic and functional information. In this paper, we introduce a computerized system for automatic detection of NPC, targeting both the primary tumor and regional nodal metastasis, on PET/CT. Methods: Candidate lesions were extracted based on the features from both PET and CT images and a priori knowledge of anatomical features and subsequently classified by a support vector machine algorithm. The system was validated with 25 PET/CT examinations from 10 patients suffering from NPC. Lesions manually contoured by experienced radiologists were used as the gold standard. Results: Results showed that the system successfully identified all 53 hypermetabolic lesions larger than 1 cm in size and excluded normal physiological uptake in brown fat, muscles, bone marrow, brain, and salivary glands. Conclusion: The system combined both imaging features and a priori clinical knowledge for classification between pathological and physiological uptake. Preliminary results showed that the system was highly accurate and promising for adoption in clinical use. © The Author(s) 2011.published_or_final_versionSpringer Open Choice, 25 May 201

Three Dimensional Nonlinear Statistical Modeling Framework for Morphological Analysis

This dissertation describes a novel three-dimensional (3D) morphometric analysis framework for building statistical shape models and identifying shape differences between populations. This research generalizes the use of anatomical atlases on more complex anatomy as in case of irregular, flat bones, and bones with deformity and irregular bone growth. The foundations for this framework are: 1) Anatomical atlases which allow the creation of homologues anatomical models across populations; 2) Statistical representation for output models in a compact form to capture both local and global shape variation across populations; 3) Shape Analysis using automated 3D landmarking and surface matching. The proposed framework has various applications in clinical, forensic and physical anthropology fields. Extensive research has been published in peer-reviewed image processing, forensic anthropology, physical anthropology, biomedical engineering, and clinical orthopedics conferences and journals. The forthcoming discussion of existing methods for morphometric analysis, including manual and semi-automatic methods, addresses the need for automation of morphometric analysis and statistical atlases. Explanations of these existing methods for the construction of statistical shape models, including benefits and limitations of each method, provide evidence of the necessity for such a novel algorithm. A novel approach was taken to achieve accurate point correspondence in case of irregular and deformed anatomy. This was achieved using a scale space approach to detect prominent scale invariant features. These features were then matched and registered using a novel multi-scale method, utilizing both coordinate data as well as shape descriptors, followed by an overall surface deformation using a new constrained free-form deformation. Applications of output statistical atlases are discussed, including forensic applications for the skull sexing, as well as physical anthropology applications, such as asymmetry in clavicles. Clinical applications in pelvis reconstruction and studying of lumbar kinematics and studying thickness of bone and soft tissue are also discussed