<figure><img src="image_1.png" /></figure>

HAL
open science

# Geometrical and morphometrical tools for the inclusion analysis of metallic alloys

Johan Debayle

► To cite this version:

Johan Debayle. Geometrical and morphometrical tools for the inclusion analysis of metallic alloys. Metallurgical Research & Technology, 2019, 116 (5), pp.508. 10.1051/metal/2019011. hal-02279921

HAL Id: hal-02279921

https://hal.science/hal-02279921

Submitted on 5 Sep 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Geometrical and morphometrical tools for the inclusion analysis of metallic alloys

Johan Debayle*

École Nationale Supérieure des Mines de Saint-Étienne, SPIN/LGF UMR CNRS 5307, 158 cours Fauriel, 42023 Saint-Étienne, France

Received: 14 November 2018 / Accepted: 12 March 2019

**Abstract.** The mechanical and use properties of metal alloys depend on several factors, including the amount and the geometry of impurities (inclusions). In this context, image analysis enables these inclusions to be studied from digital images acquired by various systems such as optical/electron microscopy or X-ray tomography. This paper therefore aims to present some geometrical and morphometrical tools of image analysis, in order to characterize inclusions in metal alloys. To achieve this quantification, many geometrical and morphometrical features are traditionally used to quantitatively describe a population of objects (inclusions). Integral geometry, via Minkowski's functionals (in 2D: area, perimeter, Euler-Poincaré number), has been particularly investigated in image analysis. Nevertheless, they are sometimes insufficient for the characterization of complex microstructures (such as aggregates/agglomerates of objects). Other quantitative parameters are then necessary in order to discriminate or group different families of objects. In particular, shape diagrams are mathematical representations in the Euclidean plane for studying the morphology (shape) of objects, regardless of their size. In addition, this representation also makes it possible to analyze the evolution from one shape to another. In conclusion, image analysis using integral geometry and shape diagrams provide efficient tools with known mathematical properties to quantitatively describe inclusions (providing separate information on size and shape). The geometrical characteristics of these inclusions could thereafter be related to the mechanical properties of the metal alloys.

**Keywords:** geometrical characterization / image analysis / integral geometry / Minkowski functionals / shape diagrams

## 1 Introduction

The mechanical and use properties of metal alloys depend on several factors, including the amount and the geometry of impurities (inclusions). In this context, image analysis enables these inclusions to be studied from digital images acquired by various systems such as electron microscopy or X-ray tomography.

This paper therefore aims to present some geometrical and morphometrical tools of image analysis, in order to characterize inclusions in metal alloys. To achieve this quantification, many geometrical and morphometrical features are traditionally used to quantitatively describe a population of objects (inclusions) [1,2].

Integral geometry [3–5], via Minkowski’s functionals [6] (in 2D: area, perimeter, Euler-Poincaré number), has been particularly investigated in image analysis. Indeed, integral geometry [3–5] has generalized convex geometry to

finite unions of convex sets. This concept well applies to digital image analysis, since the smallest element (2D pixel) is indeed a convex set [7]. Nevertheless, they are sometimes insufficient for the characterization of complex microstructures (such as aggregates/agglomerates of objects) [8]. Other quantitative parameters are then necessary in order to discriminate or group different families of objects. In particular, shape diagrams [9–11] are mathematical representations in the Euclidean plane for studying the morphology (shape) of objects, regardless of their size. In addition, this representation also makes it possible to analyze the evolution from one shape to another. It also enables a convexity discrimination to be done with a direct visualization of the shape diagrams. For very complexe structures, such as fractal objects, other descriptors could be required. Indeed, these geometrical and morphometrical descriptors could be insufficient (as it stands) to analyse such strutures but the proposed measurements can be computed on some geometrical transformations of the object to be studied, providing some measurements

<footnote>* e-mail: debayle@emse.fr</footnote>
functions and not only scalars. This kind of approach has been particularly studied in [12]. Other descriptors are based on computational geometry and stochastic geometry [13] but they will not be studied in this paper.

The following of the paper is organized as follows: Section 2 gives the fundamentals notions of the Minkowski functionals. The computational aspects as well as the properties and the limitations of these functionals are done. Section 3 presents the shape diagrams and some illustration of the properties of continuity and convexity discrimination are done. The last section concludes this work. The objective of this paper is not to give an exhaustive list of geometrical and morphometrical tools, but rather to focus on specific descriptors with a compact representation and good mathematical properties and which are easy to compute.

## 2 Geometrical tools via integral geometry

Convex geometry, based on Minkowski functionals [6], allows to geometrically characterize a convex set in $\mathbb{R}^n$. There are exactly $n+1$ Minkowski functionals, satisfying certain properties (additivity, continuity, invariance by displacements). For example on $\mathbb{R}^2$, these functionals are the area, the perimeter and the Euler number. The two first functionals are related to well-known geometrical characteristics while the Euler number is a topological characteristic corresponding to the number of connected components minus the number of holes inside the set. The Minkowski functionals provide a basis of measurements in $\mathbb{R}^n$ [14]. They are therefore linearly independent and any other functional satisfying these properties is a linear combination of these properties. Integral geometry [3–5] has generalized convex geometry to finite unions of convex sets. This concept well applies to digital image analysis, since the smallest element (2D pixel) is indeed a convex set [7]. In the following, an efficient method (based on discrete geometry) for computing these Minkowski functionals from a 2D binary image is given [12]. Thereafter, the properties of the Minkowski functionals and their limitations will be presented.

