Next Article in Journal
Nucleation Process in Explosive Boiling Phenomena of Water on Micro-Platinum Wire
Previous Article in Journal
The Ontology of Haag’s Local Quantum Physics
Previous Article in Special Issue
Multi-Channel Representation Learning Enhanced Unfolding Multi-Scale Compressed Sensing Network for High Quality Image Reconstruction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Classification of Multiple H&E Images via an Ensemble Computational Scheme

by
Leonardo H. da Costa Longo
1,*,
Guilherme F. Roberto
2,
Thaína A. A. Tosta
3,
Paulo R. de Faria
4,
Adriano M. Loyola
5,
Sérgio V. Cardoso
5,
Adriano B. Silva
6,
Marcelo Z. do Nascimento
6 and
Leandro A. Neves
1,*
1
Department of Computer Science and Statistics (DCCE), São Paulo State University (UNESP), Rua Cristóvão Colombo, 2265, São José do Rio Preto 15054-000, SP, Brazil
2
Department of Informatics Engineering, Faculty of Engineering, University of Porto, Dr. Roberto Frias, sn, 4200-465 Porto, Portugal
3
Science and Technology Institute, Federal University of São Paulo (UNIFESP), Avenida Cesare Mansueto Giulio Lattes, 1201, São José dos Campos 12247-014, SP, Brazil
4
Department of Histology and Morphology, Institute of Biomedical Science, Federal University of Uberlândia (UFU), Av. Amazonas, S/N, Uberlândia 38405-320, MG, Brazil
5
Area of Oral Pathology, School of Dentistry, Federal University of Uberlândia (UFU), R. Ceará—Umuarama, Uberlândia 38402-018, MG, Brazil
6
Faculty of Computer Science (FACOM), Federal University of Uberlândia (UFU), Avenida João Naves de Ávila 2121, Bl.B, Uberlândia 38400-902, MG, Brazil
*
Authors to whom correspondence should be addressed.
Entropy 2024, 26(1), 34; https://doi.org/10.3390/e26010034
Submission received: 27 November 2023 / Revised: 23 December 2023 / Accepted: 25 December 2023 / Published: 28 December 2023
(This article belongs to the Special Issue Information Theory in Image Processing and Pattern Recognition)

Abstract

:
In this work, a computational scheme is proposed to identify the main combinations of handcrafted descriptors and deep-learned features capable of classifying histological images stained with hematoxylin and eosin. The handcrafted descriptors were those representatives of multiscale and multidimensional fractal techniques (fractal dimension, lacunarity and percolation) applied to quantify the histological images with the corresponding representations via explainable artificial intelligence (xAI) approaches. The deep-learned features were obtained from different convolutional neural networks (DenseNet-121, EfficientNet-b2, Inception-V3, ResNet-50 and VGG-19). The descriptors were investigated through different associations. The most relevant combinations, defined through a ranking algorithm, were analyzed via a heterogeneous ensemble of classifiers with the support vector machine, naive Bayes, random forest and K-nearest neighbors algorithms. The proposed scheme was applied to histological samples representative of breast cancer, colorectal cancer, oral dysplasia and liver tissue. The best results were accuracy rates of 94.83 % to 100 % , with the identification of pattern ensembles for classifying multiple histological images. The computational scheme indicated solutions exploring a reduced number of features (a maximum of 25 descriptors) and with better performance values than those observed in the literature. The presented information in this study is useful to complement and improve the development of computer-aided diagnosis focused on histological images.

1. Introduction

Histopathology considers the study of biological tissues and their microscopic structures. Histopathologists examine tissue samples under a microscope to identify abnormal changes in cell structure and indicate possible pathological conditions [1,2]. In this process, staining techniques, such as hematoxylin and eosin (H&E), and microscopy are commonly explored in order to obtain information for the diagnosis, treatment and understanding of the progression of hyperplasias, dysplasias, metaplasias and neoplasms [1,3,4]. For example, when the H&E technique is used, the staining process aims to highlight the basophils and eosinophils in tissues [5,6]. The first dye highlights cell nuclei with a bluish color, while the second makes the cytoplasm reddish. Other cellular structures have colors derived from mixing these dyes [7]. These conditions contribute to the analysis of regions of interest, whether by specialists or automated systems aimed at classifying and recognizing cancer patterns, considered a global public health issue. Over the past decade, there has been a 20% increase in cancer incidence, and more than 25 million new cases are expected by 2030. According to [8], cancer has become an important cause of premature mortality globally and is associated with high social and economic costs. The estimated productivity losses are EUR 104.6 billion of the gross national domestic product in Europe and USD 46.3 billion of the combined gross domestic product of the BRICS countries (Brazil, Russia, India, China and South Africa).
In this context, delays in the diagnosis and treatment of cancer can increase the rates of the disease in advanced stages and consequently, mortality. On the other hand, the characterization of cellular changes and their associations with cancer are complex and challenging processes for histopathologists [3,4,9]. Methods based on computer vision and artificial intelligence techniques have provided important advances to minimize these challenges, increasing cancer diagnostic and prognostic accuracy values, especially through computer-aided diagnosis (CAD) [10,11]. In these systems, the extraction and classification of features are essential for the recognition of histological patterns commonly examined by pathologists, especially in the contexts of colorectal cancer, breast cancer, oral dysplasia and liver tissue [12,13]. It is noted, for example, that the combined use of techniques (ensemble learning), such as those based on convolutional neural networks (CNN), make it possible to characterize information at different levels and scales to verify how they relate to each other in the data space [14].
Although a CNN architecture directly performs image classifications, when the values in its internal layers (deep-learned features) were investigated separately, the results indicated more relevant performance, with new approaches in the context of medical images [15,16,17,18,19,20], including the use of transfer learning [21] or even different ensembles of descriptors [22,23,24,25,26,27]. Among the combinations with deep-learned features, the use of fractal descriptors (handcrafted), such as fractal dimension, lacunarity and percolation, deserves to be highlighted because they are capable of appropriately measuring complex shapes generally found in nature and in the context of H&E samples [12,28,29,30,31,32,33,34]. Moreover, the handcrafted features that we selected are commonly used for describing complex structures like the ones found in histopathological images. Therefore, it is noted that proposals based on ensemble learning have been indicated as one of the main research fields for the development of new models [35].
In addition, it is important to note that explainable artificial intelligence (xAI) has contributed significantly to the improvement of ensemble learning models, especially in the validation and interpretation of results [36,37], in order to ensure that accurate classifications are determined for the right reasons [38,39]. This question motivated the development of strategies based on class activation mappings (CAMs) [40], specifically, gradient-weighted class activation mappings (Grad-CAMs) and local interpretable model-agnostic explanations (LIME) [41]. These approaches can be applied to produce visualizations of the image regions that support the CNN classification process. Quantifying this type of image with fractal techniques can complement the process of classification and pattern recognition of histological images. Despite these observations, it is noted that this hypothesis has not yet been investigated in the specialized literature in the context of multiple H&E datasets such as those explored here.
The strategies highlighted previously result in a highly complex feature space, a fact that can make investigations of contexts with a reduced number of samples unfeasible [42,43], such as histological datasets commonly used to investigate colorectal cancer, breast cancer, oral dysplasia and liver tissues. This situation can be overcome by identifying the most relevant features for the classification process and, consequently, indicating more accurate and robust CAD systems. Thus, feature selection plays a critical role in identifying patterns [44], but there is no universal approach in order to define the best results for all contexts [45,46,47,48]. On the other hand, algorithms based on ranking and filters, such as ReliefF, are capable of detecting feature dependencies and provide the best solutions in different experiments [47,49,50,51]. This strategy provides sets of features with different dimensions via any desired criteria. In addition, a set of descriptors can also be highly dependent on the heuristics of the classifier used to evaluate the model [52]. This challenge can be minimized through a classification process with different heuristics (ensemble of classifiers) [53]. Techniques that are part of this set consider the so-called crowd wisdom, in which a decision is made from different perspectives and when associated, can be more accurate. The justification is that possible individual errors are compensated by the successes of the other components [53]. Therefore, investigating the most relevant combinations of descriptors and techniques for the analysis, classification and pattern recognition of H&E images remains an ongoing challenge. This makes ensemble-learning-based solutions more generalizable and robust.
Even with some initiatives observed [12,54,55,56,57], models based on ensemble learning with multiscale and multidimensional fractal descriptors, such as those investigated here, have not yet been fully explored in the literature, including the quantification with the Grad-CAM and LIME representations. In this context, some insights are still pertinent, such as whether it is possible to define standards between the techniques used to classify multiple types of H&E images; whether multiscale and multidimensional fractal descriptors indicate gains in relation to the results achieved via deep-learned features with transfer learning; whether fractal descriptors obtained from Grad-CAM and LIME representations can contribute to the performance of an ensemble learning scheme; and whether the combination of ensembles (descriptors and classifiers) indicates more competitive performance in relation to that available in the specialized literature. Some architectures such as DenseNet-121 [58], Inception-V3 [59], ResNet-50 [60] and VGG-19 [61] can still be explored to provide deep-learned features, even with the success provided in the direct classification of some types of histological images [23,62,63,64]. EfficientNet [65] has not yet been fully explored with the approaches presented here. Therefore, the strategies and conditions previously presented are useful to make knowledge comprehensible to the specialists focused on developing and improving CAD systems. Moreover, the proposed scheme to identify, select and classify the main combinations is based on aspects widely discussed in information theory, image processing and pattern recognition, especially to obtain more robust baseline schemes in this context of histological samples.
In this work, a computational scheme was defined to identify the most relevant feature ensembles in order to complement and improve the development of CAD systems focused on H&E images. The handcrafted descriptors were defined using multiscale and multidimensional fractal techniques (percolation, fractal dimension and lacunarity) to quantify the original H&E samples and the corresponding LIME and Grad-CAM representations. The deep-learned features were obtained from the DenseNet-121, EfficientNet-b2, Inception-V3, ResNet-50 and VGG-19 architectures. The descriptors were analyzed based on different ensembles, considering the ReliefF algorithm with an ensemble of classifiers (support vector machine, naive Bayes, random forest and K-nearest neighbors). The proposed methodology was applied to distinguish histological samples representative of breast cancer, colorectal cancer, oral dysplasia and liver tissue. The information and conditions obtained were detailed from each experiment. The main contributions of this work are:
  • A computational scheme capable of indicating the main ensembles of descriptors for the study of histological images, exploring the ReliefF algorithm and multiple classifiers;
  • An optimized ensemble of deep-learned features with the best results for classifying colorectal cancer, liver tissue and oral dysplasia, using a reduced number of features (up to 25 descriptors);
  • Indications of the discriminative power of ensembles based on fractal features from the LIME and CAM representations;
  • Solutions without overfitting and a more robust baseline scheme, with the necessary details for comparisons and improvements of CAD systems focused on H&E images.
In the Section 2 of this paper, the methodology is described in detail, and in Section 3, the results are presented and discussed. The conclusion is drawn in Section 4.

2. Materials and Methods

