High-Throughput Segmentation of Tiled Biological Structures using Random-Walk Distance Transforms

Abstract Various 3D imaging techniques are routinely used to examine biological materials, the results of which are usually a stack of grayscale images. In order to quantify structural aspects of the biological materials, however, they must first be extracted from the dataset in a process called segmentation. If the individual structures to be extracted are in contact or very close to each other, distance-based segmentation methods utilizing the Euclidean distance transform are commonly employed. Major disadvantages of the Euclidean distance transform, however, are its susceptibility to noise (very common in biological data), which often leads to incorrect segmentations (i.e., poor separation of objects of interest), and its limitation of being only effective for roundish objects. In the present work, we propose an alternative distance transform method, the random-walk distance transform, and demonstrate its effectiveness in high-throughput segmentation of three microCT datasets of biological tilings (i.e., structures composed of a large number of similar repeating units). In contrast to the Euclidean distance transform, the random-walk approach represents the global, rather than the local, geometric character of the objects to be segmented and, thus, is less susceptible to noise. In addition, it is directly applicable to structures with anisotropic shape characteristics. Using three case studies—tessellated cartilage from a stingray, the dermal endoskeleton of a starfish, and the prismatic layer of a bivalve mollusc shell—we provide a typical workflow for the segmentation of tiled structures, describe core image processing concepts that are underused in biological research, and show that for each study system, large amounts of biologically-relevant data can be rapidly segmented, visualized, and analyzed.

[1]  Adam P. Summers,et al.  Ontogeny of the tessellated skeleton: insight from the skeletal growth of the round stingray Urobatis halleri , 2009, Journal of anatomy.

[2]  M. Elphick,et al.  Body wall structure in the starfish Asterias rubens , 2017, Journal of anatomy.

[3]  Ullrich Köthe,et al.  Ilastik: Interactive learning and segmentation toolkit , 2011, 2011 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

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

[5]  Charles Hansen,et al.  The Visualization Handbook , 2011 .

[6]  Ullrich Köthe,et al.  An Efficient Fusion Move Algorithm for the Minimum Cost Lifted Multicut Problem , 2016, ECCV.

[7]  Laurent Najman,et al.  Geodesic Saliency of Watershed Contours and Hierarchical Segmentation , 1996, IEEE Trans. Pattern Anal. Mach. Intell..

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

[9]  James C. Weaver,et al.  Uncovering Nature's Design Strategies through Parametric Modeling, Multi‐Material 3D Printing, and Mechanical Testing   , 2017 .

[10]  Badrinath Roysam,et al.  A hybrid 3D watershed algorithm incorporating gradient cues and object models for automatic segmentation of nuclei in confocal image stacks , 2003, Cytometry. Part A : the journal of the International Society for Analytical Cytology.

[11]  A. Rack,et al.  Self-similar mesostructure evolution of the growing mollusc shell reminiscent of thermodynamically driven grain growth , 2014, Nature materials.

[12]  Serge Beucher,et al.  Watershed, Hierarchical Segmentation and Waterfall Algorithm , 1994, ISMM.

[13]  S. Weiner,et al.  Spiers Memorial Lecture. Lessons from biomineralization: comparing the growth strategies of mollusc shell prismatic and nacreous layers in Atrina rigida. , 2003, Faraday discussions.

[14]  Hans-Christian Hege,et al.  amira: A Highly Interactive System for Visual Data Analysis , 2005, The Visualization Handbook.

[15]  David Knötel,et al.  Automated segmentation of complex patterns in biological tissues: Lessons from stingray tessellated cartilage , 2017, PloS one.

[16]  Wen Yang,et al.  Flexible Dermal Armor in Nature , 2012, JOM.

[17]  H. Atkinson Overview no. 65 , 1988 .

[18]  Mats Hillert,et al.  On the theory of normal and abnormal grain growth , 1965 .

[19]  R. Lemanis,et al.  Morphological and textural evolution of the prismatic ultrastructure in mollusc shells: A comparative study of Pinnidae species. , 2019, Acta biomaterialia.

[20]  Nicu Sebe,et al.  Computer Vision – ECCV 2016 , 2016, Lecture Notes in Computer Science.

[21]  Peter Fratzl,et al.  Ultrastructural and developmental features of the tessellated endoskeleton of elasmobranchs (sharks and rays) , 2016, Journal of anatomy.

[22]  Jack Snoeyink,et al.  Computing contour trees in all dimensions , 2000, SODA '00.

[23]  S. Austad Invertebrates , 1991, Experimental Gerontology.

[24]  P. Danielsson Euclidean distance mapping , 1980 .

[25]  Jakob Andreas Bærentzen,et al.  3D distance fields: a survey of techniques and applications , 2006, IEEE Transactions on Visualization and Computer Graphics.

[26]  Herbert Edelsbrunner,et al.  Topological Persistence and Simplification , 2000, Proceedings 41st Annual Symposium on Foundations of Computer Science.

[27]  Stuart R. Stock,et al.  MicroComputed Tomography: Methodology and Applications , 2008 .

[28]  J. Dunlop,et al.  Mechanical behavior of idealized, stingray-skeleton-inspired tiled composites as a function of geometry and material properties. , 2017, Journal of the mechanical behavior of biomedical materials.

[29]  Ullrich Köthe,et al.  Probabilistic image segmentation with closedness constraints , 2011, 2011 International Conference on Computer Vision.

[30]  O. Bøggild The shell structure of the Mollusks , 1930 .

[31]  B. Sakmann,et al.  Journal of Neuroscience Methods Automated Three-dimensional Detection and Counting of Neuron Somata , 2022 .

[32]  Michael Unser,et al.  Teaching image-processing programming in Java , 2003, IEEE Signal Process. Mag..

[33]  Wilhelm Burger,et al.  Digital Image Processing - An Algorithmic Introduction using Java , 2008, Texts in Computer Science.

[34]  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).

[35]  John D. Currey,et al.  The mechanical behaviour of some molluscan hard tissues , 2009 .

[36]  Ronen Basri,et al.  Shape representation and classification using the Poisson equation , 2004, Proceedings of the 2004 IEEE Computer Society Conference on Computer Vision and Pattern Recognition, 2004. CVPR 2004..

[37]  R. Al-Raoush,et al.  Distribution of local void ratio in porous media systems from 3D X-ray microtomography images , 2006 .

[38]  A. Rack,et al.  Biomineralization as a Paradigm of Directional Solidification: A Physical Model for Molluscan Shell Ultrastructural Morphogenesis , 2018, Advanced materials.

[39]  Dieter Kotschick,et al.  The Topology and Combinatorics of Soccer Balls , 2006 .

[40]  Valerio Pascucci,et al.  Contour trees and small seed sets for isosurface traversal , 1997, SCG '97.

[41]  J J Vaquero,et al.  Applying watershed algorithms to the segmentation of clustered nuclei. , 1998, Cytometry.