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

HAL
open science

# Fuzzy Color Aura Matrices for Texture Image Segmentation

Zohra Haliche, Kamal Hammouche, Olivier Losson, Ludovic Macaire

► To cite this version:

Zohra Haliche, Kamal Hammouche, Olivier Losson, Ludovic Macaire. Fuzzy Color Aura Matrices for Texture Image Segmentation. Journal of Imaging, 2022, 8 (9), pp.244. 10.3390/jimaging8090244. hal-03773140

HAL Id: hal-03773140

https://hal.science/hal-03773140v1

Submitted on 15 May 2025

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
<figure><img src="image_3.png" /></figure>

Article

# Fuzzy Color Aura Matrices for Texture Image Segmentation

Zohra Haliche ¹ID, Kamal Hammouche ¹ID, Olivier Losson ²ID and Ludovic Macaire <fnref n="2" />,*ID

**Abstract:** Fuzzy gray-level aura matrices have been developed from fuzzy set theory and the aura concept to characterize texture images. They have proven to be powerful descriptors for color texture classification. However, using them for color texture segmentation is difficult because of their high memory and computation requirements. To overcome this problem, we propose to extend fuzzy gray-level aura matrices to fuzzy color aura matrices, which would allow us to apply them to color texture image segmentation. Unlike the marginal approach that requires one fuzzy gray-level aura matrix for each color channel, a single fuzzy color aura matrix is required to locally characterize the interactions between colors of neighboring pixels. Furthermore, all works about fuzzy gray-level aura matrices consider the same neighborhood function for each site. Another contribution of this paper is to define an adaptive neighborhood function based on information about neighboring sites provided by a pre-segmentation method. For this purpose, we propose a modified simple linear iterative clustering algorithm that incorporates a regional feature in order to partition the image into superpixels. All in all, the proposed color texture image segmentation boils down to a superpixel classification using a simple supervised classifier, each superpixel being characterized by a fuzzy color aura matrix. Experimental results on the Prague texture segmentation benchmark show that our method outperforms the classical state-of-the-art supervised segmentation methods and is similar to recent methods based on deep learning.

**Keywords:** color texture segmentation; aura matrices; fuzzy color aura matrix; SLIC superpixel; regional feature

## 1. Introduction

Color texture segmentation consists of partitioning an image into homogeneous regions with respect to color and texture properties. It is involved in various fields, such as medical image analysis [1], remote sensing [2], synthetic aperture radar [3], and fruit detection [4]. Although a wide variety of techniques have been developed, color texture segmentation remains an open and challenging problem due to the high variability of textures, combined with the great diversity of colors. Most approaches characterize each pixel by a set of texture and color features, and use them to assign one color texture class to each pixel thanks to classification algorithms [5,6]. Sections 1.1 and 1.2 discuss the state of the art of color texture features and classification algorithms used for color texture segmentation.

### 1.1. Color Texture Features

The characterization of color textures is a fundamental problem in computer vision, where a texture is generally described by some visual cues represented by statistical structures, hereafter called texture features. Pixel colors are often represented by the values or by statistic measures of color components (R, G, B) or by those derived from color spaces such as HSV, L*u*v*, or L*a*b* [7]. Various texture features are designed to characterize texture appearance. Widely used ones are based on Gabor filters, wavelet transform, Markov random field model, local binary patterns, and co-occurrence matrices [8–10].

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

### check for updates

**Citation:** Haliche, Z.; Hammouche, K.; Losson, O.; Macaire, L. Fuzzy Color Aura Matrices for Texture Image Segmentation. *J. Imaging* **2022**, *8*, 244. https://doi.org/10.3390/jimaging8090244

Academic Editor: Lucia Mancini

Received: 13 July 2022

Accepted: 2 September 2022

Published: 8 September 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ Laboratoire Vision Artificielle et Automatique des Systèmes, Université Mouloud Mammeri, Tizi-Ouzou 15000, Algeria</footnote>

<footnote n="2">² Univ. Lille, CNRS, Centrale Lille, UMR 9189 CRISTAL, F-59000 Lille, France</footnote>

<footnote>* Correspondence: ludovic.macaire@univ-lille.fr</footnote>
Color texture classification/segmentation techniques can be broadly classified into two approaches according to how they combine color texture and color features [5,8]. In the first approach, color and texture features are computed separately and are then combined by the clustering or classification process. For example, in [11], L\*, u\*, and v\* components were used as color features and the output of Gabor filters as texture features. In [12], the authors used the mean, standard deviation, and skewness of each color channel, H, S, and V, as color features, and the co-occurrence matrices computed on the intensity image were considered as texture features.

In the second approach, color and texture features are assumed to be mutually dependent. With this approach, marginal, opponent, and vector strategies have been developed to compute color texture features [5,8]. The marginal strategy assumes that texture can be separately described within each color channel, where pixels are characterized by only one color component. Texture features designed for gray-level images are then computed for each color channel and aggregated into a global feature representation. In [13], for instance, texture features were extracted from wavelet transform and co-occurrence matrices for each color channel. In [14], Leung-Malik filter banks were applied on each color channel for color texture feature extraction. The opponent strategy extends texture feature extraction to color images thanks to both within- and between-channel analyses. The within-channel part consists of computing features from each color channel separately (as in the marginal strategy), whereas between-channel (opponent) features are obtained by processing pairs of color channels jointly. This strategy was adopted by Panjwani et al. [15], who used Markov random field models to characterize a texture in terms of spatial interaction within each color channel and interaction between different color channels. The vector strategy takes advantage of color vector information. It makes it possible to analyze relationships between colors of neighboring pixels, as in the case of gray-level images. Examples of this strategy can be found in [9,16] where local binary patterns were computed from a color image. It should be noted that the vector strategy is more suitable to characterizing color texture because it is less memory-consuming and fully takes correlation between colors of neighboring pixels into account.

## 1.2. Color Texture Image Segmentation by Pixel Classification

Color texture segmentation methods can also be categorized into three ways, depending on whether the pixel classification is performed in a supervised, unsupervised, or semi-supervised manner [1,11,17]. In supervised segmentation, prior knowledge about the training samples and their class labels is needed to classify the input image pixels. Some classical supervised classifiers are used for texture segmentation [1,4,17]. Among them are the K-nearest neighbor (KNN) [18] and the Bayesian [19] classifiers, the support vector machine (SVM) [20], random forest [12], Markov random field [11], and neural networks [14]. Supervised color texture segmentation based on deep learning has been developed in the last decade [21–23]. These efficient methods generally use convolutional neural networks (CNN), such as the U-Net [24], deep visual model (DA) [25], pyramid scene parsing network (PSP-Net) [26], supervised fully convolutional network for texture (FCNT) [21], and empirical wavelet transform-based fully convolutional network for texture (EWT-FCNT) [23].

Unsupervised segmentation does not require prior knowledge and discovers different classes by clustering pixels from their features only [5,6]. For example, the popular K-means clustering algorithm is mainly employed to perform the classification of pixels [6,27].

Other color texture segmentation methods rely on graph cut techniques [28], mean shift clustering [29], estimation-maximization [1], the random Markov model [11], dictionary learning and sparse representation-based classification (DLSRC) [13], and spectral clustering [3].

Semi-supervised segmentation is suitable where only partial prior knowledge about the training samples and their class labels is available. For instance, in [17], the constrained spectral clustering algorithm was applied for semi-supervised classification of pixels characterized by color texture features.
Supervised segmentation methods, especially deep learning based methods, provide better results than unsupervised and semi-supervised approaches thanks to prior information [17]. It should also be noted that some of the above methods perform segmentation at the pixel level and others at the "superpixel" level. This is the case for the methods discussed in [1,3]. Recall that superpixels group neighboring pixels based on their spatial and color similarity, and are used as samples in order to speed up color texture segmentation.

### 1.3. Fuzzy Color Texture Features

All aforementioned methods assume that images are crisp and free from vagueness. However, in practice, color images carry more or less fuzziness because of the imprecision inherent to the discretization of both the spatial domain (sampling) and the color component levels (quantization). Therefore, boundaries separating the various image regions are not precisely defined, and pixel levels are imprecise measures of the reflectance of surfaces observed by the camera. Furthermore, the assumption that texture images are mainly represented by spatial repetitions of a pattern may not be valid any longer. Texture analysis techniques based on fuzzy concepts have then been proposed in order to take this imprecision into account.

For instance, fuzzy histograms [30], fuzzy local binary patterns [31], and local fuzzy patterns [32] are extracted from gray-level texture images. Fuzzy gray-level co-occurrence matrices (FGLCMs) are also proposed to characterize spatial interactions between gray-levels of neighboring pixels [33–37]. However, characterizing a color texture by at least three FGLCMs (one by channel) is memory expensive. To overcome this drawback, Ledoux et al. [38] proposed to extend FGLCMs to color images by defining fuzzy color sets. A color image is then represented by one single fuzzy color co-occurrence matrix (FCCM) that characterizes the local interactions between colors of neighboring pixels. Moreover, FGLCMs and FCCMs only consider spatially-invariant neighborhoods. Hammouche et al. [39] showed that adaptive neighborhoods are useful for texture analysis and provided an elegant formalism to deal with spatially-variant neighborhoods thanks to the aura concept.

In the framework based on the aura set concept, Elfadel and Picard [40] proposed a generalization of GLCMs called gray-level aura matrices (GLAMs). A GLAM quantifies the presence of a set of pixels with a specified level in the neighborhood of another set of pixels having another level. The amount of neighboring pixels with the specified level is quantified by means of the aura measure. GLAMs are used for texture representation and synthesis [41,42], image retrieval [43,44], classification [45–50], and segmentation [51–53]. A generalization of GLAMs to the fuzzy framework has been proposed by Hammouche et al. [39]. Representing each color channel by a fuzzy GLAM (FGLAM) outperforms the FGLCM representation for texture classification. Recently, FGLAMs have been used to improve the accuracy of wood species classification [54]. However, as for FGLCMs, the computation of FGLAMs for each pixel is memory and time expensive, which makes their use for color image segmentation a challenge.