### 2.1 Notions of discrete geometry

#### 2.1.1 Digital topology

Let $X$ be a two-dimensional binary image of size $l_1 \times l_2, (l_1, l_2) \in \mathbb{N}^2$. The spatial support of $X$ is denoted $D = [0, l_1] \cap \mathbb{N} \times [0, l_2] \cap \mathbb{N}$. $X$ is represented by a matrix $B = (b_x)_{x \in D}$, where $b_x$ belonging to $\{0,1\}$, is the value of the intensity of the pixel (abbreviated as “picture element”, the most small image element) of coordinates $x \in D$. By convention, the object and background of the image are represented by pixels of intensity $b_x=1$ and $b_x=0$, respectively. In particular, a set of connected pixels included in the object is a hole. The notion of neighbourhood induces topology on $D$:

- two pixels of the object (or background) with coordinates $x^1 = (x_1^1, x_2^1)$ and $x^2 = (x_1^2, x_2^2)$ are 4-adjacent if and only if $|x_1^1 - x_1^2| + |x_2^1 - x_2^2| = 1$;

<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Image composed of sixteen pixels (a) represented by crosses (x). The image is covered by intrapixel (b) or interpixel (c) cells, composed of vertices, edges and faces (vertice = ●, edge = -, face = ■). Note: image (a) is covered by sixteen intrapixel cells or nine interpixel cells [12].</figcaption></figure>

- two pixels of the object (or background) with coordinates $x^1 = (x_1^1, x_2^1)$ and $x^2 = (x_1^2, x_2^2)$ are 8-adjacents if and only if $\max(|x_1^1 - x_1^2|, |x_2^1 - x_2^2|) = 1$.

The p-neighbourhood ($p=4$ or $p=8$) of a pixel is therefore defined as the set of pixels p-adjacent to this one [Ros74]. The notion of neighbourhood makes it possible to establish the notion of connectedness [15]. Indeed, two pixels of the object (or background) are p-connected ($p=4$ or $p=8$) if there is a path of p-adjacent pixels connecting them. Similarly, a related component in a two-dimensional binary image is a set of adjacent pixels forming a path in the image. In order to satisfy Jordan's theorem, the connections must be different for the background pixels and for the object pixels. Therefore, these two-dimensional image connections are noted (4, 8) or (8, 4), the first coordinate corresponding to the object, and the second one at the bottom [16].

#### 2.1.2 Cell configuration

The “physical” spatial support of the image (i.e. the “continuous” spatial support) is covered by cells associated with pixels. A cell (size square the interpixel distance) is composed of one face, four edges and four vertices. Either a cell is centered in a pixel (cell intrapixel), or a cell is built by connecting pixels (interpixel cell). Figure 1 illustrates the overlapping of an image consisting of sixteen pixels, and its two possible representations.

4- and 8-adjacence are associated with the cell representation:

- using intrapixel cells, two pixels of the object (background) are adjacent if and only if their respective cells have a common vertex (8-adjacence);

- using interpixel cells, two pixels of the object (background) are adjacent if and only if they share the same edge (4-adjacence).

Figure 2 shows the two cell representations of the two-dimensional binary image $X$ associated with matrix $B$.

#### 2.1.3 Neighborhood configuration

