Simultaneous segmentation and beam-hardening correction in computed microtomography of rock cores

We propose a post-reconstruction correction procedure for the beam-hardening artifact that neither requires knowledge of the X-ray spectrum nor of the attenuation coefficients in multi-mineral geologic samples. The beam-hardening artifact in polychromatic X-ray computer tomography (CT) hampers segmentation of the phase assemblage in geologic samples. We show that in cylindrically shaped samples like rock cores, the X-ray attenuation value for a single phase depends mainly on the distance from the center of the cylinder. This relationship could be easily extracted from the CT data for every phase and used to infer the presence of these phases at all positions in the sample. Our new approach enables simple segmentation of the phases for further analysis. Additionally, we propose a simple mathematical exercise for parallel beam geometry to show how recent experimentally observed deviations from Beer's law could be incorporated to model the beam-hardening artifact and simultaneously be used for its removal. The beam-hardening artifact (BHA) in computer tomography depends on the shape of the scanned sample.In cylindrical samples: beam-hardening artifact is a function of spatial distance from the center of the cylinder (BHA curve).The knowledge of BHA curves of different phases are the key to segregate phases.The segmentation provides the possibility to use computer tomography data in a quantitative manner.

[1]  M. Wevers,et al.  Quantitative analysis of reservoir rocks by microfocus X-ray computerised tomography , 2000 .

[2]  Hod Lipson,et al.  Distilling Free-Form Natural Laws from Experimental Data , 2009, Science.

[3]  A. Macovski,et al.  Energy-selective reconstructions in X-ray computerised tomography , 1976, Physics in medicine and biology.

[4]  A Macovski,et al.  Energy dependent reconstruction in X-ray computerized tomography. , 1976, Computers in biology and medicine.

[5]  J. Hawkins,et al.  On Intelligence , 2004 .

[6]  Rudy Swennen,et al.  Beam hardening artifact reduction in microfocus computed tomography for improved quantitative coal characterization , 2006 .

[7]  Peter Aundal Toft,et al.  Fast curve estimation using preconditioned generalized Radon transform , 1996, IEEE Trans. Image Process..

[8]  Vladimir I Ivakhnenko A Novel Quasi-Linearization Method for CT Image Reconstruction in Scanners With a Multi-Energy Detector System , 2010, IEEE Transactions on Nuclear Science.

[9]  Matthias Franz,et al.  Reducing non-linear artifacts of multi-material objects in industrial 3D computed tomography , 2008 .

[10]  S. L. Wellington,et al.  Tomographic imaging of three‐phase flow experiments , 1987 .

[11]  P Rüegsegger,et al.  Standardization of Computed Tomography Images by Means of a Material-Selective Beam Hardening Correction , 1978, Journal of computer assisted tomography.

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

[13]  Adrian E. Raftery,et al.  Finding Curvilinear Features in Spatial Point Patterns: Principal Curve Clustering with Noise , 2000, IEEE Trans. Pattern Anal. Mach. Intell..

[14]  S. Samarasekera,et al.  Exact cone beam CT with a spiral scan. , 1998, Physics in medicine and biology.

[15]  J. Sijbers,et al.  A model-based correction method for beam hardening artefacts in X-ray microtomography , 2003 .

[16]  W. Kalender,et al.  Empirical cupping correction: a first-order raw data precorrection for cone-beam computed tomography. , 2006, Medical physics.

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

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

[19]  M. Abramowitz,et al.  Handbook of Mathematical Functions With Formulas, Graphs and Mathematical Tables (National Bureau of Standards Applied Mathematics Series No. 55) , 1965 .

[20]  R. Brooks,et al.  Beam hardening in X-ray reconstructive tomography , 1976 .

[21]  S. Krimmel,et al.  3D computed tomography using a microfocus X-ray source: Analysis of artifact formation in the reconstructed images using simulated as well as experimental projection data , 2005 .

[22]  R. Jennings,et al.  A method for comparing beam-hardening filter materials for diagnostic radiology. , 1988, Medical physics.

[23]  P. Hammersberg,et al.  Correction for beam hardening artefacts in computerised tomography. , 1998, Journal of X-ray science and technology.

[24]  J. Sijbers,et al.  An energy-based beam hardening model in tomography. , 2002, Physics in medicine and biology.

[25]  A. Macovski,et al.  Generalized image combinations in dual KVP digital radiography. , 1981, Medical physics.

[26]  Martine Wevers,et al.  Towards 3-D petrography: application of microfocus computer tomography in geological science , 2001 .

[27]  G. Herman Correction for beam hardening in computed tomography. , 1979, Physics in medicine and biology.