Assessment of prior image induced nonlocal means regularization for low-dose CT reconstruction: Change in anatomy

Prior-image-based CT reconstruction using attenuation-mismatched priors

Prior-image-based reconstruction (PIBR) methods are powerful tools for reducing radiation doses and improving the image quality of low-dose computed tomography (CT). Apart from anatomical changes, prior and current images can also have different attenuations because they originated from different scanners or from the same scanner but with different x-ray beam qualities (e.g., kVp settings, beam filters) during data acquisition. In such scenarios, with attenuation-mismatched priors, PIBR is challenging. In this work, we investigate a specific PIBR method, called statistical image reconstruction, using normal-dose image-induced nonlocal means regularization (SIR-ndiNLM), to address PIBR with such attenuation-mismatched priors and achieve quantitative low-dose CT imaging. We propose two corrective schemes for the original SIR-ndiNLM method, (1) a global histogram-matching approach and (2) a local attenuation correction approach, to account for the attenuation differences between the prior and current images in PIBR. We validate the efficacy of the proposed schemes using images acquired from dual-energy CT scanners to simulate attenuation mismatches. Meanwhile, we utilize different CT slices to simulate anatomical mismatches or changes between the prior and the current low-dose image. We observe that the original SIR-ndiNLM introduces artifacts to the reconstruction when an attenuation-mismatched prior is used. Furthermore, we find that a larger attenuation mismatch between the prior and current images results in more severe artifacts in the SIR-ndiNLM reconstruction. Our two proposed corrective schemes enable SIR-ndiNLM to effectively handle the attenuation mismatch and anatomical changes between the two images and successfully eliminate the artifacts. We demonstrate that the proposed techniques permit SIR-ndiNLM to leverage the attenuation-mismatched prior and achieve quantitative low-dose CT reconstruction from both low-flux and sparse-view data acquisitions. This work permits robust and reliable PIBR for CT data acquired using different beam settings.

Unpaired Image Denoising via Wasserstein GAN in Low-Dose CT Image with Multi-Perceptual Loss and Fidelity Loss

The use of low-dose computed tomography (LDCT) in medical practice can effectively reduce the radiation risk of patients, but it may increase noise and artefacts, which can compromise diagnostic information. The methods based on deep learning can effectively improve image quality, but most of them use a training set of aligned image pairs, which are difficult to obtain in practice. In order to solve this problem, on the basis of the Wasserstein generative adversarial network (GAN) framework, we propose a generative adversarial network combining multi-perceptual loss and fidelity loss. Multi-perceptual loss uses the high-level semantic features of the image to achieve the purpose of noise suppression by minimizing the difference between the LDCT image and the normal-dose computed tomography (NDCT) image in the feature space. In addition, L2 loss is used to calculate the loss between the generated image and the original image to constrain the difference between the denoised image and the original image, so as to ensure that the image generated by the network using the unpaired images is not distorted. Experiments show that the proposed method performs comparably to the current deep learning methods which utilize paired image for image denoising.

Anisotropic structure property based image reconstruction method for limited-angle computed tomography

Limited-angle computed tomography (CT) may appear in restricted CT scans. Since the available projection data is incomplete, the images reconstructed by filtered back-projection (FBP) or algebraic reconstruction technique (ART) often encounter shading artifacts. However, using the anisotropy property of the shading artifacts that coincide with the characteristic of limited-angle CT images can reduce the shading artifacts. Considering this concept, we combine the anisotropy property of the shading artifacts with the anisotropic structure property of an image to develop a new algorithm for image reconstruction. Specifically, we propose an image reconstruction method based on adaptive weighted anisotropic total variation (AwATV). This method, termed as AwATV method for short, is designed to preserve image structures and then remove the shading artifacts. It characterizes both of above properties. The anisotropy property of the shading artifacts accounts for reducing artifacts, and the anisotropic structure property of an image accounts for preserving structures. In order to evaluate the performance of AwATV, we use the simulation projection data of FORBILD head phantom and real CT data for image reconstruction. Experimental results show that AwATV can always reconstruct images with higher SSIM and PSNR, and smaller RMSE, which means that AwATV enables to reconstruct images with higher quality in term of artifact reduction and structure preservation.

Multi-energy computed tomography reconstruction using a nonlocal spectral similarity model

