27,310 research outputs found
Quantifying admissible undersampling for sparsity-exploiting iterative image reconstruction in X-ray CT
Iterative image reconstruction (IIR) with sparsity-exploiting methods, such
as total variation (TV) minimization, investigated in compressive sensing (CS)
claim potentially large reductions in sampling requirements. Quantifying this
claim for computed tomography (CT) is non-trivial, because both full sampling
in the discrete-to-discrete imaging model and the reduction in sampling
admitted by sparsity-exploiting methods are ill-defined. The present article
proposes definitions of full sampling by introducing four sufficient-sampling
conditions (SSCs). The SSCs are based on the condition number of the system
matrix of a linear imaging model and address invertibility and stability. In
the example application of breast CT, the SSCs are used as reference points of
full sampling for quantifying the undersampling admitted by reconstruction
through TV-minimization. In numerical simulations, factors affecting admissible
undersampling are studied. Differences between few-view and few-detector bin
reconstruction as well as a relation between object sparsity and admitted
undersampling are quantified.Comment: Revised version that was submitted to IEEE Transactions on Medical
Imaging on 8/16/201
Convex optimization problem prototyping for image reconstruction in computed tomography with the Chambolle-Pock algorithm
The primal-dual optimization algorithm developed in Chambolle and Pock (CP),
2011 is applied to various convex optimization problems of interest in computed
tomography (CT) image reconstruction. This algorithm allows for rapid
prototyping of optimization problems for the purpose of designing iterative
image reconstruction algorithms for CT. The primal-dual algorithm is briefly
summarized in the article, and its potential for prototyping is demonstrated by
explicitly deriving CP algorithm instances for many optimization problems
relevant to CT. An example application modeling breast CT with low-intensity
X-ray illumination is presented.Comment: Resubmitted to Physics in Medicine and Biology. Text has been
modified according to referee comments, and typos in the equations have been
correcte
Jump-sparse and sparse recovery using Potts functionals
We recover jump-sparse and sparse signals from blurred incomplete data
corrupted by (possibly non-Gaussian) noise using inverse Potts energy
functionals. We obtain analytical results (existence of minimizers, complexity)
on inverse Potts functionals and provide relations to sparsity problems. We
then propose a new optimization method for these functionals which is based on
dynamic programming and the alternating direction method of multipliers (ADMM).
A series of experiments shows that the proposed method yields very satisfactory
jump-sparse and sparse reconstructions, respectively. We highlight the
capability of the method by comparing it with classical and recent approaches
such as TV minimization (jump-sparse signals), orthogonal matching pursuit,
iterative hard thresholding, and iteratively reweighted minimization
(sparse signals)
- …