Optimization of image quality and acquisition time for lab-based X-ray microtomography using an iterative reconstruction algorithm

Abstract Non-invasive laboratory-based X-ray microtomography has been widely applied in many industrial and research disciplines. However, the main barrier to the use of laboratory systems compared to a synchrotron beamline is its much longer image acquisition time (hours per scan compared to seconds to minutes at a synchrotron), which results in limited application for dynamic in situ processes. Therefore, the majority of existing laboratory X-ray microtomography is limited to static imaging; relatively fast imaging (tens of minutes per scan) can only be achieved by sacrificing imaging quality, e.g. reducing exposure time or number of projections. To alleviate this barrier, we introduce an optimized implementation of a well-known iterative reconstruction algorithm that allows users to reconstruct tomographic images with reasonable image quality, but requires lower X-ray signal counts and fewer projections than conventional methods. Quantitative analysis and comparison between the iterative and the conventional filtered back-projection reconstruction algorithm was performed using a sandstone rock sample with and without liquid phases in the pore space. Overall, by implementing the iterative reconstruction algorithm, the required image acquisition time for samples such as this, with sparse object structure, can be reduced by a factor of up to 4 without measurable loss of sharpness or signal to noise ratio.

[1]  Martin J. Blunt,et al.  Pore‐by‐pore capillary pressure measurements using X‐ray microtomography at reservoir conditions: Curvature, snap‐off, and remobilization of residual CO2 , 2014 .

[2]  Nigel P. Brandon,et al.  Validation of a physically-based solid oxide fuel cell anode model combining 3D tomography and impedance spectroscopy , 2016 .

[3]  Kees Joost Batenburg,et al.  An Iterative CT Reconstruction Algorithm for Fast Fluid Flow Imaging , 2015, IEEE Transactions on Image Processing.

[4]  Ali Q. Raeini,et al.  Automatic measurement of contact angle in pore-space images , 2017 .

[5]  Jan Sijbers,et al.  Fast and flexible X-ray tomography using the ASTRA toolbox. , 2016, Optics express.

[6]  William R B Lionheart,et al.  4D-CT reconstruction with unified spatial-temporal patch-based regularization , 2015 .

[7]  Martin J. Blunt,et al.  Dynamic imaging of oil shale pyrolysis using synchrotron X‐ray microtomography , 2016 .

[8]  Jeffrey A. Fessler,et al.  Combining Ordered Subsets and Momentum for Accelerated X-Ray CT Image Reconstruction , 2015, IEEE Transactions on Medical Imaging.

[9]  Andre Phillion,et al.  Quantitative 3D Characterization of Solidification Structure and Defect Evolution in Al Alloys , 2012 .

[10]  Dorthe Wildenschild,et al.  Image processing of multiphase images obtained via X‐ray microtomography: A review , 2014 .

[11]  R. Ketcham,et al.  Acquisition, optimization and interpretation of X-ray computed tomographic imagery: applications to the geosciences , 2001 .

[12]  Martin J. Blunt,et al.  Pore-scale imaging of trapped supercritical carbon dioxide in sandstones and carbonates , 2014 .

[13]  Samuel Krevor,et al.  Dynamic fluid connectivity during steady-state multiphase flow in a sandstone , 2017, Proceedings of the National Academy of Sciences.

[14]  Martin J Blunt,et al.  X‐ray Microtomography of Intermittency in Multiphase Flow at Steady State Using a Differential Imaging Method , 2017, Water resources research.

[15]  Frederick R. Forst,et al.  On robust estimation of the location parameter , 1980 .

[16]  Nigel P. Brandon,et al.  The application of phase contrast X-ray techniques for imaging Li-ion battery electrodes , 2014 .

[17]  Martin J. Blunt,et al.  Quantification of sub-resolution porosity in carbonate rocks by applying high-salinity contrast brine using X-ray microtomography differential imaging , 2016 .

[18]  M. Vannier,et al.  Why do commercial CT scanners still employ traditional, filtered back-projection for image reconstruction? , 2009, Inverse problems.

[19]  Martin J. Blunt,et al.  Reservoir-condition pore-scale imaging of dolomite reaction with supercritical CO 2 acidified brine: Effect of pore-structure on reaction rate using velocity distribution analysis , 2018 .

[20]  Martin J Blunt,et al.  Dynamic three-dimensional pore-scale imaging of reaction in a carbonate at reservoir conditions. , 2015, Environmental science & technology.

[21]  E. Y. Sidky,et al.  How little data is enough? Phase-diagram analysis of sparsity-regularized X-ray computed tomography , 2014, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[22]  Philip J. Withers,et al.  Temporal sparsity exploiting nonlocal regularization for 4D computed tomography reconstruction , 2016, Journal of X-ray science and technology.

[23]  Y. Bresler,et al.  Sampling Requirements for Circular Cone Beam Tomography , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

[24]  David S. Eastwood,et al.  A novel high-temperature furnace for combined in situ synchrotron X-ray diffraction and infrared thermal imaging to investigate the effects of thermal gradients upon the structure of ceramic materials , 2014, Journal of synchrotron radiation.

[25]  Gengsheng Lawrence Zeng,et al.  Unmatched projector/backprojector pairs in an iterative reconstruction algorithm , 2000, IEEE Transactions on Medical Imaging.

[26]  Henrik Turbell,et al.  Cone-Beam Reconstruction Using Filtered Backprojection , 2001 .

[27]  P. Joseph An Improved Algorithm for Reprojecting Rays through Pixel Images , 1983, IEEE Transactions on Medical Imaging.

