Optimization based beam-hardening correction in CT under data integral invariant constraint

A balanced total-variation-Chambolle-Pock algorithm for EPR imaging

Total variation (TV) minimization algorithm is an effective algorithm capable of accurately reconstructing images from sparse projection data in a variety of imaging modalities including computed tomography (CT) and electron paramagnetic resonance imaging (EPRI). The data divergence constrained, TV minimization (DDcTV) model and its Chambolle-Pock (CP) solving algorithm have been proposed for CT. However, when the DDcTV-CP algorithm is applied to 3D EPRI, it suffers from slow convergence rate or divergence. We hypothesize that this is due to the magnitude imbalance between the data fidelity term and the TV regularization term. In this work, we propose a balanced TV (bTV) model incorporating a balance parameter and demonstrate its capability to avoid convergence issues for the 3D EPRI application. Simulation and real experiments show that the DDcTV-CP algorithm cannot guarantee convergence but the bTV-CP algorithm may guarantee convergence and achieve fast convergence by use of an appropriate balance parameter. Experiments also show that underweighting the balance parameter leads to slow convergence, whereas overweighting the balance parameter leads to divergence. The iteration-behavior change-law with the variation of the balance parameter is explained by use of the data tolerance ellipse and gradient descent principle. The findings and insights gained in this work may be applied to other imaging modalities and other constrained optimization problems.

Reduced iteration image reconstruction of incomplete projection CT using regularization strategy through Lp norm dictionary learning

BACKGROUND

For sparse and limited angle projection Computed Tomography (CT), the reconstructed image usually suffers from considerable artifacts due to undersampled data.

OBJECTIVE

To improve image reconstruction quality of sparse and limited angle projection CT, this study tested a novel reconstruction algorithm based on Dictionary Learning (DL) from sparse and limited projections.

METHODS

The study used signal sparse representation and feature extraction to render the DL technology, which is constrained by L2 and Lp norms, respectively. A Lp Norm Dictionary Learning term is suitable for regular term of objective function for CT image reconstruction. This is helpful for solving the objective function by combining algorithm of ART. Based on these features, the new algorithm of ART-DL-Lp is proposed for CT image reconstruction. The alternate solving strategy of the algorithm of "ART first, then adaptive DL" is provided in sequence. The impact on reconstruction results of ART-DL-Lp at different p values (0 < p < 1) is also considered.

RESULTS

For non-ideal projections with noise, the digital experiments show that ART-DL-Lp data were superior to those of ART, SART, and ART-DL-L2. Accordingly, the objective evaluation metrics for non-ideal situation of RMSE, MAE, PSNR, Residuals and SSIM are all better than those of contrasted three algorithms. The metrics curves of ART-DL-Lp algorithm are recorded as the best. In both incomplete projection situations, smaller p-value of ART-DL-Lp algorithm induces more close reconstructed images to the original form and better five objective evaluation metrics.

CONCLUSIONS

Overall, the reconstruction efficiency of the proposed ART-DL-Lp for CT imaging using the noisy incomplete projections outperforms ART, SART and ART-DL-L2 algorithms. For ART-DL-Lp algorithm, lower p-values result in better reconstruction performance.