Foetal echocardiographic segmentation

Congenital heart disease affects just under one percentage of all live births [1]. Those defects that manifest themselves as changes to the cardiac chamber volumes are the motivation for the research presented in this thesis. Blood volume measurements in vivo require delineation of the cardiac chambers and manual tracing of foetal cardiac chambers is very time consuming and operator dependent. This thesis presents a multi region based level set snake deformable model applied in both 2D and 3D which can automatically adapt to some extent towards ultrasound noise such as attenuation, speckle and partial occlusion artefacts. The algorithm presented is named Mumford Shah Sarti Collision Detection (MSSCD). The level set methods presented in this thesis have an optional shape prior term for constraining the segmentation by a template registered to the image in the presence of shadowing and heavy noise. When applied to real data in the absence of the template the MSSCD algorithm is initialised from seed primitives placed at the centre of each cardiac chamber. The voxel statistics inside the chamber is determined before evolution. The MSSCD stops at open boundaries between two chambers as the two approaching level set fronts meet. This has significance when determining volumes for all cardiac compartments since cardiac indices assume that each chamber is treated in isolation. Comparison of the segmentation results from the implemented snakes including a previous level set method in the foetal cardiac literature show that in both 2D and 3D on both real and synthetic data, the MSSCD formulation is better suited to these types of data. All the algorithms tested in this thesis are within 2mm error to manually traced segmentation of the foetal cardiac datasets. This corresponds to less than 10% of the length of a foetal heart. In addition to comparison with manual tracings all the amorphous deformable model segmentations in this thesis are validated using a physical phantom. The volume estimation of the phantom by the MSSCD segmentation is to within 13% of the physically determined volume

Rapid Segmentation Techniques for Cardiac and Neuroimage Analysis

Recent technological advances in medical imaging have allowed for the quick acquisition of highly resolved data to aid in diagnosis and characterization of diseases or to guide interventions. In order to to be integrated into a clinical work flow, accurate and robust methods of analysis must be developed which manage this increase in data. Recent improvements in in- expensive commercially available graphics hardware and General-Purpose Programming on Graphics Processing Units (GPGPU) have allowed for many large scale data analysis problems to be addressed in meaningful time and will continue to as parallel computing technology improves. In this thesis we propose methods to tackle two clinically relevant image segmentation problems: a user-guided segmentation of myocardial scar from Late-Enhancement Magnetic Resonance Images (LE-MRI) and a multi-atlas segmentation pipeline to automatically segment and partition brain tissue from multi-channel MRI. Both methods are based on recent advances in computer vision, in particular max-flow optimization that aims at solving the segmentation problem in continuous space. This allows for (approximately) globally optimal solvers to be employed in multi-region segmentation problems, without the particular drawbacks of their discrete counterparts, graph cuts, which typically present with metrication artefacts. Max-flow solvers are generally able to produce robust results, but are known for being computationally expensive, especially with large datasets, such as volume images. Additionally, we propose two new deformable registration methods based on Gauss-Newton optimization and smooth the resulting deformation fields via total-variation regularization to guarantee the problem is mathematically well-posed. We compare the performance of these two methods against four highly ranked and well-known deformable registration methods on four publicly available databases and are able to demonstrate a highly accurate performance with low run times. The best performing variant is subsequently used in a multi-atlas segmentation pipeline for the segmentation of brain tissue and facilitates fast run times for this computationally expensive approach. All proposed methods are implemented using GPGPU for a substantial increase in computational performance and so facilitate deployment into clinical work flows. We evaluate all proposed algorithms in terms of run times, accuracy, repeatability and errors arising from user interactions and we demonstrate that these methods are able to outperform established methods. The presented approaches demonstrate high performance in comparison with established methods in terms of accuracy and repeatability while largely reducing run times due to the employment of GPU hardware

Registration and Deformable Model-Based Neck Muscles Segmentation and 3D Reconstruction

