Fast binary CT using Fourier null space regularization (FNSR)

X-ray CT is increasingly being adopted in manufacturing as a non destructive inspection tool. Traditionally, industrial workflows follow a two step procedure of reconstruction followed by segmentation. Such workflows suffer from two main problems: (1) the reconstruction typically requires thousands of projections leading to increased data acquisition times. (2) The application of the segmentation process a posteriori is dependent on the quality of the original reconstruction and often does not preserve data fidelity. We present a fast iterative x-ray CT method which simultaneously reconstructs and segments an image from a limited number of projections called Fourier null space regularization (FNSR). The novelty of the approach is in the explicit updating of the image null space with values derived from a regularized image from the previous iteration, thus compensating for any missing projections and effectively regularizing the reconstruction. The speed of the method is achieved by directly applying the Fourier Slice Theorem where the non-uniform fast Fourier transform (NUFFT) is used to compute the frequency spectrum of the projections at their positions in the image k-space. At each iteration a segmented image is computed which is used to populate the null values of the image k-space effectively steering the reconstruction towards a binary solution. The effectiveness of the method to generate accurate reconstructions is demonstrated and benchmarked against other iterative reconstruction techniques using a series of numerical examples. Finally, FNSR is validated using industrial x-ray CT data where accurate reconstructions were achieved with 18 or more projections, a significant reduction from the 5000 needed by filtered back projection

Dynamical models and machine learning for supervised segmentation

This thesis is concerned with the problem of how to outline regions of interest in medical images, when the boundaries are weak or ambiguous and the region shapes are irregular. The focus on machine learning and interactivity leads to a common theme of the need to balance conflicting requirements. First, any machine learning method must strike a balance between how much it can learn and how well it generalises. Second, interactive methods must balance minimal user demand with maximal user control. To address the problem of weak boundaries,methods of supervised texture classification are investigated that do not use explicit texture features. These methods enable prior knowledge about the image to benefit any segmentation framework. A chosen dynamic contour model, based on probabilistic boundary tracking, combines these image priors with efficient modes of interaction. We show the benefits of the texture classifiers over intensity and gradient-based image models, in both classification and boundary extraction. To address the problem of irregular region shape, we devise a new type of statistical shape model (SSM) that does not use explicit boundary features or assume high-level similarity between region shapes. First, the models are used for shape discrimination, to constrain any segmentation framework by way of regularisation. Second, the SSMs are used for shape generation, allowing probabilistic segmentation frameworks to draw shapes from a prior distribution. The generative models also include novel methods to constrain shape generation according to information from both the image and user interactions. The shape models are first evaluated in terms of discrimination capability, and shown to out-perform other shape descriptors. Experiments also show that the shape models can benefit a standard type of segmentation algorithm by providing shape regularisers. We finally show how to exploit the shape models in supervised segmentation frameworks, and evaluate their benefits in user trials

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

Novel Methods for Multi-Shape Analysis