In this section, the main steps for the proposed scheme are described, exploring the combined use of deep-learned features via transfer learning with fractal descriptors obtained from original H&E images and their CAM and LIME representations. In the first step, the CNN architectures were defined, and the output layer was fine-tuned to match the classes available on each dataset. Also, in the same step, the CAM and LIME representations were generated considering the fine-tuned models. The second step defined the extraction of the deep-learned features from the selected architectures and the multiscale and multidimensional fractal features using the original H&E images and their xAI representations. In the third step, the extracted features were combined (feature vectors) in different ensembles in order to identify the most relevant information in each context. In the fourth step, the features from each ensemble were ranked and selected to define the best solutions with a reduced number of descriptors. In the fifth and last step, the discriminative capacities of the optimized ensembles were verified through a heterogeneous ensemble with four classifiers. An overview of the proposed scheme is illustrated in Figure 1, with details presented in the next sections.

2.1. Datasets

The proposed scheme was tested on four H&E-stained histological datasets, representatives of breast cancer (UCSB) from [66], colorectal cancer (CR) from [67], liver tissue (LG) from [68] and oral epithelial dysplasia (OED) from [69]. The main details about these datasets are in Table 1, with some samples presented in Figure 2 in order to illustrate each context.

2.2. Step 1—Fine-Tuning the CNN and xAI Representations

Five CNN architectures were considered in the present scheme: DenseNet-121 [58], EfficientNet-b2 [65], Inception-V3 [59], ResNet-50 [60] and VGG-19 [61]. These models were chosen considering their different image classification strategies, which provided a broader investigation of the proposed scheme with different deep-learned features. All models were obtained from the PyTorch library, with details presented in Table 2, including accuracy values achieved on the ImageNet dataset [70].
The fine-tuning step was applied to map the last layer of each CNN with the available classes, changing the final connections and weights corresponding to the total number of groups in each H&E dataset (Table 1). This strategy avoided the full network training stage and made it possible to investigate datasets with a reduced number of images. This process was performed with a k-fold cross-validation ( k = 10 ), which divided the input in k folds of approximately the same size in order to train the model from several possible samples. For each of the k iterations of the training process, a fold was used as evaluation data and the other k 1 as training data. Each iteration should train a new model, and the output selected was the model with the highest accuracy when classifying test data left out of the process’s training input. This technique was applied to reduce the possibility of overfitting.
Figure 1. An overview of the proposed scheme.
Figure 1. An overview of the proposed scheme.
Entropy 26 00034 g001
Figure 2. Examples of H&E images: breast UCSB from [66], benign (a) and malignant (b); CR from [67], benign (c) and malignant (d); LG from [68], male (e) and female (f); OED from [69], healthy (g) and severe (h).
Figure 2. Examples of H&E images: breast UCSB from [66], benign (a) and malignant (b); CR from [67], benign (c) and malignant (d); LG from [68], male (e) and female (f); OED from [69], healthy (g) and severe (h).
Entropy 26 00034 g002
In addition, each training considered 10 epochs, using the stochastic gradient descent (SGDM) strategy, and an initial learning rate l r = 0.01 with a reduction factor of 0.75 every 2 epochs; the cross-entropy function was used to calculate the adjustment on the parameters. This was repeated for each permutation of architecture and dataset. It is important to highlight that the input images were normalized considering the standard deviation and average of the ImageNet dataset’s color channel values to match the methodology used on the model’s pretraining [72]. Finally, the resulting fine-tuned models were the ones that achieved the highest accuracy for the evaluation set, independent of the epoch.

xAI Representations: LIME and Grad-CAM

To obtain each xAI representation, LIME and Grad-CAM methods were applied to every image using the fine-tuned model for the corresponding dataset. The Grad-CAM representations were defined through the last convolutional layer of each CNN architecture, using the approach of [73]. This choice was based on the idea that the deepest layers contained values related to global patterns on the input [74]. The output was a map, converted into a heatmap, with the weights indicating the contribution of each pixel to the final classification. In Figure 3, the Grad-CAM representations of some H&E images are illustrated.
For the LIME method, the results were defined via 1000 local disturbances and using the quick-shift segmentation algorithm [75]. The obtained representation indicated an image with five regions of interest, or superpixels, that were most relevant to explaining the classification result. In Figure 4, some examples of LIME representations obtained from H&E images are illustrated.

2.3. Step 2—Feature Extraction

In this study, the attributes were defined from three origins to compose the ensembles: fractal features from the H&E images; deep-learned features from the layer preceding the output on multiple CNN models; and fractal features from the LIME and CAM representations. These three groups were identified as handcrafted features, deep-learned features and xAI features, respectively.

2.3.1. Handcrafted Features: Multiscale and Multidimensional Fractal Techniques

The quantification was carried out using multiscale and multidimensional fractal techniques, specifically the fractal dimension, lacunarity and percolation approaches. The fractal dimension was based on the idea of expanding Euclidean concepts, in which measurements were contained in an n-dimensional space, with n being an integer greater than zero. The fractal dimension quantified the amount of space filled, indicating the roughness of the structure under analysis. Lacunarity is a complementary measure to the fractal dimension, quantifying the distribution and organization of pixels contained in an image. The lacunarity values represented how the patterns were organized at different observation scales. Percolation is a physical concept that can be observed in the movement and filtering of fluids through porous materials. A classic example of this phenomenon is water flowing through a glass of coffee powder. Considering the quantification process, this concept was explored to indicate the number of clusters, image porosity and cluster size [32].
Fractal techniques were calculated from probability matrices, responsible for storing the probabilities of an image containing a square region filled with each of the shapes that constituted it. This filling was verified via a distance relationship between pixel values under analysis and the size of the area in question [76]. Each matrix was obtained using the gliding-box method, which consisted of sliding a square box of side r across the entire image and checking whether the pixels were inside or outside the box [77]. Thus, given an r-sided box with a central pixel p c , a pixel p was considered when its distance d in relation to p c was less than or equal to r. This process resulted in a frequency distribution matrix N ( m , r ) , with m representing the number of pixels within a box of side r. The probability matrix P ( m , r ) was obtained through the normalization of N ( m , r ) , according to Equation (1), dividing the value of each count by n r , which represented the total number of boxes of side r contained in the image with the application of the gliding box, as seen in Equation (2).
P ( m , r ) = N ( m , r ) n r ,
n r = ( w i d t h r + 1 ) ( h e i g h t r + 1 ) .
The previously described strategy indicated a multiscale quantification due to the variation in r. In addition, as described by [78], in the proposed scheme, the Chebyshev distance or chessboard was applied to calculate d. Finally, as the application context uses colored images with RGB color space, the multidimensional strategy was applied considering a pixel via a 5-dimensional representation, such as ( x , y , r , g , b ) [78].
From these procedures, the fractal dimension F D was calculated considering the estimate of boxes of side r necessary to overlay an image, as shown in Equation (3). The lacunarity L was calculated from the first and second moments of the probability matrix, obtained according to Equations (4) and (5), respectively. These moments were combined using Equation (6) and resulted in the measurement of the lacunarity on a specific observation scale.
F D ( r ) = m = 1 r 2 P ( m , r ) m .
μ ( r ) = m = 1 r 2 m P ( m , r ) .
μ 2 ( r ) = m = 1 r 2 m 2 P ( m , r ) .
L ( r ) = μ 2 ( r ) ( μ ( r ) ) 2 ( μ ( r ) ) 2 .
Percolation measurements were also extracted using the gliding-box algorithm, following the strategy of [32]. Therefore, given a box with side r, the pixels contained in the box were considered pores and labeled using the Hoshen–Kopelman [79] algorithm. Pores with the same label were understood as part of the same cluster. This process was repeated for each box of side r. Thus, it was possible to compute the metrics C ( r ) , which indicated the average number of clusters present in each box; Q ( r ) , which defined the average coverage ratio of the largest cluster; and P ( r ) , which provided the ratio of percolating boxes.
The C ( r ) metric was obtained via an average count of the number of clusters c (present in each box i), according to Equation (7). The metric Q ( r ) was obtained considering an average of the size of the largest cluster | c m a x | in each of the i boxes, according to Equation (8). Finally, the ratio of percolating boxes was defined by dividing the number of boxes in which percolation occurred by the number of boxes computed in an image. Percolation in a box p i occurred when the ratio between the number of pores Ω i and the total number of pixels in the box r 2 exceeded the percolation threshold 0.59275 [12,32], as shown in Equation (9). Thus, P ( r ) was defined via Equation (10).
C ( r ) = i = 1 n r c i n r .
Q ( r ) = i = 1 n r | c m a x i | n r .
p i = 1 , Ω i r 2 0.59275 , 0 , Ω i r 2 < 0.59275 .
P ( r ) = i = 1 n r p i n r .
It is important to highlight that the quantification carried out in this work used matrices defined from boxes with side r, within the 3 r 41 range, according to descriptions presented by [12,78]. Also, the r parameter was set to an odd value to guarantee the existence of a central pixel in each box. Therefore, the increase in the value of r was of two units in each iteration. The probability matrix based on these parameters guaranteed a quantification on 20 different scales.
In addition, both lacunarity functions and percolation measures were also interpreted as scalar values in order to obtain representative descriptors of possible patterns existing in each observation [12,32,78,80,81,82,83]. In these proposals, the authors were able to point out how some of these curves displayed a distinct behavior for each of the classes, making them relevant to the classification process. These features were defined as:
  • Area under the curve (A): It indicates the complexity of the texture. For a discrete function consisting of N points defined in x 1 , , x n , this descriptor can be obtained via Equation (11), with a and b as the point indices that delimit the analysis range;
    A ( a , b ) = b a 2 N n = a b 1 ( f ( x n ) + f ( x n + 1 ) ) .
  • Skewness (S): it is defined via Equation (12), where N is the number of points in the function, x i is the i-th point in the function, x ¯ is the average of the function values, and a and b are the indices of the points that delimit the interval;
    S ( a , b ) = 1 N i = a b ( x i x ¯ ) 3 [ 1 N i = a b ( x i x ¯ ) 2 ] 3 .
  • Area ratio (R): From the asymmetry, the ratio between the halves of the area under the curve must also present similar values for similar classes. This descriptor was obtained through Equation (11), with a and b indicating the points that delimit the interval;
    R ( a , b ) = A ( b 2 + 1 , b ) A ( a , b 2 ) .
  • Maximum point: It indicates the value in the largest heterogeneous area of the curve. Thus, images from the same class can present similar values, for both f ( x ) and x. Totally different values are expected for different classes.
In summary, the fractal descriptors were organized into a feature vector with 116 descriptors: 20 fractal dimensions; 20 lacunarities; 20 average numbers of clusters present in each box; 20 average ratios covered by the largest cluster on each box; 20 percentages of boxes percolated; and 16 curve descriptors, 4 for each of the functions L ( r ) , C ( r ) , Q ( r ) and P ( r ) .

2.3.2. Deep-Learned Features

