Advances in dual-energy computed tomography imaging of radiological properties

Dual-energy computed tomography (DECT) has shown great potential in the reduction of uncertainties of proton ranges and low energy photon cross section estimation used in radiation therapy planning. The work presented herein investigated three contributions for advancing DECT applications. 1) A linear and separable two-parameter DECT, the basis vector model (BVM) was used to estimate proton stopping power. Compared to other nonlinear two-parameter models in the literature, the BVM model shows a comparable accuracy achieved for typical human tissues. This model outperforms other nonlinear models in estimations of linear attenuation coefficients. This is the first study to clearly illustrate the advantages of linear model not only in accurately mapping radiological quantities for radiation therapy, but also in providing a unique model for accurate linear forward projection modelling, which is needed by the statistical iterative reconstruction (SIR) and other advanced DECT reconstruction algorithms. 2) Accurate DECT requires knowledge of x-ray beam properties. Using the Birch-Marshall1 model and beam hardening correction coefficients encoded in a CT scanner’s sinogram header files, an efficient and accurate way to estimate the x-ray spectrum is proposed. The merits of the proposed technique lie in requiring no physical transmission measurement after a one-time calibration against an independently measured spectrum. This technique can also be used in monitoring the aging of x-ray CT tubes. 3) An iterative filtered back projection with anatomical constraint (iFBP-AC) algorithm was also implemented on a digital phantom to evaluate its ability in mitigating beam hardening effects and supporting accurate material decomposition for in vivo imaging of photon cross section and proton stopping power. Compared to iFBP without constraints, both algorithms demonstrate high efficiency of convergence. For an idealized digital phantom, similar accuracy was observed under a noiseless situation. With clinically achievable noise level added to the sinograms, iFBP-AC greatly outperforms iFBP in prediction of photon linear attenuation at low energy, i.e., 28 keV. The estimated mean errors of iFBP and iFBP-AC for cortical bone are 1% and 0.7%, respectively; the standard deviations are 0.6% and 5%, respectively. The achieved accuracy of iFBP-AC shows robustness versus contrast level. Similar mean errors are maintained for muscle tissue. The standard deviation achieved by iFBP-AC is 1.2%. In contrast, the standard deviation yielded by iFBP is about 20.2%. The algorithm of iFBP-AC shows potential application of quantitative measurement of DECT. The contributions in this thesis aim to improve the clinical performance of DECT

Advanced Statistical Modeling for Model-Based Iterative Reconstruction for Single-Energy and Dual-Energy X-Ray CT

Model-based iterative reconstruction (MBIR) has been increasingly broadly applied as an improvement over traditional, analytical image reconstruction methods in X-ray CT, primarily due to its significant advantage in drastic dose reduction without diagnostic loss. Early success of the method in conventional CT has encouraged the extension to a wide range of applications that includes more advanced imaging modalities, such as dual-energy X-ray CT, and more challenging imaging conditions, such as low-dose and sparse-sampling scans, each requiring refined statistical models including the data model and the prior model. In this dissertation, we developed an MBIR algorithm for dual-energy CT that included a joint data-likelihood model to account for correlated data noise. Moreover, we developed a Gaussian-Mixture Markov random filed (GM-MRF) image model that can be used as a very expressive prior model in MBIR for X-ray CT reconstruction. The GM-MRF model is formed by merging individual patch-based Gaussian-mixture models and therefore leads to an expressive MRF model with easily estimated parameters. Experimental results with phantom and clinical datasets have demonstrated the improvement in image quality due to the advanced statistical modeling

Generating spectral dental panoramic images from single energy computed tomography volumes

Purpose: To implement a framework generating synthetic spectral panoramic images from single energy CT volumes. Using the framework output to compare the synthetic images against experimental spectral panoramic images for cross-verification. Methods: A simulation framework for generating synthetic spectral panoramic images from CT volumes is described. A cone beam CT scan of an anthropomorphic head phantom is used as input. An experimental spectral panoramic image of the same phantom is acquired. Results: The output of the framework of an anthropomorphic head phantom is compared against an experimental spectral panoramic image of the same phantom. The synthetic and experimental spectral panoramic images resemble each other considerably, especially the bone features. In the soft tissue images, there are some deviations, which are a result of the differences between the experimental and synthetic processing pipelines. Conclusions: It is demonstrated that generating synthetic spectral panoramic images from single energy CT volumes is possible. The synthetic images have many similarities with the experimental results, increasing the confidence in the correctness of the information contained within experimental spectral panoramic images and indicating that the synthetic images could be useful in further research