Multi-shape analysis has the objective to recognise, classify, or quantify morphological patterns or regularities within a set of shapes of a particular object class in order to better understand the object class of interest. One important aspect of multi-shape analysis are Statistical Shape Models (SSMs), where a collection of shapes is analysed and modelled within a statistical framework. SSMs can be used as (statistical) prior that describes which shapes are more likely and which shapes are less likely to be plausible instances of the object class of interest. Assuming that the object class of interest is known, such a prior can for example be used in order to reconstruct a three-dimensional surface from only a few known surface points. One relevant application of this surface reconstruction is 3D image segmentation in medical imaging, where the anatomical structure of interest is known a-priori and the surface points are obtained (either automatically or manually) from images. Frequently, Point Distribution Models (PDMs) are used to represent the distribution of shapes, where each shape is discretised and represented as labelled point set. With that, a shape can be interpreted as an element of a vector space, the so-called shape space, and the shape distribution in shape space can be estimated from a collection of given shape samples. One crucial aspect for the creation of PDMs that is tackled in this thesis is how to establish (bijective) correspondences across the collection of training shapes. Evaluated on brain shapes, the proposed method results in an improved model quality compared to existing approaches whilst at the same time being superior with respect to runtime. The second aspect considered in this work is how to learn a low-dimensional subspace of the shape space that is close to the training shapes, where all factors spanning this subspace have local support. Compared to previous work, the proposed method models the local support regions implicitly, such that no initialisation of the size and location of these regions is necessary, which is advantageous in scenarios where this information is not available. The third topic covered in this thesis is how to use an SSM in order to reconstruct a surface from only few surface points. By using a Gaussian Mixture Model (GMM) with anisotropic covariance matrices, which are oriented according to the surface normals, a more surface-oriented fitting is achieved compared to a purely point-based fitting when using the common Iterative Closest Point (ICP) algorithm. In comparison to ICP we find that the GMM-based approach gives superior accuracy and robustness on sparse data. Furthermore, this work covers the transformation synchronisation method, which is a procedure for removing noise that accounts for transitive inconsistency in the set of pairwise linear transformations. One interesting application of this methodology that is relevant in the context of multi-shape analysis is to solve the multi-alignment problem in an unbiased/reference-free manner. Moreover, by introducing an improvement of the numerical stability, the methodology can be used to solve the (affine) multi-image registration problem from pairwise registrations. Compared to reference-based multi-image registration, the proposed approach leads to an improved registration accuracy and is unbiased/reference-free, which makes it ideal for statistical analyses

An investigation into the limitations of myocardial perfusion imaging

Myocardial Perfusion Imaging (MPI) plays a very important role in the management of patients with suspected Coronary Artery Disease and its use has grown despite the shortcomings of the technique. Significant progress has been made in identifying the causes of these shortcomings and many solutions been suggested in the literature but the clinical sensitivity and specificity of the technique is still well below optimum. Monte Carlo Simulation is a very useful tool in identifying and guiding the understanding of the existing problems in MPI and this present study utilised this method to establish the basis of the simulations to be used and the way to analyse the results so that many of the causes of the attenuation defects, when using MPI, could be identified. This was achieved by investigating the effect that the different anatomical parts of the thorax have on the attenuation defects caused. A further aspect investigated was the impact that self-absorption in the heart has on these defects. The variability of these defects were further investigated by altering the position and orientation of the heart itself within the thorax and determining the effect it has on the attenuation defects caused. Results indicate that the attenuation caused is a very complicated process, that the self-absorption of the heart plays an extremely important role and the impact of the different positions and orientation of the heart inside the thorax are also significant. The distortion caused on the images by these factors was demonstrated by the intensity losses in the basal part and an over-estimation in the apical parts, which were clearly observable on the final clinical images, with the potential to affect clinical interpretation. Attenuation correction procedures using transmission sources, have been available for some time, but have not been adopted widely, amidst concern that they introduce additional artefacts. This study determined the effectiveness of these methods by establishing the level of correction obtained and whether additional artefacts were introduced. This included the effectiveness of the compensation achieved with the use of the latest commercially available comprehensive correction techniques. The technique investigated was “Flash3D" from Siemens providing transmission based attenuation correction, depth-dependent resolution recovery and scatter correction. The comparison between the defects and intensity losses predicted by the Monte Carlo Simulations and the corrections provided by this commercial correction technique revealed that solution is compensating almost entirely for these problems and therefore do provide substantial progress in overcoming the limitations of MPI. As a result of the improvements gained from applying these commercially available techniques and the accuracy established in this study for the mentioned technique it is strongly recommended that these new techniques be embraced by the wider Nuclear Medicine community so that the limitations in MPI can be reduced in clinical environment. Non-withstanding the above gains made there remains room for improvement by overcoming the of use transmission attenuation correction techniques by replacing them with emission based techniques. In this study two new related emission based attenuation correction techniques have been suggested and investigated and provide a promising prospect of overcoming these limitations

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