In order to efficiently calculate the number of vertices $v$, edges $e$ and faces $f$ of the object (pixels of intensities equal to 1), the different neighbourhood configurations (size two
<figure><img src="image_4.png" /><figcaption>Fig. 2. Representation of intrapixel (b) and interpixel (c) cells of the binary image X associated to the matrix B (a) [12].</figcaption></figure>

by two pixels) of the original image X are determined [1,2,17–22]. Thus, to each pixel corresponds a neighborhood. In addition, sixteen configurations are possible (Tab. 1 where the values 0 and 1 represent the intensity of the pixels (1: object, 0: background).

Then, each configuration contributes for a known number of vertices, edges and faces (Eqs. (1–3)). To determine the neighborhood configurations of all pixels, an efficient algorithm consists in convolving the binary matrix B (associated with the image X) by a mask F of dimension two, whose values are powers of two, and whose origin is the pixel at the top left with value 1:

$$F = \begin{pmatrix} 1 & 4 \\ 2 & 8 \end{pmatrix}$$

Each value $\alpha \in [0, 15] \cap \mathbb{N}$ of the resulting matrix $B_*F$, where * denotes the convolution operator, corresponds to a known neighborhood configuration. Table 1 shows the possible configurations (1: object, 0: background) and the associated values $\alpha$ resulting from the convolution. The histogram $h$ of $B_*F$ gives the distribution of the neighborhood configurations of the X image. And each configuration contributes to a known number of vertices, edges and faces:

$$v = \sum_{\alpha=0}^{15} v_{\alpha} h(\alpha), \qquad (1)$$

$$e = \sum_{\alpha=0}^{15} e_{\alpha} h(\alpha), \qquad (2)$$

$$f = \sum_{\alpha=0}^{15} f_{\alpha} h(\alpha), \qquad (3)$$

where for $\alpha \in [0, 15] \cap \mathbb{N}$, $v_\alpha$, $e_\alpha$ and $f_\alpha$ are the coefficients of the linear combinations, depending on the image overlap, given in Table 2.

## 2.2 Minkowski functionals

Based on these notions of discrete geometry, the different Minkowski functionals can be computed from the number of vertices $v$, edges $e$ and faces $f$.

### 2.2.1 Euler number

The Euler number, denoted $\chi$, corresponding to the number of connected components minus the number of holes, can be computed from the image overlap by cells:

$$\chi = v - e + f, \qquad (4)$$

where $v$, $e$ and $f$ are the number of vertices, edges and faces, respectively.

It is important to note that equation (4) gives different numbers of Euler $\chi$, depending on the representation of cells. The representation of intrapixel (or interpixel) cells gives the result with (8, 4)-connexity (respectively (4, 8)-connexity). Each configuration of neighbourhoods contributes to a known number of vertices, edges and faces, and therefore for the number of Euler. Thanks to equations (1–4), the number of Euler is calculated as follows:

$$\chi^{(8,4)} = \sum_{\alpha=0}^{15} \chi_{\alpha}^{(8,4)} h(\alpha), \qquad (5)$$

$$\chi^{(4,8)} = \sum_{\alpha=0}^{15} \chi_{\alpha}^{(4,8)} h(\alpha), \qquad (6)$$

where $\chi_\alpha = v_\alpha - e_\alpha + f_\alpha$ is given in the Table 2. The Euler number is then efficiently computed by using the coefficients associated to each configuration, depending on the image overlap.

### 2.2.2 Perimeter

Two methods, giving different results, allow to estimate the perimeter of an object. Using intrapixel cells, the perimeter, denoted $P$ is determined by counting the number of edges common to two pixels of different intensities, thanks to the following formula:

$$P = -4f + 2e. \qquad (7)$$

The discretization implies that this perimeter is underestimated. Following the equations (1–3), $P$ is computed as:

$$P = \sum_{\alpha=0}^{15} P_{\alpha} h(\alpha), \qquad (8)$$

where $P_\alpha = v_\alpha - e_\alpha + f_\alpha$ is given in the Table 2.

The Crofton perimeter, denoted $P$ independent of the representation of the cells used, gives a better estimation. It associates the length of a curve with the number of times a “random” line intersects the set [23]. Let $\gamma$ be a planar curve, $l$ a line oriented in direction $\varphi$ and of length $p$, $n_\gamma(l)$ the number of points at which $\gamma$ and $l$ intersect. Crofton’s formula expresses the arc length of the $\gamma$ curve in terms of the spatial integral of all oriented lines:

$$\text{length}(\gamma) = \frac{1}{4} \iint n_{\gamma}(\varphi, p) d\varphi dp,$$
<table><caption>**Table 1.** Neighborhood configurations with the associated values $\alpha$.</caption><thead><tr><th>α</th><th>0</th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th><th>11</th><th>12</th><th>13</th><th>14</th><th>15</th></tr></thead><tbody><tr><th>0</th><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>1</td></tr><tr><th>0</th><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td><td>0</td><td>1</td><td>1</td><td>1</td><td>0</td><td>1</td><td>1</td></tr><tr><th>α</th><td>8</td><td>9</td><td>10</td><td>11</td><td>12</td><td>13</td><td>14</td><td>15</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>1</th><td>0</td><td>1</td><td>0</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>1</td><td>1</td></tr><tr><th>0</th><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td><td>0</td><td>1</td><td>1</td><td>1</td><td>0</td><td>1</td><td>1</td></tr></tbody></table>

<table><caption>**Table 2.** Coefficients of the linear combinations, depending or not on the cell representation [12]. The lines highlighted with the gray color correspond to the direct computation of the Minkowski functionals.</caption><thead><tr><th colspan="15">Intrapixel cells</th></tr><tr><th>α</th><th>0</th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th><th>11</th><th>12</th><th>13</th><th>14</th><th>15</th></tr></thead><tbody><tr><th>f<sub>α</sub></th><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td></tr><tr><th>ε<sub>α</sub></th><td>0</td><td>2</td><td>1</td><td>2</td><td>1</td><td>2</td><td>2</td><td>2</td><td>0</td><td>2</td><td>1</td><td>2</td><td>1</td><td>2</td><td>2</td><td>2</td></tr><tr><th>ν<sub>α</sub></th><td>0</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><th>χ<sub>α</sub><sup>(8,4)</sup></th><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>-1</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>-1</td><td>0</td></tr><tr><th>P<sub>α</sub></th><td>0</td><td>0</td><td>2</td><td>0</td><td>2</td><td>0</td><td>4</td><td>0</td><td>0</td><td>0</td><td>2</td><td>0</td><td>2</td><td>0</td><td>4</td><td>0</td></tr></tbody><thead><tr><th colspan="15">Interpixel cells</th></tr><tr><th>α</th><th>0</th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th><th>11</th><th>12</th><th>13</th><th>14</th><th>15</th></tr></thead><tbody><tr><th>f<sub>α</sub></th><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td></tr><tr><th>ε<sub>α</sub></th><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>2</td></tr><tr><th>ν<sub>α</sub></th><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>1</td></tr><tr><th>χ<sub>α</sub><sup>(4,8)</sup></th><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>-1</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr></tbody><thead><tr><th colspan="15">Intrapixel or interpixel cells</th></tr><tr><th>α</th><th>0</th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th><th>11</th><th>12</th><th>13</th><th>14</th><th>15</th></tr></thead><tbody><tr><th>i<sup>(1)</sup><sub>α</sub></th><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td><td>0</td><td>0</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><th>i<sup>(2π/4)</sup><sub>α</sub></th><td>0</td><td>0</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><th>i<sup>(2π/2)</sup><sub>α</sub></th><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><th>i<sup>(3π/4)</sup><sub>α</sub></th><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><th>P<sup>(4,8)</sup><sub>α</sub></th><td>0</td><td>π/2</td><td>0</td><td>0</td><td>0</td><td>π/2</td><td>0</td><td>0</td><td>π/2</td><td>π</td><td>0</td><td>0</td><td>π/2</td><td>π</td><td>0</td><td>0</td><td>0</td></tr><tr><th>P<sup>(8,4)</sup><sub>α</sub></th><td>0</td><td>π/4 + π/4√2</td><td>π/4√2</td><td>π/2√2</td><td>0</td><td>π/4 + π/4√2</td><td>0</td><td>π/4√2</td><td>π/4</td><td>π/2</td><td>π/4√2</td><td>π/4√2</td><td>π/4</td><td>π/2</td><td>0</td><td>0</td><td>0</td></tr><tr><th>A<sub>α</sub></th><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0</td><td>1</td><td>1</td></tr></tbody></table>

In the discrete case, only the angle orientations $0, \frac{\pi}{4}, \frac{\pi}{2}$, and $\frac{3\pi}{4}$ are considered; they are chosen according to the connectivity. The number of intercepts for each of these lines is counted, and denoted respectively $i^0, i^{\pi/4}, i^{\pi/2}$ and $i^{3\pi/4}$. They are then standardized by 1 or $\sqrt{2}$ depending on the orientation. The average is finally calculated by multiplying by $\pi$, a value representing the integral of the orientations, and by dividing by the number of orientations used. The Crofton perimeter is then calculated as follows:

$$ \bar{P}_{\alpha}^{(4,8)} = \frac{\pi}{2} (i^0 + i^{\pi/2}) \quad \text{for the (4,8) - connectivity, (9)} $$

$$ \bar{P}_{\alpha}^{(8,4)} = \frac{\pi}{4} \left( i^0 + \frac{i^{\pi/4}}{\sqrt{2}} + i^{\pi/2} + \frac{i^{3\pi/4}}{\sqrt{2}} \right) \quad \text{for the (8,4) - connectivity.} \quad (10) $$

This perimeter is computed thanks to the following linear combinations that provide coefficients to the different numbers of intercepts.

$$ \bar{P}_{\alpha}^{(4,8)} = \sum_{\alpha=0}^{15} \bar{P}_{\alpha}^{(4,8)} h(\alpha), \quad (11) $$

$$ \bar{P}_{\alpha}^{(8,4)} = \sum_{\alpha=0}^{15} \bar{P}_{\alpha}^{(8,4)} h(\alpha), \quad (12) $$

where the coefficients $\bar{P}_{\alpha}^{(4,8)}$ and $\bar{P}_{\alpha}^{(8,4)}$ are given in Table 2.

For wired objects, considered as objects of dimension one, the two previously defined perimeters do not give robust results, since they double the length of a line.

### 2.2.3 Area

The area (or surface), denoted $A$, of the set is generally defined by the number of pixels of the set. According to the representation of the cells, it is then equal to the number of vertices $v$ (interpixel cells) or the number of faces $f$ (intrapixel cells), the result being always the same. Generally, this overestimates the area. Thanks to equations (1–3), the area $A$ is calculated as follows:

$$ A = \sum_{\alpha=0}^{15} A_{\alpha} h(\alpha), \quad (13) $$

where the coefficients $A_{\alpha}$ are given in Table 2.
<figure><img src="image_5.png" /><figcaption>Fig. 3. Four binary images (300 × 300 pixels or 1 × 1 mm) of borosilicate glass visually different but having the same Minkowski functionals [8]: area 0.437 mm², perimeter 30.28 mm, number of Euler 11 or 16 depending on the connectivity).</figcaption></figure>

