Image processing of multiphase images obtained via X‐ray microtomography: A review

Easier access to X-ray microtomography (μCT) facilities has provided much new insight from high-resolution imaging for various problems in porous media research. Pore space analysis with respect to functional properties usually requires segmentation of the intensity data into different classes. Image segmentation is a nontrivial problem that may have a profound impact on all subsequent image analyses. This review deals with two issues that are neglected in most of the recent studies on image segmentation: (i) focus on multiclass segmentation and (ii) detailed descriptions as to why a specific method may fail together with strategies for preventing the failure by applying suitable image enhancement prior to segmentation. In this way, the presented algorithms become very robust and are less prone to operator bias. Three different test images are examined: a synthetic image with ground-truth information, a synchrotron image of precision beads with three different fluids residing in the pore space, and a μCT image of a soil sample containing macropores, rocks, organic matter, and the soil matrix. Image blur is identified as the major cause for poor segmentation results. Other impairments of the raw data like noise, ring artifacts, and intensity variation can be removed with current image enhancement methods. Bayesian Markov random field segmentation, watershed segmentation, and converging active contours are well suited for multiclass segmentation, yet with different success to correct for partial volume effects and conserve small image features simultaneously.

[1]  Guillermo Sapiro,et al.  Segmentation-free skeletonization of gray-scale images via PDEs , 2000, Proceedings 2000 International Conference on Image Processing (Cat. No.00CH37101).

[2]  M. Tuller,et al.  Segmentation of X‐ray computed tomography images of porous materials: A crucial step for characterization and quantitative analysis of pore structures , 2009 .

[3]  Nikhil R. Pal,et al.  Image thresholding: Some new techniques , 1993, Signal Process..

[4]  Trond Varslot,et al.  IMAGE REGISTRATION: ENHANCING AND CALIBRATING X-RAY MICRO-CT IMAGING , 2008 .

[5]  Pau-Choo Chung,et al.  A Fast Algorithm for Multilevel Thresholding , 2001, J. Inf. Sci. Eng..

[6]  Raoul Kopelman,et al.  Percolation and cluster distribution. I. Cluster multiple labeling technique and critical concentration algorithm , 1976 .

[7]  Hans-Jörg Vogel,et al.  Soil‐structure development including seasonal dynamics in a long‐term fertilization experiment , 2011 .

[8]  M. L. Porter,et al.  Image analysis algorithms for estimating porous media multiphase flow variables from computed microtomography data: a validation study , 2010 .

[9]  W. E. Blanz,et al.  Image segmentation by pixel classification , 1981, Pattern Recognit..

[10]  Richard A. Ketcham,et al.  New algorithms for ring artifact removal , 2006, SPIE Optics + Photonics.

[11]  N. Otsu A threshold selection method from gray level histograms , 1979 .

[12]  Rosemary Knight,et al.  Numerical modeling of microscopic fluid distribution in porous media , 1990 .

[13]  P. Lions,et al.  Image selective smoothing and edge detection by nonlinear diffusion. II , 1992 .

[14]  Wolfgang Fink,et al.  Segmentation of X-Ray CT Data of Porous Materials: A Review of Global and Locally Adaptive Algorithms , 2015 .

[15]  Zhou Delong,et al.  Minimum Error Thresholding Based on Two Dimensional Histogram , 2009, 2009 WRI World Congress on Computer Science and Information Engineering.

[16]  W. Guitang,et al.  A new method for image segmentation , 2009, 2009 Asia-Pacific Conference on Computational Intelligence and Industrial Applications (PACIIA).

[17]  Adrian Sheppard,et al.  Techniques for image enhancement and segmentation of tomographic images of porous materials , 2004 .

[18]  Philippe Gouze,et al.  Permeability impairment of a limestone reservoir triggered by heterogeneous dissolution and particles migration during CO2‐rich injection , 2013 .

[19]  Cédric Gommes,et al.  Practical methods for measuring the tortuosity of porous materials from binary or gray‐tone tomographic reconstructions , 2009 .

[20]  Serge Beucher,et al.  Use of watersheds in contour detection , 1979 .

