Distance, granulometry, skeleton

In this chapter, we present a series of concepts and operators based on the notion of distance. As often with mathematical morphology, there exists more than one way to present ideas, that are simultaneously equivalent and complementary. Here, our problem is to present methods to characterize sets of points based on metric, geometry and topology considerations. An important concept is that of the skeleton, which is of fundamental importance in pattern recognition, and has many practical application

Skeletonization and Partitioning of Digital Images Using Discrete Morse Theory

We show how discrete Morse theory provides a rigorous and unifying foundation for defining skeletons and partitions of grayscale digital images. We model a grayscale image as a cubical complex with a real-valued function defined on its vertices (the voxel values). This function is extended to a discrete gradient vector field using the algorithm presented in Robins, Wood, Sheppard TPAMI 33:1646 (2011). In the current paper we define basins (the building blocks of a partition) and segments of the skeleton using the stable and unstable sets associated with critical cells. The natural connection between Morse theory and homology allows us to prove the topological validity of these constructions; for example, that the skeleton is homotopic to the initial object. We simplify the basins and skeletons via Morse-theoretic cancellation of critical cells in the discrete gradient vector field using a strategy informed by persistent homology. Simple working Python code for our algorithms for efficient vector field traversal is included. Example data are taken from micro-CT images of porous materials, an application area where accurate topological models of pore connectivity are vital for fluid-flow modelling

Removing interference components in time frequency representations using morphological operators

Time-frequency representations have been of great interest in the analysis and classification of non-stationary signals. The use of highly selective transformation techniques is a valuable tool for obtaining accurate information for studies of this type. The Wigner-Ville distribution has high time and frequency selectivity in addition to meeting some interesting mathematical properties. However, due to the bi-linearity of the transform, interference terms emerge when the transform is applied over multi-component signals. In this paper, we propose a technique to remove cross-components from the Wigner-Ville transform using image processing algorithms. The proposed method exploits the advantages of non-linear morphological filters, using a spectrogram to obtain an adequate marker for the morphological processing of the Wigner-Ville transform. Unlike traditional smoothing techniques, this algorithm provides cross-term attenuations while preserving time-frequency resolutions. Moreover, it could also be applied to distributions with different interference geometries. The method has been applied to a set of different time-frequency transforms, with promising results. © 2011 Elsevier Inc. All rights reserved.This work was supported by the National R&D Program under Grant TEC2008-02975 (Spain), FEDER programme and Generalitat Valenciana CMAP 340.Gómez García, S.; Naranjo Ornedo, V.; Miralles Ricós, R. (2011). Removing interference components in time frequency representations using morphological operators. Journal of Visual Communication and Image Representation. 22(1):401-410. doi:10.1016/j.jvcir.2011.03.007S40141022

Mathematical Methods for the Quantification of Actin-Filaments in Microscopic Images

In cell biology confocal laser scanning microscopic images of the actin filament of human osteoblasts are produced to assess the cell development. This thesis aims at an advanced approach for accurate quantitative measurements about the morphology of the bright-ridge set of these microscopic images and thus about the actin filament. Therefore automatic preprocessing, tagging and quantification interplay to approximate the capabilities of the human observer to intuitively recognize the filaments correctly. Numerical experiments with random models confirm the accuracy of this approach

Skeletonization methods for image and volume inpainting

Image and shape restoration techniques are increasingly important in computer graphics. Many types of restoration techniques have been proposed in the 2D image-processing and according to our knowledge only one to volumetric data. Well-known examples of such techniques include digital inpainting, denoising, and morphological gap filling. However efficient and effective, such methods have several limitations with respect to the shape, size, distribution, and nature of the defects they can find and eliminate. We start by studying the use of 2D skeletons for the restoration of two-dimensional images. To this end, we show that skeletons are useful and efficient for volumetric data reconstruction. To explore our hypothesis in the 3D case, we first overview the existing state-of-the-art in 3D skeletonization methods, and conclude that no such method provides us with the features required by efficient and effective practical usage. We next propose a novel method for 3D skeletonization, and show how it complies with our desired quality requirements, which makes it thereby suitable for volumetric data reconstruction context. The joint results of our study show that skeletons are indeed effective tools to design a variety of shape restoration methods. Separately, our results show that suitable algorithms and implementations can be conceived to yield high end-to-end performance and quality of skeleton-based restoration methods. Finally, our practical applications can generate competitive results when compared to application areas such as digital hair removal and wire artifact removal

Computing the Component-Labeling and the Adjacency Tree of a Binary Digital Image in Near Logarithmic-Time

Connected component labeling (CCL) of binary images is one of the fundamental operations in real time applications. The adjacency tree (AdjT) of the connected components offers a region-based representation where each node represents a region which is surrounded by another region of the opposite color. In this paper, a fully parallel algorithm for computing the CCL and AdjT of a binary digital image is described and implemented, without the need of using any geometric information. The time complexity order for an image of m × n pixels under the assumption that a processing element exists for each pixel is near O(log(m+ n)). Results for a multicore processor show a very good scalability until the so-called memory bandwidth bottleneck is reached. The inherent parallelism of our approach points to the direction that even better results will be obtained in other less classical computing architectures.Ministerio de Economía y Competitividad MTM2016-81030-PMinisterio de Economía y Competitividad TEC2012-37868-C04-0