In this study, the deep-learned features were extracted from five CNN architectures: DenseNet-121 [58], EfficientNet-b2 [65], Inception-V3 [59], ResNet-50 [60] and VGG-19 [61]. The results were five vectors of deep-learned features. It is important to highlight that in general, the values in the initial layers of a CNN define the quantification of local patterns, such as shape, edge and color. The deeper layers are useful for identifying global patterns, such as texture and semantics [74]. Therefore, in order to explore the global patterns, the normalization layer after the last dense block was chosen from the DenseNet-121 architecture. This layer provided 1024 values. From the EfficientNet-b2, Inception-V3 and ResNet-50 architectures, each final average pooling layer contributed 1408, 2048 and 2048 values, respectively. From VGG-19, the extraction occurred in the last fully connected layer before the output, providing a vector with 4096 values.

2.4. Step 3—Feature Ensemble

This step was defined to organize the handcrafted and deep-learned descriptors with the corresponding ensembles, concatenating the features in order to analyze their discriminative capabilities in each of the H&E datasets [84]. It is important to emphasize that each representation obtained through xAI techniques was quantified with fractal approaches (Section 2.3.1), completing the handcrafted set. For each of the five CNN architectures, two feature vectors of xAI representations were defined via fractal techniques, with one from the Grad-CAM images and the other from the LIME images. Each vector considered 116 descriptors (described in Section 2.3.1): 20 fractal dimensions; 20 lacunarity values; 20 average numbers of clusters present in each box; 20 average ratios covered by the largest cluster on each box; 20 percentages of boxes that percolated; and 16 curve descriptors, 4 for each of the functions L ( r ) , C ( r ) , Q ( r ) and P ( r ) .
In this context, the descriptors explored here were organized into 55 distinct compositions, with 16 individual compositions according to their origins and the number of available descriptors. Among the vectors, three groups of distinct origins were defined: handcrafted, deep-learned and xAI. The handcrafted vectors were obtained through the fractal techniques (F), indicating 116 descriptors. The deep-learned vectors were obtained via DenseNet-121 (D) with 1024 features, EfficientNet-b2 (E) with 1408 features, Inception-V3 (I) with 2048 features, ResNet-50 (R) with 2048 features and VGG-19 (V) with 4096 features. Finally, the xAI vectors were composed of handcrafted descriptors obtained through the application of fractal techniques. In this case, the xAI vectors were defined with 116 descriptors, organized into two groups: those resulting from Grad-CAM explanations, with Grad-CAM DenseNet-121 ( D C A M ), Grad-CAM EfficientNet-b2 ( E C A M ), Grad-CAM Inception-V3 ( I C A M ), Grad-CAM ResNet-50 ( R C A M ) and Grad-CAM VGG-19 ( V C A M ); and those resulting from the LIME explanations, with LIME DenseNet-121 ( D L I M E ), LIME EfficientNet-b2 ( E L I M E ), LIME Inception-V3 ( I L I M E ), LIME ResNet-50 ( R L I M E ) and LIME VGG-19 ( V L I M E ). The 16 vector compositions with their origins and descriptor numbers are indicated in Table 3.
The remaining 39 feature vectors were created with ensembles of these individual vectors by aggregation, as illustrated in Figure 5. The ensembles were divided into three groups: an ensemble of handcrafted and deep-learned features; an ensemble of deep-learned features; and an ensemble of xAI features. The first group was composed of six vectors, five of them associating the fractal vector with each of the deep-learned features from a CNN, and the other was composed via fractal vector with all the deep-learned features from all CNN models, as presented in Table 4. The ensemble group of deep-learned features was composed of 11 vectors, with 10 permutations of deep-learned features from 2 distinct architectures and one with all available deep-learned features, as shown in Table 5. The last group contained 22 vectors with ensembles of xAI representation descriptors, organized according to the methodology applied to deep-learned ensembles (see Table 6).

2.5. Step 4—Feature Selection

The results achieved in Step 3 were vectors with high dimensions, which can lead to classifications with overfitting [43]. In order to avoid this problem, the solution was to apply a dimensionality reduction based on the ranking of the ReliefF technique [85]. This technique was considered to control the number of descriptors under analysis and due to its success found in other studies [45,47,86]. It is important to highlight that the analyzed vectors had dimensions between 116 and 10,740. In this context, considering that the total samples available in our experiments ranged from 58 to 265 images (Table 2), the tests were carried out with totals between 5 and 25 descriptors, starting the analysis process with the maximum value of descriptors and exploring decrements of 5 descriptors, until reaching the smallest dimension [83].

2.6. Step 5—Classifier Ensemble and Evaluation Metrics

The last step consisted of carrying out the classifications based on the vectors with an ensemble of heterogeneous classifiers, representative of different categories: support vector machine (SVM) [87], based on functions; naive bayes [88], based on probabilities; random forest [89], based on decision trees; and instance-based K-nearest neighbors [90]. The classifications were combined based on the average of probabilities. In this strategy, the average of all probabilities resulting from each of the classifiers was calculated, so the class with the highest average was selected as the answer (assigned class).
The results were analyzed considering the accuracy metric capable of indicating the global performance of the model (among all the classifications, how many the model classified correctly) [91]. Accuracy was determined through Equation (14), where true positives ( T P ) are the positive values that the model correctly classified as positive; true negatives ( T N ) are the negative values that the model correctly classified as negative; false positives ( F P ) are the negative values that the model incorrectly classified as positive; and false negatives ( F N ) are the positive values that the model incorrectly classified as negative.
A c c u r a c y = T P + V N T P + F P + T N + F N .
In addition, the classification process also considered the k-fold cross-validation approach, with k = 10 . It is important to highlight that the descriptor selection process was applied to each training set of each fold. The selected descriptors were used to perform the classification on the corresponding test set. This technique was applied to reduce the possibility of overfitting occurring, since the result was obtained through the classification average [92]. Figure 6 illustrates this cross-validation process with the selection of descriptors.

2.7. Software Packages and Execution Environment

In the proposed scheme, the Pytorch 1.9.0 [71] machine learning library was used to define the convolutional network models. Grad-CAM explanations were obtained via the pytorch-grad-cam 1.3.1 library [93]. The LIME explanations were obtained using the lime 0.2.0.1 library [94], implemented in python [41]. The application of each model was carried out in a remote environment on Google Colab, which provided a free console for executing codes in Python3 and a GPU for processing with 12 GB of available memory. The fractal descriptors were implemented and executed in parts in Matlab R2019b [95] and Python3. Feature selection and classification were processes carried out using the Weka v3.8.5 platform [96]. For this last part, the executions were carried out using an Intel® Core™ i5-8265U 1.60 GHz processor and 8 GB of RAM (Intel, Santa Clara, CA, USA).

3. Results and Discussion

The proposed scheme was tested on datasets of histological images (Section 2.1), with comparisons of benign versus malignant (UCSB and CR datasets), healthy versus severe (OED dataset), and male versus female (LG dataset). A performance overview is presented in Section 3.2, with results through the main associations against those from CNNs applied directly to classify the samples and from studies available in the specialized literature.

3.1. Feature Ensemble Performance

The combinations were tested according to details presented in Section 2.4, and the 10 highest accuracy values using the smallest number of descriptors were computed to define the average performance in each category. The average performance values are in Table 7 with the highest values in bold.
From the average values, it is noted that the ensemble of deep-learned features provided the highest values in all datasets, with an emphasis on the CR dataset (average accuracy of 100 % ). In the other datasets, this type of solution indicated accuracy values of 99.66% (LG), 97.23% (OED) and 92.93% (UCSB). Moreover, these facts confirm the feasibility of using the transfer learning strategy in order to obtain relevant compositions for the analysis of histological images. In addition, the ensemble of handcrafted with deep-learned features was the second association capable of indicating relevant accuracy values for the CR (99.76%), LG (99.02%) and OED (96.49%) datasets. Classifications using descriptors based on the xAI representations presented less expressive results, with average accuracy values between 78% and 89% (approximately), but with clear indications about the discriminative potential that can be explored on other research fronts. These facts are relevant contributions to support solutions without overfitting and more robust baseline schemes commonly explored for the improvement of CAD systems focused on H&E images.
These results were analyzed with the Friedman test in order to verify whether there were statistically significant differences between the compositions. The p-value was contrasted against α = 0.05 . Thus, if p-value < α , the difference was considered significant. The test was based on the averages of the first 10 results present in each descriptor category. Thus, taking into account a comparison among the six types of feature vectors, the p-value was 0.0033 , indicating statistically significant differences, with an emphasis on the ensemble of deep-learned features against handcrafted, xAI and xAI ensemble. Each pairwise p-value is displayed in Table 8. The differences between the obtained results via an ensemble of deep-learned features and those from an ensemble of handcrafted with deep-learned features were not statistically significant. In addition, the Friedman test provided an average ranking between the performance values of the main associations tested here. The best-ranked combination was the ensemble of deep-learned features, with average accuracy values ranging from 92.93% (UCSB) to 100% (CR).

3.1.1. Details of the Top 10 Solutions

The best combination was defined based on the highest accuracy value using the smallest number of descriptors. Considering this criterion, the obtained rankings are displayed in Table 9, Table 10, Table 11 and Table 12 with the first 10 solutions for the CR, LG, OED and UCSB datasets, respectively. It is noted that the average ranking of the solution validates its position.
In relation to the data collected from the CR dataset (Table 9), it is possible to observe that the 10 best results indicated an accuracy of 100 % . It is noted that the descriptors from the DenseNet-121 (D) and EfficientNet-b2 (E) models were those that contributed the most to these results. In these cases, the vectors were defined with a maximum of 20 descriptors, minimizing overfitting.
When the LG dataset was considered (Table 10), it was observed that only one combination indicated an accuracy value of 100%: the ensemble of deep-learned features from the DenseNet-121 (D) and ResNet-50 (R), exploring only 25 descriptors. The contribution of descriptors from the DenseNet-121 was a highlight, present in all vectors of the best results.
For the OED dataset (see Table 11), the highest accuracy was 97.97 % with an ensemble of descriptors from the Inception-V3 (I) and VGG-19 (V) architectures, using a reduced number of attributes (20 features). Another combination that achieved the same performance considered the descriptors from the EfficientNet-b2 (E) and Inception-V3 (I) models. However, this last combination involved 25 attributes. In addition, Inception-V3 (I) was the architecture that contributed the most to the best results, followed by DenseNet-121 (D). The main vectors were also defined with up to 25 descriptors.
Finally, considering the UCSB dataset, which is shown in Table 12, it is noted that the highest accuracy value was 94.83 % . This performance was defined via the association of deep-learned features from DenseNet-121 (D) with EfficientNet-b2 (E) exploring a total of 25 descriptors. From the top 10 solutions, it is observed that nine compositions were defined based on the DenseNet-121 or ResNet-50 architectures. The total number of descriptors present in the vectors follows the previously found pattern consisting of solutions with a maximum of 25 attributes.

3.1.2. Feature Summary

