An Open Source GPU Accelerated Framework for Flexible Algebraic Reconstruction at Synchrotron Light Sources

The recent developments in detector technology made possible 4D (3D + time) X-ray microtomography with high spatial and time resolutions. The resolution and duration of such exper- iments is currently limited by destructive X-ray radiation . Algebraic reconstruction technique (ART) can incorporate a priori knowledge into a reconstruction model that will allow us to apply some ap- proaches to reduce an imaging dose and keep a good enough reconstruction quality. However, these techniques are very computationally demanding. In this paper we present a framework for ART reconstruction based on OpenCL technology. Our approach treats an algebraic method as a compo- sition of interacting blocks which perform different tasks , such as projection selection, minimization, projecting and regularization. These tasks are realised us ing multiple algorithms differing in perfor- mance, the quality of reconstruction, and the area of applicability. Our framework allows to freely combine algorithms to build the reconstruction chain. All algorithms are implemented with OpenCL and are able to run on a wide range of parallel hardware. As well the framework is easily scalable to clustered environment with MPI. We will describe the architecture of ART framework and evaluate the quality and performance on latest generation of GPU hardware from NVIDIA and AMD.

[1]  John Banhart,et al.  Early pore formation in aluminium foams studied by synchrotron-based microtomography and 3-D image analysis , 2009 .

[2]  E. Sidky,et al.  Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization , 2008, Physics in medicine and biology.

[3]  Kees Joost Batenburg,et al.  DART: A Practical Reconstruction Algorithm for Discrete Tomography , 2011, IEEE Transactions on Image Processing.

[4]  A. Rack,et al.  Synchrotron-based radioscopy employing spatio-temporal micro-resolution for studying fast phenomena in liquid metal foams , 2009, Journal of synchrotron radiation.

[5]  B. De Man,et al.  Distance-driven projection and backprojection , 2002, 2002 IEEE Nuclear Science Symposium Conference Record.

[6]  Faycal Kharfi,et al.  Mathematics and Physics of Computed Tomography (CT): Demonstrations and Practical Examples , 2013 .

[7]  Simon Rit,et al.  The Reconstruction Toolkit (RTK), an open-source cone-beam CT reconstruction toolkit based on the Insight Toolkit (ITK) , 2014 .

[8]  Emmanuel Brun,et al.  PyHST2: an hybrid distributed code for high speed tomographic reconstruction with iterative reconstruction and a priori knowledge capabilities , 2013, ArXiv.

[9]  P. Joseph An Improved Algorithm for Reprojecting Rays through Pixel Images , 1982 .

[10]  Marc Teboulle,et al.  A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems , 2009, SIAM J. Imaging Sci..

[11]  K. Mueller,et al.  Fast and accurate three-dimensional reconstruction from cone-beam projection data using algebraic methods , 1998 .

[12]  Jeffrey A. Fessler,et al.  3D Forward and Back-Projection for X-Ray CT Using Separable Footprints , 2010, IEEE Transactions on Medical Imaging.

[13]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[14]  J. Boone,et al.  Comprar Essential Physics Of Medical Imaging, International Edition | Jerrold T. Bushberg | 9781451118100 | Lippincott Williams & Wilkins , 2012 .

[15]  P. Gilbert Iterative methods for the three-dimensional reconstruction of an object from projections. , 1972, Journal of theoretical biology.

[16]  A. Kak,et al.  Simultaneous Algebraic Reconstruction Technique (SART): A Superior Implementation of the Art Algorithm , 1984, Ultrasonic imaging.

[17]  J. Strzelczyk The Essential Physics of Medical Imaging , 2003 .

[18]  Emmanuel Brun,et al.  A Dictionary Learning Approach with Overlap for the Low Dose Computed Tomography Reconstruction and Its Vectorial Application to Differential Phase Tomography , 2014, PloS one.

[19]  K J Batenburg,et al.  Performance improvements for iterative electron tomography reconstruction using graphics processing units (GPUs). , 2011, Journal of structural biology.

[20]  Tilo Baumbach,et al.  X-ray phase-contrast in vivo microtomography probes new aspects of Xenopus gastrulation , 2013, Nature.

[21]  Andreas Kopmann,et al.  UFO: A Scalable GPU-based Image Processing Framework for On-line Monitoring , 2012, 2012 IEEE 14th International Conference on High Performance Computing and Communication & 2012 IEEE 9th International Conference on Embedded Software and Systems.

[22]  Marc Teboulle,et al.  Fast Gradient-Based Algorithms for Constrained Total Variation Image Denoising and Deblurring Problems , 2009, IEEE Transactions on Image Processing.

[23]  Tom Goldstein,et al.  The Split Bregman Method for L1-Regularized Problems , 2009, SIAM J. Imaging Sci..

[24]  I. Daubechies,et al.  An iterative thresholding algorithm for linear inverse problems with a sparsity constraint , 2003, math/0307152.

[25]  O Betz,et al.  On the possibilities of hard X-ray imaging with high spatio-temporal resolution using polychromatic synchrotron radiation. , 2010, Journal of X-ray science and technology.

[26]  G. Herman,et al.  Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and x-ray photography. , 1970, Journal of theoretical biology.

[27]  James D. Murray Method of steepest descents , 1984 .

[28]  Y. X. Wang,et al.  Nuclear Instruments and Methods in Physics Research Section B : Beam Interactions with Materials and Atoms , 2018 .