6 research outputs found
Fast algorithms and efficient GPU implementations for the Radon transform and the back-projection operator represented as convolution operators
The Radon transform and its adjoint, the back-projection operator, can both
be expressed as convolutions in log-polar coordinates. Hence, fast algorithms
for the application of the operators can be constructed by using FFT, if data
is resampled at log-polar coordinates. Radon data is typically measured on an
equally spaced grid in polar coordinates, and reconstructions are represented
(as images) in Cartesian coordinates. Therefore, in addition to FFT, several
steps of interpolation have to be conducted in order to apply the Radon
transform and the back-projection operator by means of convolutions.
Both the interpolation and the FFT operations can be efficiently implemented
on Graphical Processor Units (GPUs). For the interpolation, it is possible to
make use of the fact that linear interpolation is hard-wired on GPUs, meaning
that it has the same computational cost as direct memory access. Cubic order
interpolation schemes can be constructed by combining linear interpolation
steps which provides important computation speedup.
We provide details about how the Radon transform and the back-projection can
be implemented efficiently as convolution operators on GPUs. For large data
sizes, speedups of about 10 times are obtained in relation to the computational
times of other software packages based on GPU implementations of the Radon
transform and the back-projection operator. Moreover, speedups of more than a
1000 times are obtained against the CPU-implementations provided in the MATLAB
image processing toolbox
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
Distributed optimization for nonrigid nano-tomography
Resolution level and reconstruction quality in nano-computed tomography
(nano-CT) are in part limited by the stability of microscopes, because the
magnitude of mechanical vibrations during scanning becomes comparable to the
imaging resolution, and the ability of the samples to resist beam damage during
data acquisition. In such cases, there is no incentive in recovering the sample
state at different time steps like in time-resolved reconstruction methods, but
instead the goal is to retrieve a single reconstruction at the highest possible
spatial resolution and without any imaging artifacts. Here we propose a joint
solver for imaging samples at the nanoscale with projection alignment,
unwarping and regularization. Projection data consistency is regulated by dense
optical flow estimated by Farneback's algorithm, leading to sharp sample
reconstructions with less artifacts. Synthetic data tests show robustness of
the method to Poisson and low-frequency background noise. Applicability of the
method is demonstrated on two large-scale nano-imaging experimental data sets.Comment: Manuscript and supplementary materia