[21]  Mark L. Rivers,et al.  Using X-ray computed tomography in hydrology: systems, resolutions, and limitations , 2002 .

[22]  Ryan T. Armstrong,et al.  Linking pore-scale interfacial curvature to column-scale capillary pressure , 2012 .

[23]  Jean-Michel Morel,et al.  Nonlocal Image and Movie Denoising , 2008, International Journal of Computer Vision.

[24]  John C. Russ,et al.  The Image Processing Handbook , 2016, Microscopy and Microanalysis.

[25]  Sankar K. Pal,et al.  A review on image segmentation techniques , 1993, Pattern Recognit..

[26]  Josiane Zerubia,et al.  Bayesian image classification using Markov random fields , 1996, Image Vis. Comput..

[27]  Mark L. Rivers,et al.  Comparison of image segmentation methods in simulated 2D and 3D microtomographic images of soil aggregates , 2011 .

[28]  P.K Sahoo,et al.  A survey of thresholding techniques , 1988, Comput. Vis. Graph. Image Process..

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

[30]  Alasdair N. Houston,et al.  Effect of scanning and image reconstruction settings in X-ray computed microtomography on quality and segmentation of 3D soil images , 2013 .

[31]  D. Wildenschild,et al.  X-ray imaging and analysis techniques for quantifying pore-scale structure and processes in subsurface porous medium systems , 2013 .

[32]  Bülent Sankur,et al.  Survey over image thresholding techniques and quantitative performance evaluation , 2004, J. Electronic Imaging.

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

[34]  W. Brent Lindquist,et al.  Image Thresholding by Indicator Kriging , 1999, IEEE Trans. Pattern Anal. Mach. Intell..

[35]  Nabih N. Abdelmalek,et al.  Maximum likelihood thresholding based on population mixture models , 1992, Pattern Recognit..

[36]  M. Costanza-Robinson,et al.  X-ray microtomography determination of air-water interfacial area-water saturation relationships in sandy porous media. , 2008, Environmental science & technology.

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

[38]  Hans-Jörg Vogel,et al.  Quantification of soil structure based on Minkowski functions , 2010, Comput. Geosci..

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

[40]  Jean-Michel Morel,et al.  A non-local algorithm for image denoising , 2005, 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'05).

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

[42]  Anil K. Jain,et al.  Goal-Directed Evaluation of Binarization Methods , 1995, IEEE Trans. Pattern Anal. Mach. Intell..

[43]  Jitendra Malik,et al.  Scale-Space and Edge Detection Using Anisotropic Diffusion , 1990, IEEE Trans. Pattern Anal. Mach. Intell..

[44]  Patrick Pérez,et al.  Poisson image editing , 2003, ACM Trans. Graph..

[45]  S. D. Yanowitz,et al.  A new method for image segmentation , 1988, [1988 Proceedings] 9th International Conference on Pattern Recognition.

[46]  A Sheppard,et al.  On the challenges of measuring interfacial characteristics of three‐phase fluid flow with x‐ray microtomography , 2014, Journal of microscopy.

[47]  Andrew K. C. Wong,et al.  A new method for gray-level picture thresholding using the entropy of the histogram , 1985, Comput. Vis. Graph. Image Process..

[48]  John F. Canny,et al.  A Computational Approach to Edge Detection , 1986, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[49]  C. V. Jawahar,et al.  Investigations on fuzzy thresholding based on fuzzy clustering , 1997, Pattern Recognit..

[50]  Hans-Jörg Vogel,et al.  Segmentation of X-ray microtomography images of soil using gradient masks , 2010, Comput. Geosci..

[51]  Carsten Raven,et al.  Numerical removal of ring artifacts in microtomography , 1998 .

[52]  Márcio Portes de Albuquerque,et al.  Image thresholding using Tsallis entropy , 2004, Pattern Recognit. Lett..

[53]  William G. Gray,et al.  Approximation of interfacial properties in multiphase porous medium systems , 2007 .