It is also possible to count the number of faces obtained with the representation of the interpixel cells. In this case, the area is generally underestimated, but the wired objects have no area (they are therefore considered as objects of dimension one), this result being more robust.

## 2.3 Properties

The $n+1$ Minkowski functionals $MF$ satisfy the following properties, for all compact convex sets $X$ and $Y$ in $\mathbb{R}^n$:

- increasingness: $X \subseteq Y \Rightarrow MF(X) \leq MF(Y)$;

- invariance to rigid transformations $t$: $MF(X) = MF(X_t)$ where $X_t$ denotes the transformation of $X$ by the rigid operator $t$ (translation and rotation);

- homogeneity: $MF(\lambda X) = \lambda^d MF(X)$ for $\lambda \in \mathbb{R}$ and $d$ the dimensionnality of $MF$;

- additivity: $MF(X) + MF(Y) = MF(X \cup Y) + MF(X \cap Y)$;

- continuity:

$$ \lim_{i \to +\infty} d_H(X_i, X) = 0 \Rightarrow \lim_{i \to +\infty} MF(X_i)MF(X) \text{ where} $$

$d_H$ denotes the Hausdorff distance [DD06] and $(X_i)_i$ denotes a family of convex sets.

Among these five properties, three are still valid for the finite unions of convex sets: invariance by rigid transformation, homogeneity and additivity.

