Generation of ground truth images to validate micro-CT image-processing pipelines

Abstract

Digital rock technology and pore-scale physics have become increasingly relevant topics in a wide range of porous media with important applications in subsurface engineering. This technology relies heavily on images of pore space and pore-level fluid distribution determined by X-ray microcomputed tomography (micro-CT). Digital images of pore space (or pore-scale fluid distribution) are typically obtained as gray-level images that first need to be processed and segmented to obtain the binary images that uniquely represent rock and pore (including fluid phases). This processing step is not trivial. Rock complexity, image quality, noise, and other artifacts prohibit the use of a standard processing workflow. Instead, an array of strategies of increasing sophistication has been developed. Typical processing pipelines consist of filtering, segmentation, and postprocessing steps. For each step, various choices and different options exist. This makes selection and validation of an optimum processing pipeline difficult. Using Darcy-scale quantities as a benchmark is not a good option because of rock heterogeneity and different scales of observation. Here, we present a conceptual workflow where noisy images are derived from a ground truth by systematically including typical image artifacts and noise. Artifacts and noise are not simply added to the images. Instead, tomographic forward projection and reconstruction steps are used to incorporate the artifacts in a physically correct way. A proof of concept of this workflow is demonstrated by comparing seven different image-segmentation pipelines ranging from absolute thresholding to a machine-learning approach (Trainable Weka Segmentation). The Trainable Weka Segmentation showed the best performance of the tested methods.