To circumvent these constraints, we adopted a vector strategy and propose to extend FGLAMs to fuzzy color aura matrices (FCAMs). An FCAM makes it possible to locally characterize the interactions between colors of neighboring pixels. While one FGLAM must be computed for each color channel, a single low-dimensionnal FCAM is required to describe the color texture.

### 1.4. FCAM for Image Segmentation by Superpixel Classification

In this study, we applied FCAM to color texture image segmentation. As FCAM is based on a locally-adaptive neighborhood, it can characterize the texture represented by small connected pixel subsets with different shapes, i.e., by superpixels. We use the simple linear iterative clustering (SLIC) scheme [55] to generate superpixels from a color image. These are then classified using a simple supervised classifier to segment the color texture image.
The remainder of this paper is organized as follows. In Section 2, we first give an overview of basic SLIC algorithm; then we propose a modified version that incorporates the regional information. Section 3 introduces some definitions of the fuzzy color aura concept, and explains how to characterize the color texture of each superpixel by an FCAM. In Section 4, we first present the datasets used in the experiments; then we give details about the proposed color texture segmentation method based on FCAMs. Next, we assess the regional SLIC algorithm and discuss its parameter settings. Last, we compare the segmentation results achieved by our supervised segmentation approach with those obtained by fuzzy texture features and by several other state-of-the-art color texture segmentation methods. Concluding remarks about the contribution of this paper are given in Section 5.

## 2. Superpixel

The proposed color texture image segmentation is based on the classification of superpixels. A superpixel is a compact set of connected sites with similar properties. It is usually adopted to replace the pixel grid in an image in order to reduce the computational burden of subsequent processing. Superpixels are generated by over-segmentation algorithms using color, spatial, and/or texture information. Among superpixel generation algorithms, simple linear iterative clustering (SLIC) is widely used due to its simplicity, speed, and ability to adhere to image boundaries. In this section, we briefly give an overview of SLIC algorithm used to generate superpixels. Then, we propose a modified version that incorporates the regional information.

### 2.1. Basic SLIC

The SLIC algorithm [55] generates a desired number of regular superpixels, with la ow computational overhead, by clustering sites based on their spatial and color features. Let $I$ be an RGB image defined on a lattice $S$, such that each site $s \in S$ is characterized by three color components: $I(s) = (I^R(s), I^G(s), I^B(s))^{\top}$. The RGB color components are transformed into the $La^*b^*$ color space, so that SLIC represents each site $s(x_s, y_s)$ by a five-dimensional feature vector: $(I^L(s), I^a(s), I^b(s), x_s, y_s)^{\top}$. SLIC follows a $k$-means clustering strategy but searches the nearest cluster center according to the distance $D(s, s') = \sqrt{d_c^2(s, s') + m^2 \cdot d_s^2(s, s')}/S^2$ between two sites $s$ and $s'$, where $d_c^2(s, s') = (I^L(s) - I^L(s'))^2 + (I^a(s) - I^a(s'))^2 + (I^b(s) - I^b(s'))^2$ and $d_s^2(s, s') = (x_s - x_{s'})^2 + (y_s - y_{s'})^2$ measure the color and spatial proximity. The compactness parameter $m$ is set to its default value of 1, and the maximum spatial distance within a cluster is defined as the sampling step $S = \sqrt{|S|/P}$, where $|S|$ is the total number of sites in the image and $P$ is the number of superpixels. Figure 1b shows the results of the segmentation achieved by SLIC on the color image of Figure 1a. It should be emphasized that SLIC takes a post-processing step to enforce connectivity by merging small isolated superpixels with nearby larger ones. Therefore, the actual number of superpixels produced by SLIC can be slightly lower than the desired number of superpixels.

<figure><img src="image_6.png" /><figcaption>**Figure 1.** SLIC pre-segmentation results ($P$ = 400 superpixels). $f_R$ in (c) is rescaled to $[0, 255]$. SLIC: simple linear iterative clustering.</figcaption></figure>
## 2.2. Regional SLIC

The basic SLIC algorithm achieves good pre-segmentation on color images, but may fail to find boundaries of textures. Indeed, it successfully detects homogeneous regions but somewhat fails to separate areas with different textures (see Figure 1a and arrows on Figure 1b). To improve this pre-segmentation, we propose a modified version of the algorithm, referred to as regional SLIC, that takes color, spatial, and regional information into account. From the basic SLIC pre-segmentation result, we compute the regional feature $f_R$ of each superpixel R as [56]:

$$f_R = [p_R - p_R \cdot \log(p_R)] + \frac{1}{|\mathcal{N}_R|} \sum_{R' \in \mathcal{N}_R} [p_{R'} - p_{R'} \cdot \log(p_{R'})], \quad (1)$$

where $p_R = |R|/|S|$ is the area ratio of R to the whole image, and $\mathcal{N}_R$ is the set of superpixels that are adjacent to R. The regional feature $f_R$ is a sum of two terms: the first one reflects the size of the superpixel R and monotonously increases with its size (because $p_R \in [0, 1]$), and the second one reflects the superpixel context by taking the influence of adjacent superpixels into account. A small superpixel surrounded by small superpixels provides a low value of $f_R$, and conversely, if both R and adjacent superpixels are large, $f_R$ is high (see Figure 1c). At each site s of superpixel R, we replace the luminance component $I^L(s)$ by the regional feature $f_R$ and we apply the SLIC algorithm again with $(f_R, I^a(s), I^b(s), x_s, y_s)^T$ as the feature vector (see Algorithm 1).