[28]  Francesco De Carlo,et al.  TomoPy: a framework for the analysis of synchrotron tomographic data , 2014, Journal of synchrotron radiation.

[29]  Christoph Rau,et al.  Dynamics of snap-off and pore-filling events during two-phase fluid flow in permeable media , 2017, Scientific Reports.

[30]  Y. Nesterov A method for solving the convex programming problem with convergence rate O(1/k^2) , 1983 .

[31]  P. J. Huber Robust Estimation of a Location Parameter , 1964 .

[32]  Veerle Cnudde,et al.  High-resolution X-ray computed tomography in geosciences: A review of the current technology and applications , 2013 .

[33]  Benoit Recur,et al.  Bayesian approach to time-resolved tomography. , 2015, Optics express.

[34]  R. Guedouar,et al.  A comparative study between matched and mis-matched projection/back projection pairs used with ASIRT reconstruction method , 2010 .

[35]  Christoph H. Arns,et al.  Image-based relative permeability upscaling from the pore scale , 2016 .

[36]  Martin J Blunt,et al.  Reaction Rates in Chemically Heterogeneous Rock: Coupled Impact of Structure and Flow Properties Studied by X-ray Microtomography. , 2017, Environmental science & technology.

[37]  H. Malcolm Hudson,et al.  Accelerated image reconstruction using ordered subsets of projection data , 1994, IEEE Trans. Medical Imaging.

[38]  Frieder Enzmann,et al.  Real-time 3D imaging of Haines jumps in porous media flow , 2013, Proceedings of the National Academy of Sciences.

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

[40]  R. Armstrong,et al.  Critical capillary number: Desaturation studied with fast X‐ray computed microtomography , 2014 .

[41]  J. Fessler Statistical Image Reconstruction Methods for Transmission Tomography , 2000 .

[42]  P. Withers,et al.  Quantitative X-ray tomography , 2014 .

[43]  J. C. Elliott,et al.  X‐ray microtomography , 1982, Journal of microscopy.

[44]  William R B Lionheart,et al.  SparseBeads data: benchmarking sparsity-regularized computed tomography , 2017 .

[45]  S. Stock Recent advances in X-ray microtomography applied to materials , 2008 .

[46]  M. Blunt,et al.  Multi-scale multi-dimensional microstructure imaging of oil shale pyrolysis using X-ray micro-tomography, automated ultra-high resolution SEM, MAPS Mineralogy and FIB-SEM , 2017 .

[47]  Tapan Mukerji,et al.  Digital rock physics benchmarks - Part I: Imaging and segmentation , 2013, Comput. Geosci..

[48]  A. Kingston,et al.  Dynamic tomography with a priori information. , 2011, Applied optics.

[49]  Kees Joost Batenburg,et al.  TVR-DART: A More Robust Algorithm for Discrete Tomography From Limited Projection Data With Automated Gray Value Estimation , 2016, IEEE Transactions on Image Processing.

[50]  Christoph H. Arns,et al.  Pore Scale Characterization of Carbonates Using X-Ray Microtomography , 2005 .

[51]  Philippe Gouze,et al.  X-ray microtomography characterization of porosity, permeability and reactive surface changes during dissolution. , 2011, Journal of contaminant hydrology.

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

[53]  Peter M. Joseph,et al.  View sampling requirements in fan beam computed tomography. , 1980 .

[54]  M. Blunt,et al.  Pore-scale imaging and modelling , 2013 .

[55]  P M Joseph,et al.  View sampling requirements in fan beam computed tomography. , 1980, Medical physics.

[56]  Veerle Cnudde,et al.  Fast laboratory-based micro-computed tomography for pore-scale research: Illustrative experiments and perspectives on the future , 2016 .

[57]  F. Boas,et al.  CT artifacts: Causes and reduction techniques , 2012 .

[58]  Emil Y. Sidky,et al.  Quantifying Admissible Undersampling for Sparsity-Exploiting Iterative Image Reconstruction in X-Ray CT , 2011, IEEE Transactions on Medical Imaging.

[59]  Arash Aghaei,et al.  Direct pore-to-core up-scaling of displacement processes: Dynamic pore network modeling and experimentation , 2015 .

[60]  Kees Joost Batenburg,et al.  Fast Tomographic Reconstruction From Limited Data Using Artificial Neural Networks , 2013, IEEE Transactions on Image Processing.

[61]  Martin J. Blunt,et al.  Pore-scale contact angle measurements at reservoir conditions using X-ray microtomography , 2014 .

[62]  Philip J. Withers,et al.  Mapping fibre failure in situ in carbon fibre reinforced polymers by fast synchrotron X-ray computed tomography , 2017 .

[63]  Martin J. Blunt,et al.  Visualization and quantification of capillary drainage in the pore space of laminated sandstone by a porous plate method using differential imaging X‐ray microtomography , 2017 .

[64]  Sally J. Marshall,et al.  The X-ray tomographic microscope: Three-dimensional perspectives of evolving microstructures , 1994 .

[65]  Martin J. Blunt,et al.  Multiphase Flow in Permeable Media: A Pore-Scale Perspective , 2017 .

[66]  Martin J. Blunt,et al.  The Imaging of Dynamic Multiphase Fluid Flow Using Synchrotron-Based X-ray Microtomography at Reservoir Conditions , 2015, Transport in Porous Media.

[67]  W. B. Lindquist,et al.  Pore and throat size distributions measured from synchrotron X-ray tomographic images of Fontaineble , 2000 .

[68]  Stephen J. Neethling,et al.  Multi-scale quantification of leaching performance using X-ray tomography , 2016 .

[69]  L. Feldkamp,et al.  Practical cone-beam algorithm , 1984 .