Regarding Hadwiger's theorem [14], one other main important property is that any homogeneous and continuous functional that is invariant under rigid motions can be represented as a sum of the $n+1$ Minkowski functionals. In $\mathbb{R}^2$, it means that such a functional $\mu$ computed on a set $X$ can be expressed as:

$$ \mu(X) = c_0\chi(X) + c_1P(X) + c_2A(X), \quad (14) $$

where $c_0$, $c_1$ and $c_2$ are real-valued coefficients.

## 2.4 Limitations

Integral geometry provide useful and efficient tools for the geometrical characterization of sets via the Minkowski functionals. Nevertheless, they are not enough to discriminate complex spatial structures [8], as shown in Figure 3.

It is then necessary to provide other characteristics to be able to discriminate these spatial structures. In the literature, there exists a lot of other characteristics. In

the following, morphometrical functionals, and more particularly shape diagrams, will be presented. The advantage is that they have a compact representation, some good mathematical properties and a low computational cost.

# 3 Morphometrical tools via shape diagrams

Shape diagrams [9–11] are representations in the Euclidean plane introduced to study the morphology of 2D connected compact sets. Such a set is represented by a point within a shape diagram whose coordinate axes are morphometrical functionals defined as normalized ratios of geometrical functionals. In this section, the construction of the shape diagrams will be firstly given. The second part will be focused on some properties of continuity and convexity discrimination.

## 3.1 Morphometrical functionals

In addition to area $A$ and perimeter $P$, other geometrical functionals have been studied such as the radii $r$ and $R$ of the inscribed and circumscribed circles respectively, and the minimum and maximum Feret diameters and $d$ respectively [24]. The Feret diameter is a measure of an object size along a specified direction. In general, it can be defined as the distance between the two parallel planes restricting the object perpendicular to that direction. It is therefore also called the caliper diameter, referring to the measurement of the object size with a caliper. $\omega$ and $d$ are the minimal and maximal diameter on the different possible orientations. Figure 4 illustrates these four additional geometrical functionals.

In practice, the Feret diameters are easily determined by firstly computing the convex hull of the object. Indeed, the Feret diameters for an object and for its convex hull are equal. The convex hull being a polygon, one can use the “rotating callipers” algorithm [25] to directly determine the largest and smallest projections.

For a connected compact set, the relationships between these geometrical functionals are constrained by geometric inequalities [26,27]. These geometric inequalities link geometrical functionals by pairs. Furthermore, they allow to determine the morphometrical functionals (Tab. 3).
<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Geometrical functionals: radii of inscribed (r) and circumscribed (R) circles, minimum (ω) and maximum (d) Feret diameters.</figcaption></figure>

The morphometrical functionals are invariant under similitude transformations (consequently, they do not depend on the global size of the set) and are defined as ratios between two geometrical functionals. In these ratios, the units of the numerator and the denominator are dimensionally homogeneous and the result has therefore no unit. Moreover, a normalization by a constant value (scalar multiplication) allows to have a ratio that ranges within [0,1]. For each morphometrical functional, the scalar value depends directly on the associated geometric inequality [28]. In total, there are fifteen morphometrical functionals for a connected compact set. $4\pi A/P^2$, $r/R$, $\omega/d$ and $4R/P$ are four examples of these morphometrical functionals. Their concrete meanings are the roundness, the circularity, the diameter constancy and the thinness, respectively.

## 3.2 Shape diagrams

Shape diagrams can be defined by these fifteen morphometrical functionals. Each shape diagram enable to represent the morphometry of any connected compact set from two morphometrical functionals (that is to say from three geometrical functionals because the two denominators use the same geometrical functionals).

Let be any triplet of the considered six geometrical functionals (A, P, r, R,ω, d) and (M₁, M₂) be some particular morphometrical functionals valued in [0.1]². A shape diagram D is represented in the plane domain [0.1]² (whose axis coordinates are the morphometrical functionals M₁ and M₂) where any connected compact set S is mapped onto a point (x, y). Mathematically, a shape diagram D is obtained from the following mapping:

$$D: \begin{cases} K(E^2) \rightarrow [0, 1]^2 \\ S \mapsto (x, y) \end{cases},$$

