Impact of the non‐negativity constraint in model‐based iterative reconstruction from CT data

Abstract

Purpose Model‐based iterative reconstruction is a promising approach to achieve dose reduction without affecting image quality in diagnostic x‐ray computed tomography (CT). In the problem formulation, it is common to enforce non‐negative values to accommodate the physical non‐negativity of x‐ray attenuation. Using this a priori information is believed to be beneficial in terms of image quality and convergence speed. However, enforcing non‐negativity imposes limitations on the problem formulation and the choice of optimization algorithm. For these reasons, it is critical to understand the value of the non‐negativity constraint. In this work, we present an investigation that sheds light on the impact of this constraint. Methods We primarily focus our investigation on the examination of properties of the converged solution. To avoid any possibly confounding bias, the reconstructions are all performed using a provably converging algorithm started from a zero volume. To keep the computational cost manageable, an axial CT scanning geometry with narrow collimation is employed. The investigation is divided into five experimental studies that challenge the non‐negativity constraint in various ways, including noise, beam hardening, parametric choices, truncation, and photon starvation. These studies are complemented by a sixth one that examines the effect of using ordered subsets to obtain a satisfactory approximate result within 50 iterations. All studies are based on real data, which come from three phantom scans and one clinical patient scan. The reconstructions with and without the non‐negativity constraint are compared in terms of image similarity and convergence speed. In select cases, the image similarity evaluation is augmented with quantitative image quality metrics such as the noise power spectrum and closeness to a known ground truth. Results For cases with moderate inconsistencies in the data, associated with noise and bone‐induced beam hardening, our results show that the non‐negativity constraint offers little benefit. By varying the regularization parameters in one of the studies, we observed that sufficient edge‐preserving regularization tends to dilute the value of the constraint. For cases with strong data inconsistencies, the results are mixed: the constraint can be both beneficial and deleterious; in either case, however, the difference between using the constraint or not is small relative to the overall level of error in the image. The results with ordered subsets are encouraging in that they show similar observations. In terms of convergence speed, we only observed one major effect, in the study with data truncation; this effect favored the use of the constraint, but had no impact on our ability to obtain the converged solution without constraint. Conclusions Our results did not highlight the non‐negativity constraint as being strongly beneficial for diagnostic CT imaging. Altogether, we thus conclude that in some imaging scenarios, the non‐negativity constraint could be disregarded to simplify the op...

Frequency Response and Distortion Properties of Reconstruction Algorithms in Computed Tomography

“…The intention of using the non-negative constraint in the SIRT algorithm was not primarily to reduce noise ( 28 ) . Thus, the attenuation of air was set as the constraint level, and noise originating from CT numbers above air may not be directly reduced by this type of mechanism.…”

“…Thus, the attenuation of air was set as the constraint level, and noise originating from CT numbers above air may not be directly reduced by this type of mechanism. However, it may reduce the convergence time ( 28 ) and due to its mechanism, by cutting off the CT number at a certain level, it may help identify other nonlinear distortion effects such as those arising from metal artefacts. Metal artefacts may occur due to cut-offs at the maximum CT number of the CT system ( 14 ) .…”

Empirical unbiased zero-count correction in photon counting CT

“…1,2 Correction schemes that have been developed to address non-positive measurements in EID-CT imaging include the naïve replacement of non-positive measurements with artificial positive values, 1,3,4 the filtering of the raw counts, 1, 5-8 and iterative reconstruction methods with imposed non-negativity constraints or appropriate statistical modeling. 4,[9][10][11] However, these preexisitng correction strategies may introduce CT number bias by replacing zeroes with artificial posiitive values, do not address the statistical bias associated with the non-linear log transformation, and may degrade spatial resolution or introduce correlations.…”

Kernel-Free Quadratic Surface Support Vector Regression with Non-Negative Constraints

“…At present, several types of constraints have been utilized, including non-negative constraints [1][2][3][4], monotonicity constraints [5][6][7][8], smoothing constraints [9][10][11], etc. Powell et al [12] proposed a Bayesian hierarchical model for estimating constraints conditional random fields to analyze the relationship between air pollution and health.…”