Multi-phase classification by a least-squares support vector machine approach in tomography images of geological samples

Abstract. Image processing of X-ray-computed polychromatic cone-beam micro-tomography (μXCT) data of geological samples mainly involves artefact reduction and phase segmentation. For the former, the main beam-hardening (BH) artefact is removed by applying a best-fit quadratic surface algorithm to a given image data set (reconstructed slice), which minimizes the BH offsets of the attenuation data points from that surface. A Matlab code for this approach is provided in the Appendix. The final BH-corrected image is extracted from the residual data or from the difference between the surface elevation values and the original grey-scale values. For the segmentation, we propose a novel least-squares support vector machine (LS-SVM, an algorithm for pixel-based multi-phase classification) approach. A receiver operating characteristic (ROC) analysis was performed on BH-corrected and uncorrected samples to show that BH correction is in fact an important prerequisite for accurate multi-phase classification. The combination of the two approaches was thus used to classify successfully three different more or less complex multi-phase rock core samples.

[1]  Johan A. K. Suykens,et al.  Multiclass LS-SVMs: Moderated Outputs and Coding-Decoding Schemes , 2002, Neural Processing Letters.

[2]  Manuel Dierick,et al.  A novel beam hardening correction method requiring no prior knowledge, incorporated in an iterative reconstruction algorithm , 2012 .

[3]  Florian Fusseis,et al.  A brief guide to synchrotron radiation-based microtomography in (structural) geology and rock mechanics , 2014 .

[4]  K. Steiner,et al.  SPECTRO-MICROSCOPIC CHARACTERIZATION OF BIOGEOCHEMICAL INTERFACES IN SOIL 3 D simulation of the permeability tensor in a soil aggregate on basis of nanotomographic imaging and LBE solver , 2011 .

[5]  Zuleima T. Karpyn,et al.  Pore-Scale Lattice Boltzmann Modeling and 4D X-ray Computed Microtomography Imaging of Fracture-Matrix Fluid Transfer , 2013, Transport in Porous Media.

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

[7]  Johan A. K. Suykens,et al.  Least Squares Support Vector Machine Classifiers , 1999, Neural Processing Letters.

[8]  Alexander G. Schwing,et al.  Onset of Oil Mobilization and Nonwetting-Phase Cluster-Size Distribution , 2015 .

[9]  Anders Kaestner,et al.  Imaging and image processing in porous media research , 2008 .

[10]  D. Selvathi,et al.  Brain MRI Slices Classification Using Least Squares Support Vector Machine , 2007 .

[11]  JovanovićZoran,et al.  Simultaneous segmentation and beam-hardening correction in computed microtomography of rock cores , 2013 .

[12]  Wolfram Rühaak,et al.  Processing of rock core microtomography images: Using seven different machine learning algorithms , 2016, Comput. Geosci..

[13]  Vladimir Vapnik,et al.  The Nature of Statistical Learning , 1995 .

[14]  Christopher D. Brown,et al.  Receiver operating characteristics curves and related decision measures: A tutorial , 2006 .

[15]  E. Spangenberg,et al.  Microtomographic quantification of hydraulic clay mineral displacement effects during a CO2 sequestration experiment with saline aquifer sandstone. , 2013, Environmental science & technology.

[16]  Tom Fawcett,et al.  An introduction to ROC analysis , 2006, Pattern Recognit. Lett..

[17]  J. Schwarz,et al.  Simulation of Fluid Flow on Fractures and Implications for Reactive Transport Simulations , 2013, Transport in Porous Media.

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

[19]  Dorthe Wildenschild,et al.  Efficiently engineering pore-scale processes: The role of force dominance and topology during nonwetting phase trapping in porous media , 2015 .

[20]  Mohammad Hossein Fatemi,et al.  Classification of central nervous system agents by least squares support vector machine based on their structural descriptors: A comparative study , 2012 .

[21]  Patrick Haffner,et al.  Support vector machines for histogram-based image classification , 1999, IEEE Trans. Neural Networks.

[22]  V. Cnudde,et al.  Software tools for quantification of X-ray microtomography at the UGCT , 2007 .

[23]  Ethem Alpaydin,et al.  Introduction to machine learning , 2004, Adaptive computation and machine learning.

[24]  Vladimir N. Vapnik,et al.  The Nature of Statistical Learning Theory , 2000, Statistics for Engineering and Information Science.

[25]  Johan A. K. Suykens,et al.  Least Squares Support Vector Machines , 2002 .

[26]  N. Aronszajn.  Theory of Reproducing Kernels. , 1950 .

[27]  Sabine Van Huffel,et al.  Weighted LS-SVM for function estimation applied to artifact removal in bio-signal processing , 2010, 2010 Annual International Conference of the IEEE Engineering in Medicine and Biology.

[28]  J. Hanley,et al.  The meaning and use of the area under a receiver operating characteristic (ROC) curve. , 1982, Radiology.

[29]  Sotiris B. Kotsiantis,et al.  Supervised Machine Learning: A Review of Classification Techniques , 2007, Informatica.

[30]  János Urai,et al.  Multi-scale characterization of porosity in Boom Clay (HADES-level, Mol, Belgium) using a combination of X-ray μ-CT, 2D BIB-SEM and FIB-SEM tomography , 2015 .

[31]  J. Suykens,et al.  A tutorial on support vector machine-based methods for classification problems in chemometrics. , 2010, Analytica chimica acta.

[32]  R. Armstrong,et al.  Fast X-ray Micro-Tomography of Multiphase Flow in Berea Sandstone: A Sensitivity Study on Image Processing , 2014, Transport in Porous Media.

[33]  C. Steefel,et al.  Pore-scale controls on calcite dissolution rates from flow-through laboratory and numerical experiments. , 2014, Environmental science & technology.

[34]  G. Wise,et al.  A theoretical analysis of the properties of median filters , 1981 .

[35]  Markus Tuller,et al.  Application of Segmentation for Correction of Intensity Bias in X‐Ray Computed Tomography Images , 2010 .

[36]  William G. Gray,et al.  Interfacial area measurements for unsaturated flow through a porous medium , 2004 .

[37]  Frieder Enzmann,et al.  Simultaneous segmentation and beam-hardening correction in computed microtomography of rock cores , 2013, Comput. Geosci..

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

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

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

[41]  H. Vogel,et al.  Water Flow and Solute Transport in Unsaturated Sand—A Comprehensive Experimental Approach , 2015 .

[42]  Frieder Enzmann,et al.  Natural micro-scale heterogeneity induced solute and nanoparticle retardation in fractured crystalline rock. , 2012, Journal of contaminant hydrology.

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