<table><caption>**Table 3.** Shape functionals for simply connected compact sets. A, P, r, R, ω, d, denote the area, perimeter, radii of the inscribed and circumscribed circles, minimum and maximum Feret diameters [24], respectively.</caption><thead><tr><th>Geometrical functionals</th><th>Geometric inequalities</th><th>Morphological functionals</th><th>Extremal sets</th></tr></thead><tbody><tr><td>r, R</td><td>r ≤ R</td><td>r/R</td><td>C</td></tr><tr><td>ω, R</td><td>ω ≤ 2R</td><td>ω/2R</td><td>C</td></tr><tr><td>A, R</td><td>A ≤ πR<sup>2</sup></td><td>A/πR<sup>2</sup></td><td>C</td></tr><tr><td>d, R</td><td>d ≤ 2R</td><td>d/2R</td><td>Y</td></tr><tr><td>r, d</td><td>2r ≤ d</td><td>2r/d</td><td>C</td></tr><tr><td>ω, d</td><td>ω ≤ d</td><td>ω/d</td><td>W</td></tr><tr><td>A, d</td><td>4A ≤ πd<sup>2</sup></td><td>4A/πd<sup>2</sup></td><td>C</td></tr><tr><td>R, d</td><td>√3R ≤ d</td><td>√3R/d</td><td>Z</td></tr><tr><td>r, P</td><td>2πr ≤ P</td><td>2πr/P</td><td>C</td></tr><tr><td>ω, P</td><td>πω ≤ P</td><td>πω/P</td><td>W</td></tr><tr><td>A, P</td><td>4πA ≤ P<sup>2</sup></td><td>4πA/P<sup>2</sup></td><td>C</td></tr><tr><td>d, P</td><td>2d ≤ P</td><td>2d/P</td><td>L</td></tr><tr><td>R, P</td><td>4R ≤ P</td><td>4R/P</td><td>L</td></tr><tr><td>r, A</td><td>πr<sup>2</sup> ≤ A</td><td>πr<sup>2</sup>/A</td><td>C</td></tr><tr><td>r, ω</td><td>2r ≤ ω</td><td>2r/ω</td><td>X</td></tr></tbody></table>

Extremal sets are the sets for which an inequality becomes an equality.

C: the disks; W: the constant width compact convex sets; L: the line segments; X: many compact convex sets; Y: many simply connected compact sets; Z every compact convex set of diameter d containing an equilateral triangle of side-length d.

where $K(E^2)$ denotes the compact sets of the Euclidean 2D plane. Using all the fifteen morphometrical functionals, twenty-two shape diagrams are defined, denoted $(D_k)_{k \in \{1,30\} \cup \{7,10\} \cup \{17,20\}}$, respectively (Tab. 4).

Some geometric inequalities, and consequently some shape diagrams, are restricted to convex shapes [29], that are not considered in this paper. A detailed comparative study has been performed in order to analyze the representation relevance of these twenty-two shape diagrams [29–31].
For example, the shape diagram $D_{24} : (A, R, P)$ is obtained from the following mapping.

$$D_{24}(A, R, P) : \begin{cases} K(\mathbb{E}^2) \to [0, 1]^2 \\ S \to \left(\frac{4\pi A}{P^2}, \frac{4R}{P}\right), \end{cases}$$

The concrete meanings of the morphometrical functionals $4\pi A/P^2$ and $4R/P$ are the roundness and the thinness, respectively.

Figure 5 illustrates three shape diagrams where some elementary sets are mapped.

## 3.3 Properties

The properties of the different shape diagrams have been studied in [29–31]. In the following, the properties of continuity and convexity discrimination are presented in order to show the better relevance of some shape diagrams.

### 3.3.1 Continuity

Let two 2D connected compact sets and the connected compact set class allowing to switch from one to the other using a degree of freedom of the set. For example, the semi-circle goes to the semi-disk through semi-rings whose the cavity radius decreases (Fig. 6).

For each pair $i$ of sets, a class $C_i$ of connected compact sets with one degree of freedom could be defined. Thus, a curve denoted $C_{i,k}$ from each simply connected compact set class $C_i$ is created in each shape diagram $D_k$. These curves are related to the continuity of the shape diagram with respect to a small deformation of the set. This process can be used for elementary some pairs of connected compact sets.

In the following Figure 7, eight classes of connected compact sets are illustrated on three shape diagrams. It shows the location of the continuity curves from one shape to another.

A detailed study of this property can be found in [29–31].

### 3.3.2 Convexity discrimination

The convexity discrimination first requires the definition of the shape convexity. A set is convex when the line segment which joins any two points in it lies totally within the set. In other terms, the shape convexity could be quantified with the probability that two points in the set lies totally within it. Convexity parameters are commonly used in the analysis of shapes. The measurement value of the shape convexity of any set ranges between 0 and 1 (it is a probability). A convex set gives the value 1. Furthermore, the less the parameter value is high, the less the shape is convex. The convexity measurement can be computed, for instance, by the ratio $c = A/A_c$ where $A_c$ is the area of the convex hull of the set. This convexity parameter $c$ has the following desirable properties:

- its value is always a number within $[0,1]$;

- its value is 1 if and only if the set is convex;

