48,223 research outputs found
Automatic alignment for three-dimensional tomographic reconstruction
In tomographic reconstruction, the goal is to reconstruct an unknown object
from a collection of line integrals. Given a complete sampling of such line
integrals for various angles and directions, explicit inverse formulas exist to
reconstruct the object. Given noisy and incomplete measurements, the inverse
problem is typically solved through a regularized least-squares approach. A
challenge for both approaches is that in practice the exact directions and
offsets of the x-rays are only known approximately due to, e.g. calibration
errors. Such errors lead to artifacts in the reconstructed image. In the case
of sufficient sampling and geometrically simple misalignment, the measurements
can be corrected by exploiting so-called consistency conditions. In other
cases, such conditions may not apply and we have to solve an additional inverse
problem to retrieve the angles and shifts. In this paper we propose a general
algorithmic framework for retrieving these parameters in conjunction with an
algebraic reconstruction technique. The proposed approach is illustrated by
numerical examples for both simulated data and an electron tomography dataset
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
A Tensor-Based Dictionary Learning Approach to Tomographic Image Reconstruction
We consider tomographic reconstruction using priors in the form of a
dictionary learned from training images. The reconstruction has two stages:
first we construct a tensor dictionary prior from our training data, and then
we pose the reconstruction problem in terms of recovering the expansion
coefficients in that dictionary. Our approach differs from past approaches in
that a) we use a third-order tensor representation for our images and b) we
recast the reconstruction problem using the tensor formulation. The dictionary
learning problem is presented as a non-negative tensor factorization problem
with sparsity constraints. The reconstruction problem is formulated in a convex
optimization framework by looking for a solution with a sparse representation
in the tensor dictionary. Numerical results show that our tensor formulation
leads to very sparse representations of both the training images and the
reconstructions due to the ability of representing repeated features compactly
in the dictionary.Comment: 29 page
- …