Multi-energy computed tomography (MECT) is able to acquire simultaneous multi-energy measurements from one scan. In addition, it allows material differentiation and quantification effectively. However, due to the limited energy bin width, the number of photons detected in an energy-specific channel is smaller than that in traditional CT, which results in image quality degradation. To address this issue, in this work, we develop a statistical iterative reconstruction algorithm to acquire high-quality MECT images and high-accuracy material-specific images. Specifically, this algorithm fully incorporates redundant self-similarities within nonlocal regions in the MECT image at one bin and rich spectral similarities among MECT images at all bins. For simplicity, the presented algorithm is referred to as 'MECT-NSS'. Moreover, an efficient optimization algorithm is developed to solve the MECT-NSS objective function. Then, a comprehensive evaluation of parameter selection for the MECT-NSS algorithm is conducted. In the experiment, the datasets include images from three phantoms and one patient to validate and evaluate the MECT-NSS reconstruction performance. The qualitative and quantitative results demonstrate that the presented MECT-NSS can successfully yield better MECT image quality and more accurate material estimation than the competing algorithms.

Tensor decomposition and non-local means based spectral CT image denoising

BACKGROUNDAs one type of the state-of-the-art detectors, photon counting detectors are used in spectral computed tomography (CT) to classify the received photons into several energy channels and generate multichannel projections simultaneously. However, FBP reconstructed images contain severe noise due to the low photon counts in each energy channel.OBJECTIVEA spectral CT image denoising method based on tensor-decomposition and non-local means (TDNLM) is proposed.METHODSIn a CT image, it is widely accepted that there exists self-similarity over the spatial domain. In addition, because a multichannel CT image is obtained from the same object at different energies, images among different channels are highly correlated. Motivated by these two characteristics of the spectral CT images, tensor decomposition and non-local means are employed to recover fine structures in spectral CT images. Moreover, images in all energy channels are added together to form a high signal-to-noise ratio image, which is applied to encourage the signal preservation of the TDNLM. The combination of TD, NLM and the guidance of a high-quality image enhances the low-dose spectral CT, and a parameter selection strategy is designed to achieve the optimal image quality.RESULTSThe effectiveness of the developed algorithm is validated on both numerical simulations and realistic preclinical applications. The root mean square error (RMSE) and the structural similarity (SSIM) are used to quantitatively assess the image quality. The proposed method successfully restored high-quality images (average RMSE=0.0217 cm-1 and SSIM=0.987) from noisy spectral CT images (average RMSE=0.225 cm-1 and SSIM=0.633). In addition, RMSE of each decomposed material component is also remarkably reduced. Compared to the state-of-the-art iterative spectral CT reconstruction algorithms, the proposed method achieves comparable performance with dramatically reduced computational cost, resulting in a speedup of >50.CONCLUSIONSThe outstanding denoising performance, the high computational efficiency and the adaptive parameter selection strategy make the proposed method practical for spectral CT applications.

Regularization Analysis and Design for Prior-Image-Based X-Ray CT Reconstruction

Prior-image-based reconstruction (PIBR) methods have demonstrated great potential for radiation dose reduction in computed tomography applications. PIBR methods take advantage of shared anatomical information between sequential scans by incorporating a patient-specific prior image into the reconstruction objective function, often as a form of regularization. However, one major challenge with PIBR methods is how to optimally determine the prior image regularization strength which balances anatomical information from the prior image with data fitting to the current measurements. Too little prior information yields limited improvements over traditional model-based iterative reconstruction, while too much prior information can force anatomical features from the prior image not supported by the measurement data, concealing true anatomical changes. In this paper, we develop quantitative measures of the bias associated with PIBR. This bias exhibits as a fractional reconstructed contrast of the difference between the prior image and current anatomy, which is quite different from traditional reconstruction biases that are typically quantified in terms of spatial resolution or artifacts. We have derived an analytical relationship between the PIBR bias and prior image regularization strength and illustrated how this relationship can be used as a predictive tool to prospectively determine prior image regularization strength to admit specific kinds of anatomical change in the reconstruction. Because bias is dependent on local statistics, we further generalized shift-variant prior image penalties that permit uniform (shift invariant) admission of anatomical changes across the imaging field of view. We validated the mathematical framework in phantom studies and compared bias predictions with estimates based on brute force exhaustive evaluation using numerous iterative reconstructions across regularization values. The experimental results demonstrate that the proposed analytical approach can predict the bias-regularization relationship accurately, allowing for prospective determination of the prior image regularization strength in PIBR. Thus, the proposed approach provides an important tool for controlling image quality of PIBR methods in a reliable, robust, and efficient fashion.