Thanks to this regional SLIC method, the lattice S is partitioned into P superpixels $\{P^p\}_{p=1}^P$, i.e., $\cup_{p=1}^P P^p = S$ and $P^p \cap P^{p'} = \emptyset$, for $p \neq p'$. A superpixel P^p is defined as a set of connected sites, i.e., any two sites s and s' in P^p are connected by at least one path composed of sites in P^p. Figure 1d shows the final superpixels extracted from the color texture image of Figure 1a. Note (see arrows) that areas with different textures are better delineated with the regional SLIC than with its basic version, despite the regional feature describing the superpixel contextual information in the pre-segmented image and not being able to be considered as a texture feature.

**Algorithm 1** Regional SLIC.

**Input:** RGB image I

1. Convert I from RGB to $La^*b^*$ color space;

2. Apply SLIC to S where each site s is characterized by $(I^L(s), I^a(s), I^b(s), x_s, y_s)^T$;

3. Extract the regional feature $f_R$ from each superpixel R of the SLIC map using Equation (1);

4. Apply SLIC to S where each site s is characterized by $(f_R, I^a(s), I^b(s), x_s, y_s)^T$.

**Output:** Partition of I into superpixels $\{P^p\}_{p=1}^P$

## 3. Fuzzy Color Aura

In this section, we extend the aura concept introduced by Elfadel and Picard for gray-level images [40] to color images. We explain how it is used to characterize a superpixel thanks to a matrix of aura cardinals. As this matrix is huge when all RGB colors are considered, we propose to extend the color aura concept to fuzzy colors whose number is reduced.

### 3.1. Color Aura Set

Let I be an RGB image defined on a lattice S, such that each site s ∈ S is characterized by its color I(s). For a given color x ∈ RGB, we define the set of sites with this color as $S_x = \{s \in S, I(s) = x\}$. Then, $\{S_x\}_{x \in RGB}$ is a partition of S—i.e., $\cup_{x \in RGB} S_x = S$ and $S_x \cap S_{x'} = \emptyset$ for $x \neq x'$.
Given two colors $(x, x') \in \text{RGB}^2$, we define the color aura set $A_{S_{x'}}(S_x)$ of $S_x$ with respect to $S_{x'}$ as:

$$A_{S_{x'}}(S_x) = \bigcup_{s \in S} (N_{S_x}^s \cap S_{x'}), \qquad (2)$$

where $N_{S_x}^s$ is the neighboring site set of each site $s \in S_x$:

$$N_{S_x}^s = \{ \mathbf{r} \in S, (\mathbf{s} \in S_x) \wedge (\mathbf{r} \in N^s) \}. \qquad (3)$$

$N^s = \{ \mathbf{r} \in S, \|\mathbf{r} - \mathbf{s}\|_{\infty} \le d \}$ is the neighborhood of a site $\mathbf{s}$. In this paper, $d = 1$ such that $N^s$ is only composed of the eight closest sites from $\mathbf{s}$.

$A_{S_{x'}}(S_x)$ is the subset of $S_{x'}$ composed of the sites that are present in the neighborhood of those of $S_x$. It provides an interpretation of the presence of $S_{x'}$ in the neighborhood of $S_x$. The color aura set is a generalization of the gray-level aura set introduced by Elfadel and Picard [40] to the color case.

Figure 2a shows a color image, defined on a lattice $S$ of $7 \times 7$ pixels, and three color sets, $S_R$, $S_G$, and $S_B$. For the neighborhood $N^s$, the aura set $A_{S_G}(S_G)$ of $S_B$ with respect to $S_G$ is composed of the green sites marked as circles. Note that it differs from the aura set $A_{S_B}(S_G)$ of $S_G$ with respect to $S_B$, composed of the blue sites marked as diamonds.

<figure><img src="image_7.png" /><figcaption>Figure 2. Examples of two aura sets in a color image. (a) Original image composed of 3 color site sets, $S_R$, $S_G$, and $S_B$, and aura sets $A_{S_G}(S_G)$ (circles) and $A_{S_B}(S_G)$ (diamonds). (b) Superpixels $P^1$ (left) and $P^2$ (right), and aura sets $A_{S_G}^1(S_G)$ (①) and $A_{S_B}^1(S_G)$ (②). $A_{S_G}^2(S_G)$ and $A_{S_B}^2(S_G)$ are empty. (c) Superpixels $P^1$ (left) and $P^2$ (right), and aura sets $A_{S_G}^2(S_G)$ (②) and $A_{S_B}^2(S_G)$ (③). $A_{S_G}^1(S_G)$ and $A_{S_B}^1(S_G)$ are empty.</figcaption></figure>

### 3.2. Color Aura Set in a Superpixel

Rather than building an aura set from all image sites, we focus on a superpixel $P^p$ and define the color aura set $A_{S_{x'}}^p(S_x)$ of $S_x$ with respect to $S_{x'}$ for the superpixel $P^p$ as:

$$A_{S_{x'}}^p(S_x) = \bigcup_{s \in S} (N_{S_x}^{p,s} \cap S_{x'}), \qquad (4)$$

where $N_{S_x}^{p,s}$ is the neighboring site set of each site $s \in S_x$ within the superpixel $P^p$:

$$N_{S_x}^{p,s} = \{ \mathbf{r} \in S, (\mathbf{s} \in S_x \cap P^p) \wedge (\mathbf{r} \in N^s \cap P^p) \}. \qquad (5)$$

The image of Figure 2a contains two textures, one represented by red and green vertical stripes on the left three columns, and another represented by red and blue horizontal stripes on the right three columns. These two textures are separated by the fourth column composed of red, green, and blue sites. Figure 2b shows the partitioning of this image into two superpixels delimited by a vertical black line. $P^1$ covers the four left columns and $P^2$ the three right ones. The aura sets $A_{S_G}^1(S_G)$ and $A_{S_B}^1(S_G)$ are not empty, since $P^1$ contains neighboring green and blue sites, but $A_{S_G}^2(S_G)$ and $A_{S_B}^2(S_G)$ are empty because $P^2$ contains
no green site. Figure 2c shows another partition where $P^1$ covers the three left columns and $P^2$ the four right ones. In that case, $A_{S_G}^1(S_B)$ and $A_{S_B}^1(S_G)$ are empty because $P^1$ contains no blue site. This example illustrates that color aura sets depend on superpixel edges. Indeed, along the fourth column that separates the two textures, only one green site (on first row) and one blue site (second row) belong to the aura sets of superpixels in both Figure 2b,c.

### 3.3. Color Aura Cardinal

The aura measure was introduced by Elfadel and Picard [40] to characterize an aura set by a number that expresses the amount of mixing between neighboring site sets. We use here a simpler measure [39] and quantify the color aura set of a color site set $S_x$ with respect to another color site set $S_{x'}$ within the superpixel $P^p$ thanks to its cardinal defined as:

$$m^p(\mathbf{x}, \mathbf{x}') = |A_{S_{x'}}^p(S_x)|. \qquad (6)$$

The aura cardinal measures for all the possible pairs of color sets, $(S_x, S_{x'})$, $(\mathbf{x}, \mathbf{x}') \in \text{RGB}^2$, are gathered in a matrix. This color aura cardinal matrix $[m^p(\mathbf{x}, \mathbf{x}')]$ can be then considered as the texture feature of the superpixel $P^p$. However, when color components of the image are defined on 256 levels, as is classically the case, the number of possible colors reaches $256^3$ and $m^p$ would be of size $256^3 \times 256^3$. As such, a huge memory requirement is needed in practice, so we propose to decrease the number of analyzed colors by introducing fuzzy colors.

### 3.4. Fuzzy Color

To form a small subset $\mathcal{C}$ of C colors among the $256^3$ possible ones, we use the uniform quantization technique for its simplicity of implementation. For each color component $k \in \{R, G, B\}$, the full level range $[0, 255]$ is divided into $C^k$ disjoint intervals $[0, L^k - 1]$, $[L^k, 2L^k - 1]$, ..., $[(C^k - 1)L^k, 255]$, of respective width $L^k \doteq 256/C^k$. The $C^k$ centers $\{[(L^k - 1)/2], [(3L^k - 1)/2], \dots, [256 - (L^k + 1)/2]\}$ of the intervals define the k-th component of the colors in $\mathcal{C}$. The numbers $C^R, C^G$, and $C^B$ are chosen so that the number of colors $C = C^R \cdot C^G \cdot C^B$ is much lower than $256^3$.

In the fuzzy framework [38], a fuzzy color $\tilde{c}$ is characterized by its membership function $\mu_{\tilde{c}}: \text{RGB} \rightarrow [0, 1]$. The membership degree $\mu_{\tilde{c}}(\mathbf{x})$ of any color $\mathbf{x} \in \text{RGB}$ is defined using its infinity norm or Euclidean distance to the crisp counterpart $\mathbf{c} \in \mathcal{C}$ of $\tilde{c}$, thanks to either:

* the crisp membership function:

$$\mu_{\tilde{c}}(\mathbf{x}) = \begin{cases} 1 & \text{if } \|\mathbf{x} - \mathbf{c}\|_{\infty} \le \lfloor \frac{L^k}{2} \rfloor, \\ 0 & \text{otherwise,} \end{cases} \qquad (7)$$

* the symmetrical Gaussian function:

$$\mu_{\tilde{c}}(\mathbf{x}) = \exp\left(-\frac{\|\mathbf{x} - \mathbf{c}\|_2^2}{2\alpha^2}\right), \qquad (8)$$

* the triangular function:

$$\mu_{\tilde{c}}(\mathbf{x}) = \max\left(1 - \frac{\|\mathbf{x} - \mathbf{c}\|_2}{\beta}, 0\right), \qquad (9)$$

* or the fuzzy C-means (FCM) membership function:

$$\mu_{\tilde{c}}(\mathbf{x}) = \frac{1}{\sum_{\mathbf{c}' \in \mathcal{C}} \left( \frac{\|\mathbf{x} - \mathbf{c}\|_2}{\|\mathbf{x} - \mathbf{c}'\|_2} \right)^2} \cdot \qquad (10)$$
Here, $\alpha$ and $\beta$ are real positive constants used to control the span of the fuzzy color, and $\zeta$ is any real number greater than 1. In this paper, we set $\beta = (L^R, L^G, L^B)$, $\alpha = \beta / \sqrt{2\ln(2)}$, and $\zeta = 2$. The parameters $\alpha$ and $\beta$ are chosen to ensure that $\mu_{\tilde{\epsilon}}(x) = 0.5$ at the bounds of each color domain of width $(L^R, L^G, L^B)^T$ centered at c. Figure 3 shows the shapes of the four membership functions computed at c = (192, 64, 64), for plane (R, G).

<figure><img src="image_8.png" /><figcaption>Figure 3. Membership functions $\mu_{\tilde{\epsilon}}$ at c = (192, 64, 64) with $C^R = C^G = C^B = 2$.</figcaption></figure>

### 3.5. Fuzzy Color Aura Set in a Superpixel

A fuzzy color site set $S_c, c \in \mathcal{C}$, is defined by its membership function $\mu_{S_c}$. The membership degree $\mu_{S_c}(s)$ of each site $s \in S$ to $S_c$ is the membership degree $\mu_{\tilde{\epsilon}}(\mathbf{I}(s))$ of its color $\mathbf{I}(s) \in \text{RGB}$ to the fuzzy color $\tilde{\epsilon}$ [38]:

$$ \mu_{S_c}(\mathbf{s}) = \mu_{\tilde{\epsilon}}(\mathbf{I}(\mathbf{s})). \qquad (11) $$

From there, we define the fuzzy color aura set $A_{S_c'}^p(S_c)$ of $S_c$ with respect to $S_c'$, for any color pair $(c, c') \in \mathcal{C}^2$, as the fuzzy site set with the following membership degree at each site $\mathbf{r} \in S$:

$$ \mu_{A_{S_c'}^p(S_c)}(\mathbf{r}) = \min\left\{\sup_{\mathbf{s} \in S}[\min(\mu_{S_c}(\mathbf{s}), n^p(\mathbf{r}, \mathbf{s}))], \mu_{S_c'}(\mathbf{r})\right\}. \qquad (12) $$

The neighborhood function $n^p(\mathbf{r}, \mathbf{s})$ expresses the membership degree of $\mathbf{r}$ to the neighborhood of $\mathbf{s}$ within the superpixel $P^p, p \in [1, P]$. This function may have any support size and shape, and may take any real value between 0 and 1. In the simplest binary case, we design it as:

$$ n^p(\mathbf{r}, \mathbf{s}) = \begin{cases} 1 & \text{if } (\mathbf{s} \in P^p) \land (\mathbf{r} \in N^s \cap P^p), \\ 0 & \text{otherwise.} \end{cases} \qquad (13) $$

As a justification of Equation (12), we consider the fuzzy color aura set by analogy with the crisp case [39]. The crisp set union operator $\cup$ of Equation (4) is transcribed by
the fuzzy operator sup to get the membership degree of the fuzzy color aura set S<sub>c</sub> with respect to S<sub>c'</sub> at each site r ∈ S:

$$ \mu_{A_{S_c'}^p(S_c)}(\mathbf{r}) = \sup_{\mathbf{s} \in S} \left( \mu_{N_{S_c}^{p,s} \cap S_{c'}}(\mathbf{r}) \right). \quad (14) $$

The crisp set intersection operator ∩ is transcribed by the fuzzy operator min, such that the fuzzy counterpart of the color site set N<sub>S<sub>x</sub></sub><sup>p,s</sup> ∩ S<sub>x'</sub> of Equation (4) is defined by its membership function given by:

$$ \mu_{N_{S_c}^{p,s} \cap S_{c'}}(\mathbf{r}) = \min \left[ \mu_{N_{S_c}^{p,s}}(\mathbf{r}), \mu_{S_{c'}}(\mathbf{r}) \right]. \quad (15) $$

The fuzzy counterpart of N<sub>S<sub>x</sub></sub><sup>p,s</sup> (see Equation (5)) is defined for the fuzzy color site set S<sub>c</sub>, c ∈ C, by the following membership degree at each site r ∈ S:

$$ \mu_{N_{S_c}^{p,s}}(\mathbf{r}) = \min(\mu_{S_c}(\mathbf{s}), n^p(\mathbf{r}, \mathbf{s})). \quad (16) $$

Plugging Equation (16) into (15) and the result into (14) provides the definition (12) of the fuzzy aura set after swapping the first two operators.

### 3.6. Fuzzy Color Aura Cardinal

The fuzzy color aura cardinal of S<sub>c</sub> with respect to S<sub>c'</sub> within P<sup>p</sup> is exactly defined as in the crisp case (see Equation (6)) and directly follows from Equation (12):

$$ \tilde{m}^p(\mathbf{c}, \mathbf{c}') = \sum_{\mathbf{r} \in P^p} \mu_{A_{S_c'}^p(S_c)}(\mathbf{r}) = \sum_{\mathbf{r} \in P^p} \min \left\{ \sup_{\mathbf{s} \in P^p} \left[ \min(\mu_{S_c}(\mathbf{s}), n^p(\mathbf{r}, \mathbf{s})) \right], \mu_{S_{c'}}(\mathbf{r}) \right\}. \quad (17) $$

Following the crisp scheme, we define a fuzzy color aura matrix (FCAM) as the collection of all fuzzy color aura cardinals of a superpixel for the C<sup>2</sup> color pairs (c, c') ∈ C<sup>2</sup>. Note that to define the FCAM, we only consider a few of the 256<sup>3</sup> colors on which the RGB image is defined—namely, the set C of C colors such that C ≪ 256<sup>3</sup>. This provides a compact FCAM of small size C × C that is a suitable color texture descriptor of a superpixel for image segmentation purpose.

To make superpixels of different sizes comparable, their FCAMs are normalized element-wise to sum up to one:

$$ \tilde{m}^p(\mathbf{c}, \mathbf{c}') = \frac{\tilde{m}^p(\mathbf{c}, \mathbf{c}')}{\sum_{(\mathbf{y}, \mathbf{y}') \in \mathcal{C}^2} \tilde{m}^p(\mathbf{y}, \mathbf{y}')}. \quad (18) $$

Finally, each superpixel P<sup>p</sup> of an image I is characterized by C<sup>2</sup> features that are the elements of its normalized FCAM m<sup>p</sup>.

## 4. Experiments

In this section, we first present the experimental dataset and explain how FCAM features are used to segment its color texture images. To evaluate the performance of the proposed approach, we then assess the regional SLIC algorithm, discuss parameter settings, and study the relevance of FCAM features. We finally compare the segmentation results achieved by our supervised segmentation approach with those obtained by several state-of-the-art color texture segmentation methods.

### 4.1. Experimental Setup

#### 4.1.1. Dataset

As the experimental dataset, we used the challenging Prague texture segmentation benchmark [57]. It contains 20 color texture mosaics to be segmented (input test images), some of which are shown in Figure 4 (top) with their corresponding segmentation ground
truth images. Each of these images represents from 3 to 12 classes and has been synthetically generated from the same number of original color texture images. The original Prague dataset is composed of 89 images (one image per texture class) grouped into 10 categories of natural and artificial textures. Figure 4 (bottom) shows some of these images that form the training dataset. All images of Prague database are of size 512 × 512 pixels.

<figure><img src="image_9.png" /><figcaption>**Figure 4.** Examples of test images (**top** two rows) and training images (from “flowers” and “made” categories, **bottom** two rows) from the Prague dataset. Note that the framed test image represents $K_I = 6$ classes from the sole “flowers” category.</figcaption></figure>

### 4.1.2. Color Texture Image Segmentation Based on FCAMs

Our method of color texture image segmentation based on FCAMs is a supervised superpixel classification procedure. Each test image **I** represents $K_I$ classes whose pixels should be retrieved by the segmentation. To this end, we used a simple artificial neural network known as the extreme learning machine (ELM) that comes with a very fast learning algorithm [58]. It consists of three fully-connected neuron layers: an input layer of $C^2$ neurons that receives FCAM elements, a single hidden layer, and an output layer with $K_I$ neurons. The number of neurons in the hidden layer was empirically set to $100 \cdot K_I$. The initial weights of hidden neurons were set to random values, and the output weights were determined according to a least-square solution [58].

The proposed color texture segmentation of any test image **I** based on FCAMs follows the two successive stages outlined in Figure 5. The ELM is first trained with a set of $K_I \cdot T$ training samples, where $T$ is the number of prototype sites per class that are randomly selected on each of the $K_I$ training images. Each training sample is the FCAM feature computed over a square patch $W^t$, $t \in \{1, ..., T\}$, centered at a prototype site and of size $(2W + 1)^2$. The trained ELM is then used to segment **I** as follows. First, **I** is segmented into superpixels using the regional SLIC method. Then, the FCAM of each superpixel is fed as input into the trained ELM, whose output provides the estimated texture class. All sites in the superpixel are finally assigned to this class. A refinement procedure can be applied to further improve segmentation accuracy. This consists of reassigning superpixels smaller than 0.5% of the image size with the label of the largest adjacent superpixel [21]. All steps of the proposed color texture image segmentation are summarized in Algorithm 2.
<figure><img src="image_10.png" /><figcaption>**Figure 5.** Flowchart of the proposed color texture image segmentation (SLIC: simple linear iterative clustering, FCAM: fuzzy color aura maxtrix, ELM: extreme learning machine).</figcaption></figure>

**Algorithm 2** Color texture image segmentation.

**Input:** Test image I, $K_I$ training images

**Parameters:** Number $T$ of prototypes per class, number $P$ of superpixels, number $C$ of fuzzy colors, patch half width $W$, membership function $\mu_{\tilde{\epsilon}}$.

**Step 1: Training stage**

1. In each training image, randomly select $T$ prototype sites.

2. At each prototype site, compute the normalized FCAM $\tilde{m}^t$ over a square patch $W^t$ of size $(2W+1)^2$ using Equation (18).

3. Train the ELM classifier with the $K_I \cdot T$ normalized FCAMs of prototypes.

**Step 2: Segmentation stage**

1. Run regional SLIC (Algorithm 1) on I to provide $P$ superpixels $\{P^p\}_{p=1}^P$.

2. Compute the normalized FCAM $\tilde{m}^p$ of each superpixel $P^p$ using Equation (18).

3. Feed $\tilde{m}^p$ into the trained ELM and assign each pixel in $P^p$ to the ELM output class.

**Step 3: Refinement (optional)**

**Output:** Segmented image

#### 4.1.3. Regional SLIC—Preliminary Assessment

To demonstrate the relevance of the proposed regional SLIC algorithm, we compare its results with those achieved by the basic SLIC algorithm thanks to four standard metrics. The achievable segmentation accuracy (ASA) quantifies the segmentation performance achievable by assigning each superpixel to the ground truth region with the highest overlap [55,59]. The boundary recall (BR) assesses the boundary adherence with respect to the ground truth boundaries [60]. The under-segmentation error (UE) evaluates the segmentation boundary accuracy as the overlap between superpixels and ground truth regions [55,59]. The compactness (COM) measures the compactness of superpixels [60]. Higher BR, ASA, and COM, and smaller UE, indicate better pre-segmentation.

Both the basic and regional SLIC algorithms require one to set the number of super-pixels. This number must be small enough to get superpixels that are large enough to contain homogeneous textures, and small enough for the superpixels to finely fit boundaries between textures. We empirically found that $P = 400$ is a good trade-off for the Prague dataset. Table 1 shows the average and standard deviation of the four metrics obtained by the basic and modified SLIC algorithms over the 20 images of Prague dataset. From this table, we can see that the results achieved by the regional SLIC are better than those obtained by the basic SLIC according to all metrics.
<table><caption>**Table 1.** Basic and regional SLIC performances on the Prague dataset. Arrows denote the metric direction (the higher ↑ or smaller ↓, the better), and best results are in bold.</caption><thead><tr><th>Metric</th><th>↑BR</th><th>↓UE</th><th>↑ASA</th><th>↑COM</th></tr></thead><tbody><tr><td>Basic SLIC</td><td>0.46 ± 0.15</td><td>0.36 ± 0.08</td><td>0.79 ± 0.05</td><td>0.56 ± 0.05</td></tr><tr><td>Regional SLIC</td><td><b>0.52 ± 0.13</b></td><td><b>0.32 ± 0.08</b></td><td><b>0.82 ± 0.05</b></td><td><b>0.63 ± 0.04</b></td></tr></tbody></table>

### 4.1.4. Parameter Settings

The proposed color texture segmentation requires one to set several parameters. For the training stage, the number of prototypes was set to $T = 1000$, since CNN-based methods use 1000 training texture mosaic images for each image to segment [23]. The size $(2W + 1)^2$ of the patch centered at each prototype site and used for FCAM computation during the training stage was set according to the image size $N$ and the number $P$ of superpixels such that $(2W + 1)^2 \approx N/P$. As $N = 512 \times 512$ and $P = 400$, the patch size was set to $27 \times 27$ pixels.

To characterize a site by an FCAM in the segmentation stage, we have to set both the number $C$ of fuzzy colors and the membership function $\mu_{\tilde{c}}$ that defines the membership degree $\mu_{\tilde{c}}(s)$ of each site $s$ to any fuzzy color $\tilde{c}$ according to its color $I(s)$. We consider very few colors in order to evaluate how the memory cost of FCAMs can be reduced while preserving their relevance as texture features. Specifically, the number of fuzzy colors one of these values each time: $C = 2 \times 2 \times 1$, $C = 2 \times 2 \times 2$, $C = 2 \times 3 \times 2$, $C = 2 \times 4 \times 2$, $C = 3 \times 3 \times 2$, $C = 3 \times 4 \times 2$, or $C = 4 \times 4 \times 2$, so that the FCAM size was $4 \times 4$, $8 \times 8$, $12 \times 12$, $16 \times 16$, $18 \times 18$, $24 \times 24$, or $32 \times 32$. Note that we privileged the G color component over $R$ and $B$ because it is similar to luminance; other combinations provide close classification accuracy results. Regarding the membership function, we had to choose one among the four functions (crisp, Gaussian, triangular, or FCM) presented in Section 3.4.

We retained a given number of fuzzy colors and membership function on the grounds of segmentation accuracy. To make this choice independent of the classification algorithm, we used the nearest neighbor classifier (1NN) instead of the ELM, and with no refinement step. Figure 6 shows the average accuracy obtained with each membership function according to the number $C$ of fuzzy colors over the 20 images. From this figure, we can see that the fuzzy membership functions (Gaussian, triangular, and FCM) largely outperform the crisp one, especially when the number of fuzzy colors is low. However, when $C \ge 12$, all fuzzy membership functions lead to similar segmentation accuracies, even if the Gaussian membership function seems to perform slightly better. Moreover, accuracies do not vary significantly beyond $C = 16$. In the following, all experiments were therefore performed with $C = 16$ fuzzy colors and the Gaussian membership function.

<figure><img src="image_11.png" /><figcaption>**Figure 6.** Average accuracy (%) vs. number of colors ($C$) for different membership functions.</figcaption></figure>
## 4.2. Comparison with Other Fuzzy Texture Features

We here compare the relevance of FCAMs with respect to other fuzzy texture features, namely, the fuzzy gray-level co-occurrence matrix (FGLCM) [33], the fuzzy color co-occurrence matrix (FCCM) [38], and the fuzzy gray-level aura cardinal matrix (FGLAM) [39]. Note that FGLCM and FCCM are similar to a fuzzy color aura matrix computed with a fuzzy aura local measure, as shown in [39]. FGLCM and FGLAM features are formed by the concatenation of the marginal features of the three *R*, *G*, and *B* components; and a single FCCM captures the interactions among neighboring pixel colors. For a fair comparison, all these matrices were computed with the same parameters as the FCAM (Gaussian membership function, $C^R = C^G = C^B = 16$ for FGLCM and FGLAM, and $C = 16$ for FCCM).

The comparison is based on the segmentation accuracy as the main performance measure and on the computational time required by each feature extraction method. Experiments were carried out using the 1NN classifier instead of the ELM, on a computer with an Intel Core i7 3.60 GHz CPU and 8 GB of RAM. The results over the 20 images of the Prague dataset are summarized in Table 2. They clearly show that FCAM is more relevant, three times faster to compute, and requires less memory than marginal features (FGLCM and FGLAM). FCAM is also more efficient and slightly faster to compute than FCCM.

<table><caption>**Table 2.** Average accuracy (%) and computation time (s) for different fuzzy features on the Prague dataset. “nr” means no segmentation refinement, and “wr” means with segmentation refinement.</caption><thead><tr><td>Feature</td><td>Size</td><td>Accuracy nr</td><td>Accuracy wr</td><td>Comp. Time</td></tr></thead><tbody><tr><td>FGLCM</td><td>3 × 16<sup>2</sup> = 768</td><td>86.69</td><td>90.97</td><td>123.53</td></tr><tr><td>FCCM</td><td>16<sup>2</sup> = 256</td><td>86.52</td><td>90.77</td><td>37.62</td></tr><tr><td>FGLAM</td><td>3 × 16<sup>2</sup> = 768</td><td>87.24</td><td>91.28</td><td>112.32</td></tr><tr><td>FCAM</td><td>16<sup>2</sup> = 256</td><td><b>87.82</b></td><td><b>92.00</b></td><td><b>33.12</b></td></tr></tbody></table>

## 4.3. Comparison with State-of-the-Art Supervised Segmentation Methods

In this section, we compare the results obtained by the proposed segmentation method with the results of several state-of-the-art segmentation methods, on the Prague dataset. The assessment of segmentation performance was based on the conventional measures provided on the Prague texture segmentation website [57], which include: (1) region-based criteria: correct segmentation (CS), over-segmentation (OS), under-segmentation (US), missed error (ME), and noise error (NE); (2) pixel-wise based criteria: omission error (O), commission error (C), class accuracy (CA), recall (CO), precision (CC), type I error (I.), type II error (II.), mean class accuracy estimate (EA), mapping score (MS), root mean square proportion estimation error (RM), and comparison index (CI); (3) consistency-error criteria: global consistency error (GCE) and local consistency error (LCE).

The methods involved in the comparison were: (1) the MRF algorithm based on a Markov random field pixel classification model [11], (2) the COF algorithm that uses the co-occurrence features and the 1NN classifier [57], (3) the Con-Col algorithm [57], (4) the supervised fully convolutional network for texture (FCNT) algorithm [23] without refinement or (5) with refinement, (6) the empirical-wavelet-transform-based fully convolutional network for texture (EWT-FCNT) [23] that combines the empirical wavelet transform with FCNT, (7) U-Net [24], (8) the deep visual model (DA) [25], (9) the pyramid scene parsing network (PSP-Net) [26], and (10) our proposed method without refinement (FCAM nr) or (11) with refinement (FCAM wr).

The segmentation results of EWT-FCNT, FCNT, MRF, COF, and Con-Col were taken from the Prague benchmark website [57]; those of U-Net, DA, and PSP-Net were taken from [23]. All these methods except MRF, COF, and Con-Col, involve a refinement step to improve performance. Regarding the number of prototypes, no information is available for MRF, COF, and Con-Col. In contrast, the CNN-based methods used a training set of 1000 texture mosaic images of size 512 × 512 pixels, specifically created from the original Prague dataset, for each image to segment.
From Table 3, we can see that our proposed method, with and even without refinement, outperforms the classical supervised methods, namely, MRF, COF, and Con-Col. In contrast, EWT-FCNT and other CNN-based methods provided better segmentation results than our method. However, by carefully analyzing these results, we found that—except the EWT-FCNT method, which provided exceptional results—our method provided similar results (to U-Net, FCNT, DA, and PSP-Net). For example, the gaps between its accuracy (CO = 95.20% with refinement) and those of U-Net, FCNT, DA, and PSP-Net were lower than 1.7%. We can also see that our FCAM-based method did not suffer from under-segmentation nor over-segmentation (US and OS measures are equal to 0) compared to deep learning methods. It is also noticeable that our method without refinement outperformed FCNT without refinement according to most of the criteria.

<table><caption>**Table 3.** Results of supervised methods on the Prague dataset (20 test images). Arrows ↑, ↓ denote the required criterion direction; “nr” means no segmentation refinement, and “wr” means with segmentation refinement. Best results are in bold, and the second best ones in italic.</caption><thead><tr><th>Criteria</th><th>MRF</th><th>COF</th><th>Con-Col</th><th>FCNTnr</th><th>FCNTwr</th><th>EWT-FCNT</th><th>U-Net</th><th>DA</th><th>PSP-Net</th><th>FCAMnr</th><th>FCAMwr</th></tr></thead><tbody><tr><td>↑ CS</td><td>46.11</td><td>52.48</td><td>84.57</td><td>87.52</td><td>96.01</td><td><strong>98.45</strong></td><td>96.71</td><td>94.18</td><td>96.45</td><td>84.24</td><td>91.27</td></tr><tr><td>↓ OS</td><td>0.81</td><td><strong>0.00</strong></td><td><strong>0.00</strong></td><td><strong>0.00</strong></td><td>1.56</td><td><strong>0.00</strong></td><td>1.71</td><td><strong>0.00</strong></td><td><strong>0.17</strong></td><td><strong>0.00</strong></td><td><strong>0.00</strong></td></tr><tr><td>↓ US</td><td>4.18</td><td>1.94</td><td>1.70</td><td><strong>0.00</strong></td><td>1.20</td><td><strong>0.00</strong></td><td><strong>0.00</strong></td><td>1.18</td><td><strong>0.41</strong></td><td><strong>0.00</strong></td><td><strong>0.00</strong></td></tr><tr><td>↓ ME</td><td>44.82</td><td>41.55</td><td>9.50</td><td>6.70</td><td>0.78</td><td><strong>0.37</strong></td><td><strong>0.68</strong></td><td>3.42</td><td>1.23</td><td>11.45</td><td>5.93</td></tr><tr><td>↓ NE</td><td>45.29</td><td>40.97</td><td>10.22</td><td>6.90</td><td>0.89</td><td><strong>0.46</strong></td><td><strong>0.48</strong></td><td>3.24</td><td>1.12</td><td>11.39</td><td>5.36</td></tr><tr><td>↓ O</td><td>14.52</td><td>20.74</td><td>7.00</td><td>7.46</td><td>2.72</td><td><strong>0.93</strong></td><td><strong>0.72</strong></td><td>3.13</td><td>2.75</td><td>4.91</td><td>2.96</td></tr><tr><td>↓ C</td><td>16.77</td><td>22.10</td><td>5.34</td><td>6.16</td><td>2.29</td><td><strong>1.04</strong></td><td><strong>0.70</strong></td><td>1.32</td><td>2.39</td><td>5.89</td><td>2.72</td></tr><tr><td>↑ CA</td><td>65.42</td><td>67.01</td><td>86.21</td><td>87.08</td><td>93.95</td><td><strong>97.67</strong></td><td>95.86</td><td>94.53</td><td>93.89</td><td>87.28</td><td>91.54</td></tr><tr><td>↑ CO</td><td>76.19</td><td>77.86</td><td>92.02</td><td>92.61</td><td>96.73</td><td><strong>98.78</strong></td><td>96.91</td><td>96.23</td><td>96.06</td><td>92.60</td><td>95.20</td></tr><tr><td>↑ CC</td><td>80.30</td><td>78.34</td><td>92.68</td><td>93.26</td><td>97.02</td><td><strong>98.81</strong></td><td>97.38</td><td>97.01</td><td>96.41</td><td>93.65</td><td>95.96</td></tr><tr><td>↓ I.</td><td>23.81</td><td>22.14</td><td>7.98</td><td>7.39</td><td>3.27</td><td><strong>1.22</strong></td><td>3.09</td><td>3.77</td><td>3.94</td><td>7.40</td><td>4.80</td></tr><tr><td>↓ II.</td><td>4.82</td><td>4.40</td><td>1.70</td><td>1.49</td><td>0.68</td><td><strong>0.25</strong></td><td><strong>0.41</strong></td><td>0.58</td><td>0.69</td><td>1.32</td><td>0.87</td></tr><tr><td>↑ EA</td><td>75.40</td><td>76.21</td><td>91.72</td><td>92.68</td><td>96.68</td><td><strong>98.77</strong></td><td>97.01</td><td>96.24</td><td>96.08</td><td>92.58</td><td>95.13</td></tr><tr><td>↑ MS</td><td>64.29</td><td>66.79</td><td>88.03</td><td>88.92</td><td>95.10</td><td><strong>98.17</strong></td><td>95.37</td><td>94.35</td><td>94.08</td><td>88.90</td><td>92.80</td></tr><tr><td>↓ RM</td><td>6.43</td><td>4.47</td><td>2.08</td><td>1.38</td><td>0.86</td><td><strong>0.24</strong></td><td><strong>0.61</strong></td><td>1.07</td><td>0.70</td><td>1.64</td><td>1.24</td></tr><tr><td>↑ CI</td><td>76.69</td><td>77.05</td><td>92.02</td><td>92.81</td><td>96.77</td><td><strong>98.78</strong></td><td>97.08</td><td>96.41</td><td>96.15</td><td>92.84</td><td>95.35</td></tr><tr><td>↓ GCE</td><td>25.79</td><td>23.94</td><td>11.76</td><td>12.54</td><td>5.55</td><td>2.33</td><td><strong>2.13</strong></td><td>3.50</td><td>4.67</td><td>11.60</td><td>7.45</td></tr><tr><td>↓ LCE</td><td>20.68</td><td>19.69</td><td>8.61</td><td>9.94</td><td>3.75</td><td><strong>1.68</strong></td><td><strong>1.46</strong></td><td>2.47</td><td>3.52</td><td>8.76</td><td>5.31</td></tr></tbody></table>

To thoroughly analyze segmentation results, Table 4 presents the accuracy obtained on each individual test image by our proposed FCAM method; the handcrafted supervised methods MRF, COF, and Con-Col; and the two CNN-based methods, FCNT (with refinement) and EWT-FCNT. Only the segmentation results of individual test images obtained by the compared methods (EWT-FCNT, FCNT, MRF, COF, and Con-Col) are available on the Prague [57] benchmark website. Table 4 shows that accuracy obtained by our FCAM method was higher than 96% for half of the 20 test images and outperformed MRF, COF and Con-Col for almost all images. For 6 of the 20 tested images, our method also outperformed FCNT and ranked second behind EWT-FCNT. Even better, it provided the best accuracy for image 11. For the other images, the accuracy obtained by our method is close to that of FCNT, though for images 7, 8, 10, and 13, the obtained accuracy is poor (below 91%). This explains why the average accuracy dropped to 95.20%.

Figure 7 shows a visual comparison of two handcrafted methods (COF and Con-Col) and two deep learning-based ones (FCNT wr and EWT-FCNT) with the proposed method. The results of these methods are publicly available online [57]. Overall, FCAM provided satisfactory visual segmentation results that are close to the ground truth. Unlike other methods, FCAM is not prone to over-segmentation, and classification errors almost only occurred at region boundaries. This is mainly due to our pre-segmentation based on regional SLIC that, despite its superiority over the basic SLIC, lacks the accuracy to correctly determine the boundaries between two or more different texture regions.
<table><caption>**Table 4.** Accuracy (CO) for each test image of the Prague dataset. “nr” means no segmentation refinement, and “wr” means with segmentation refinement. Best results are in bold, and the second best ones in italic.</caption><thead><tr><th>Image</th><th>MRF</th><th>COF</th><th>Con-Col</th><th>FCNTwr</th><th>EWT-FCNT</th><th>FCAMnr</th><th>FCAMwr</th></tr></thead><tbody><tr><td>01</td><td>99.79</td><td>96.19</td><td>96.92</td><td>99.12</td><td><strong>99.91</strong></td><td>99.54</td><td>99.72</td></tr><tr><td>02</td><td>77.36</td><td>93.02</td><td>92.70</td><td>97.71</td><td><strong>99.51</strong></td><td>92.66</td><td>96.79</td></tr><tr><td>03</td><td>95.60</td><td>96.56</td><td>92.33</td><td>97.47</td><td><strong>99.21</strong></td><td>98.66</td><td>99.03</td></tr><tr><td>04</td><td>73.20</td><td>68.63</td><td>93.73</td><td>98.41</td><td><strong>98.77</strong></td><td>96.38</td><td>98.58</td></tr><tr><td>05</td><td>89.72</td><td>89.67</td><td>91.64</td><td>96.64</td><td><strong>98.95</strong></td><td>92.95</td><td>93.83</td></tr><tr><td>06</td><td>84.00</td><td>57.59</td><td>95.78</td><td>97.30</td><td><strong>99.52</strong></td><td>96.97</td><td>97.64</td></tr><tr><td>07</td><td>67.74</td><td>58.08</td><td>89.71</td><td>96.09</td><td><strong>96.71</strong></td><td>79.65</td><td>84.70</td></tr><tr><td>08</td><td>58.12</td><td>71.27</td><td>90.10</td><td>95.67</td><td><strong>98.80</strong></td><td>84.16</td><td>90.21</td></tr><tr><td>09</td><td>72.95</td><td>61.29</td><td>93.23</td><td>96.70</td><td><strong>99.35</strong></td><td>94.90</td><td>95.47</td></tr><tr><td>10</td><td>71.52</td><td>58.97</td><td>85.72</td><td>92.52</td><td><strong>96.69</strong></td><td>86.36</td><td>89.51</td></tr><tr><td>11</td><td>61.28</td><td>80.88</td><td>96.82</td><td>96.72</td><td>95.32</td><td>94.44</td><td><strong>97.12</strong></td></tr><tr><td>12</td><td>81.55</td><td>63.01</td><td>87.20</td><td>96.03</td><td><strong>99.51</strong></td><td>92.74</td><td>95.84</td></tr><tr><td>13</td><td>74.35</td><td>84.32</td><td>77.02</td><td>96.10</td><td><strong>98.48</strong></td><td>87.99</td><td>90.11</td></tr><tr><td>14</td><td>91.29</td><td>90.32</td><td>96.51</td><td>97.75</td><td><strong>99.17</strong></td><td>94.25</td><td>96.39</td></tr><tr><td>15</td><td>57.77</td><td>79.61</td><td>96.26</td><td>97.70</td><td><strong>99.56</strong></td><td>95.31</td><td>98.61</td></tr><tr><td>16</td><td>61.33</td><td>60.63</td><td>91.19</td><td>94.31</td><td><strong>99.46</strong></td><td>86.23</td><td>93.92</td></tr><tr><td>17</td><td>62.74</td><td>72.31</td><td>91.92</td><td>96.96</td><td><strong>99.38</strong></td><td>91.21</td><td>95.10</td></tr><tr><td>18</td><td>77.81</td><td>92.96</td><td>96.16</td><td>98.24</td><td><strong>99.58</strong></td><td>97.18</td><td>98.54</td></tr><tr><td>19</td><td>76.07</td><td>94.98</td><td>97.02</td><td>98.55</td><td><strong>99.78</strong></td><td>97.64</td><td>99.39</td></tr><tr><td>20</td><td>89.68</td><td>86.87</td><td>88.48</td><td>94.66</td><td><strong>97.91</strong></td><td>92.82</td><td>93.44</td></tr><tr><td>Average</td><td>76.19</td><td>77.86</td><td>92.02</td><td>96.73</td><td><strong>98.78</strong></td><td>92.60</td><td>95.20</td></tr></tbody></table>

However, it is important to underline that deep learning-based methods (including EWT-FCNT) need to create a large training set (1000 images for each image to segment) from the original training database. Moreover, they require an expensive learning step to extract features and classify pixels. In order to get an overview on the computational costs of our approach compared to CNN-based methods, the following section gives the computational times measured during both training and segmentation phases.

### 4.4. *Processing Time*

In Section 4.2, we estimated the FCAM computational time only at superpixels. In order to provide an overview of the computational requirements of the entire FCAM-based color texture segmentation method, we estimated its overall runtime over the Prague dataset on an 3.60 GHz Intel Core i7 computer with 8 GB RAM. Table 5 displays the average computing time of our proposed method and those of CNN-based methods (EWT-FCNT, FCNT, U-Net, DA, and PSP-Net). These processing times are separated into two parts: the training time and the segmentation time. The segmentation times for CNNs were taken from [23]. They were measured on a laptop computer with 2.5 GHz quad-core Intel Core i7 processor and 16 GB memory, equipped with GTX 1080Ti external GPU with 11 GB memory. The training times of these CNNs were unfortunately not provided. From this table, we can see that the processing time required by our method during the segmentation stage is the highest. It is about 41.14 s for a 512 × 512 image. FCAM computation for superpixels consumes most of this time (33.12 s). The remaining time is mainly shared among the regional SLIC computation (8 s) and the classification of superpixels by ELM (0.02 s). It should be noted that the reported computing time for our method was obtained with MATLAB code without any GPU acceleration, and consequently, our method could be implemented for real-time applications with this enhancement.
<figure><img src="image_12.png" /><figcaption>**Figure 7.** Sample segmentation results on Prague dataset, from left to right: image 01, image 03, image 06, image 15, and image 19.</figcaption></figure>

In contrast, the total duration of the training stage for our method remained lower than 5 min and can be detailed as 260 s to compute the FCAM at the prototypes, and 0.9 s to train the ELM. The overall training time (260.9 s) was much lower than with CNNs, which typically require several hours on powerful computers equipped with large-memory GPUs.
<table><caption>**Table 5.** Average computation time per image of FCNT, EWT-FCNT, U-Net, DA, PSP-Net, and FCAM on the Prague dataset.</caption><thead><tr><th>Stage</th><th>FCNT</th><th>EWT-FCNT</th><th>U-Net</th><th>DA</th><th>PSP-Net</th><th>FCAM</th></tr></thead><tbody><tr><td>Segmentation</td><td>3.61 ms</td><td>1.830 s</td><td>4.98 ms</td><td>3.80 ms</td><td>14.39 ms</td><td>41.14 s</td></tr><tr><td>Training</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>260.9 s</td></tr></tbody></table>

## 5. Conclusions

In this paper, we introduced fuzzy color aura cardinal matrices (FCAMs) to locally characterize colors and textures, and applied them for color texture image segmentation. The FCAM feature makes it possible to locally characterize the interactions between colors of neighboring sites. A single low-dimensional FCAM is required to describe the color texture at each site, unlike in the marginal approach where a fuzzy gray-level aura matrix (FGLAM) must be computed for each color channel.

The proposed color texture image segmentation is based on the classification of super-pixels, generated from a modified version of the SLIC algorithm to incorporate regional information. An FCAM is then computed for each superpixel thanks to a locally-adaptive neighborhood function. The superpixels are finally classified using a simple supervised ELM classifier. Experiments on the Prague texture segmentation benchmark showed that the proposed color texture segmentation based on FCAMs outperforms the classical state-of-the-art segmentation methods and is competitive with recent methods based on deep learning. However, unlike CNN-based approaches that require an expensive learning procedure and a large training set of segmented texture images, whose construction is time-consuming, our method is applied straightforwardly from a much smaller database using ELM-based classification.

Despite its respectable performance, our segmentation method remains sensitive to pre-segmentation results. Although the regional SLIC improves segmentation results in comparison with the basic SLIC, the detection of color texture boundaries is still not accurate enough. In future work, we intend to use a texture-aware superpixel procedure that takes into account the properties of the available color spaces. The proposed segmentation method was developed in a supervised context; we also plan to adapt it to unsupervised color texture segmentation.

**Author Contributions:** Conceptualization, Z.H., K.H., O.L. and L.M.; methodology, Z.H., K.H., O.L. and L.M.; software, Z.H.; writing—review and editing, Z.H., K.H., O.L. and L.M. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Borovec, J.; Švihlík, J.; Kybic, J.; Habart, D. Supervised and unsupervised segmentation using superpixels, model estimation, and graph cut. *J. Electron. Imaging* **2017**, *26*, 061610. [CrossRef]

2. Mikeš, S.; Haindl, M.; Scarpa, G.; Gaetano, R. Benchmarking of remote sensing segmentation methods. *IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens.* **2015**, *8*, 2240–2248. [CrossRef]

3. Akbarizadeh, G.; Rahmani, M. Efficient combination of texture and color features in a new spectral clustering method for PolSAR image segmentation. *Natl. Acad. Sci. Lett.* **2017**, *40*, 117–120. [CrossRef]

4. Zhang, C.; Zou, K.; Pan, Y. A method of apple image segmentation based on color-texture fusion feature and machine learning. *Agronomy* **2020**, *10*, 972. [CrossRef]

5. Ilea, D.E.; Whelan, P.F. Image segmentation based on the integration of colour–texture descriptors–A review. *Pattern Recognit.* **2011**, *44*, 2479–2501. [CrossRef]
6. Haindl, M.; Mikeš, S. A competition in unsupervised color image segmentation. *Pattern Recognit.* **2016**, *57*, 136–151. [CrossRef]

7. Garcia-Lamont, F.; Cervantes, J.; López, A.; Rodriguez, L. Segmentation of images by color features: A survey. *Neurocomputing* **2018**, *292*, 1–27. [CrossRef]

8. Bianconi, F.; Harvey, R.W.; Southam, P.; Fernández, A. Theoretical and experimental comparison of different approaches for color texture classification. *J. Electron. Imaging* **2011**, *20*, 043006. [CrossRef]

9. Ledoux, A.; Losson, O.; Macaire, L. Color local binary patterns: Compact descriptors for texture classification. *J. Electron. Imaging* **2016**, *25*, 061404. [CrossRef]

10. Liu, L.; Chen, J.; Fieguth, P.; Zhao, G.; Chellappa, R.; Pietikäinen, M. From BoW to CNN: Two decades of texture representation for texture classification. *Int. J. Comput. Vis.* **2019**, *127*, 74–109. [CrossRef]

11. Kato, Z.; Pong, T.C.; Lee, J.C.M. Color image segmentation and parameter estimation in a Markovian framework. *Pattern Recognit. Lett.* **2001**, *22*, 309–321. [CrossRef]

12. Nakyoung, O.; Choi, J.; Kim, D.; Kim, C. Supervised classification and segmentation of textured scene images. In Proceedings of the International Conference on Consumer Electronics (ICCE 2015), Las Vegas, NV, USA, 9–12 January 2015; IEEE: Piscataway, NJ, USA, 2015; pp. 473–476. [CrossRef]

13. Yang, S.; Lv, Y.; Ren, Y.; Yang, L.; Jiao, L. Unsupervised images segmentation via incremental dictionary learning based sparse representation. *Inf. Sci.* **2014**, *269*, 48–59. [CrossRef]

14. Laboreiro, V.R.S.; de Araujo, T.P.; Bessa Maia, J.E. A texture analysis approach to supervised face segmentation. In Proceedings of the IEEE Symposium on Computers and Communications (ISCC 2014), Funchal, Portugal, 23–26 June 2014; IEEE: Piscataway, NJ, USA, 2014; pp. 1–6. [CrossRef]

15. Panjwani, D.; Healey, G. Markov random field models for unsupervised segmentation of textured color images. *IEEE Trans. Pattern Anal. Mach. Intell.* **1995**, *17*, 939–954. [CrossRef]

16. Chen, K.M.; Chen, S.Y. Color texture segmentation using feature distributions. *Pattern Recognit. Lett.* **2002**, *23*, 755–771. [CrossRef]

17. Salmi, A.; Hammouche, K.; Macaire, L. Constrained feature selection for semisupervised color-texture image segmentation using spectral clustering. *J. Electron. Imaging* **2021**, *30*, 013014. [CrossRef]

18. Jenicka, S. Supervised Texture-Based Segmentation Using Basic Texture Models. In *Land Cover Classification of Remotely Sensed Images*; Springer: Berlin/Heidelberg, Germany, 2021; pp. 73–88.

19. Al-Kadi, O.S. Supervised texture segmentation: A comparative study. In Proceedings of the IEEE Jordan Conference on Applied Electrical Engineering and Computing Technologies 2011 (AEECT 2011), Amman, Jordan, 6–8 December 2011; IEEE: Piscataway, NJ, USA, 2011; pp. 1–5. [CrossRef]

20. Yang, H.Y.; Wang, X.Y.; Wang, Q.Y.; Zhang, X.J. LS-SVM based image segmentation using color and texture information. *J. Vis. Commun. Image Represent.* **2012**, *23*, 1095–1112. [CrossRef]

21. Andrearczyk, V.; Whelan, P.F. Texture segmentation with fully convolutional networks. *arXiv* **2017**, arXiv:1703.05230.

22. Garcia-Garcia, A.; Orts-Escolano, S.; Oprea, S.; Villena-Martinez, V.; Martinez-Gonzalez, P.; Garcia-Rodriguez, J. A survey on deep learning techniques for image and video semantic segmentation. *Appl. Soft Comput.* **2018**, *70*, 41–65. [CrossRef]

23. Huang, Y.; Zhou, F.; Gilles, J. Empirical curvelet based fully convolutional network for supervised texture image segmentation. *Neurocomputing* **2019**, *349*, 31–43. [CrossRef]

24. Ronneberger, O.; Fischer, P.; Brox, T. U-Net: Convolutional networks for biomedical image segmentation. In Proceedings of the International Conference on Medical Image Computing and Computer-Assisted Intervention (MICCAI 2015), Munich, Germany, 5–9 October 2015; Springer: Cham, Switzerland, 2015; pp. 234–241. [CrossRef]

25. Wang, W.; Shen, J. Deep visual attention prediction. *IEEE Trans. Image Process.* **2017**, *27*, 2368–2378. [CrossRef]

26. Zhao, H.; Shi, J.; Qi, X.; Wang, X.; Jia, J. Pyramid Scene Parsing Network. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR 2017), Honolulu, HI, USA, 21–26 July 2017; IEEE: Piscataway, NJ, USA, 2017; pp. 6230–6239. [CrossRef]

27. Yuan, J.; Wang, D.; Cheriyadat, A.M. Factorization-based texture segmentation. *IEEE Trans. Image Process.* **2015**, *24*, 3488–3497. [CrossRef]

28. Yang, Y.; Han, S.; Wang, T.; Tao, W.; Tai, X.C. Multilayer graph cuts based unsupervised color-texture image segmentation using multivariate mixed Student's t-distribution and regional credibility merging. *Pattern Recognit.* **2013**, *46*, 1101–1124. [CrossRef]

29. Akbulut, Y.; Guo, Y.; Şengür, A.; Aslan, M. An effective color texture image segmentation algorithm based on hermite transform. *Appl. Soft Comput.* **2018**, *67*, 494–504. [CrossRef]

30. Barcelo, A.; Montseny, E.; Sobrevilla, P. Fuzzy texture unit and fuzzy texture spectrum for texture characterization. *Fuzzy Sets Syst.* **2007**, *158*, 239–252. [CrossRef]

31. Keramidas, E.; Iakovidis, D.; Maroulis, D. Fuzzy binary patterns for uncertainty-aware texture representation. *Electron. Lett. Comput. Vis. Image Anal.* **2011**, *10*, 63–78. [CrossRef]

32. Vieira, R.T.; de Oliveira Chierici, C.E.; Ferraz, C.T.; Gonzaga, A. Local Fuzzy Pattern: A New Way for Micro-pattern Analysis. In Proceedings of the 13th International Conference on Intelligent Data Engineering and Automated Learning (IDEAL 2012), Natal, Brazil, 29–31 August 2012; Lecture Notes in Computer Science; Springer: Cham, Switzerland, 2012; Volume 7435, pp. 602–611. [CrossRef]

33. Jawahar, C.; Ray, A. Incorporation of gray-level imprecision in representation and processing of digital images. *Pattern Recognit. Lett.* **1996**, *17*, 541–546. [CrossRef]
34. Cheng, H.D.; Chen, C.H.; Chiu, H.H. Image segmentation using fuzzy homogeneity criterion. *Inf. Sci.* **1997**, *98*, 237–262. [CrossRef]

35. Sen, D.; Pal, S.K. Image segmentation using global and local fuzzy statistics. In Proceedings of the IEEE India Council International Conference (INDICON 2006), New Delhi, India, 15–17 September 2006; IEEE: Piscataway, NJ, USA, 2006; pp. 1–6. [CrossRef]

36. Munklang, Y.; Auephanwiriyakul, S.; Theera-Umpon, N. A novel fuzzy co-occurrence matrix for texture feature extraction. In Proceedings of the 13th International Conference on Computational Science and its Applications (ICCSA 2013), Ho Chi Minh City, Vietnam, 24–27 June 2013; Lecture Notes in Computer Science; Springer: Berlin/Heidelberg, Germany, 2013; Volume 7973, pp. 246–257. [CrossRef]

37. Khaldi, B.; Kherfi, M.L. Modified integrative color intensity co-occurrence matrix for texture image representation. *J. Electron. Imaging* **2016**, *25*, 053007. [CrossRef]

38. Ledoux, A.; Losson, O.; Macaire, L. Texture classification with fuzzy color co-occurrence matrices. In Proceedings of the IEEE International Conference on Image Processing (ICIP 2015), Québec City, QC, Canada, 27–30 September 2015; IEEE: Piscataway, NJ, USA, 2015; pp. 1429–1433. [CrossRef]

39. Hammouche, K.; Losson, O.; Macaire, L. Fuzzy aura matrices for texture classification. *Pattern Recognit.* **2016**, *53*, 212–228. [CrossRef]

40. Elfadel, I.M.; Picard, R.W. Gibbs random fields, cooccurrences, and texture modeling. *IEEE Trans. Pattern Anal. Mach. Intell.* **1994**, *16*, 24–37. [CrossRef]

41. Qin, X.; Yang, Y.H. Basic gray level aura matrices: Theory and its application to texture synthesis. In Proceedings of the 10th International Conference on Computer Vision (ICCV’05), Beijing, China, 17–20 October 2005; IEEE: Piscataway, NJ, USA, 2005; Volume 1, pp. 128–135. [CrossRef]

42. Qin, X.; Yang, Y.H. Aura 3D textures. *IEEE Trans. Vis. Comput. Graph.* **2007**, *13*, 379–389. [CrossRef] [PubMed]

43. Qin, X.; Yang, Y.H. Similarity measure and learning with gray level aura matrices (GLAM) for texture image retrieval. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Washington, DC, USA, 27 June–2 July 2004; IEEE: Piscataway, NJ, USA, 2004; Volume 1, pp. I-326–I-333. [CrossRef]

44. Wiesmüller, S.; Chandy, D.A. Content based mammogram retrieval using gray level aura matrix. *Int. J. Comput. Commun. Inf. Syst.* **2010**, *2*, 217–223.

45. Liao, S.; Chung, A.C.S. Texture classification by using advanced local binary patterns and spatial distribution of dominant patterns. In Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP’07), Honolulu, HI, USA, 15–20 April 2007; IEEE: Piscataway, NJ, USA, 2007; pp. 1221–1224. [CrossRef]

46. Wang, Y.; Gao, X.; Fu, R.; Jian, Y. Dayside corona aurora classification based on X-gray level aura matrices. In Proceedings of the 9th ACM International Conference on Image and Video Retrieval (CIVR 2010), Xi’an, China, 5–7 July 2010; ACM: New York, NY, USA, 2010; pp. 282–287. [CrossRef]

47. Hannan, M.A.; Arebey, M.; Begum, R.A.; Basri, H. Gray level aura matrix: An image processing approach for waste bin level detection. In Proceedings of the World Congress on Sustainable Technologies (WCST 2011), London, UK, 7–10 November 2011; IEEE: Piscataway, NJ, USA, 2011; pp. 77–82. [CrossRef]

48. Khalid, M.; Yusof, R.; Khairuddin, A.S.M. Improved tropical wood species recognition system based on multi-feature extractor and classifier. *Int. J. Electr. Comput. Eng.* **2011**, *5*, 1495–1501.

49. Hannan, M.A.; Arebey, M.; Begum, R.A.; Basri, H. An automated solid waste bin level detection system using a gray level aura matrix. *Waste Manag.* **2012**, *32*, 2229–2238. [CrossRef] [PubMed]

50. Yusof, R.; Khalid, M.; Khairuddin, A.S.M. Application of kernel-genetic algorithm as nonlinear feature selection in tropical wood species recognition system. *Comput. Electron. Agric.* **2013**, *93*, 68–77. [CrossRef]

51. Haliche, Z.; Hammouche, K. The gray level aura matrices for textured image segmentation. *Analog. Integr. Circuits Signal Process.* **2011**, *69*, 29–38. [CrossRef]

52. Haliche, Z.; Hammouche, K.; Postaire, J.G. Texture image segmentation based on the elements of gray level aura matrices. In Proceedings of the Global Summit on Computer & Information Technology (GSCIT 2014), Sousse, Tunisia, 14–16 June 2014; IEEE: Piscataway, NJ, USA, 2014; pp. 1–6. [CrossRef]

53. Haliche, Z.; Hammouche, K.; Postaire, J.G. A fast algorithm for texture feature extraction from gray level aura matrices. *Int. J. Circuits Syst. Signal Process.* **2015**, *9*, 54–61.

54. Han, J.C.; Zhao, P.; Wang, C.K. Wood species recognition through FGLAM textural and spectral feature fusion. *Wood Sci. Technol.* **2021**, *55*, 535–552. [CrossRef]

55. Achanta, R.; Shaji, A.; Smith, K.; Lucchi, A.; Fua, P.; Süsstrunk, S. SLIC superpixels compared to state-of-the-art superpixel methods. *IEEE Trans. Pattern Anal. Mach. Intell.* **2012**, *34*, 2274–2282. [CrossRef]

56. Chen, X.; Zheng, C.; Yao, H.; Wang, B. Image segmentation using a unified Markov random field model. *IET Image Process.* **2017**, *11*, 860–869. [CrossRef]

57. Mikes, S.; Haindl, M. Texture Segmentation Benchmark. *IEEE Trans. Pattern Anal. Mach. Intell.* **2021**, *44*, 5647–5663. [CrossRef] [PubMed]

58. Huang, G.B.; Zhu, Q.Y.; Siew, C.K. Extreme learning machine: Theory and applications. *Neurocomputing* **2006**, *70*, 489–501. [CrossRef]
59. Liu, M.Y.; Tuzel, O.; Ramalingam, S.; Chellappa, R. Entropy rate superpixel segmentation. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR'11), Colorado Springs, CO, USA, 20-25 June 2011; IEEE: Piscataway, NJ, USA, 2011; pp. 2097-2104. [CrossRef]

60. Stutz, D.; Hermans, A.; Leibe, B. Superpixels: An evaluation of the state-of-the-art. *Comput. Vis. Image Underst.* **2018**, *166*, 1-27.
[CrossRef]