The best result was defined through the ensembles of deep-learned features, considering the top 10 solutions for each H&E dataset. For the CR dataset (Figure 7a), it is observed that the compositions were defined mainly based on the EfficientNet-b2 (E) model, representing from 60% to 90% of the features present in 7 of the top 10 solutions. Regarding the LG dataset (Figure 7b), the obtained vectors from the DenseNet-121 (D) architecture predominated (9 of the 10 compositions), indicating from 48% to 100% of the total features in each ensemble. When the OED dataset (Figure 7c) was considered, the highest frequency was obtained through Inception-V3 (I), present in 7 of the 10 compositions, but with a smaller presence of features in the ensembles (from 10% to 40%). Finally, the compositions for the UCSB dataset (Figure 7d) were based on more homogeneous compositions among DenseNet-121 (D), EfficientNet-b2 (E) and ResNet-50 (R). However, in each ensemble, there was a higher incidence of features from the ResNet-50 (R) architecture (from 73% to 100% of the total descriptors).
Based on the details presented previously, some patterns and/or behaviors were observed. The ensembles explored here effectively contributed to a distinction between the H&E images, surpassing the results provided via single-source descriptors. The ensemble of deep-learned features supported the main solutions, with highlights, in terms of occurrence, for the descriptors extracted from the DenseNet-121 (present in 28 of the 40 vectors) and EfficientNet-b2 (present in 19 of the 40 vectors) models. Finally, considering the dimensionalities of the feature vectors, compositions involving 10 to 25 deep-learned features were sufficient to determine the top 10 solutions.

3.2. Proposed Scheme versus Fine-Tuned CNN Classifications

In order to verify possible gains of the proposed scheme in relation to CNN models applied directly to each H&E dataset, the performance values of the architectures (DenseNet-121 [58], EfficientNet-b2 [65], Inception -V3 [59], ResNet-50 [60] and VGG-19 [61]) were collected after the fine-tuning process. The methodological details were presented in Section 2.2. In addition, Figure 8, Figure 9, Figure 10 and Figure 11 show the best accuracy rates via CNNs against those of the proposed scheme. It is important to highlight that other experiments could be defined, involving different conditions and limits, but those carried out here provided sufficient information to support the contributions of our investigation.
When image classifications were given directly by the CNN architectures, the DenseNet-121 and EfficientNet-b2 models indicated the lowest performance values, ranging from 67.81 % to 79.26 % and from 54.35 % to 80.35 % , respectively. The VGG-19 and ResNet-50 networks provided the best performance: LG ( 88.29 % ) and OED ( 97.59 % ) classified via VGG-19; CR ( 91.26 % ) and UCSB ( 88.52 % ) through the ResNet-50. However, these accuracy values were lower than those achieved through the proposed scheme, exploring an ensemble of deep-learned features with four heterogeneous classifiers. Also, through the information detailed here, it was noted that the classification performance via a CNN model was not directly related to that achieved through deep-learned features used in an independent process. Most of the best results explored features from the DenseNet-121 and EfficientNet-b2 networks, indicating the lowest performance. These results confirm the contributions obtained in this study.

3.3. Performance Overview in Relation to the Literature

Different techniques are available in the specialized literature in order to investigate patterns in histological images, such as for the CR, LG, OED and UCSB datasets [12,13,19,55,69,97,98,99,100]. Therefore, an illustrative overview is important to show the quality of our proposal, indicated in Table 13, Table 14, Table 15 and Table 16 for each of the datasets.
From this illustrative overview, it is possible to conclude that the proposed scheme provided solutions that surpassed a single type of descriptor or even other relevant associations [19,55,69,97,99,100,106,112]. Moreover, the computational scheme was capable of indicating optimized ensembles with the best results for classifying colorectal cancer, liver tissue and oral dysplasia, considering a maximum of 25 descriptors. Finally, these solutions without overfitting and with more robust baseline schemes are for improving CAD systems focused on H&E images. These contributions complement the proposal presented here in this illustrative overview.

4. Conclusions

In this work, a computational scheme was developed in order to define the main ensembles of descriptors for the study of histological images, exploring their ranking based on the ReliefF algorithm with a robust ensemble of classifiers (four heterogeneous algorithms). The handcrafted descriptors were established from multiscale and multidimensional fractal techniques (fractal dimension, lacunarity and percolation) and applied to quantify H&E images and their Grad-CAM and LIME representations. The deep-learned features were obtained from multiple CNN architectures, considering the transfer learning strategy. The experiments were carried out on H&E images, representative of breast cancer, colorectal cancer, oral dysplasia and liver tissue.
From the results, the ensemble of deep-learned features provided the highest values in all datasets, with accuracy rates of 94.83 % (UCSB), 97.97 % (OED) and 100 % (CR and LG), exploring a reduced number of features (up to 25 attributes). The descriptors were mainly obtained from the DenseNet-121 and EfficientNet-b2 architectures. In addition, the proposed scheme also indicated that handcrafted ensembles with deep-learned features provided expressive distinctions in the contexts of multiple histological images, with accuracy rates of 99.76% (CR), 99.02% (LG) and 96.49% (OED). This type of composition indicated better performance values than those achieved through individualized analyses, commonly observed in the specialized literature, whether exploring only deep learning [102,103,106,112] or handcrafted techniques [33,69,74,107]. In both categories of ensembles, this study provided useful details and conditions for the community interested in the development and improvement of models for classifying patterns in H&E samples. In relation to the experiments exploring xAI representations, the results were less expressive, with average accuracy values ranging from 78% to 89%. On the other hand, this type of composition achieved results similar to those of directly applied networks, responsible for providing the xAI representations. Thus, we believe that there are still several study avenues to understand the full information capacity present in this type of representation, seeking to improve CAD system designs focused on histological images.
The best solutions were analyzed in relation to the results obtained from consolidated machine-learning techniques, directly applying CNN models to classify the histological datasets. This process considered the DenseNet-121, EfficientNet-b2, Inception-V3, VGG-19 and ResNet-50 architectures. The results were accuracy values from 54.35 to 97.59 . The VGG-19 and ResNet-50 networks indicated the best rates: LG ( 88.29 % ) and OED ( 97.59 % ) via VGG-19; CR ( 91.26 % ) and UCSB ( 88.52 % ) through ResNet-50. These performance values were lower than those achieved through the best solutions with the proposed scheme. When an illustrative overview was considered in relation to the specialized literature, it was possible to conclude that the proposed scheme provided solutions that surpassed a single type of descriptor or even other relevant associations [55,69,97,106,112]. Moreover, the computational scheme was capable of indicating optimized ensembles with the best results for classifying colorectal cancer, liver tissue and oral dysplasia. Therefore, these conditions highlight the ability of the proposed scheme to present solutions without overfitting and a more robust baseline scheme, with the necessary details for the analysis and testing of CAD systems, focused on H&E samples. Regarding classifications involving representative samples of breast cancer (UCSB dataset), the proposed scheme provided a lower performance, indicating a possible limit of the main solution in this context.
In future works, it is intended to (1) expand the number of handcrafted techniques to quantify H&E images and their representations, especially to define the possible limits involving combinations via xAI; (2) carry out new tests after applying adjustments to the parameters of the CNN architectures in order to verify their impacts on the xAI representations and corresponding quantification; (3) explore multiview learning approaches to complement multiple representations, including investigations into possible gains after applying learning enrichment strategies with fractal techniques.

Author Contributions

