While model-based reconstruction methods have been successfully applied to flat-panel cone-beam CT (FP-CBCT) systems, typical implementations ignore both spatial correlations in the projection data as well as system blurs due to the detector and focal spot in the x-ray source. In this work, we develop a forward model for flat-panel-based systems that includes blur and noise correlation associated with finite focal spot size and an indirect detector (e.g., scintillator). This forward model is used to develop a staged reconstruction framework where projection data are deconvolved and log-transformed, followed by a generalized least-squares reconstruction that utilizes a non-diagonal statistical weighting to account for the correlation that arises from the acquisition and data processing chain. We investigate the performance of this novel reconstruction approach in both simulated data and in CBCT test-bench data. In comparison to traditional filtered backprojection and model-based methods that ignore noise correlation, the proposed approach yields a superior noise-resolution tradeoff. For example, for a system with 0.34 mm FWHM scintillator blur and 0.70 FWHM focal spot blur, using the a correlated noise model instead of an uncorrelated noise model increased resolution by 42% (with variance matched at 6.9 × 10−8 mm−2). While this advantage holds across a wide range of systems with differing blur characteristics, the improvements are greatest for systems where source blur is larger than detector blur.

Abstract-We present a novel reconstruction algorithm based on a general cone-beam CT forward model which is capable of incorporating the blur and noise correlations that are exhibited in flat-panel CBCT measurement data. Specifically, the proposed model may include scintillator blur, focal-spot blur, and noise correlations due to light spread in the scintillator. The proposed algorithm (GPL-BC) uses a Gaussian PenalizedLikelihood objective function which incorporates models of Blur and Correlated noise. In a simulation study, GPL-BC was able to achieve lower bias as compared to deblurring followed by FDK as well as a model-based reconstruction method without integration of measurement blur. In the same study, GPL-BC was able to achieve better line-pair reconstructions (in terms of segmented-image accuracy) as compared to deblurring followed by FDK, a model based method without blur, and a model based method with blur but not noise correlations. A prototype extremities quantitative cone-beam CT test bench was used to image a physical sample of human trabecular bone. These data were used to compare reconstructions using the proposed method and model based methods without blur and/or correlation to a registered µCT image of the same bone sample. The GPL-BC reconstructions resulted in more accurate trabecular bone segmentation. Multiple trabecular bone metrics, including Trabecular Thickness (Tb.Th.) were computed for each reconstruction approach as well as the µCT volume. The GPL-BC reconstruction provided the most accurate Tb.Th. measurement, 0.255 mm, as compared to the µCT derived value of 0.193 mm, followed by the GPL-B reconstruction, the GPL-I reconstruction, and then the FDK reconstruction (0.271 mm, 0.309 mm, and 0.335 mm, respectively).

Spectral information in CT maybe used for material decomposition to produce accurate reconstructions of material density and to separate materials with similar overall attenuation. Traditional methods separate the reconstruction and decomposition steps, often resulting in undesirable trade-offs (e.g. sampling constraints, a simplified spectral model). In this work, we present a model-based material decomposition algorithm which performs the reconstruction and decomposition simultaneously using a multienergy forward model. In a kV-switching simulation study, the presented method is capable of reconstructing iodine at 0.5 mg ml −1 with a contrast-tonoise ratio greater than two, as compared to 3.0 mg ml −1 for image domain decomposition. The presented method also enables novel acquisition methods, which was demonstrated in this work with a combined kV-switching/split-filter acquisition explored in simulation and physical test bench studies. This novel design used four spectral channels to decompose three materials: water, iodine, and gadolinium. In simulation, the presented method accurately reconstructed concentration value estimates with RMSE values of 4.86 mg ml −1 for water, 0.108 mg ml −1 for iodine and 0.170 mg ml −1 for gadolinium. In test-bench data, the RMSE values were 134 mg ml −1 ,5.26 mg ml −1 and 1.85 mg ml −1 , respectively. These studies demonstrate the ability of modelbased material decomposition to produce accurate concentration estimates in challenging spatial/ spectral sampling acquisitions.

Model-based iterative reconstruction techniques for CT that include a description of the noise statistics and a physical forward model of the image formation process have proven to increase image quality for many applications. Specifically, including models of the system blur into the physical forward model and thus implicitly performing a deconvolution of the projections during tomographic reconstruction, could demonstrate distinct improvements, especially in terms of resolution. However, the results strongly rely on an exact characterization of all components contributing to the system blur. Such characterizations can be laborious and even a slight mismatch can diminish image quality significantly. Therefore, we introduce a novel objective function, which enables us to jointly estimate system blur parameters during tomographic reconstruction. Conventional objective functions are biased in terms of blur and can yield lowest cost to blurred reconstructions with low noise levels. A key feature of our objective function is a new normalized sparsity measure for CT based on total-variation regularization, constructed to be less biased in terms of blur. We outline a solving strategy for jointly recovering low-dimensional blur parameters during tomographic reconstruction. We perform an extensive simulation study, evaluating the performance of the regularization and the dependency of the different parts of the objective function on the blur parameters. Scenarios with different regularization strengths and system blurs are investigated, demonstrating that we can recover the blur parameter used for the simulations. The proposed strategy is validated and the dependency of the objective function with the number of iterations is analyzed. Finally, our approach is experimentally validated on test-bench data of a human wrist phantom, where the estimated blur parameter coincides well with visual inspection. Our findings are not restricted to attenuation-based CT and may facilitate the recovery of more complex imaging model parameters.

Material decomposition for imaging multiple contrast agents in a single acquisition has been made possible by spectral CT: a modality which incorporates multiple photon energy spectral sensitivities into a single data collection. This work presents an investigation of a new approach to spectral CT which does not rely on energy-discriminating detectors or multiple x-ray sources. Instead, a tiled pattern of K-edge filters are placed in front of the x-ray to create spatially encoded spectra data. For improved sampling, the spatial-spectral filter is moved continuously with respect to the source. A model-based material decomposition algorithm is adopted to directly reconstruct multiple material densities from projection data that is sparse in each spectral channel. Physical effects associated with the x-ray focal spot size and motion blur for the moving filter are expected to impact overall performance. In this work, those physical effects are modeled and a performance analysis is conducted. Specifically, experiments are presented with simulated focal spot widths between 0.2 mm and 4.0 mm. Additionally, filter motion blur is simulated for a linear translation speeds between 50 mm/s and 450 mm/s. The performance differential between a 0.2 mm and a 1.0 mm focal spot is less than 15% suggesting feasibility of the approach with realistic x-ray tubes. Moreover, for reasonable filter actuation speeds, higher speeds are shown to decrease error (due to improved sampling) despite motion-based spectral blur.