[54]  Luc Vincent,et al.  Watersheds in Digital Spaces: An Efficient Algorithm Based on Immersion Simulations , 1991, IEEE Trans. Pattern Anal. Mach. Intell..

[55]  John D. Austin,et al.  Adaptive histogram equalization and its variations , 1987 .

[56]  P. Sivakumar,et al.  A REVIEW ON IMAGE SEGMENTATION TECHNIQUES , 2016 .

[57]  D. Wildenschild,et al.  Quantitative Analysis of Flow Processes in a Sand Using Synchrotron‐Based X‐ray Microtomography , 2005, Vadose Zone Journal.

[58]  David L. Silverstein,et al.  Prediction of Air-Water Interfacial Area in Wet Unsaturated Porous Media , 2000 .

[59]  Richard A. Ketcham,et al.  Computational methods for quantitative analysis of three-dimensional features in geological specimens , 2005 .

[60]  Manfred Krafczyk,et al.  Tomographical Imaging and Mathematical Description of Porous Media Used for the Prediction of Fluid Distribution , 2006 .

[61]  Jan Sijbers,et al.  Reduction of ring artefacts in high resolution micro-CT reconstructions. , 2004, Physics in medicine and biology.

[62]  Wei Wang,et al.  Observer-dependent variability of the thresholding step in the quantitative analysis of soil images and X-ray microtomography data , 2010 .

[63]  Hans-Jörg Vogel,et al.  Topological characterization of pore space in soil — sample preparation and digital image-processing , 1996 .

[64]  Jos B. T. M. Roerdink,et al.  The Watershed Transform: Definitions, Algorithms and Parallelization Strategies , 2000, Fundam. Informaticae.

[65]  J. Alex Stark,et al.  Adaptive image contrast enhancement using generalizations of histogram equalization , 2000, IEEE Trans. Image Process..

[66]  T. W. Ridler,et al.  Picture thresholding using an iterative selection method. , 1978 .

[67]  Paul Fieguth,et al.  Statistical fusion of two-scale images of porous media , 2009 .

[68]  Sungzoon Cho,et al.  Improvement of kittler and illingworth's minimum error thresholding , 1989, Pattern Recognit..

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

[70]  Emmanuelle Gouillart,et al.  scikit-image: image processing in Python , 2014, PeerJ.

[71]  Roger L. Nielsen,et al.  TRACE. FOR: a program for the calculation of combined major and trace-element liquid lines of descent for natural magmatic systems , 1988 .

[72]  Paul L. Rosin Unimodal thresholding , 2001, Pattern Recognit..

[73]  Christoph H. Arns,et al.  Pore scale imaging and modelling of coal properties , 2015 .

[74]  Thomas S. Huang,et al.  A fast two-dimensional median filtering algorithm , 1979 .

[75]  Ole Wendroth,et al.  Discrimination of Soil Phases by Dual Energy X‐ray Tomography , 1999 .

[76]  Wolfgang Fink,et al.  Three‐Dimensional Multiphase Segmentation of X‐Ray CT Data of Porous Materials Using a Bayesian Markov Random Field Framework , 2012 .

[77]  T. Fort,et al.  Determination of the air-water interfacial area in wet unsaturated porous media , 1996 .

[78]  B. Münch,et al.  Stripe and ring artifact removal with combined wavelet--Fourier filtering. , 2009, Optics express.

[79]  Du-Ming Tsai,et al.  A fast thresholding selection procedure for multimodal and unimodal histograms , 1995, Pattern Recognit. Lett..

[80]  홍기상 A pseudo-distance map for the segmentation-free skeletonization of gray-scale images , 2001 .

[81]  Alasdair N. Houston,et al.  Adaptive-window indicator kriging: A thresholding method for computed tomography images of porous media , 2013, Comput. Geosci..

[82]  Philippe Renard,et al.  Connectivity metrics for subsurface flow and transport , 2013 .

[83]  Azriel Rosenfeld,et al.  Image Segmentation by Pixel Classification in (Gray Level, Edge Value) Space , 1978, IEEE Transactions on Computers.

[84]  J. Besag On the Statistical Analysis of Dirty Pictures , 1986 .