Conceptualization, methodology, validation, formal analysis, investigation, writing—original draft preparation, writing—review and editing, L.H.d.C.L., G.F.R., T.A.A.T., P.R.d.F., A.M.L., S.V.C., A.B.S., M.Z.d.N. and L.A.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the: Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (CAPES)—Finance Code 001; National Council for Scientific and Technological Development CNPq (grants #132940/2019-1, #313643/2021-0 and #311404/2021-9); the State of Minas Gerais Research Foundation—FAPEMIG (grant #APQ-00578-18); São Paulo Research Foundation—FAPESP (grant #2022/03020-1).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
H&EHematoxylin and eosin
CADComputer-aided diagnosis
CNNConvolutional neural network
xAIExplainable artificial intelligence
CAMClass activation mapping
Grad-CAMGradient-weighted class activation mapping
LIMELocal interpretable model-agnostic explanations
UCSBBreast cancer dataset
CRColorectal cancer dataset
LGLiver tissue dataset
OEDOral epithelial dysplasia dataset
SGDMStochastic gradient descent
FFractal techniques
DDenseNet-121
EEfficientNet-b2
IInception-V3
RResNet-50
VVGG-19
D C A M Grad-CAM representation via DenseNet-121
E C A M Grad-CAM representation via EfficientNet-b2
R C A M Grad-CAM representation via ResNet-50
I C A M Grad-CAM representation via Inception-V3
V C A M Grad-CAM representation via VGG-19
D L I M E LIME representation via DenseNet-121
E L I M E LIME representation via EfficientNet-b2
I L I M E LIME representation via Inception-V3
R L I M E LIME representation via ResNet-50
V L I M E LIME representation via VGG-19
VPTrue positives
VNTrue negatives
FPFalse positives
FNFalse negatives

References

  1. Gurcan, M.N.; Boucheron, L.E.; Can, A.; Madabhushi, A.; Rajpoot, N.M.; Yener, B. Histopathological image analysis: A review. IEEE Rev. Biomed. Eng. 2009, 2, 147–171. [Google Scholar] [CrossRef] [PubMed]
  2. Titford, M. A short history of histopathology technique. J. Histotechnol. 2006, 29, 99–110. [Google Scholar] [CrossRef]
  3. Angel Arul Jothi, J.; Mary Anita Rajam, V. A survey on automated cancer diagnosis from histopathology images. Artif. Intell. Rev. 2017, 48, 31–81. [Google Scholar] [CrossRef]
  4. Shmatko, A.; Ghaffari Laleh, N.; Gerstung, M.; Kather, J.N. Artificial intelligence in histopathology: Enhancing cancer research and clinical oncology. Nat. Cancer 2022, 3, 1026–1038. [Google Scholar] [CrossRef] [PubMed]
  5. Carleton, H.M.; Drury, R.A.B.; Wallington, E.A. Carleton’s Histological Technique; Oxford University Press: New York, NY, USA, 1980. [Google Scholar]
  6. Titford, M. The long history of hematoxylin. Biotech. Histochem. 2005, 80, 73–78. [Google Scholar] [CrossRef] [PubMed]
  7. Chan, J.K.C. The Wonderful Colors of the Hematoxylin–Eosin Stain in Diagnostic Surgical Pathology. Int. J. Surg. Pathol. 2014, 22, 12–32. [Google Scholar] [CrossRef] [PubMed]
  8. Frick, C.; Rumgay, H.; Vignat, J.; Ginsburg, O.; Nolte, E.; Bray, F.; Soerjomataram, I. Quantitative estimates of preventable and treatable deaths from 36 cancers worldwide: A population-based study. Lancet Glob. Health 2023, 11, e1700–e1712. [Google Scholar] [CrossRef]
  9. Higgins, C. Applications and challenges of digital pathology and whole slide imaging. Biotech. Histochem. 2015, 90, 341–347. [Google Scholar] [CrossRef]
  10. Hu, W.; Li, X.; Li, C.; Li, R.; Jiang, T.; Sun, H.; Huang, X.; Grzegorzek, M.; Li, X. A state-of-the-art survey of artificial neural networks for whole-slide image analysis: From popular convolutional neural networks to potential visual transformers. Comput. Biol. Med. 2023, 161, 107034. [Google Scholar] [CrossRef]
  11. He, L.; Long, L.R.; Antani, S.; Thoma, G.R. Histology image analysis for carcinoma detection and grading. Comput. Methods Programs Biomed. 2012, 107, 538–556. [Google Scholar] [CrossRef]
  12. Roberto, G.F.; Lumini, A.; Neves, L.A.; do Nascimento, M.Z. Fractal Neural Network: A new ensemble of fractal geometry and convolutional neural networks for the classification of histology images. Expert Syst. Appl. 2021, 166, 114103. [Google Scholar] [CrossRef]
  13. Nanni, L.; Brahnam, S.; Ghidoni, S.; Maguolo, G. General purpose (GenP) bioimage ensemble of handcrafted and learned features with data augmentation. arXiv 2019, arXiv:1904.08084. [Google Scholar]
  14. Wang, W.; Yang, Y.; Wang, X.; Wang, W.; Li, J. Development of convolutional neural network and its application in image classification: A survey. Opt. Eng. 2019, 58, 040901. [Google Scholar] [CrossRef]
  15. Wang, Z.; Li, M.; Wang, H.; Jiang, H.; Yao, Y.; Zhang, H.; Xin, J. Breast Cancer Detection Using Extreme Learning Machine Based on Feature Fusion With CNN Deep Features. IEEE Access 2019, 7, 105146–105158. [Google Scholar] [CrossRef]
  16. Mahbod, A.; Schaefer, G.; Ellinger, I.; Ecker, R.; Pitiot, A.; Wang, C. Fusing fine-tuned deep features for skin lesion classification. Comput. Med. Imaging Graph. 2019, 71, 19–29. [Google Scholar] [CrossRef] [PubMed]
  17. Alinsaif, S.; Lang, J. Histological Image Classification using Deep Features and Transfer Learning. In Proceedings of the 2020 17th Conference on Computer and Robot Vision (CRV), Ottawa, ON, Canada, 13–15 May 2020; pp. 101–108. [Google Scholar] [CrossRef]
  18. Zerouaoui, H.; Idri, A.; El Alaoui, O. A new approach for histological classification of breast cancer using deep hybrid heterogenous ensemble. Data Technol. Appl. 2023, 57, 245–278. [Google Scholar] [CrossRef]
  19. de Oliveira, C.I.; do Nascimento, M.Z.; Roberto, G.F.; Tosta, T.A.; Martins, A.S.; Neves, L.A. Hybrid models for classifying histological images: An association of deep features by transfer learning with ensemble classifier. Multimed. Tools Appl. 2023, 1–24. [Google Scholar] [CrossRef]
  20. Taino, D.F.; Ribeiro, M.G.; Roberto, G.F.; Zafalon, G.F.; do Nascimento, M.Z.; Tosta, T.A.; Martins, A.S.; Neves, L.A. Analysis of cancer in histological images: Employing an approach based on genetic algorithm. Pattern Anal. Appl. 2021, 24, 483–496. [Google Scholar] [CrossRef]
  21. Torrey, L.; Shavlik, J. Transfer learning. In Handbook of Research on Machine Learning Applications and Trends: Algorithms, Methods, and Techniques; IGI Global: Hershey, PA, USA, 2010; pp. 242–264. [Google Scholar]
  22. Nanni, L.; Ghidoni, S.; Brahnam, S. Handcrafted vs. non-handcrafted features for computer vision classification. Pattern Recognit. 2017, 71, 158–172. [Google Scholar] [CrossRef]
  23. Hagerty, J.R.; Stanley, R.J.; Almubarak, H.A.; Lama, N.; Kasmi, R.; Guo, P.; Drugge, R.J.; Rabinovitz, H.S.; Oliviero, M.; Stoecker, W.V. Deep learning and handcrafted method fusion: Higher diagnostic accuracy for melanoma dermoscopy images. IEEE J. Biomed. Health Inform. 2019, 23, 1385–1391. [Google Scholar] [CrossRef]
  24. Wei, L.; Su, R.; Wang, B.; Li, X.; Zou, Q.; Gao, X. Integration of deep feature representations and handcrafted features to improve the prediction of N6-methyladenosine sites. Neurocomputing 2019, 324, 3–9. [Google Scholar] [CrossRef]
  25. Hasan, A.M.; Jalab, H.A.; Meziane, F.; Kahtan, H.; Al-Ahmad, A.S. Combining Deep and Handcrafted Image Features for MRI Brain Scan Classification. IEEE Access 2019, 7, 79959–79967. [Google Scholar] [CrossRef]
  26. Li, S.; Xu, P.; Li, B.; Chen, L.; Zhou, Z.; Hao, H.; Duan, Y.; Folkert, M.; Ma, J.; Huang, S.; et al. Predicting lung nodule malignancies by combining deep convolutional neural network and handcrafted features. Phys. Med. Biol. 2019, 64, 175012. [Google Scholar] [CrossRef] [PubMed]
  27. Tripathi, S.; Singh, S.K. Ensembling handcrafted features with deep features: An analytical study for classification of routine colon cancer histopathological nuclei images. Multimed. Tools Appl. 2020, 79, 34931–34954. [Google Scholar] [CrossRef]
  28. Fortin, C.S.; Kumaresan, R.; Ohley, W.J.; Hoefer, S. Fractal dimension in the analysis of medical images. IEEE Eng. Med. Biol. Mag. 1992, 11, 65–71. [Google Scholar] [CrossRef]
  29. Aralica, G.; Ivelj, M.Š.; Pačić, A.; Baković, J.; Periša, M.M.; Krištić, A.; Konjevoda, P. Prognostic Significance of Lacunarity in Preoperative Biopsy of Colorectal Cancer. Pathol. Oncol. Res. 2020, 26, 2567–2576. [Google Scholar] [CrossRef]
  30. Jain, A.K.; Duin, R.P.W.; Mao, J. Statistical pattern recognition: A review. IEEE Trans. Pattern Anal. Mach. Intell. 2000, 22, 4–37. [Google Scholar] [CrossRef]
  31. Qin, J.; Puckett, L.; Qian, X. Image Based Fractal Analysis for Detection of Cancer Cells. In Proceedings of the 2020 IEEE International Conference on Bioinformatics and Biomedicine (BIBM), Seoul, Republic of Korea, 16–19 December 2020; pp. 1482–1486. [Google Scholar] [CrossRef]
  32. Roberto, G.F.; Neves, L.A.; Nascimento, M.Z.; Tosta, T.A.; Longo, L.C.; Martins, A.S.; Faria, P.R. Features based on the percolation theory for quantification of non-Hodgkin lymphomas. Comput. Biol. Med. 2017, 91, 135–147. [Google Scholar] [CrossRef]
  33. Roberto, G.F.; Nascimento, M.Z.; Martins, A.S.; Tosta, T.A.; Faria, P.R.; Neves, L.A. Classification of breast and colorectal tumors based on percolation of color normalized images. Comput. Graph. 2019, 84, 134–143. [Google Scholar] [CrossRef]
  34. Candelero, D.; Roberto, G.F.; Do Nascimento, M.Z.; Rozendo, G.B.; Neves, L.A. Selection of cnn, haralick and fractal features based on evolutionary algorithms for classification of histological images. In Proceedings of the 2020 IEEE International Conference on Bioinformatics and Biomedicine (BIBM), Seoul, Republic of Korea, 16–19 December 2020; pp. 2709–2716. [Google Scholar]
  35. Khan, A.; Sohail, A.; Zahoora, U.; Qureshi, A.S. A survey of the recent architectures of deep convolutional neural networks. Artif. Intell. Rev. 2020, 53, 5455–5516. [Google Scholar] [CrossRef]
  36. Khedkar, S.; Subramanian, V.; Shinde, G.; Gandhi, P. Explainable AI in healthcare. In Proceedings of the 2nd International Conference on Advances in Science & Technology (ICAST), Mumbai, India, 8–9 April 2019. [Google Scholar]
  37. Wells, L.; Bednarz, T. Explainable ai and reinforcement learning—A systematic review of current approaches and trends. Front. Artif. Intell. 2021, 4, 550030. [Google Scholar] [CrossRef] [PubMed]
  38. Samek, W.; Müller, K.R. Towards explainable artificial intelligence. In Explainable AI: Interpreting, Explaining and Visualizing Deep Learning; Springer: Berlin/Heidelberg, Germany, 2019; pp. 5–22. [Google Scholar]
  39. Gunning, D.; Aha, D. DARPA’s explainable artificial intelligence (XAI) program. AI Mag. 2019, 40, 44–58. [Google Scholar]
  40. Zhou, B.; Khosla, A.; Lapedriza, A.; Oliva, A.; Torralba, A. Learning deep features for discriminative localization. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Las Vegas, NV, USA, 27–30 June 2016; pp. 2921–2929. [Google Scholar]
  41. Ribeiro, M.T.; Singh, S.; Guestrin, C. “Why should i trust you?” Explaining the predictions of any classifier. In Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, San Francisco, CA, USA, 13–17 August 2016; pp. 1135–1144. [Google Scholar]
  42. Bishop, C.M. Pattern Recognition and Machine Learning; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  43. Watanabe, S. Pattern Recognition: Human and Mechanical; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 1985. [Google Scholar]
  44. Bolón-Canedo, V.; Sánchez-Marono, N.; Alonso-Betanzos, A.; Benítez, J.M.; Herrera, F. A review of microarray datasets and applied feature selection methods. Inf. Sci. 2014, 282, 111–135. [Google Scholar] [CrossRef]
  45. Urbanowicz, R.J.; Meeker, M.; La Cava, W.; Olson, R.S.; Moore, J.H. Relief-based feature selection: Introduction and review. J. Biomed. Inform. 2018, 85, 189–203. [Google Scholar] [CrossRef]
  46. Manhrawy, I.I.; Qaraad, M.; El-Kafrawy, P. Hybrid feature selection model based on relief-based algorithms and regulizer algorithms for cancer classification. Concurr. Comput. Pract. Exp. 2021, 33, e6200. [Google Scholar] [CrossRef]
  47. Ghosh, P.; Azam, S.; Jonkman, M.; Karim, A.; Shamrat, F.J.M.; Ignatious, E.; Shultana, S.; Beeravolu, A.R.; De Boer, F. Efficient prediction of cardiovascular disease using machine learning algorithms with relief and LASSO feature selection techniques. IEEE Access 2021, 9, 19304–19326. [Google Scholar] [CrossRef]
  48. Zebari, R.; Abdulazeez, A.; Zeebaree, D.; Zebari, D.; Saeed, J. A comprehensive review of dimensionality reduction techniques for feature selection and feature extraction. J. Appl. Sci. Technol. Trends 2020, 1, 56–70. [Google Scholar] [CrossRef]
  49. Li, M.; Ma, X.; Chen, C.; Yuan, Y.; Zhang, S.; Yan, Z.; Chen, C.; Chen, F.; Bai, Y.; Zhou, P.; et al. Research on the auxiliary classification and diagnosis of lung cancer subtypes based on histopathological images. IEEE Access 2021, 9, 53687–53707. [Google Scholar] [CrossRef]
  50. Burçak, K.C.; Uğuz, H. A New Hybrid Breast Cancer Diagnosis Model Using Deep Learning Model and ReliefF. Trait. Signal 2022, 39, 521–529. [Google Scholar] [CrossRef]
  51. Silva, A.B.; De Oliveira, C.I.; Pereira, D.C.; Tosta, T.A.; Martins, A.S.; Loyola, A.M.; Cardoso, S.V.; De Faria, P.R.; Neves, L.A.; Do Nascimento, M.Z. Assessment of the association of deep features with a polynomial algorithm for automated oral epithelial dysplasia grading. In Proceedings of the 2022 35th SIBGRAPI Conference on Graphics, Patterns and Images (SIBGRAPI), Natal, Brazil, 24–27 October 2022; Volume 1, pp. 264–269. [Google Scholar]
  52. Dietterich, T.G. Ensemble methods in machine learning. In Proceedings of the International Workshop on Multiple Classifier Systems, Cagliari, Italy, 21–23 June 2000; Springer: Berlin/Heidelberg, Germany, 2000; pp. 1–15. [Google Scholar]
  53. Sagi, O.; Rokach, L. Ensemble learning: A survey. WIREs Data Min. Knowl. Discov. 2018, 8, e1249. [Google Scholar] [CrossRef]
  54. Lumini, A.; Nanni, L. Convolutional neural networks for ATC classification. Curr. Pharm. Des. 2018, 24, 4007–4012. [Google Scholar] [CrossRef] [PubMed]
  55. Kausar, T.; Wang, M.; Idrees, M.; Lu, Y. HWDCNN: Multi-class recognition in breast histopathology with Haar wavelet decomposed image based convolution neural network. Biocybern. Biomed. Eng. 2019, 39, 967–982. [Google Scholar] [CrossRef]
  56. Toğaçar, M.; Cömert, Z.; Ergen, B. Enhancing of dataset using DeepDream, fuzzy color image enhancement and hypercolumn techniques to detection of the Alzheimer’s disease stages by deep learning model. Neural Comput. Appl. 2021, 33, 9877–9889. [Google Scholar] [CrossRef]
  57. Maia, B.M.S.; de Assis, M.C.F.R.; de Lima, L.M.; Rocha, M.B.; Calente, H.G.; Correa, M.L.A.; Camisasca, D.R.; Krohling, R.A. Transformers, convolutional neural networks, and few-shot learning for classification of histopathological images of oral cancer. Expert Syst. Appl. 2023, 241, 122418. [Google Scholar] [CrossRef]
  58. Huang, G.; Liu, Z.; Van Der Maaten, L.; Weinberger, K.Q. Densely connected convolutional networks. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Honolulu, HI, USA, 21–26 July 2017; pp. 4700–4708. [Google Scholar]
  59. Szegedy, C.; Vanhoucke, V.; Ioffe, S.; Shlens, J.; Wojna, Z. Rethinking the inception architecture for computer vision. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Las Vegas, NV, USA, 27–30 June 2016; pp. 2818–2826. [Google Scholar]
  60. He, K.; Zhang, X.; Ren, S.; Sun, J. Deep residual learning for image recognition. In Proceedings of the IEEE conference on Computer Vision and Pattern Recognition, Las Vegas, NV, USA, 27–30 June 2016; pp. 770–778. [Google Scholar]
  61. Simonyan, K.; Zisserman, A. Very deep convolutional networks for large-scale image recognition. arXiv 2014, arXiv:1409.1556. [Google Scholar]
  62. Reddy, A.S.B.; Juliet, D.S. Transfer Learning with ResNet-50 for Malaria Cell-Image Classification. In Proceedings of the 2019 International Conference on Communication and Signal Processing (ICCSP), Tamil Nadu, India, 4–6 April 2019; pp. 0945–0949. [Google Scholar] [CrossRef]
  63. Ganguly, A.; Das, R.; Setua, S.K. Histopathological Image and Lymphoma Image Classification using customized Deep Learning models and different optimization algorithms. In Proceedings of the 2020 11th International Conference on Computing, Communication and Networking Technologies (ICCCNT), Kharagpur, India, 1–3 July 2020; pp. 1–7. [Google Scholar] [CrossRef]
  64. Al-Haija, Q.A.; Adebanjo, A. Breast Cancer Diagnosis in Histopathological Images Using ResNet-50 Convolutional Neural Network. In Proceedings of the 2020 IEEE International IOT, Electronics and Mechatronics Conference (IEMTRONICS), Vancouver, BC, Canada, 9–12 September 2020; pp. 1–7. [Google Scholar] [CrossRef]
  65. Tan, M.; Le, Q. Efficientnet: Rethinking model scaling for convolutional neural networks. In Proceedings of the International Conference on Machine Learning, PMLR, Long Beach, CA, USA, 9–15 June 2019; pp. 6105–6114. [Google Scholar]
  66. Gelasca, E.D.; Byun, J.; Obara, B.; Manjunath, B. Evaluation and Benchmark for Biological Image Segmentation. In Proceedings of the IEEE International Conference on Image Processing, San Diego, CA, USA, 12–15 October 2008. [Google Scholar]
  67. Sirinukunwattana, K.; Pluim, J.P.; Chen, H.; Qi, X.; Heng, P.A.; Guo, Y.B.; Wang, L.Y.; Matuszewski, B.J.; Bruni, E.; Sanchez, U.; et al. Gland segmentation in colon histology images: The glas challenge contest. Med. Image Anal. 2017, 35, 489–502. [Google Scholar] [CrossRef] [PubMed]
  68. AGEMAP. The Atlas of Gene Expression in Mouse Aging Project (AGEMAP). 2020.
  69. Silva, A.B.; Martins, A.S.; Tosta, T.A.A.; Neves, L.A.; Servato, J.P.S.; de Araújo, M.S.; de Faria, P.R.; do Nascimento, M.Z. Computational analysis of histological images from hematoxylin and eosin-stained oral epithelial dysplasia tissue sections. Expert Syst. Appl. 2022, 193, 116456. [Google Scholar] [CrossRef]
  70. Deng, J.; Dong, W.; Socher, R.; Li, L.; Li, K.; Fei-Fei, L. ImageNet: A large-scale hierarchical image database. In Proceedings of the 2009 IEEE Conference on Computer Vision and Pattern Recognition, Miami, FL, USA, 20–25 June 2009; pp. 248–255. [Google Scholar] [CrossRef]
  71. Paszke, A.; Gross, S.; Massa, F.; Lerer, A.; Bradbury, J.; Chanan, G.; Killeen, T.; Lin, Z.; Gimelshein, N.; Antiga, L.; et al. PyTorch: An Imperative Style, High-Performance Deep Learning Library. In Advances in Neural Information Processing Systems 32; Curran Associates, Inc.: New York, NY, USA, 2019; pp. 8024–8035. [Google Scholar]
  72. torchvision.models. 2021. Available online: https://pytorch.org/vision/stable/models.html (accessed on 24 December 2023).
  73. Selvaraju, R.R.; Cogswell, M.; Das, A.; Vedantam, R.; Parikh, D.; Batra, D. Grad-cam: Visual explanations from deep networks via gradient-based localization. In Proceedings of the IEEE International Conference on Computer Vision, Venice, Italy, 22–29 October 2017; pp. 618–626. [Google Scholar]
  74. dos Santos, F.P.; Ponti, M.A. Alignment of Local and Global Features from Multiple Layers of Convolutional Neural Network for Image Classification. In Proceedings of the 2019 32nd SIBGRAPI Conference on Graphics, Patterns and Images (SIBGRAPI), Rio de Janeiro, Brazil, 28–31 October 2019; pp. 241–248. [Google Scholar]
  75. Van der Walt, S.; Schönberger, J.L.; Nunez-Iglesias, J.; Boulogne, F.; Warner, J.D.; Yager, N.; Gouillart, E.; Yu, T. scikit-image: Image processing in Python. PeerJ 2014, 2, e453. [Google Scholar] [CrossRef]
  76. Cantrell, C.D. Modern Mathematical Methods for Physicists and Engineers; Cambridge University Press: Cambridge, MA, USA, 2000. [Google Scholar]
  77. Allain, C.; Cloitre, M. Characterizing the lacunarity of random and deterministic fractal sets. Phys. Rev. A 1991, 44, 3552. [Google Scholar] [CrossRef]
  78. Ivanovici, M.; Richard, N. Fractal dimension of color fractal images. IEEE Trans. Image Process. 2010, 20, 227–235. [Google Scholar] [CrossRef]
  79. Hoshen, J.; Kopelman, R. Percolation and cluster distribution. I. Cluster multiple labeling technique and critical concentration algorithm. Phys. Rev. B 1976, 14, 3438. [Google Scholar] [CrossRef]
  80. Căliman, A.; Ivanovici, M. Psoriasis image analysis using color lacunarity. In Proceedings of the 2012 13th International Conference on Optimization of Electrical and Electronic Equipment (OPTIM), Brasov, Romania, 24–26 May 2012; pp. 1401–1406. [Google Scholar]
  81. Segato dos Santos, L.F.; Neves, L.A.; Rozendo, G.B.; Ribeiro, M.G.; Zanchetta do Nascimento, M.; Azevedo Tosta, T.A. Multidimensional and fuzzy sample entropy (SampEnMF) for quantifying H&E histological images of colorectal cancer. Comput. Biol. Med. 2018, 103, 148–160. [Google Scholar] [CrossRef] [PubMed]
  82. Tosta, T.A.; Bruno, D.O.; Longo, L.C.; do Nascimento, M.Z. Colour Feature Extraction and Polynomial Algorithm for Classification of Lymphoma Images. In Proceedings of the Progress in Pattern Recognition, Image Analysis, Computer Vision, and Applications: 24th Iberoamerican Congress, CIARP 2019, Havana, Cuba, 28–31 October 2019; Springer Nature: Berlin/Heidelberg, Germany, 2019; Volume 11896, p. 262. [Google Scholar]
  83. Ribeiro, M.G.; Neves, L.A.; do Nascimento, M.Z.; Roberto, G.F.; Martins, A.S.; Azevedo Tosta, T.A. Classification of colorectal cancer based on the association of multidimensional and multiresolution features. Expert Syst. Appl. 2019, 120, 262–278. [Google Scholar] [CrossRef]
  84. Dasigi, V.; Mann, R.C.; Protopopescu, V.A. Information fusion for text classification—An experimental comparison. Pattern Recognit. 2001, 34, 2413–2425. [Google Scholar] [CrossRef]
  85. Kononenko, I.; Šimec, E.; Robnik-Šikonja, M. Overcoming the myopia of inductive learning algorithms with RELIEFF. Appl. Intell. 1997, 7, 39–55. [Google Scholar] [CrossRef]
  86. Demir, S.; Key, S.; Baygin, M.; Tuncer, T.; Dogan, S.; Belhaouari, S.B.; Poyraz, A.K.; Gurger, M. Automated knee ligament injuries classification method based on exemplar pyramid local binary pattern feature extraction and hybrid iterative feature selection. Biomed. Signal Process. Control 2022, 71, 103191. [Google Scholar] [CrossRef]
  87. Cortes, C.; Vapnik, V. Support-vector networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar] [CrossRef]
  88. Pedrini, H.; Schwartz, W.R. Análise de Imagens Digitais: Princípios, Algoritmos e Aplicações; Thomson Learning: Chicago, IL, USA, 2008. [Google Scholar]
  89. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef]
  90. Cunningham, P.; Delany, S.J. k-Nearest Neighbour Classifiers. arXiv 2020, arXiv:2004.04523. [Google Scholar]
  91. Martinez, E.; Louzada, F.; Pereira, B. A curva ROC para testes diagnósticos. Cad Saúde Coletiva 2003, 11, 7–31. [Google Scholar]
  92. Kohavi, R. A study of cross-validation and bootstrap for accuracy estimation and model selection. In Proceedings of the International Joint Conferences on Artificial Intelligence, Montreal, QC, Canada, 20–25 August 1995; Volume 14, pp. 1137–1145. [Google Scholar]
  93. Gildenblat, J. PyTorch Library for CAM Methods. 2021. Available online: https://github.com/jacobgil/pytorch-grad-cam (accessed on 24 December 2023).
  94. Ribeiro, M.T. Lime. 2016. Available online: https://github.com/marcotcr/lime (accessed on 20 December 2023).
  95. MATLAB, version 9.7.0 (R2019b); The MathWorks Inc.: Natick, MA, USA, 2019.
  96. Frank, E.; Hall, M.A.; Holmes, G.; Kirkby, R.; Pfahringer, B.; Witten, I.H. Weka: A machine learning workbench for data mining. In Data Mining and Knowledge Discovery Handbook: A Complete Guide for Practitioners and Researchers; Maimon, O., Rokach, L., Eds.; Springer: Berlin/Heidelberg, Germany, 2005; pp. 1305–1314. [Google Scholar]
  97. Tavolara, T.E.; Niazi, M.K.K.; Arole, V.; Chen, W.; Frankel, W.; Gurcan, M.N. A modular cGAN classification framework: Application to colorectal tumor detection. Sci. Rep. 2019, 9, 18969. [Google Scholar] [CrossRef] [PubMed]
  98. Yu, C.; Chen, H.; Li, Y.; Peng, Y.; Li, J.; Yang, F. Breast cancer classification in pathological images based on hybrid features. Multimed. Tools Appl. 2019, 78, 21325–21345. [Google Scholar] [CrossRef]
  99. Zhang, R.; Zhu, J.; Yang, S.; Hosseini, M.S.; Genovese, A.; Chen, L.; Rowsell, C.; Damaskinos, S.; Varma, S.; Plataniotis, K.N. HistoKT: Cross Knowledge Transfer in Computational Pathology. In Proceedings of the ICASSP 2022—2022 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Singapore, 22–27 May 2022; pp. 1276–1280. [Google Scholar]
  100. Azarmehr, N.; Shephard, A.; Mahmood, H.; Rajpoot, N.; Khurram, S.A. A neural architecture search based framework for segmentation of epithelium, nuclei and oral epithelial dysplasia grading. In Proceedings of the Medical Imaging with Deep Learning, Zürich, Switzerland, 6–8 July 2022. [Google Scholar]
  101. Dabass, M.; Vig, R.; Vashisth, S. Five-grade cancer classification of colon histology images via deep learning. In Communication and Computing Systems; CRC Press: Boca Raton, FL, USA, 2019; pp. 18–24. [Google Scholar]
  102. Sena, P.; Fioresi, R.; Faglioni, F.; Losi, L.; Faglioni, G.; Roncucci, L. Deep learning techniques for detecting preneoplastic and neoplastic lesions in human colorectal histological images. Oncol. Lett. 2019, 18, 6101–6107. [Google Scholar] [CrossRef] [PubMed]
  103. Bentaieb, A.; Hamarneh, G. Adversarial Stain Transfer for Histopathology Image Analysis. IEEE Trans. Med. Imaging 2018, 37, 792–802. [Google Scholar] [CrossRef] [PubMed]
  104. Awan, R.; Al-Maadeed, S.; Al-Saady, R.; Bouridane, A. Glandular structure-guided classification of microscopic colorectal images using deep learning. Comput. Electr. Eng. 2020, 85, 106450. [Google Scholar] [CrossRef]
  105. Di Ruberto, C.; Putzu, L.; Arabnia, H.; Quoc-Nam, T. A feature learning framework for histology images classification. In Emerging Trends in Applications and Infrastructures for Computational Biology, Bioinformatics, and Systems Biology: Systems and Applications; Elsevier B.V.: Amsterdam, The Netherlands, 2016; pp. 37–48. [Google Scholar]
  106. Andrearczyk, V.; Whelan, P.F. Deep learning for biomedical texture image analysis. In Proceedings of the Irish Machine Vision & Image Processing Conference, Kildare, Ireland, 30 August–1 September 2017; Irish Pattern Recognition & Classification Society (IPRCS): Dublin, Ireland, 2017. [Google Scholar]
  107. Watanabe, K.; Kobayashi, T.; Wada, T. Semi-supervised feature transformation for tissue image classification. PLoS ONE 2016, 11, e0166413. [Google Scholar] [CrossRef] [PubMed]
  108. Adel, D.; Mounir, J.; El-Shafey, M.; Eldin, Y.A.; El Masry, N.; AbdelRaouf, A.; Abd Elhamid, I.S. Oral epithelial dysplasia computer aided diagnostic approach. In Proceedings of the 2018 13th International Conference on Computer Engineering and Systems (ICCES), Cairo, Egypt, 18–19 December 2018; pp. 313–318. [Google Scholar]
  109. Krishnan, M.M.R.; Shah, P.; Choudhary, A.; Chakraborty, C.; Paul, R.R.; Ray, A.K. Textural characterization of histopathological images for oral sub-mucous fibrosis detection. Tissue Cell 2011, 43, 318–330. [Google Scholar] [CrossRef]
  110. Li, Y.; Xie, X.; Shen, L.; Liu, S. Reverse active learning based atrous DenseNet for pathological image classification. BMC Bioinform. 2019, 20, 445. [Google Scholar] [CrossRef]
  111. Feng, Y.; Zhang, L.; Yi, Z. Breast cancer cell nuclei classification in histopathology images using deep neural networks. Int. J. Comput. Assist. Radiol. Surg. 2018, 13, 179–191. [Google Scholar] [CrossRef]
  112. Papastergiou, T.; Zacharaki, E.I.; Megalooikonomou, V. Tensor decomposition for multiple-instance classification of high-order medical data. Complexity 2018, 2018, 8651930. [Google Scholar] [CrossRef]
  113. Araújo, T.; Aresta, G.; Castro, E.; Rouco, J.; Aguiar, P.; Eloy, C.; Polónia, A.; Campilho, A. Classification of breast cancer histology images using convolutional neural networks. PLoS ONE 2017, 12, e0177544. [Google Scholar] [CrossRef] [PubMed]
