In this study, we scanned the core of a cylindrical soil sample (60mm diameter and 100mm height) by X-ray Computed Tomography (CT) producing 300 consecutive 2D digital images with 16-bit gray level depth and a resolution of 32 microns (image size 676 × 676 pixels). The aim of this work was to determine the geometry and spatial distribution of the elements in a sample, related in this case to pore, solid and gravel, inside each 2D image for the latter reconstruction of the corresponding 3D approximation of the elements using the total set of 300 soil images. Therefore, it was possible to determine the relative percentage of each element present in each 2D image and, correspondingly, the structure and total percentage in the 3D reconstruction. The identification of elements in the 2D image slices was very well accomplished using three standard segmentation algorithms: k-Means, Fuzzy c-Means and Otsu multilevel. In order to compare and evaluate the quality of results, a non-uniformity (NU) measure was applied such that low values were indicative of homogeneous regions. Due to the depth of the greyscale of the images, the results were very similar with comparable statistics and homogeneity (NU values) among the detected materials of the three algorithms. That suggests that the pore, solid and gravel spaces were very well identified, and this is reflected through their connectivity in the 3D reconstruction. Additionally, the gray level depth was reduced to 8 bits and the same study was undertaken. In this case, the quality of results was comparable to the previous ones, as the number of elements and NU values were very close. However, this also depends largely on the high resolution of the images. Thereby, the soil sample of this work was very well characterized using the simplest and most common algorithms for image segmentation thanks to the high contrast and resolution, and regardless the depth of the grey-level.
Read full abstract