<table><caption>**Table 4.** The twenty two shape diagrams axes coordinates for 2D simply connected compact sets.</caption><thead><tr><th>Shape diagrams</th><th colspan="2">Axes coordinates</th></tr></thead><tbody><tr><td>$D_1 : (\omega, r, R)$</td><td>$x = \omega/2R$</td><td>$y = r/R$</td></tr><tr><td>$D_2 : (\omega, A, R)$</td><td>$x = \omega/2R$</td><td>$y = A/\pi R^2$</td></tr><tr><td>$D_3 : (r, A, R)$</td><td>$x = r/R$</td><td>$y = A/\pi R^2$</td></tr><tr><td>$D_4 : (A, d, R)$</td><td>$x = A/\pi R^2$</td><td>$y = d/2R$</td></tr><tr><td>$D_5 : (\omega, d, R)$</td><td>$x = \omega/2R$</td><td>$y = d/2R$</td></tr><tr><td>$D_6 : (r, d, R)$</td><td>$x = r/R$</td><td>$y = d/2R$</td></tr><tr><td>$D_{11} : (\omega, r, d)$</td><td>$x = \omega/d$</td><td>$y = 2r/d$</td></tr><tr><td>$D_{12} : (\omega, A, d)$</td><td>$x = \omega/d$</td><td>$y = 4A/\pi d^2$</td></tr><tr><td>$D_{13} : (r, A, d)$</td><td>$x = 2r/d$</td><td>$y = 4A/\pi d^2$</td></tr><tr><td>$D_{14} : (A, R, d)$</td><td>$x = 4A/\pi d^2$</td><td>$y = \sqrt{3}R/d$</td></tr><tr><td>$D_{15} : (\omega, R, d)$</td><td>$x = \omega/d$</td><td>$y = \sqrt{3}R/d$</td></tr><tr><td>$D_{16} : (r, R, d)$</td><td>$x = 2r/d$</td><td>$y = \sqrt{3}R/d$</td></tr><tr><td>$D_{21} : (\omega, r; P)$</td><td>$x = \pi\omega/P$</td><td>$y = 2\pi r/P$</td></tr><tr><td>$D_{22} : (\omega, A, P)$</td><td>$x = \pi\omega/P$</td><td>$y = 4\pi A/P^2$</td></tr><tr><td>$D_{23} : (r, A, P)$</td><td>$x = 2\pi r/P$</td><td>$y = 4\pi A/P^2$</td></tr><tr><td>$D_{24} : (A, R, P)$</td><td>$x = 4\pi A/P^2$</td><td>$y = 4R/P$</td></tr><tr><td>$D_{25} : (\omega, R, P)$</td><td>$x = \pi\omega/P$</td><td>$y = 4R/P$</td></tr><tr><td>$D_{26} : (r, R, P)$</td><td>$x = 2\pi r/P$</td><td>$y = 4R/P$</td></tr><tr><td>$D_{27} : (A, d, P)$</td><td>$x = 4\pi A/P^2$</td><td>$y = 2d/P$</td></tr><tr><td>$D_{28} : (\omega, d, P)$</td><td>$x = \pi\omega/P$</td><td>$y = 2d/P$</td></tr><tr><td>$D_{29} : (r, d, P)$</td><td>$x = 2\pi r/P$</td><td>$y = 2d/P$</td></tr><tr><td>$D_{30} : (d, R, P)$</td><td>$x = 2d/P$</td><td>$y = 4R/P$</td></tr></tbody></table>

- it is invariant under similitude transformations;

- there is a simple and fast computing algorithm.

Figure 8 shows this convexity parameter for some specific sets. The colormap corresponds to the different values of this parameter from 0 (black) to 1 (dark red), i.e from non-convex sets to convex sets.
<figure><img src="image_7.png" /><figcaption>**Fig. 5.** Illustration of three specific shape diagrams where 24 elementary sets are mapped. It shows that some diagrams (using some specific morphometrical functionals) are not suitable for shape discrimination and classification. For example, different shapes can be located on the same place within the diagram $D_{12}$, highlighting the non-unicity of such a representation.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 6.** Illustration of the sets (represented by a curve on the shape diagrams) coming from a transformation (using a degree of freedom of the set) from one set to another [12].</figcaption></figure>

Figure 9 illustrates three shape diagram where 1370 various discrete sets are mapped, and on which the boundaries of the convex domain are superimposed (the domain where all compact convex sets are located [11,28]).

One can see that the shape diagram $D_{24}$ presents a good convexity discrimination [31]. The sets mapped at the bottom left (sets with both low roundness values and low thinness values) are strongly concave, whereas the sets mapped at the top (sets with high thinness values) are stongly convex. In the two other shape diagrams, it would be more difficult to discriminate the convexity of the sets since the different color points are much more mixed.

A more detailed study, including other properties of the shape diagrams, can be found in [29–31]. Among the different shape diagrams, it has been shown that the shape diagram $D_{24}:(A, R, P)$ is a good candidate for shape analysis regarding its representation relevance and discrimination power.

# 4 Conclusions

The purpose of this paper was to propose geometrical and morphometrical tools for characterizing inclusions in metal alloys from digital images acquired by various systems such as electron microscopy or X-ray tomography. The objective
<figure><img src="image_9.png" /><figcaption>**Fig. 7.** Illustration of eight continuity curves (representing the locations of the transformed sets from one set to another) onto three shape diagrams [12].</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 8.** Illustration of the convexity parameter for some specific sets. The colormap corresponds to the different values of this parameter from 0 (black) to 1 (dark red), i.e. from non-convex sets to convex sets [12].</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 9.** Family of 1370 discrete sets which are mapped onto three shape diagrams. The colormap corresponds to the different values of the convexity parameter from 0 (black) to 1 (dark red) of the sets, i.e. from non-convex sets to convex sets [12].</figcaption></figure>

was not to give an exhaustive list of geometrical and morphometrical tools, but rather to focus on specific descriptors. In this context, the Minkowski functionals (geometrical tools) and the shape diagrams (morphometrical tools) have been presented. It has been shown that these descriptors have a compact representation, are computationally effective and have some good mathematical properties. For sure, they could be insufficient (as it

stands) to analyse very complex structures but the proposed measurements can be computed on some geometrical transformations of the object to be studied, providing some measurements functions and not only scalars. This kind of approach has been particularly studied in [12]. In conclusion, image analysis using integral geometry and shape diagrams provide efficient tools to quantitatively describe inclusions (providing separate information on size
and shape). The geometrical characteristics of these inclusions could thereafter be related to the mechanical properties of the metal alloys.