Proceedings of the second "international Traveling Workshop on Interactions between Sparse models and Technology" (iTWIST'14)

The implicit objective of the biennial "international - Traveling Workshop on Interactions between Sparse models and Technology" (iTWIST) is to foster collaboration between international scientific teams by disseminating ideas through both specific oral/poster presentations and free discussions. For its second edition, the iTWIST workshop took place in the medieval and picturesque town of Namur in Belgium, from Wednesday August 27th till Friday August 29th, 2014. The workshop was conveniently located in "The Arsenal" building within walking distance of both hotels and town center. iTWIST'14 has gathered about 70 international participants and has featured 9 invited talks, 10 oral presentations, and 14 posters on the following themes, all related to the theory, application and generalization of the "sparsity paradigm": Sparsity-driven data sensing and processing; Union of low dimensional subspaces; Beyond linear and convex inverse problem; Matrix/manifold/graph sensing/processing; Blind inverse problems and dictionary learning; Sparsity and computational neuroscience; Information theory, geometry and randomness; Complexity/accuracy tradeoffs in numerical methods; Sparsity? What's next?; Sparse machine learning and inference.Comment: 69 pages, 24 extended abstracts, iTWIST'14 website: http://sites.google.com/site/itwist1

Algorithms for enhanced artifact reduction and material recognition in computed tomography

Computed tomography (CT) imaging provides a non-destructive means to examine the interior of an object which is a valuable tool in medical and security applications. The variety of materials seen in the security applications is higher than in the medical applications. Factors such as clutter, presence of dense objects, and closely placed items in a bag or a parcel add to the difficulty of the material recognition in security applications. Metal and dense objects create image artifacts which degrade the image quality and deteriorate the recognition accuracy. Conventional CT machines scan the object using single source or dual source spectra and reconstruct the effective linear attenuation coefficient of voxels in the image which may not provide the sufficient information to identify the occupying materials. In this dissertation, we provide algorithmic solutions to enhance CT material recognition. We provide a set of algorithms to accommodate different classes of CT machines. First, we provide a metal artifact reduction algorithm for conventional CT machines which perform the measurements using single X-ray source spectrum. Compared to previous methods, our algorithm is robust to severe metal artifacts and accurately reconstructs the regions that are in proximity to metal. Second, we propose a novel joint segmentation and classification algorithm for dual-energy CT machines which extends prior work to capture spatial correlation in material X-ray attenuation properties. We show that the classification performance of our method surpasses the prior work's result. Third, we propose a new framework for reconstruction and classification using a new class of CT machines known as spectral CT which has been recently developed. Spectral CT uses multiple energy windows to scan the object, thus it captures data across higher energy dimensions per detector. Our reconstruction algorithm extracts essential features from the measured data by using spectral decomposition. We explore the effect of using different transforms in performing the measurement decomposition and we develop a new basis transform which encapsulates the sufficient information of the data and provides high classification accuracy. Furthermore, we extend our framework to perform the task of explosive detection. We show that our framework achieves high detection accuracy and it is robust to noise and variations. Lastly, we propose a combined algorithm for spectral CT, which jointly reconstructs images and labels each region in the image. We offer a tractable optimization method to solve the proposed discrete tomography problem. We show that our method outperforms the prior work in terms of both reconstruction quality and classification accuracy

4-D Tomographic Inference: Application to SPECT and MR-driven PET

Emission tomographic imaging is framed in the Bayesian and information theoretic framework. The first part of the thesis is inspired by the new possibilities offered by PET-MR systems, formulating models and algorithms for 4-D tomography and for the integration of information from multiple imaging modalities. The second part of the thesis extends the models described in the first part, focusing on the imaging hardware. Three key aspects for the design of new imaging systems are investigated: criteria and efficient algorithms for the optimisation and real-time adaptation of the parameters of the imaging hardware; learning the characteristics of the imaging hardware; exploiting the rich information provided by depthof- interaction (DOI) and energy resolving devices. The document concludes with the description of the NiftyRec software toolkit, developed to enable 4-D multi-modal tomographic inference

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