Figure 3. Examples of xAI representations based on the Grad-CAM technique, with (a) the original image, (b) the weights mapped to indicate the contribution of each pixel, (c) the mapping transformed into a heatmap and (d) the heatmap overlaying the original image.
Figure 3. Examples of xAI representations based on the Grad-CAM technique, with (a) the original image, (b) the weights mapped to indicate the contribution of each pixel, (c) the mapping transformed into a heatmap and (d) the heatmap overlaying the original image.
Entropy 26 00034 g003
Figure 4. Examples of LIME representations with (a) indicating the original image and (b) the five selected superpixels for the explanation.
Figure 4. Examples of LIME representations with (a) indicating the original image and (b) the five selected superpixels for the explanation.
Entropy 26 00034 g004
Figure 5. Illustration of the ensemble of features by aggregation.
Figure 5. Illustration of the ensemble of features by aggregation.
Entropy 26 00034 g005
Figure 6. Illustration of the k-fold cross-validation strategy applied in this step.
Figure 6. Illustration of the k-fold cross-validation strategy applied in this step.
Entropy 26 00034 g006
Figure 7. Proportion of features from each CNN in the ensembles, observing the 10 best accuracy values for the CR (a), LG (b), OED (c) and UCSB (d) datasets.
Figure 7. Proportion of features from each CNN in the ensembles, observing the 10 best accuracy values for the CR (a), LG (b), OED (c) and UCSB (d) datasets.
Entropy 26 00034 g007
Figure 8. Ranking of accuracy values (%) provided by different approaches in classifying the CR dataset.
Figure 8. Ranking of accuracy values (%) provided by different approaches in classifying the CR dataset.
Entropy 26 00034 g008
Figure 9. Ranking of accuracy values (%) provided by different approaches in classifying the LG dataset.
Figure 9. Ranking of accuracy values (%) provided by different approaches in classifying the LG dataset.
Entropy 26 00034 g009
Figure 10. Ranking of accuracy values (%) provided by different approaches in classifying the OED dataset.
Figure 10. Ranking of accuracy values (%) provided by different approaches in classifying the OED dataset.
Entropy 26 00034 g010
Figure 11. Ranking of accuracy values (%) provided by different approaches in classifying the UCSB dataset.
Figure 11. Ranking of accuracy values (%) provided by different approaches in classifying the UCSB dataset.
Entropy 26 00034 g011
Table 1. Main details of the H&E-stained histological datasets.
Table 1. Main details of the H&E-stained histological datasets.
DatasetImage TypeNumber of ClassesClassesNumber of SamplesResolution
UCSB [66]Breast cancer2Malignant and benign58 (32/26) 896 × 768
CR [67]Colorectal cancer2Malignant and benign165 (74/91)from 567 × 430 to 775 × 552
LG [68]Liver tissue2Male and female265 (150/115) 417 × 312
OED [69]Oral epithelial dysplasia2Healthy and severe148 (74/74) 450 × 250
Table 2. Details of the pretrained architectures and respective accuracies (data from [71]) explored to define the deep-learned features and xAI representations.
Table 2. Details of the pretrained architectures and respective accuracies (data from [71]) explored to define the deep-learned features and xAI representations.
ArchitectureParameterLayersAccuracy (ImageNet)
DenseNet-121 [58] 8 × 10 6 12191.97%
EfficientNet-b2 [65] 9.1 × 10 6 32495.31%
Inception-V3 [59] 2.7 × 10 7 4893.45%
ResNet-50 [60] 2.6 × 10 7 5092.86%
VGG-19 [61] 1.4 × 10 8 1990.87%
Table 3. Feature vectors obtained via different techniques.
Table 3. Feature vectors obtained via different techniques.
OriginCompositionNumber of Features
HandcraftedFractals (F)116
Deep learnedDenseNet-121 (D)1024
EfficientNet-b2 (E)1408
Inception-V3 (I)2048
ResNet-50 (R)2048
VGG-19 (V)4096
xAIGrad-CAM DenseNet-121 ( D C A M )116
Grad-CAM EfficientNet-b2 ( E C A M )116
Grad-CAM Inception-V3 ( I C A M )116
Grad-CAM ResNet-50 ( R C A M )116
Grad-CAM VGG-19 ( V C A M )116
LIME DenseNet-121 ( D L I M E )116
LIME EfficientNet-b2 ( E L I M E )116
LIME Inception-V3 ( I L I M E )116
LIME ResNet-50 ( R L I M E )116
LIME VGG-19 ( V L I M E )116
Table 4. Ensembles of handcrafted and deep-learned features.
Table 4. Ensembles of handcrafted and deep-learned features.
CompositionNumber of Features
F + D1140
F + E1524
F + I2164
F + R2164
F + V4212
F + D + E + I + R + V10,740
Table 5. Ensembles of deep-learned features.
Table 5. Ensembles of deep-learned features.
CompositionNumber of Features
D + E2432
D + I3072
D + R3072
D + V5120
E + I3456
E + R3456
E + V5504
I + R4096
I + V6144
R + V6144
D + E + I + R + V10,624
Table 6. Ensembles of xAI features.
Table 6. Ensembles of xAI features.
CompositionNumber of Features
D C A M + E C A M 232
D C A M + I C A M 232
D C A M + R C A M 232
D C A M + V C A M 232
E C A M + I C A M 232
E C A M + R C A M 232
E C A M + V C A M 232
I C A M + R C A M 232
I C A M + V C A M 232
R C A M + V C A M 232
D C A M + E C A M + I C A M + R C A M + V C A M 580
D L I M E + E L I M E 232
D L I M E + I L I M E 232
D L I M E + R L I M E 232
D L I M E + V L I M E 232
E L I M E + I L I M E 232
E L I M E + R L I M E 232
E L I M E + V L I M E 232
I L I M E + R L I M E 232
I L I M E + V L I M E 232
R L I M E + V L I M E 232
D L I M E + E L I M E + I L I M E + R L I M E + V L I M E 580
Table 7. Average accuracy values (%) computed from the 10 best results in each combination and for each type of H&E image.
Table 7. Average accuracy values (%) computed from the 10 best results in each combination and for each type of H&E image.
CRLGOEDUCSB
Handcrafted84.48% ± 2.0590.42% ± 5.2387.03% ± 1.7872.41% ± 1.89
Deep learned99.27% ± 0.7698.49% ± 0.6596.28% ± 0.5591.38% ± 1.54
xAI83.82% ± 0.9086.98% ± 2.5480.88% ± 0.4378.28% ± 1.14
Ensemble of handcrafted
and deep learned
99.76% ± 0.3099.02% ± 0.5196.49% ± 0.6690.52% ± 0.86
Ensemble of deep learned100%99.66% ± 0.1197.23% ± 0.4792.93% ± 0.93
Ensemble of xAI86.18% ± 1.4589.62% ± 0.4278.85% ± 0.6178.45% ± 1.16
Table 8. Ranking of the best associations with their respective pairwise p-values, according to average accuracy values and the Friedman test.
Table 8. Ranking of the best associations with their respective pairwise p-values, according to average accuracy values and the Friedman test.
p-ValueEnsemble of
Deep Learned
Ensemble of Handcrafted
and Deep Learned
Deep LearnedHandcraftedEnsemble of xAIxAIAverage Ranking
Ensemble of deep learned-0.39550.23960.01910.01910.00661
Ensemble of handcrafted
and deep learned
0.3955-0.73130.10060.10060.03802.25
Deep learned0.23960.7313-0.18190.18190.07352.75
Handcrafted0.01910.10060.1819-1.00000.60734.75
Ensemble of xAI0.01910.10060.18191.0000-0.60734.75
xAI0.00660.03800.07350.60730.6073-5.5
Table 9. Top 10 results for the classification of the CR dataset with feature vectors composed of ensembles of deep-learned features.
Table 9. Top 10 results for the classification of the CR dataset with feature vectors composed of ensembles of deep-learned features.
Feature VectorSizeAccuracy (%)F-1 Score
D + E101001.000
E + V101001.000
D + E + I + R + V101001.000
D + E151001.000
D + I151001.000
E + V151001.000
D + E + I + R + V151001.000
D + E201001.000
D + I201001.000
D + V201001.000
Table 10. LG dataset: top 10 results exploring the ensembles of deep-learned features.
Table 10. LG dataset: top 10 results exploring the ensembles of deep-learned features.
Feature VectorSizeAccuracy (%)F-1 Score
D + R251001.000
D + E1099.620.996
D + E1599.620.996
D + I1599.620.996
D + R1599.620.996
D + E2099.620.996
D + I2099.620.996
D + R2099.620.996
D + E2599.620.996
D + I2599.620.996
Table 11. OED dataset: top 10 results via the ensembles of deep-learned features.
Table 11. OED dataset: top 10 results via the ensembles of deep-learned features.
Feature VectorSizeAccuracy (%)F1-Score
I + V2097.970.980
E + I2597.970.980
D + E2097.300.973
D + I2097.300.973
I + R2097.300.973
D + R2597.300.973
I + V2597.300.973
D + I1096.620.966
D + E1596.620.966
I + V1596.620.966
Table 12. UCSB dataset: Top 10 results exploring the ensembles of deep-learned features.
Table 12. UCSB dataset: Top 10 results exploring the ensembles of deep-learned features.
Feature VectorSizeAccuracy (%)F1-Score
D + E2594.830.948
D + E1593.100.931
E + R1593.100.931
I + R1593.100.931
D + E2093.100.931
D + I2593.100.931
E + I2593.100.931
I + R2593.100.931
E + R1091.380.914
D + R1591.380.914
Table 13. Classification of colorectal samples: accuracy values (%) provided by different approaches.
Table 13. Classification of colorectal samples: accuracy values (%) provided by different approaches.
MethodApproachAccuracy (%)
ProposedDenseNet-121 and EfficientNet-b2 (ensemble of deep-learned features)100%
Roberto et al. [12]ResNet-50, fractal dimension,
lacunarity and percolation (ensemble of handcrafted and deep-learned features)
99.39%
Dabass et al. [101]31-layer CNN (deep learning)96.97%
de Oliveira et al. [19]ResNet50 (activation_48_relu layer), ReliefF and 35
deep-learned features
98.00%
Tavolara et al. [97]GAN and U-Net (deep learning)94.02%
Sena et al. [102]12-layer CNN (deep learning)93.28%
Segato dos Santos et al. [81]Sample entropy and fuzzy logic (handcrafted)91.39%
Roberto et al. [33]Percolation (handcrafted)90.90%
Bentaieb and Hamarneh [103]U-Net and AlexNet (deep learning)87.50%
Zhang et al. [99]ResNet deep-tuning (DL)86.67%
Awan et al. [104]Color normalization, U-Net and GoogLeNet (deep learning)85.00%
Table 14. Classification of liver tissue: accuracy values (%) via different approaches.
Table 14. Classification of liver tissue: accuracy values (%) via different approaches.
MethodApproachAccuracy (%)
ProposedDenseNet-121 and ResNet-50 (ensemble of deep-learned features)100%
Di Ruberto et al. [105]Statistical analysis and texture features (handcrafted)100%
Nanni et al. [13]6 CNNs and handcrafted
features (ensemble of handcrafted and deep-learned features)
100%
Roberto et al. [12]ResNet-50, fractal dimension, lacunarity
and percolation (ensemble of handcrafted and deep-learned features)
99.62%
de Oliveira et al. [19]ResNet50 (activation_48_relu layer), ReliefF and 5
deep-learned features
99.32%
Andrearczyk and Whelan [106]Texture CNN (deep learning)99.10%
Watanabe et al. [107]GIST descriptor, PCA and LDA (handcrafted)93.70%
Table 15. Classification of oral dysplasia: accuracy rates (%) provided by different methods.
Table 15. Classification of oral dysplasia: accuracy rates (%) provided by different methods.
MethodApproachAccuracy (%)
ProposedInception-V3 and VGG-19 (ensemble of deep-learned features)97.97%
Adel et al. [108]SIFT, SURF, ORB (handcrafted)92.80%
Azarmehr et al. [100]Neural architecture search and handcrafted descriptors
(morphological and nonmorphological)
95.20%
Silva et al. [69]Morphological and nonmorphological features (handcrafted)92.40%
Maia et al. [57]Densenet12191.91%
Krishnan et al. [109]Fractal dimension, wavelet, Brownian movement
and Gabor filters (handcrafted)
88.38%
Table 16. Classification of breast cancer: accuracy values (%) trough different strategies.
Table 16. Classification of breast cancer: accuracy values (%) trough different strategies.
MethodApproachAccuracy (%)
Li et al. [110]RefineNet and Atrous DenseNet (deep learning)97.63%
Yu et al. [98]CNN, LBP, SURF, GLCM and other
handcrafted features (ensemble handcrafted and deep-learned features)
96.67%
ProposedDenseNet-121 and EfficientNet-b2 (ensemble of deep-learned features)94.83%
Feng et al. [111]Stacked denoising autoencoder (deep learning)94.41%
Kausar et al. [55]Color normalization, Haar wavelet decomposition and
16-layer CNN (deep learning)
91.00%
Roberto et al. [12]ResNet-50, fractal dimension, lacunarity and
percolation (ensemble of handcrafted and deep-learned features)
89.66%
Roberto et al. [33]Percolation (handcrafted)86.20%
Papastergiou et al. [112]Spacial decomposition and tensors (deep learning)84.67%
Araújo et al. [113]Color normalization, 13-layer CNN and SVM (deep learning)83.30%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Longo, L.H.d.C.; Roberto, G.F.; Tosta, T.A.A.; de Faria, P.R.; Loyola, A.M.; Cardoso, S.V.; Silva, A.B.; do Nascimento, M.Z.; Neves, L.A. Classification of Multiple H&E Images via an Ensemble Computational Scheme. Entropy 2024, 26, 34. https://doi.org/10.3390/e26010034

AMA Style

Longo LHdC, Roberto GF, Tosta TAA, de Faria PR, Loyola AM, Cardoso SV, Silva AB, do Nascimento MZ, Neves LA. Classification of Multiple H&E Images via an Ensemble Computational Scheme. Entropy. 2024; 26(1):34. https://doi.org/10.3390/e26010034

Chicago/Turabian Style

Longo, Leonardo H. da Costa, Guilherme F. Roberto, Thaína A. A. Tosta, Paulo R. de Faria, Adriano M. Loyola, Sérgio V. Cardoso, Adriano B. Silva, Marcelo Z. do Nascimento, and Leandro A. Neves. 2024. "Classification of Multiple H&E Images via an Ensemble Computational Scheme" Entropy 26, no. 1: 34. https://doi.org/10.3390/e26010034

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop