Regularization by Denoising Sub-sampled Newton Method for Spectral CT Multi-Material Decomposition
arXiv:2103.13909 · doi:10.1098/rsta.2020.0191
Abstract
Spectral Computed Tomography (CT) is an emerging technology that enables to estimate the concentration of basis materials within a scanned object by exploiting different photon energy spectra. In this work, we aim at efficiently solving a model-based maximum-a-posterior problem to reconstruct multi-materials images with application to spectral CT. In particular, we propose to solve a regularized optimization problem based on a plug-in image-denoising function using a randomized second order method. By approximating the Newton step using a sketching of the Hessian of the likelihood function, it is possible to reduce the complexity while retaining the complex prior structure given by the data-driven regularizer. We exploit a non-uniform block sub-sampling of the Hessian with inexact but efficient Conjugate gradient updates that require only Jacobian-vector products for denoising term. Finally, we show numerical and experimental results for spectral CT materials decomposition.
arXiv commentAccepted in Philosophical Transactions A, issue "Synergistic tomographic image reconstruction (Part 1)"
References in corpus (10)
- Beyond a Gaussian Denoiser: Residual Learning of Deep CNN for Image Denoising
- Sketching as a Tool for Numerical Linear Algebra
- Mini-Batch Semi-Stochastic Gradient Descent in the Proximal Setting
- An Online Plug-and-Play Algorithm for Regularized Image Reconstruction
- An algorithm for constrained one-step inversion of spectral CT data
- Sub-sampled Newton Methods with Non-uniform Sampling
- Sub-Sampled Newton Methods I: Globally Convergent Algorithms
- Sub-Sampled Newton Methods II: Local Convergence Rates
- Polyquant CT: direct electron and mass density reconstruction from a single polyenergetic source
- Regularization by Denoising Sub-sampled Newton Method for Spectral CT Multi-Material Decomposition