PREDICTION AND CONTROL OF IMAGE PROPERTIES IN ADVANCED COMPUTED TOMOGRAPHY

Computed Tomography (CT) is an important technique that is in widespread use for disease diagnosis, monitoring, and interventional procedures. There are many varieties of CT including cone-beam CT (CBCT) that has exceptional high spatial resolution and spectral CT that incorporates energy-dependent measurements for advanced material discrimination. The goal of this research is to quantify image properties using a prospective prediction framework for advanced reconstruction in CBCT and spectral CT systems. These predictors analyze the dependencies of image properties on system configuration, acquisition strategy, and reconstruction regularization design. The prospective estimation of image properties facilitates novel system and acquisition design, adaptive and task-driven imaging, and tuning of regularization for robust and reliable performance. The proposed research quantifies the image properties of model-based iterative reconstruction (MBIR) in CBCT and model-based material decomposition (MBMD) in spectral CT, including spatial resolution, the generalized response to local perturbations, and noise correlation. Predictions are derived with a realistic system model including physical blur, noise correlation, and a poly-energetic model that applies to a variety of spectral CT protocols. Reconstruction methods combining data statistical fidelity and various advanced regularization designs are explored. Prediction accuracy is validated with measured image properties in both simulation and physical experiments. The theoretical understanding is applied to applications with prospective reconstruction regularization design

Model-Based Iterative Reconstruction in Cone-Beam Computed Tomography: Advanced Models of Imaging Physics and Prior Information

Cone-beam computed tomography (CBCT) represents a rapidly developing imaging modality that provides three-dimensional (3D) volumetric images with sub-millimeter spatial resolution and soft-tissue visibility from a single gantry rotation. CBCT tends to have small footprint, mechanical simplicity, open geometry, and low cost compared to conventional diagnostic CT. Because of these features, CBCT has been used in a variety of specialty diagnostic applications, image-guided radiation therapy (on-board CBCT), and surgical guidance (e.g., C-arm based CBCT). However, the current generation of CBCT systems face major challenges in low-contrast, soft-tissue image quality – for example, in the detection of acute intracranial hemorrhage (ICH), which requires a fairly high level of image uniformity, spatial resolution, and contrast resolution. Moreover, conventional approaches in both diagnostic and image-guided interventions that involve a series of imaging studies fail to leverage the wealth of patient-specific anatomical information available from previous scans. Leveraging the knowledge gained from prior images holds the potential for major gains in image quality and dose reduction. Model-based iterative reconstruction (MBIR) attempts to make more efficient use of the measurement data by incorporating a forward model of physical detection processes. Moreover, MBIR allows incorporation of various forms of prior information into image reconstruction, ranging from image smoothness and sharpness to patient-specific anatomical information. By leveraging such advantages, MBIR has demonstrated improved tradeoffs between image quality and radiation dose in multi-detector CT in the past decade and has recently shown similar promise in CBCT. However, the full potential of MBIR in CBCT is yet to be realized. This dissertation explores the capabilities of MBIR in improving image quality (especially low-contrast, soft-tissue image quality) and reducing radiation dose in CBCT. The presented work encompasses new MBIR methods that: i) modify the noise model in MBIR to compensate for noise amplification from artifact correction; ii) design regularization by explicitly incorporating task-based imaging performance as the objective; iii) mitigate truncation effects in a computationally efficient manner; iv) leverage a wealth of patient-specific anatomical information from a previously acquired image; and v) prospectively estimate the optimal amount of prior image information for accurate admission of specific anatomical changes. Specific clinical challenges are investigated in the detection of acute ICH and surveillance of lung nodules. The results show that MBIR can substantially improve low-contrast, soft-tissue image quality in CBCT and enable dose reduction techniques in sequential imaging studies. The thesis demonstrates that novel MBIR methods hold strong potential to overcome conventional barriers to CBCT image quality and open new clinical applications that would benefit from high-quality 3D imaging

Statistical Sinogram Restoration in Dual-Energy CT for PET Attenuation Correction

Dual-energy (DE) X-ray computed tomography (CT) has been found useful in various applications. In medical imaging, one promising application is using low-dose DECT for attenuation correction in positron emission tomography (PET). Existing approaches to sinogram material decomposition ignore noise characteristics and are based on logarithmic transforms, producing noisy component sinogram estimates for low-dose DECT. In this paper, we propose two novel sinogram restoration methods based on statistical models: penalized weighted least square (PWLS) and penalized likelihood (PL), yielding less noisy component sinogram estimates for low-dose DECT than classical methods. The proposed methods consequently provide more precise attenuation correction of the PET emission images than do previous methods for sinogram material decomposition with DECT. We report simulations that compare the proposed techniques and existing approaches.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85900/1/Fessler11.pd

Modeling and Development of Iterative Reconstruction Algorithms in Emerging X-ray Imaging Technologies