Whiplash is a very common ailment encountered in clinical practice that is usually a result of vehicle accidents but also domestic activities and sports injuries. It is normally caused when neck organs (specifically muscles) are impaired. Whiplash-associated disorders include acute headaches, neck pain, stiffness, arm dislocation, abnormal sensations, and auditory and optic problems, the persistence of which may be chronic or acute. Insurance companies compensate almost fifty percent of claims lodged due to whiplash injury through compulsory third party motor insurance. The morphological structures of neck muscles undergo hypertrophy or atrophy following damage caused to them by accidents. Before any medical treatment is applied , any such change needs to be known which requires 3D visualization of the neck muscles through a proper segmentation of them because the neck contains many other sensitive organs such as nerves, blood vessels, the spinal cord and trachea. The segmentation of neck muscles in medical images is a more challenging task than those of other muscles and organs due to their similar densities and compactness, low resolutions and contrast in medical images, anatomical variabilities among individuals, noise, inhomogeneity of medical images and false boundaries created by intra-muscular fat. Traditional segmentation algorithms, such as those used in thresholding and clustering-based methods, are not applicable in this project and also not suitable for medical images. Although there are some techniques available in clinical research for segmenting muscles such as thigh, tongue, leg, hip and pectoral ones, to the best of author's knowledge, there are no methods available for segmenting neck muscles due to the challenges described above. In the first part of this dissertation, an atlas-based method for segmenting MR images, which uses linear and non-linear registration frameworks, is proposed, with output from the registration process further refined by a novel parametric deformable model. The proposed method is tested on real clinical data of both healthy and non-healthy individuals. During the last few decades, registration- and deformable model-based segmentation methods have been very popular for medical image segmentation due to their incorporation of prior information. While registration-based segmentation techniques can preserve topologies of objects in an image, accuracy of atlas-based segmentation depends mainly on an effective registration process. In this study, the registration framework is designed in a novel way in which images are initially registered by a distinct 3D affine transformation and then aligned by a local elastic geometrical transformation based on discrete cosines and registered firstly slice-wise and then block-wise. The numbers of motion parameters are changed in three different steps per frame. This proposed registration framework can handle anatomical variabilities and pathologies by confining its parameters in local regions. Also, as warping of the framework relies on number of motion parameters, similarities between two images, gradients of floating image and coordinate mesh grid values, it can easily manage pathological and anatomical variabilities using a hierarchical parameter scheme. The labels transferred from atlas can be improved by deformable model-based segmentation. Although geometric deformable models have been widely used in many biomedical applications over recent years, they cannot work in the context of neck muscles segmentation due to noise, background clutter and similar objects touching each other. Another important drawback of geometric deformable models is that they are many times slower than parametric deformable ones. Therefore, the segmentation results produced by the registration process are ameliorated using a multiple-object parametric deformable model which is discussed in detail in the second part of this thesis. This algorithm uses a novel Gaussian potential energy distribution which can adapt to topological changes and does not require re-parameterization. Also, it incorporates a new overlap removal technique which ensures that there are no overlaps or gaps inside an object. Furthermore, stopping criteria of vertices are designed so that difference between boundaries of the deformable model and actual object is minimal. The multiple-object parametric deformable model is also applied in a template contours propagation-based segmentation technique, as discussed in the third part of this dissertation. This method is semi-automatic, whereby a manual delineation of middle image in a MRI data set is required. It can handle anatomical variabilities more easily than atlas-based segmentation because it can segment any individual's data irrespective of his/her age, weight and height with low computational complexity and it does not depend on other data as it operates semi-automatically. In it, initial model contour resides close to the object's boundary, with degree of closeness dependent on slice thicknesses and gaps between the slices

Automated Analysis of 3D Stress Echocardiography

__Abstract__ The human circulatory system consists of the heart, blood, arteries, veins and capillaries. The heart is the muscular organ which pumps the blood through the human body (Fig. 1.1,1.2). Deoxygenated blood flows through the right atrium into the right ventricle, which pumps the blood into the pulmonary arteries. The blood is carried to the lungs, where it passes through a capillary network that enables the release of carbon dioxide and the uptake of oxygen. Oxygenated blood then returns to the heart via the pulmonary veins and flows from the left atrium into the left ventricle. The left ventricle then pumps the blood through the aorta, the major artery which supplies blood to the rest of the body [Drake et a!., 2005; Guyton and Halt 1996]. Therefore, it is vital that the cardiovascular system remains healthy. Disease of the cardiovascular system, if untreated, ultimately leads to the failure of other organs and death

Medical image analysis methods for anatomical surface reconstruction using tracked 3D ultrasound

The thesis focuses on a study of techniques for acquisition and reconstruction of surface data from anatomical objects by means of tracked 3D ultrasound. In the context of the work two experimental scanning systems are developed and tested on both artificial objects and biological tissues. The first system is based on the freehand ultrasound principle and utilizes a conventional 2D ultrasound transducer coupled with an electromechanical 3D position tracker. The main properties and the basic features of this system are discussed. A number of experiments show that its accuracy in the close to ideal conditions reaches 1.2 mm RMS. The second proposed system implements the sequential triggered scanning approach. The system consists of an ultrasound machine, a workstation and a scanning body (a moving tank filled with liquid and a transducer fixation block) that performs transducer positioning and tracking functions. The system is tested on artificial and real bones. The performed experiments illustrate that it provides significantly better accuracy than the freehand ultrasound (about 0.2 mm RMS) and allows acquiring regular data with a good precision. This makes such a system a promising tool for orthopaedic and trauma surgeons during contactless X-ray-free examinations of injured extremities. The second major subject of the thesis concerns development of medical image analysis methods for 3D surface reconstruction and 2D object detection. We introduce a method based on mesh-growing surface reconstruction that is designed for noisy and sparse data received from 3D tracked ultrasound scanners. A series of experiments on synthetic and ultrasound data show an appropriate reconstruction accuracy. The reconstruction error is measured as the averaged distance between the faces of the mesh and the points from the cloud. Dependently on the initial settings of the method the error varies in range 0.04 - 0.2% for artificial data and 0.3 - 0.7 mm for ultrasound bone data. The reconstructed surfaces correctly interpolate the original point clouds and demonstrate proper smoothness. The next significant problem considered in the work is 2D object detection. Although medical object detection is not integrated into the developed scanning systems, it can be used as a possible further extension of the systems for automatic detection of specific anatomical structures. We analyse the existent object detection methods and introduce a modification of the one based on the popular Generalized Hough Transform (GHT). Unlike the original GHT, the developed method is invariant to rotation and uniform scaling, and uses an intuitive two-point parametrization. We propose several implementations of the feature-to-vote conversion function with the corresponding vote analysis principles. Special attention is devoted to a study of the hierarchical vote analysis and its probabilistic properties. We introduce a parameter space subdivision strategy that reduces the probability of vote peak omission, and show that it can be efficiently implemented in practice using the Gumbel probability distribution

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

