2,387 research outputs found
Fast Mojette Transform for Discrete Tomography
A new algorithm for reconstructing a two dimensional object from a set of one
dimensional projected views is presented that is both computationally exact and
experimentally practical. The algorithm has a computational complexity of O(n
log2 n) with n = N^2 for an NxN image, is robust in the presence of noise and
produces no artefacts in the reconstruction process, as is the case with
conventional tomographic methods. The reconstruction process is approximation
free because the object is assumed to be discrete and utilizes fully discrete
Radon transforms. Noise in the projection data can be suppressed further by
introducing redundancy in the reconstruction. The number of projections
required for exact reconstruction and the response to noise can be controlled
without comprising the digital nature of the algorithm. The digital projections
are those of the Mojette Transform, a form of discrete linogram. A simple
analytical mapping is developed that compacts these projections exactly into
symmetric periodic slices within the Discrete Fourier Transform. A new digital
angle set is constructed that allows the periodic slices to completely fill all
of the objects Discrete Fourier space. Techniques are proposed to acquire these
digital projections experimentally to enable fast and robust two dimensional
reconstructions.Comment: 22 pages, 13 figures, Submitted to Elsevier Signal Processin
The Discrete radon transform: A more efficient approach to image reconstruction
The Radon transform and its inversion are the mathematical keys that enable tomography. Radon transforms are defined for continuous objects with continuous projections at all angles in [0,π). In practice, however, we pre-filter discrete projections take
A parallel butterfly algorithm
The butterfly algorithm is a fast algorithm which approximately evaluates a
discrete analogue of the integral transform \int K(x,y) g(y) dy at large
numbers of target points when the kernel, K(x,y), is approximately low-rank
when restricted to subdomains satisfying a certain simple geometric condition.
In d dimensions with O(N^d) quasi-uniformly distributed source and target
points, when each appropriate submatrix of K is approximately rank-r, the
running time of the algorithm is at most O(r^2 N^d log N). A parallelization of
the butterfly algorithm is introduced which, assuming a message latency of
\alpha and per-process inverse bandwidth of \beta, executes in at most O(r^2
N^d/p log N + \beta r N^d/p + \alpha)log p) time using p processes. This
parallel algorithm was then instantiated in the form of the open-source
DistButterfly library for the special case where K(x,y)=exp(i \Phi(x,y)), where
\Phi(x,y) is a black-box, sufficiently smooth, real-valued phase function.
Experiments on Blue Gene/Q demonstrate impressive strong-scaling results for
important classes of phase functions. Using quasi-uniform sources, hyperbolic
Radon transforms and an analogue of a 3D generalized Radon transform were
respectively observed to strong-scale from 1-node/16-cores up to
1024-nodes/16,384-cores with greater than 90% and 82% efficiency, respectively.Comment: To appear in SIAM Journal on Scientific Computin
Full field inversion in photoacoustic tomography with variable sound speed
Recently, a novel measurement setup has been introduced to photoacoustic
tomography, that collects data in the form of projections of the full 3D
acoustic pressure distribution at a certain time instant. Existing imaging
algorithms for this kind of data assume a constant speed of sound. This
assumption is not always met in practice and thus leads to erroneous
reconstructions. In this paper, we present a two-step reconstruction method for
full field detection photoacoustic tomography that takes variable speed of
sound into account. In the first step, by applying the inverse Radon transform,
the pressure distribution at the measurement time is reconstructed point-wise
from the projection data. In the second step, one solves a final time wave
inversion problem where the initial pressure distribution is recovered from the
known pressure distribution at the measurement time. For the latter problem, we
derive an iterative solution approach, compute the required adjoint operator,
and show its uniqueness and stability
Four-dimensional tomographic reconstruction by time domain decomposition
Since the beginnings of tomography, the requirement that the sample does not
change during the acquisition of one tomographic rotation is unchanged. We
derived and successfully implemented a tomographic reconstruction method which
relaxes this decades-old requirement of static samples. In the presented
method, dynamic tomographic data sets are decomposed in the temporal domain
using basis functions and deploying an L1 regularization technique where the
penalty factor is taken for spatial and temporal derivatives. We implemented
the iterative algorithm for solving the regularization problem on modern GPU
systems to demonstrate its practical use
- …