Many new promising X-ray-based biomedical imaging technologies have emerged over the last two decades. Five different novel X-ray based imaging technologies are discussed in this dissertation: differential phase-contrast tomography (DPCT), grating-based phase-contrast tomography (GB-PCT), spectral-CT (K-edge imaging), cone-beam computed tomography (CBCT), and in-line X-ray phase contrast (XPC) tomosynthesis. For each imaging modality, one or more specific problems prevent them being effectively or efficiently employed in clinical applications have been discussed. Firstly, to mitigate the long data-acquisition times and large radiation doses associated with use of analytic reconstruction methods in DPCT, we analyze the numerical and statistical properties of two classes of discrete imaging models that form the basis for iterative image reconstruction. Secondly, to improve image quality in grating-based phase-contrast tomography, we incorporate 2nd order statistical properties of the object property sinograms, including correlations between them, into the formulation of an advanced multi-channel (MC) image reconstruction algorithm, which reconstructs three object properties simultaneously. We developed an advanced algorithm based on the proximal point algorithm and the augmented Lagrangian method to rapidly solve the MC reconstruction problem. Thirdly, to mitigate image artifacts that arise from reduced-view and/or noisy decomposed sinogram data in K-edge imaging, we exploited the inherent sparseness of typical K-edge objects and incorporated the statistical properties of the decomposed sinograms to formulate two penalized weighted least square problems with a total variation (TV) penalty and a weighted sum of a TV penalty and an l1-norm penalty with a wavelet sparsifying transform. We employed a fast iterative shrinkage/thresholding algorithm (FISTA) and splitting-based FISTA algorithm to solve these two PWLS problems. Fourthly, to enable advanced iterative algorithms to obtain better diagnostic images and accurate patient positioning information in image-guided radiation therapy for CBCT in a few minutes, two accelerated variants of the FISTA for PLS-based image reconstruction are proposed. The algorithm acceleration is obtained by replacing the original gradient-descent step by a sub-problem that is solved by use of the ordered subset concept (OS-SART). In addition, we also present efficient numerical implementations of the proposed algorithms that exploit the massive data parallelism of multiple graphics processing units (GPUs). Finally, we employed our developed accelerated version of FISTA for dealing with the incomplete (and often noisy) data inherent to in-line XPC tomosynthesis which combines the concepts of tomosynthesis and in-line XPC imaging to utilize the advantages of both for biological imaging applications. We also investigate the depth resolution properties of XPC tomosynthesis and demonstrate that the z-resolution properties of XPC tomosynthesis is superior to that of conventional absorption-based tomosynthesis. To investigate all these proposed novel strategies and new algorithms in these different imaging modalities, we conducted computer simulation studies and real experimental data studies. The proposed reconstruction methods will facilitate the clinical or preclinical translation of these emerging imaging methods

IMPROVED IMAGE QUALITY IN CONE-BEAM COMPUTED TOMOGRAPHY FOR IMAGE-GUIDED INTERVENTIONS

In the past few decades, cone-beam computed tomography (CBCT) emerged as a rapidly developing imaging modality that provides single rotation 3D volumetric reconstruction with sub-millimeter spatial resolution. Compared to the conventional multi-detector CT (MDCT), CBCT exhibited a number of characteristics that are well suited to applications in image-guided interventions, including improved mechanical simplicity, higher portability, and lower cost. Although the current generation of CBCT has shown strong promise for high-resolution and high-contrast imaging (e.g., visualization of bone structures and surgical instrumentation), it is often believed that CBCT yields inferior contrast resolution compared to MDCT and is not suitable for soft-tissue imaging. Aiming at expanding the utility of CBCT in image-guided interventions, this dissertation concerns the development of advanced imaging systems and algorithms to tackle the challenges of soft-tissue contrast resolution. The presented material includes work encompassing: (i) a comprehensive simulation platform to generate realistic CBCT projections (e.g., as training data for deep learning approaches); (ii) a new projection domain statistical noise model to improve the noise-resolution tradeoff in model-based iterative reconstruction (MBIR); (iii) a novel method to avoid CBCT metal artifacts by optimization of the source-detector orbit; (iv) an integrated software pipeline to correct various forms of CBCT artifacts (i.e., lag, glare, scatter, beam hardening, patient motion, and truncation); (v) a new 3D reconstruction method that only reconstructs the difference image from the image prior for use in CBCT neuro-angiography; and (vi) a novel method for 3D image reconstruction (DL-Recon) that combines deep learning (DL)-based image synthesis network with physics-based models based on Bayesian estimation of the statical uncertainty of the neural network. Specific clinical challenges were investigated in monitoring patients in the neurological critical care unit (NCCU) and advancing intraoperative soft-tissue imaging capability in image-guided spinal and intracranial neurosurgery. The results show that the methods proposed in this work substantially improved soft-tissue contrast in CBCT. The thesis demonstrates that advanced imaging approaches based on accurate system models, novel artifact reduction methods, and emerging 3D image reconstruction algorithms can effectively tackle current challenges in soft-tissue contrast resolution and expand the application of CBCT in image-guided interventions