## References

1. J. Ohser, F. Mücklich, Statistical analysis of microstructures in materials science, John Wiley and Sons, New York, USA, 2000

2. J. Ohser, K. Schladitz, Image processing and analysis, Clarendon Press Oxford, Oxford, UK, 2006

3. W. Blaschke, Vorlesungen über integralgeometrie, VEB, Berlin, 1955

4. K.R. Mecke, Integral geometry in statistical physics, Int. J. Mod. Phys. B **12**(9), 861–899 (1998)

5. L.A. Santalo, Integral geometry and geometric probability, Cambridge University Press, Cambridge, UK, 2004

6. H. Minkowski, Volumen und oberfläche, Math. Ann. **57**, 447–495 (2003)

7. K.R. Michielsen, H. De Raedt, Integral-geometry morphological image analysis, Phys. Rep. **347**, 461–538 (2001)

8. V. Schulz, Description and reconstruction of microscopic random heterogenous media in order to estimate macroscopic hydraulic functions, PhD thesis, University of Heidelberg, 2003

9. W. Blaschke, Konvexe Bereiche gegebener konstanter Breite und kleinsten Inhalts, Math. Ann. **76**, 504–513 (1915)

10. W. Blaschke, Eine frage über konvexe körper, Jahresbericht Deutsch, Math.-Verein. **25**, 121–125 (1916)

11. L.A. Santalo, Sobre los sistemas completos de desigualdades entre tres elementos de una figura plana, Math. Notae. **17**, 82–104 (1961)

12. S. Rivollier, Analyse d'image géométrique et morphométrique par diagrammes de forme et voisinages adaptatifs généraux, PhD thesis, École Nationale Supérieure des Mines de Saint-Étienne, France, 2010

13. S.N. Chiu, D. Stoyan, W. Kendall, J. Mecke, Stochastic geometry and its applications, John Wiley & Sons, Chester, UK, 2013

14. H. Hadwiger, Vorlesungen über inhalt, oberfläsche und isoperimetrie, Springer-Verlag, Berlin Heidelberg, Germany, 1957

15. A. Rosenfeld, Connectivity in digital pictures, J. ACM **17**(1), 146–160 (1970)

16. A. Rosenfeld, A converse to the jordan curve theorem for digital curves, Inf. Control **29**, 292–293 (1975)

17. W. Nagel, J. Ohser, K. Pischang, An integral-geometric approach for the euler-poincaré characteristic of spatial images, J. Microsc. **198**, 54–62 (2000)

18. C. Lang, J. Ohser, R. Hilfer, On the analysis of spatial binary images, J. Microsc. **203**(3), 303–313 (2001)

19. J. Ohser, W. Nagel, K. Schladitz, The Euler number of discretised sets – On the choice of adjacency in homogeneous lattices, in: Morphology of condensed matter, Springer-Verlag, Berlin Heidelberg, Germany, 2002, pp. 275–798

20. J. Ohser, W. Nagel, K. Schladitz, The euler number of discretised sets – Surprising results in three dimensions, Imag. Anal. Stereol. **22**, 11–19 (2003)

21. K. Schladitz, J. Ohser, W. Nagel, Measuring intrinsic volumes in digital 3d images, Discret. Geom. Comp. Imag. **4245**, 247–258 (2006)

22. K. Sandfort, J. Ohser, Labeling of n-dimensional images with choosable adjacency of the pixels, Imag. Anal. Stereol. **28**, 45–61 (2009)

23. M.W. Crofton, On the theory of local probability, applied to straight lines drawn at random in a plane, Philos. Trans. R. Soc. Lond. **158**, 181–199 (1868)

24. L.R. Feret, La grosseur des grains des matières pulvérulentes, Premières Communications de la Nouvelle Association Internationale pour l'Essai des Matériaux, Groupe D, 1930, pp. 428–436

25. M.I. Shamos, Computational geometry, PhD thesis, Yale University, 1978

26. Y.D. Burago, V.A. Zalgaller, Geometric inequalities, Springer-Verlag, Berlin Heidelberg, Germany, 1988

27. P.R. Scott, P.W. Awyong, Inequalities for convex sets, Inequal. Pure Appl. Math. **1**(1–6), 1–6 (2000)

28. M.A.H. Cifre, G. Salinas, S.S. Gomis, Complete systems of inequalities, Inequal. Pure Appl. Math. **2**(1), 1–12 (2001)

29. S. Rivollier, J. Debayle, J.C. Pinoli, Shape diagrams for 2D compact sets – Part I: analytic convex sets, Aust. J. Math. Anal. Appl. **7**(2–3), 1–27 (2010)

30. S. Rivollier, J. Debayle, J.C. Pinoli, Shape diagrams for 2D compact sets – Part II: analytic simply connected sets, Aust. J. Math. Anal. Appl. **7**(2–4), 1–21 (2010)

31. S. Rivollier, J. Debayle, J.C. Pinoli, Shape diagrams for 2D compact sets – Part III: convexity discrimination for analytic and discretized simply connected sets, Aust. J. Math. Anal. Appl. **7**(2–5), 1–18 (2010)

Cite this article as: Johan Debayle, Geometrical and morphometrical tools for the inclusion analysis of metallic alloys, Metall. Res. Technol. **116**, 508 (2019)