3D shape extraction segmentation and representation of soil microstructures using generalized cylinders

This paper focuses on the modeling of soil microstructures using generalized cylinders, with a specific application to pore space. The geometric modeling of these microstructures is a recent area of study, made possible by the improved performance of computed tomography techniques. X-scanners provide very-high-resolution 3D volume images ([email protected]) of soil samples in which pore spaces can be extracted by thresholding. However, in most cases, the pore space defines a complex volume shape that cannot be approximated using simple analytical functions. We propose representing this shape using a compact, stable, and robust piecewise approximation by means of generalized cylinders. This intrinsic shape representation conserves its topological and geometric properties. Our algorithm includes three main processing stages. The first stage consists in describing the volume shape using a minimum number of balls included within the shape, such that their union recovers the shape skeleton. The second stage involves the optimum extraction of simply connected chains of balls. The final stage copes with the approximation of each simply optimal chain using generalized cylinders: circular generalized cylinders, tori, cylinders, and truncated cones. This technique was applied to several data sets formed by real volume computed tomography soil samples. It was possible to demonstrate that our geometric representation supplied a good approximation of the pore space. We also stress the compactness and robustness of this method with respect to any changes affecting the initial data, as well as its coherence with the intuitive notion of pores. During future studies, this geometric pore space representation will be used to simulate biological dynamics.

[1]  Leonardo Flórez-Valencia.  Modèle d'état de cylindre généralisé et la quantification de sténoses artérielles en imagerie 3D , 2006 .

[2]  M. Blunt.  Flow in porous media — pore-network models and multiphase flow , 2001 .

[3]  Dimitris N. Metaxas,et al.  Dynamic 3D models with local and global deformations: deformable superquadrics , 1990, [1990] Proceedings Third International Conference on Computer Vision.

[4]  Edith Perrier,et al.  DXSoil, a library for 3D image analysis in soil science , 2002 .

[5]  Ivan Bricault,et al.  From Volume Medical Images to Quadratic Surface Patches , 1997, Comput. Vis. Image Underst..

[6]  Jean-Daniel Boissonnat,et al.  Stability and Computation of Medial Axes - a State-of-the-Art Report , 2009, Mathematical Foundations of Scientific Visualization, Computer Graphics, and Massive Data Exploration.

[7]  O. Monga,et al.  3D geometric structures and biological activity: Application to microbial soil organic matter decomposition in pore space , 2008 .

[8]  F. Chazal,et al.  The λ-medial axis , 2005 .

[9]  A. Ringrose-Voase.  Measurement of soil macropore geometry by image analysis of sections through impregnated soil , 1996, Plant and Soil.

[10]  A. McBratney,et al.  Acquisition and analysis of three‐component digital images of soil pore structure. II. Application to seed beds in a fallow management trial , 1992 .

[11]  Olivier Monga,et al.  Representing geometric structures in 3D tomography soil images: Application to pore-space modeling , 2007, Comput. Geosci..

[12]  Frédéric Chazal,et al.  The "lambda-medial axis" , 2005, Graph. Model..

[13]  W. B. Lindquist,et al.  Investigating 3D geometry of porous media from high resolution images , 1999 .

[14]  A. McBratney,et al.  Acquisition and analysis of three–component digital images of soil pore structure. I. Method , 1992 .

[15]  Olivier Monga,et al.  An Optimal Region Growing Algorithm for Image Segmentation , 1987, Int. J. Pattern Recognit. Artif. Intell..

[16]  Mariette Yvinec,et al.  Algorithmic geometry , 1998 .

[17]  Baba C. Vemuri,et al.  Shape Modeling with Front Propagation: A Level Set Approach , 1995, IEEE Trans. Pattern Anal. Mach. Intell..

[18]  Rainer Horn,et al.  Three-dimensional quantification of intra-aggregate pore-space features using synchrotron-radiation-based microtomography , 2008 .

[19]  Y. Capowiez.  Differences in burrowing behaviour and spatial interaction between the two earthworm species Aporrectodea nocturna and Allolobophora chlorotica , 2000, Biology and Fertility of Soils.

[20]  Steven L. Bryant,et al.  Network model evaluation of permeability and spatial correlation in a real random sphere packing , 1993 .

[21]  Christian Barillot,et al.  Robust 3D Segmentation of Anatomical Structures with Level Sets , 2000, MICCAI.

[22]  Guillaume Caumon,et al.  Visualization of grids conforming to geological structures: a topological approach , 2005, Comput. Geosci..

[23]  Matthew T. Balhoff,et al.  Coupling pore-scale networks to continuum-scale models of porous media , 2007, Comput. Geosci..

[24]  Bryant,et al.  Prediction of relative permeability in simple porous media. , 1992, Physical review. A, Atomic, molecular, and optical physics.

[25]  James F. Blinn,et al.  A generalization of algebraic surface drawing , 1982, SIGGRAPH.

[26]  J. Bouma,et al.  Parameters for describing soil macroporosity derived from staining patterns , 1998 .

[27]  Ivan Neulander,et al.  Rendering generalized cylinders using the A-buffer , 1997 .

[28]  Pascal Frey,et al.  MEDIT : An interactive Mesh visualization Software , 2001 .

[29]  Stig Bakke,et al.  A Predictive Network Model for Three-Phase Flow in Porous Media , 2000 .