Statistical CT reconstruction using region-aware texture preserving regularization learning from prior normal-dose CT image

In some clinical applications, prior normal-dose CT (NdCT) images are available, and the valuable textures and structure features in them may be used to promote follow-up low-dose CT (LdCT) reconstruction. This study aims to learn texture information from the NdCT images and leverage it for follow-up LdCT image reconstruction to preserve textures and structure features. Specifically, the proposed reconstruction method first learns the texture information from those patches with similar structures in NdCT image, and the similar patches can be clustered by searching context features efficiently from the surroundings of the current patch. Then it utilizes redundant texture information from the similar patches as a priori knowledge to describe specific regions in the LdCT image. The advanced region-aware texture preserving prior is termed as 'RATP'. The main advantage of the PATP prior is that it can properly learn the texture features from available NdCT images and adaptively characterize the region-specific structures in the LdCT image. The experiments using patient data were performed to evaluate the performance of the proposed method. The proposed RATP method demonstrated superior performance in LdCT imaging compared to the filtered back projection (FBP) and statistical iterative reconstruction (SIR) methods using Gaussian regularization, Huber regularization and the original texture preserving regularization.

Regularization strategies in statistical image reconstruction of low-dose x-ray CT: A review

Statistical image reconstruction (SIR) methods have shown potential to substantially improve the image quality of low-dose x-ray computed tomography (CT) as compared to the conventional filtered back-projection (FBP) method. According to the maximum a posteriori (MAP) estimation, the SIR methods are typically formulated by an objective function consisting of two terms: (a) a data-fidelity term that models imaging geometry and physical detection processes in projection data acquisition, and (b) a regularization term that reflects prior knowledge or expectations of the characteristics of the to-be-reconstructed image. SIR desires accurate system modeling of data acquisition, while the regularization term also has a strong influence on the quality of reconstructed images. A variety of regularization strategies have been proposed for SIR in the past decades, based on different assumptions, models, and prior knowledge. In this paper, we review the conceptual and mathematical bases of these regularization strategies and briefly illustrate their efficacies in SIR of low-dose CT.

Prospective Image Quality Analysis and Control for Prior-Image-Based Reconstruction of Low-Dose CT

Purpose

Prior-image-based reconstruction (PIBR) is a powerful tool for low-dose CT, however, the nonlinear behavior of such approaches are generally difficult to predict and control. Similarly, traditional image quality metrics do not capture potential biases exhibited in PIBR images. In this work, we identify a new bias metric and construct an analytical framework for prospectively predicting and controlling the relationship between prior image regularization strength and this bias in a reliable and quantitative fashion.

Methods

Bias associated with prior image regularization in PIBR can be described as the fraction of actual contrast change (between the prior image and current anatomy) that appears in the reconstruction. Using local approximation of the nonlinear PIBR objective, we develop an analytical relationship between local regularization, fractional contrast reconstructed, and true contrast change. This analytic tool allows prediction bias properties in a reconstructed PIBR image and includes the dependencies on the data acquisition, patient anatomy and change, and reconstruction parameters. Predictions are leveraged to provide reliable and repeatable image properties for varying data fidelity in simulation and physical cadaver experiments.

Results

The proposed analytical approach permits accurate prediction of reconstructed contrast relative to a gold standard based on exhaustive search based on numerous iterative reconstructions. The framework is used to control regularization parameters to enforce consistent change reconstructions over varying fluence levels and varying numbers of projection angles - enabling bias properties that are less location- and acquisition-dependent.

Conclusions

While PIBR methods have demonstrated a substantial ability for dose reduction, image properties associated with those images have been difficult to express and quantify using traditional metrics. The novel framework presented in this work not only quantifies this bias in an intuitive fashion, but it gives a way to predict and control the bias. Reliable and predictable reconstruction methods are a requirement for clinical imaging systems and the proposed framework is an important step translating PIBR methods to clinical application.