Generalized Titarenko's algorithm for ring artefacts reduction

Abstract

A fast algorithm for ring artefact reduction in high-resolution micro-tomography with synchrotron radiation is presented. The new method is a generalization of the one proposed by Titarenko and collaborators, with a complete sinogram restoration prior to reconstruction with classical algorithms. The generalized algorithm can be performed in linear time and is easy to implement. Compared with the original approach, with an explicit solution, this approach is fast through the use of the conjugate gradient method. Also, low/high-resolution sinograms can be restored using higher/lower-order derivatives of the projections. Using different order for the derivative is an advantage over the classical Titarenko's approach. Several numerical results using the proposed method are provided, supporting our claims.

High Areal Energy 3D‐Interdigitated Micro‐Supercapacitors in Aqueous and Ionic Liquid Electrolytes

“…Pixel was ≈20 nm large but the true spatial resolution was 60 nm (60 nm features with at least 3 pixels can be seen). Reconstruction was made with Tomopy using the algorithm from the ASTRA toolbox and accurately approximating algebraic tomographic reconstruction by filtered back projection. The final 3D object was obtained using FIJI and AMIRA softwares.…”

Dissolution‐Assisted Pattern Formation During Olivine Carbonation

“…Total 720 projection images were acquired in 180°rotation angle range (i.e., 0.25°per angular step) with 2 s exposure time for each image. Three-dimensional reconstructions were performed with Tomopy, an open source collaborative framework for the analysis of synchrotron tomographic data Gürsoy et al, 2014;Pelt & Batenburg, 2015) using the generalized Titarenko's algorithm for ring removal (Miqueles et al, 2014). Movie S1 is a 3-D rendering visualization proposing a walk through the reconstructed slices, highlighting the segmentation of the wormhole shape porosity inside olivine grains.…”

Fast Algorithms and Efficient GPU Implementations for the Radon Transform and the Back-Projection Operator Represented as Convolution Operators

“…In some situations it is preferable to use an iterative method for doing reconstruction from tomographic data. This could for instance be because of missing data, e.g., that data for some angles are missing; suppression of artifacts [24]; or that additional information about the noise contamination can be used to improve the reconstruction results compared to direction filtered back-projections (2). Iterative reconstruction methods rely on applying the forward and back-projection operators several times.…”