Machine learning-based automated segmentation with a feedback loop for 3D synchrotron micro-CT

Die Entwicklung von Synchrotronlichtquellen der dritten Generation hat die Grundlage für die Untersuchung der 3D-Struktur opaker Proben mit einer Auflösung im Mikrometerbereich und höher geschaffen. Dies führte zur Entwicklung der Röntgen-Synchrotron-Mikro-Computertomographie, welche die Schaffung von Bildgebungseinrichtungen zur Untersuchung von Proben verschiedenster Art förderte, z.B. von Modellorganismen, um die Physiologie komplexer lebender Systeme besser zu verstehen. Die Entwicklung moderner Steuerungssysteme und Robotik ermöglichte die vollständige Automatisierung der Röntgenbildgebungsexperimente und die Kalibrierung der Parameter des Versuchsaufbaus während des Betriebs. Die Weiterentwicklung der digitalen Detektorsysteme führte zu Verbesserungen der Auflösung, des Dynamikbereichs, der Empfindlichkeit und anderer wesentlicher Eigenschaften. Diese Verbesserungen führten zu einer beträchtlichen Steigerung des Durchsatzes des Bildgebungsprozesses, aber auf der anderen Seite begannen die Experimente eine wesentlich größere Datenmenge von bis zu Dutzenden von Terabyte zu generieren, welche anschließend manuell verarbeitet wurden. Somit ebneten diese technischen Fortschritte den Weg für die Durchführung effizienterer Hochdurchsatzexperimente zur Untersuchung einer großen Anzahl von Proben, welche Datensätze von besserer Qualität produzierten. In der wissenschaftlichen Gemeinschaft besteht daher ein hoher Bedarf an einem effizienten, automatisierten Workflow für die Röntgendatenanalyse, welcher eine solche Datenlast bewältigen und wertvolle Erkenntnisse für die Fachexperten liefern kann. Die bestehenden Lösungen für einen solchen Workflow sind nicht direkt auf Hochdurchsatzexperimente anwendbar, da sie für Ad-hoc-Szenarien im Bereich der medizinischen Bildgebung entwickelt wurden. Daher sind sie nicht für Hochdurchsatzdatenströme optimiert und auch nicht in der Lage, die hierarchische Beschaffenheit von Proben zu nutzen. Die wichtigsten Beiträge der vorliegenden Arbeit sind ein neuer automatisierter Analyse-Workflow, der für die effiziente Verarbeitung heterogener Röntgendatensätze hierarchischer Natur geeignet ist. Der entwickelte Workflow basiert auf verbesserten Methoden zur Datenvorverarbeitung, Registrierung, Lokalisierung und Segmentierung. Jede Phase eines Arbeitsablaufs, die eine Trainingsphase beinhaltet, kann automatisch feinabgestimmt werden, um die besten Hyperparameter für den spezifischen Datensatz zu finden. Für die Analyse von Faserstrukturen in Proben wurde eine neue, hochgradig parallelisierbare 3D-Orientierungsanalysemethode entwickelt, die auf einem neuartigen Konzept der emittierenden Strahlen basiert und eine präzisere morphologische Analyse ermöglicht. Alle entwickelten Methoden wurden gründlich an synthetischen Datensätzen validiert, um ihre Anwendbarkeit unter verschiedenen Abbildungsbedingungen quantitativ zu bewerten. Es wurde gezeigt, dass der Workflow in der Lage ist, eine Reihe von Datensätzen ähnlicher Art zu verarbeiten. Darüber hinaus werden die effizienten CPU/GPU-Implementierungen des entwickelten Workflows und der Methoden vorgestellt und der Gemeinschaft als Module für die Sprache Python zur Verfügung gestellt. Der entwickelte automatisierte Analyse-Workflow wurde erfolgreich für Mikro-CT-Datensätze angewandt, die in Hochdurchsatzröntgenexperimenten im Bereich der Entwicklungsbiologie und Materialwissenschaft gewonnen wurden. Insbesondere wurde dieser Arbeitsablauf für die Analyse der Medaka-Fisch-Datensätze angewandt, was eine automatisierte Segmentierung und anschließende morphologische Analyse von Gehirn, Leber, Kopfnephronen und Herz ermöglichte. Darüber hinaus wurde die entwickelte Methode der 3D-Orientierungsanalyse bei der morphologischen Analyse von Polymergerüst-Datensätzen eingesetzt, um einen Herstellungsprozess in Richtung wünschenswerter Eigenschaften zu lenken