Next Article in Journal
The Impact of Total Tumor Diameter on Lymph Node Metastasis and Tumor Recurrence in Papillary Thyroid Carcinomas
Next Article in Special Issue
Advancements in Artificial Intelligence for Medical Computer-Aided Diagnosis
Previous Article in Journal
Confounders of Ultrasound Attenuation Imaging in a Linear Probe Using the Canon Aplio i800 System: A Phantom Study
Previous Article in Special Issue
Machine Learning-Based Predictive Models for Detection of Cardiovascular Diseases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Dimensional Feature Selection for Automatic Classification of Coronary Stenosis Using an Evolutionary Algorithm

by
Miguel-Angel Gil-Rios
1,
Ivan Cruz-Aceves
2,*,
Arturo Hernandez-Aguirre
3,
Ernesto Moya-Albor
4,
Jorge Brieva
4,
Martha-Alicia Hernandez-Gonzalez
5 and
Sergio-Eduardo Solorio-Meza
6
1
Tecnologías de Información, Universidad Tecnológica de León, Blvd. Universidad Tecnológica 225, Col. San Carlos, León 37670, Mexico
2
CONACYT, Centro de Investigación en Matemáticas (CIMAT), A.C., Jalisco S/N, Col. Valenciana, Guanajuato 36000, Mexico
3
Departamento de Computación, Centro de Investigación en Matemáticas (CIMAT), A.C., Jalisco S/N, Col. Valenciana, Guanajuato 36000, Mexico
4
Facultad de Ingeniería, Universidad Panamericana, Augusto Rodin 498, Ciudad de México 03920, Mexico
5
Unidad Médica de Alta Especialidad (UMAE), Hospital de Especialidades No. 1. Centro Médico Nacional del Bajio, IMSS, Blvd. Adolfo López Mateos esquina Paseo de los Insurgentes S/N, Col. Los Paraisos, León 37320, Mexico
6
División Ciencias de la Salud, Universidad Tecnológica de México, Campus León, Blvd. Juan Alonso de Torres 1041, Col. San José del Consuelo, León 37200, Mexico
*
Author to whom correspondence should be addressed.
Diagnostics 2024, 14(3), 268; https://doi.org/10.3390/diagnostics14030268
Submission received: 18 December 2023 / Revised: 11 January 2024 / Accepted: 23 January 2024 / Published: 26 January 2024
(This article belongs to the Special Issue Artificial Intelligence Advances for Medical Computer-Aided Diagnosis)

Abstract

:
In this paper, a novel strategy to perform high-dimensional feature selection using an evolutionary algorithm for the automatic classification of coronary stenosis is introduced. The method involves a feature extraction stage to form a bank of 473 features considering different types such as intensity, texture and shape. The feature selection task is carried out on a high-dimensional feature bank, where the search space is denoted by O ( 2 n ) and n = 473 . The proposed evolutionary search strategy was compared in terms of the Jaccard coefficient and accuracy classification with different state-of-the-art methods. The highest feature selection rate, along with the best classification performance, was obtained with a subset of four features, representing a 99 % discrimination rate. In the last stage, the feature subset was used as input to train a support vector machine using an independent testing set. The classification of coronary stenosis cases involves a binary classification type by considering positive and negative classes. The highest classification performance was obtained with the four-feature subset in terms of accuracy ( 0.86 ) and Jaccard coefficient ( 0.75 ) metrics. In addition, a second dataset containing 2788 instances was formed from a public image database, obtaining an accuracy of 0.89 and a Jaccard Coefficient of 0.80 . Finally, based on the performance achieved with the four-feature subset, they can be suitable for use in a clinical decision support system.

1. Introduction

Coronary artery disease (CAD) stands as a leading cause of mortality in the majority of developed countries [1]. According to the British Heart Foundation (BHF) [2], coronary heart disease was the main cause of death in the year 2021 around the world. In Figure 1, a comparative chart of different death causes with data extracted from the BHF is illustrated.
According to Figure 1, coronary heart disease presents the highest rate of cases, with 9.2 million registered cases around the world.
In coronary artery disease, atherosclerosis leads to the development of coronary stenosis at various locations [3]. Nowadays, X-ray coronary angiograms are the gold standard for the detection of coronary stenosis in clinical practice. Consequently, a cardiology specialist must exhaustively examine the entire angiogram, and according to their expertise, all regions in which coronary stenosis cases can occur are labeled by hand. Figure 2 presents X-ray coronary angiogram samples with stenosis regions labeled by a specialist in cardiology.
The coronary stenosis problem has been studied in digital image processing, in which several challenging issues must be addressed, such as the presence of noise and weak contrast. In the literature, the method proposed by Saad [4] requires a previous vessel segmentation of a coronary angiogram to identify atherosclerosis using a vessel-width variation measure. A coronary stenosis measure grading method was proposed by Kishore and Jayanthi [5], using the pixel intensities of a previously enhanced image. Alternative methods, such as Brieva et al. [6], used a Hessian-based approach to extract different texture and shape vessel features to classify positive and negative stenosis cases.
Moreover, machine learning-based techniques have been used to address the stenosis classification problem. A naive Bayes classifier was used by Taki et al. [7] to classify calcified and non-calcified coronary artery plaques. The method of Cruz–Aceves et al. [8] uses a 3 D intensity feature vector computed from the histogram of an image to detect a specific type of stenosis using a Bayesian-based classifier. Giannoglou et al. [9] proposed a fuzzy criterion for a feature selection process in atherosclerotic plaques. Chen et al. [10] use a 6 D vector of shape features for the detection of coronary artery disease.
The main disadvantage of previous methods is the use of a fixed threshold measurement in order to classify coronary stenosis cases. Since image datasets of coronary angiograms have different size, contrast, and noise levels, the artery feature values present considerable variations. In consequence, the obtained results are highly dependent on the applied vessel enhancement method. In the proposed method, a number of spatial and frequency domain filters have been adopted for vessel enhancement in order to capture relevant information from different domains.
The method proposed by Antczak and Liberadzki [11], generates synthetic coronary stenosis and non-stenosis patches aiming to enhance the performance of a convolutional neural network (CNN) from scratch. Data augmentation techniques [12] are also used to generate a large number of instances that are used in the training and testing steps of a CNN. In addition, using the explainable AI concept, it is difficult to identify what features are really useful for correct classification and what they represent, which is the main drawback of the CNN [13]. Convolutional neural networks are still considered “black box” systems that do not offer any insight or explanations on how the decision is obtained, limiting the clinical applicability due to their lack of transparency and interpretability. Since the proposed method works with identifiable features, it is possible to know their relevance in the coronary stenosis classification by performance employing statistical analysis.
In the present paper, a novel strategy to perform high-dimensional feature selection using an evolutionary algorithm for the classification of coronary stenosis is proposed. The method involves a binary classification in order to identify positive and negative coronary stenosis cases. An automatic feature selection step is driven by a hybrid-evolutionary algorithm over a high-dimensional bank of features. The accuracy metric was used as a fitness function, while the discrimination rate was also maximized. Shape feature values such as vessel length, bifurcation points and others are dependent on a vessel enhancement technique. Therefore, different enhancement methods were applied to the original images. This strategy allows the formation of a feature bank of 473 different features. To obtain an optimal subset of features, an automatic feature selection stage is performed using a K-nearest neighbor classifier (KNN) to select relevant features. The classification of positive and negative stenosis cases is performed adequately in terms of the accuracy and Jaccard coefficients, which are useful to evaluate the rate of true-positive coronary stenosis cases and avoid the rate of true-negative cases. Since the problem involves a high-dimensional search space, which can be expressed as O ( 2 473 ) , a hybrid-evolutionary algorithm is appropriate for addressing the feature selection optimization problem. An evolutionary algorithm is a high-dimensional optimization technique for working in discrete and continuous domains using the Darwinian theory about the evolution of biological organisms throughout various generations. The hybrid-evolutionary strategy achieved a discrimination rate of 0.99 , obtaining a subset of four features from the extracted bank of 473. In the experiments, a database containing 608 images was used for training of the proposed method (508 for training and 100 for testing). In addition, a second public domain database [11], of 2788 coronary patch images was used for testing.
The rest of this paper is structured as follows. The background methods of vessel enhancement methods, along with intensity, texture, and shape features, are presented in Section 2. In Section 3, the bank of 473 features, the hybrid-evolutionary algorithm, and the performance metrics of the proposed method are introduced. Experimental results are described in Section 4, and conclusions are presented in Section 5.

2. Materials and Methods

2.1. Experiment Materials

For the experiments, two distinct banks of images were used. The first image database was formed from a bank of images provided by the Mexican Institute for Social Healthcare and authorized by an institutional review board only for research purposes under the reference R-2019-1001-078. The bank contains 180 digital images of coronary angiograms, which are 512 × 512 pixels and have a grayscale color scheme. It also included the corresponding ground-truth images with the coronary stenosis cases labeled by a cardiology specialist. From the provided image bank, 304 patches of size 64 × 64 pixels containing coronary stenosis cases were extracted. In addition, 304 additional patches with non-stenosis cases were extracted in order to form a balanced database. A total set of 608 image patches were extracted to form the first database. The size of the patches corresponds to the area enclosing the stenosis cases and was labeled and validated by the specialist.
In order to assess the obtained results, a second image database was formed from the Antczak [11] image database, which is in the public domain. Each image corresponds to a coronary patch of size 32 × 32 pixels in a grayscale color scheme. The original bank is formed by 122 natural coronary patches containing a stenosis case. In addition, it contains 1394 natural coronary patches with non-coronary stenosis cases. Since the proportion of positive and negative stenosis cases is unbalanced, 1272 additional synthetic patches with a coronary stenosis case were added, which were taken from the same image bank. Consequently, the second image database consists of 2788 images with a balance of positive and negative coronary stenosis cases.

2.2. Feature Extraction

The description and measurement of objects of interest, properties of an image, or a specific region is commonly known as feature extraction [14]. It is possible to extract distinct feature types, as described in the literature [15,16]. According to their nature, features can be classified into texture, intensity, and shape.

2.2.1. Intensity Features

Intensity features are relevant in digital image processing because they are related to the corresponding value for each pixel in the image. Here, the five minimum, maximum, median, average, and standard deviation statistical measures of the pixel intensity have been extracted.

2.2.2. Texture Features

Texture features are relevant in different cardiovascular problems [17,18]. One of the most used approaches for the extraction of texture-related features over images is the gray-level co-occurrence matrix (GLCM) [19]. The GLCM computes the frequency of variation between the intensity levels of a pixel. It is expressed as a matrix whose rows and columns correspond to the entire image’s pixel intensities. The frequencies of intensity variations are computed in a specific spatial relationship denoted by ( Δ x , Δ y ) between two different pixels with intensity levels i and j as follows:
C Δ x , Δ y ( i , j ) = x = 1 n y = 1 m 1 , if I ( x , y ) = i and I ( x + Δ x , y + Δ y ) = j 0 , otherwise ,
where C Δ x , Δ y ( i , j ) is the frequency in which two pixels with intensities i and j at a specific offset ( Δ x , Δ y ) occur, and n and m represent the height and width of the image.
Additionally, it is also possible to extract texture features from the output of alternative representation methods, such as the Radon transform [20]. The Radon transform is the projection of the image intensity along with a radial line oriented at some specific angle, which can be computed as follows:
R ( ρ , θ ) = f ( x , y ) δ ( ρ x cos θ y sin θ ) d x d y ,
where R ( ρ , θ ) is the Radon Transform of a function f ( x , y ) at an angle θ , δ ( r ) is the Dirac delta function, and δ ( ρ x cos θ y sin θ ) forces the integration of f ( x , y ) along the line ρ x cos θ y sin θ = 0 .

2.2.3. Shape Features

Shape-based features enable the extraction of quantifiable information concerning various aspects related to the artery shape; for example, a segment length, the tortuosity level present in a determined arterial section, or the number of bifurcations present on it and the vessel width. Nevertheless, to extract shape-based features, it is necessary to perform a previous vessel enhancement process over the original image in which useless information such as noise and background are identified. Consequently, the shape measure values are highly dependent on the applied filtering method. The use of the Hessian matrix and eigenvalues methodology [21] has proved to be adequate for vessel enhancement. However, an automatic thresholding strategy such as the Otsu method [22] is necessary to separate vessel and non-vessel pixels.
The vessel skeleton has been used to extract vessel-shape-related information. To obtain the corresponding vessel skeleton from a binary segmented image, the medial-axis transform method has been commonly applied [23].

2.3. Vessel Enhancement Methods

Vessel enhancement methods are useful for discriminating irrelevant non-vessel information on coronary angiograms. In the literature, spatial and frequency domain filters have been applied to the vessel enhancement problem, achieving suitable results. In the experiments, eight state-of-the-art vessel enhancement methods have been adopted, which are mentioned below.
  • Spatial domain filters
    (a)
    Hessian-based methods.
    • Vesselness measure [24]. The Frangi method computes a vesselness measure using the eigenvalues of a Hessian matrix.
    • Hessian matrix and clustering [25]. The Hessian matrix is also used by Salem et al. for vessel enhancement by computing the largest eigenvalue and vessel orientation over all scales.
    (b)
    Morphological top-hat filter [26]. In mathematical morphology, the top-hat filter is useful to enhance images with non-uniform illumination. Because of this property, the top-hat operator has been used to enhance vessel-like structures [27,28].
    (c)
    Multi-scale line detection [29]. An alternative approach that has been used for artery enhancement is the linear matched filter. This method works under the assumption that blood vessels can be modeled by linear segments that share the same orientation and length.
    (d)
    Gaussian matched filter
    • Single-scale Gaussian filter (GMF). In this approach, a gray-scale template is formed from a Gaussian distribution, which is convolved with the input image.
    • Multi-scale Gaussian filter. The main limitation of GMF is the use of a fixed vessel diameter represented by the σ parameter in which non-corresponding vessel diameters will be distinguished. In order to overcome this disadvantage, a multi-scale Gaussian matched filter was proposed by Cruz–Aceves et al. [30] considering different vessel width scales.
  • Frequency domain filters
    (a)
    Gabor filter
    • Single-scale Gabor filter [31]. The Gabor filter is a Gaussian curve modulated by a sinusoidal function, which is useful for the detection of directional features. In addition, Rangayyan et al. [32] simplified the matching template equation so it is governed by only two parameters.
    • Multi-scale Gabor filter [33]. Similar to the GMF, the use of a fixed vessel diameter represented by the τ parameter will only detect the main artery tree and, as a consequence, discriminate vessels with diameters lower than τ . In order to overcome this disadvantage, Rangayyan et al. proposed a multi-scale Gabor filter for retinal vessels.

2.4. Metaheuristics

The classification techniques learn and predict by classifying instances defined by their features. Therefore, the classification accuracy performance is highly dependent on the used feature set since not all of the used features could be relevant for the classification process. In this context, a feature selection task is necessary after the feature extraction stage is performed. However, the feature selection task is turned into a high-dimensional complexity problem when the number of involved features is elevated because the number of different combinations that are required to find the most suitable feature subset is denoted by 2 n , where n is the number of involved features. Consequently, the use of high-dimensional optimization algorithms is appropriate to address the feature selection problem.

2.4.1. Simulated Annealing

Simulated annealing (SA) is a metaheuristic that was abstracted from an industrial process. In the annealing process, the material is exposed to a certain high temperature, and after, a controlled cooling process is performed. Simultaneously, care must be taken in order to preserve certain molecular alignments in the material to ensure their quality. This process was adapted as a computational search technique by Kirpatrick et al. [34] to solve combinatorial and continuous optimization problems. The algorithm is governed by the T m i n , T m a x and T s t e p parameters, which refer to the minimum, maximum, and changing-step temperatures, respectively. At each iteration step, a new solution is generated based on a computed probability that involves the current temperature and the decreasing parameter Δ E , which represents the objective function response. The probability is calculated from the Boltzmann distribution as follows:
P ( Δ E , T ) = f ( s ) f ( s ) T ,
where P ( Δ E , T ) is the probability computed from the Boltzmann distribution, and f ( s ) and f ( s ) denote the objective function value obtained with the current and the previous SA solution, respectively.

2.4.2. Boltzmann Univariate Marginal Distribution Algorithm (BUMDA)

BUMDA [35] is a population-based method that uses the estimation of distribution to generate new individuals. The main idea of BUMDA is the use of a distribution probability computed from the best solutions of the current generation in order to generate the new one [36]. The Boltzmann probability distribution used by BUMDA is calculated as follows:
μ = j W ( X j ) x j , where W ( X j ) = g ( X j ) X j g ( X j ) ,
ν = j W ( X j ) ( X j μ ) 2 , where W ( X j ) = g ( X j ) X j g ( X j ) + 1 ,
where μ is the objective function average, ν is the objective function variance that was obtained from the population. g ( X j ) corresponds to the value of the objective function obtained by the individual j t h , which is an individual of the population X. Consequently, in order to generate the next population, a fraction of the current one that contains the best individuals is used to produce the new generation, as follows ( n p o p is the population size):
θ t + 1 = f ( x npop ) if t = 1 , f ( x npop 2 ) if f ( x npop 2 ) > = θ t , f ( x i ) when f ( x i ) > = θ t | i = npop 2 + 1 npop ,

2.5. Machine Learning-Based Classifiers

Classifiers are useful for deciding if a specific instance belongs to part of one class or another. For the coronary stenosis classification problems, they are useful for determining if an image or a region over it corresponds to a positive stenosis case or a negative one.

2.5.1. K-Nearest Neighbor

K-nearest neighbor (KNN) represents a fast classification method that was first proposed by Evelyn Fix and Joseph Hodges in 1952 [37]. Later, in 1967, Thomas Cover and Peter E. Hart expanded the initial proposal by introducing the concept of nearest neighbor [38]. The KNN is governed only by the k parameter, which is a positive integer and indicates the number of associated nearest neighbors that a new instance will have in order to measure its probability of membership to different classes. The KNN inputs are labeled vectors in a multidimensional feature space. In the first stage, the training of the KNN model consists only of the storage of the feature vectors and their corresponding class. In the second stage, the KNN classifies new instances by measuring their frequency among the k nearest instances to determine the label of the new instance. In addition, the “nearest” term is associated with the similarity concept in which measurement is commonly based on a distance metric (commonly, Euclidean distance).

2.5.2. Support Vector Machine

Originally conceived as a linear separator for binary classification in supervised learning, the Support Vector Machine (SVM) faces challenges when dealing with instances characterized by significant data overlaps, making linear separability unattainable [39]. To address this issue, SVM has the capability of projecting instances from their original representation space into higher-dimensional orders, enabling successful classification [40]. To execute these projections, SVM leverages instances situated on both sides of the separation boundary, whether it be a line, plane, or hyperplane. The SVM is then formulated as follows [41]:
f ( x ) = W T ϕ ( X ) + b ,
where W is the weight vector and normal to the hyperplane, ϕ is the projection function or kernel, b is the bias or threshold, and X is the data point to be classified.

3. Proposed Method

The proposed strategy consists of three stages. The first stage corresponds to the feature extraction in order to form a bank of 473 features involving intensity, shape, and texture types. For shape features, 8 vessel enhancement methods from the state-of-the-art were used. In the second stage, a subset of features is selected (feature selection) by using a hybrid-evolutionary algorithm in order to maximize the classification accuracy in training data while minimizing the number of features. In the final stage, the selected subset is tested for the classification of coronary stenosis using an independent test set of angiograms. In Figure 3, the steps of the proposed strategy are illustrated.
The bank of 473 features is formed as follows. Intensity-based features, such as minimum, maximum, median, mean, and standard deviation of the pixel intensities, are present in the original image. In addition, these features were also extracted from the responses of the different enhancement methods. In texture features, the Haralik [42] methodology was applied in order to obtain 14 distinct texture-related features. Consequently, 50 different shape features were extracted, including those used by Welikala [43]. Since shape feature values depend on the previous image enhancement process used, different methods were applied, such as those by Frangi et al. and Salem et al., as well as single and multi-scale Gaussian matched filtering, single and multi-scale Gabor filtering, linear multi-scale, and a multi-scale top-hat operator. All extracted features are described below.
  • Intensity-based Features
The intensity features correspond to the statistical measures of standard deviation, minimum, maximum, average, and median of the pixel intensities. Those features were computed from the original image and from the filter response of the 8 enhancement methods described previously. Consequently, 45 intensity-based features were extracted.
  • Texture Features
From the original image, the set of 14 texture features proposed by Haralik [42] were computed. In addition, this set of features was also applied to the Radon transform response. The total number of texture features is 28.
  • Shape Features
For shape features, 14 of them were extracted as described in the Welikala methodology [43], as follows:
  • The total number of vessel pixels.
  • The total number of vessel segments.
  • Vessel density.
  • Tortuosity.
  • The minimum vessel length.
  • The maximum vessel length.
  • The median vessel length.
  • The mean vessel length.
  • The standard deviation length.
  • The number of bifurcation points.
  • Gray level coefficient of variation.
  • Gradient mean.
  • Gradient coefficient of variation.
In addition, in the study by Gil et al. [23], 5 shape features were extracted considering continuous arterial sections and their corresponding segments delimited by tortuosity as follows:
  • The minimum standard deviation of the segments in length pixels considering all arterial sections. Since each arterial section is composed of continuous segments, it is possible to measure the length of each segment and compute the standard deviation for each section. Therefore, if several arterial sections are present in the image, it is possible to obtain statistical measures over the arterial sections.
  • The maximum standard deviation of segments in length pixels considering all arterial sections.
  • The median standard deviation of the segments in length pixels considering all arterial sections.
  • The average of standard deviations of the segments in length pixels considering all arterial sections.
  • The variance of the standard deviations of the segments in length pixels considering all arterial sections.
In addition, 25 shape features were computed as following:
  • Minimum perimeter. The perimeter of an arterial section is the length of its boundary.
  • Maximum perimeter.
  • Median perimeter.
  • Mean perimeter.
  • Standard deviation of the perimeters.
  • Minimum compactness. It can be computed as follows:
    Compactness = Perimeter 2 Area .
  • Maximum compactness.
  • Median compactness.
  • Mean compactness.
  • Standard deviation of compactness.
  • Minimum circularity ratio. It can be computed as follows:
    Circularity Ratio = 4 · π · Area Perimeter 2 .
    Similar to previous measures, for images containing several arterial sections, it is possible to compute circularity for each section and obtain statistical measurements.
  • Maximum circularity ratio.
  • Median circularity ratio.
  • Mean circularity ratio.
  • Standard deviations of the circularity ratios.
  • Minimum rectangularity. It can be computed as follows:
    Rectangularity = Area of Arterial Region Area of Bounding Rectangle of Arterial Region .
  • Maximum rectangularity.
  • Median rectangularity.
  • Mean rectangularity.
  • Standard deviation of rectangularities.
  • Minimum elongatedness. It can be computed as follows:
    Elongatedness = l w ,
    where l is the arterial section length in pixels, and w represents the vessel width in pixels.
  • Maximum elongatedness.
  • Median elongatedness.
  • Mean elongatedness.
  • Standard deviation of elongatedness.
Finally, 6 shape-density features were also extracted.
  • Minimum vessel pixel density of all arterial sections present in the patch.
  • Maximum vessel pixel density.
  • Median vessel pixel density.
  • Mean vessel pixel density.
  • Standard deviation of the vessel pixel densities.
  • Sum of the vessel pixel densities of all arterial sections.
A set with 50 distinct shape-based features was described previously. Since shape feature values are dependent on the applied enhancement method, they were extracted from the 8 responses corresponding to each applied enhancement method in order to extract a set with 400 shape-related features.
After the feature extraction process is concluded, a numeric feature dataset is generated and partitioned randomly into training and testing instances in a balanced manner. The feature selection task is performed on the training dataset, and it is turned into a search process that is conducted by the hybrid-evolutionary algorithm involving the BUMDA and SA metaheuristics.
Since the total number of extracted features is 473, the identification of an optimal feature subset using an exhaustive search process involves a computational cost of O( 2 473 ), which is highly difficult to perform. By involving a single search evolutionary method, the problem can be solved partially. However, due to the high-dimensional complexity of the problem, it is possible to improve the solution achieved by the evolutionary method at each iteration, applying a refined search. This will lead to the use of a hybrid-evolutionary method in which the main goal of feature selection is to identify an optimal subset that improves classification accuracy and reduces model complexity.
The proposed method for the automatic feature selection task is formed by the BUMDA and the SA strategies. The use of BUMDA is adequate because the use of the Boltzmann distribution produces spread populations in comparison with the UMDA and other related techniques [35], which also decreases the risk of falling into local-optima solutions. In addition, since the SA algorithm is a single-solution search method, it is suitable to improve the best solution produced by the BUMDA at each iteration in a refined search step. The combination of these metaheuristics produces a hybrid-evolutionary method focused on the search for the best suitable feature subset by considering the minimization of its size and the maximization of the accuracy classification performance.
Since the feature selection task involves the evaluation of each feature subset produced by the search techniques, a classification model must be trained. By considering that BUMDA is governed by the population size ( p s ) and the number of generations ( n g ), and the SA is governed by the number of iterations ( n i ) computed from the initial, final, and step temperature, the total number of classification models to train is computed as p s × n g × n i . For instance, if p s = 100 , n g = 1000 , and n i = 1000 , the total number of different classifiers to be evaluated is 1 × 10 8 . Consequently, a fast-convergent classification method such as the K-nearest neighbor is suitable in this stage. Subsequently, when the feature selection stage is completed, a more complex classification technique, such as the SVM, can be used in a testing stage to improve the classification performance.
Figure 4 illustrates the flowchart corresponding to the proposed hybrid- evolutionary method.
To evaluate the feature selection performance of the proposed method, the feature decreasing rate (FDR) metric was used and it can be computed as follows:
FDR = 1 Number of Selected Features Total Number of Features .
The FDR metric is used as a way to measure the feature selection performance. When the number of selected features decreases, the FDR increases, allowing the hybrid-evolutionary method to conduct the search for an optimal feature subset by maximizing the FDR.
The final step of the proposed method is the classification of an independent test set using the subset of selected features. In the experiments, the first database of 608 images was divided into training (508 images) and testing (100 images) sets. Both datasets use a balance of positive and negative cases. Each patch is in a grayscale color scheme and their corresponding size is 64 × 64 pixels. In Figure 5, a subset of sample images, along with its filter response, are presented.
In order to evaluate the obtained results, a dataset was formed from the second image database described in Section 2.1. The dataset was divided into the training set with ≈65% (1828 instances) and testing set (960 instances). In Figure 6, sample patches corresponding to the second dataset, are presented.
To determine the classification performance, the accuracy and the Jaccard Coefficient metrics were used. The accuracy metric can be calculated as follows:
A c c = T P + T N T P + T N + F P + F N ,
where A c c is the classification accuracy, T P is the fraction of positive cases classified correctly, T N is the fraction of negative cases classified correctly, F P is the fraction of negative cases classified as positive, and F N represents the fraction of positive cases classified as negative. Moreover, the Jaccard coefficient ( J C ) is calculated as follows:
J C = T P T P + F P + F N ,
The accuracy and the J C metrics are convenient for measuring binary classification performance. The accuracy metric allows us to know how well the classification of positive and negative stenosis cases is performed. However, in coronary stenosis classification, it is important to know the performance of positive coronary stenosis classification, which is evaluated by the J C metric. These two metrics are the most commonly used in image binary classification problems.

4. Results and Discussion

In this section, the results achieved in each stage of the proposed method are presented and discussed. Moreover, they are compared with other methods from the literature. All the experiments were implemented in MATLAB software version 2018 and executed on a computer with an Intel Core i7 processor and 8 GB of RAM.
In the first stage, a bank of 473 features was extracted. This stage involved the original images as well as the responses of the different vessel enhancement methods applied to them. In Table 1, a summary of the distinct feature types extracted is described.
After the feature extraction step was concluded, an optimal subset of features was obtained by applying a hybrid-evolutionary algorithm. The feature subset was evaluated in terms of the training accuracy and FDR metrics. Furthermore, the obtained results were compared with other search methods from the literature. In Table 2, different search methods are described with their corresponding parameter settings, including a statistical analysis of the FDR for each of them.
By considering the results described in Table 2, the highest FDR was obtained using the hybrid-evolutionary strategy, which uses the smallest feature subset of four shape features, as follows:
  • Mean Intensity extracted from the Frangi method response.
  • Average standard deviation of the segments in length pixels for all arterial sections, extracted from the Frangi filter response.
  • Gradient mean extracted from the linear multi-scale method response.
  • Gradient coefficient of variation extracted from the top-hat method response.
In addition, the training classification accuracy performance was also considered. For almost all techniques, the FDR variation was small, which means that the search behavior was stable. By involving the FDR and the training classification accuracy, a performance evaluation was conducted. Figure 7 illustrates the performance chart of the best trial for each compared strategy considering the FDR and the training accuracy.
The classification performance of each method was measured using the testing set. In Table 3, the performance obtained using the testing dataset, which corresponds to the first database, is described in terms of the accuracy and Jaccard coefficient metrics. In addition to the KNN-based classifier, an SVM was also trained in order to measure the performance of the selected features under different classifiers. It was established that there was a maximum number of 1000 iterations for KNN and the SVM and a cross-validation with k = 10 . The values for those parameters were chosen as a tradeoff between an optimal result and the computational time.
The results described in Table 3 shows that the highest rate, in terms of the accuracy and Jaccard coefficient was obtained using the subset of 4 features selected by the hybrid-evolutionary method along with the KNN and the SVM classifiers. Correspondingly, using only 4 of the 473 features represents a discrimination rate of 0.99 % over the initial bank of features. Moreover, the obtained results are similar to those achieved using the Deep Learning CNN-16C and GLNet architectures, which exhibits the robustness of the proposed strategy.
To evaluate the efficiency of the previously discussed results, an additional dataset formed from the publically available database [11], was used for training and testing a KNN and an SVM with the bank of 4 features described in Table 4.
Table 4 shows that the highest classification rate in terms of the accuracy and Jaccard coefficient metrics was also achieved by the proposed method, along with the SVM-based classifier. This result is relevant because it shows that the SVM classification performance is competitive using the bank of four features. Similarly, results obtained using the KNN-based classifier were close to those achieved by the SVM, also using the 4 D feature vector. Consequently, the accuracy and Jaccard coefficient rates show how competitive the feature subset was in classifying stenosis cases. An analysis of the frequency for each of the selected features is described in Figure 8.
According to Figure 8, almost all selected features by the hybrid-evolutionary method have a presence in the other contrasted search methods, which provides evidence of their relevance in the classification process.
Based on Figure 9, it is relevant that two of the four features that were selected by the hybrid-evolutionary strategy are in the group of features with the highest selection probabilities considering all trials of all compared techniques. This finding is important because it statistically validates the relevance of the final selected 4 D feature vector for the automatic classification of coronary stenosis cases. In addition, the mean time required for the classification of a single testing instance was ≈0.02 s.
Even when the achieved results were competitive, the weak vessel contrast present in almost all coronary angiograms and the continuous heart movement decreased the response accuracy of vessel enhancement methods, leading to classification errors. Since the human body is composed of a variety of molecular substances, the reaction of the contrast medium to the X-rays is not enough to produce an optimal artery distinction in the majority of the coronary angiograms. However, although technological innovations have been produced to improve the generation of coronary angiograms [47], the elevated cost of medical imaging devices makes it difficult to adopt in the short term.
One of the strengths of the proposed methods is the identification of a subset with relevant features, which allows performing a positive and negative classification of coronary stenosis cases with a high accuracy rate. However, the conformation of the feature bank can be a limitation in the proposed method and also a robust discrete optimization technique for high-dimensional problems.
After the testing stage was finished, the obtained results with the proposed method achieved a high classification performance in terms of the accuracy and the JC metrics using only a 4 D feature vector. For instance, using the proposed method with the first testing set, accuracy and JC rates of 0.87 and 0.76 were achieved, respectively. It is relevant that the obtained accuracy was slightly higher than that achieved by the deep learning model proposed by Antczak et al., which was 0.86 . In addition, the hybrid-evolutionary strategy results also overcome to those obtained using single search metaheuristics. The single BUMDA and SA methods achieved a classification accuracy performance of 0.82 and 0.85 in their best result. Accordingly, the single SA technique was close to the best accuracy result. However, the single SA found a 15 D feature vector with respect to the 4 D feature vector found using the hybrid-evolutionary strategy, which shows how the hybrid approach was relevant to producing an optimal feature vector in terms of size and classification performance. Furthermore, when the second testing dataset was used to measure the classification performance, accuracy and JC rates of 0.89 and 0.80 were achieved, respectively, which surpassed those obtained by all the compared methods, including the results obtained with deep learning techniques ( 0.86 and 0.74 for accuracy and JC metrics) and single metaheuristics.

5. Conclusions

In this paper, a novel strategy consisting of three stages was presented for the automatic classification of positive and negative coronary stenosis cases. In the first step, a bank of 473 features was formed, which represents a computational search complexity of O ( 2 473 ) , which was explored using different search strategies from the literature. Results achieved in this stage were relevant since only 4 of the 473 features (decreasing rate of 0.99 ) were selected using the hybrid-evolutionary method across all of the other methods compared. It was also relevant that three of the four features are of shape type: Frangi filter—average of the standard deviations of the segments length across all arterial sections, the linear multi-scale filter—gradient mean, and the tophat—gradient coefficient of variation. The remaining feature is intensity type and corresponds to the mean intensity, which is also extracted from the response of the Frangi enhancement method. It is relevant that none of the selected features were obtained from the original image directly. In the testing stage, the proposed method achieved the highest classification performance in terms of accuracy and JC metrics. For the first testing set, the highest accuracy and JC rates were 0.87 and 0.86 , respectively. Correspondingly, with the second dataset, the highest performance in terms of accuracy and JC was 0.84 and 0.89 , respectively. Additionally, taking into account the computational time and classification accuracy of the proposed method based on four selected features, it can be a potential method as part of a computer-aided diagnosis system in cardiology. Finally, future work can be conducted in fast convergence of feature selection for high-dimensional spaces in order to improve the computational time without decreasing the classification performance rate.

Author Contributions

Conceptualization, M.-A.G.-R., I.C.-A., M.-A.H.-G. and S.-E.S.-M.; formal analysis, M.-A.H.-G.; investigation, S.-E.S.-M.; methodology, M.-A.G.-R., I.C.-A., E.M.-A. and J.B.; project administration, I.C.-A. and A.H.-A.; software, M.-A.G.-R. and A.H.-A.; validation, M.-A.G.-R., I.C.-A., M.-A.H.-G. and S.-E.S.-M.; visualization, E.M.-A. and J.B.; writing—original draft, M.-A.G.-R. and I.C.-A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by CONACyT under Project IxM-CONACyT No. 3150-3097. Ernesto Moya-Albor and Jorge Brieva would like to thank Facultad de Ingeniería of Universidad Panamericana for supporting this work. The APC was partially funded by Facultad de Ingeniería, Universidad Panamericana.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The first database presented in this article are not readily available because ethic and legal restrictions. The second database is not our authorship. It is publicly available at: https://github.com/KarolAntczak/DeepStenosisDetection.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BUMDABoltzmann univariate marginal distribution algorithm
FDRfeature decreasing rate
GAgenetic algorithm
KNNK-nearest neighbor
NSFnumber of selected features
SAsimulated annealing
SVMsupport vector machine
UMDAunivariate marginal distribution algorithm

References

  1. Duggan, J.P.; Peters, A.S.; Trachiotis, G.D.; Antevil, J.L. Epidemiology of Coronary Artery Disease. Surg. Clin. 2022, 102, 499–516. [Google Scholar] [CrossRef] [PubMed]
  2. British-Heart-Foundation. Global Hearth and Cirsculatory Diseases Factsheet. Available online: https://www.bhf.org.uk/-/media/files/for-professionals/research/heart-statistics/bhf-cvd-statistics-global-factsheet.pdf?rev=f323972183254ca0a1043683a9707a01&hash=5AA21565EEE5D85691D37157B31E4AAA (accessed on 8 January 2024).
  3. Frąk, W.; Wojtasińska, A.; Lisińska, W.; Młynarska, E.; Franczyk, B.; Rysz, J. Pathophysiology of Cardiovascular Diseases: New Insights into Molecular Mechanisms of Atherosclerosis, Arterial Hypertension, and Coronary Artery Disease. Biomedicines 2022, 10, 1938. [Google Scholar] [CrossRef] [PubMed]
  4. Saad, I.A. Segmentation of Coronary Artery Images and Detection of Atherosclerosis. J. Eng. Appl. Sci. 2018, 13, 7381–7387. [Google Scholar] [CrossRef]
  5. Kishore, A.N.; Jayanthi, V. Automatic stenosis grading system for diagnosing coronary artery disease using coronary angiogram. Int. J. Biomed. Eng. Technol. 2019, 31, 260–277. [Google Scholar] [CrossRef]
  6. Brieva, J.; Gálvez, M.; Toumoulin, C. Coronary extraction and stenosis quantification in X-ray angiographic imaging. In Proceedings of the The 26th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, San Francisco, CA, USA, 1–5 September 2004; Volume 26, pp. 1714–1717. [Google Scholar]
  7. Taki, A.; Roodaki, A.; Setahredan, S.K.; Zoroofi, R.A.; Konig, A.; Navab, N. Automatic segmentation of calcified plaques and vessel borders in IVUS images. Int. J. Comput. Assist. Radiol. Surg. 2008, 2008, 347–354. [Google Scholar] [CrossRef]
  8. Cruz-Aceves, I.; Cervantes-Sanchez, F.; Hernandez-Aguirre, A. Automatic Detection of Coronary Artery Stenosis Using Bayesian Classification and Gaussian Filters Based on Differential Evolution. In Hybrid Intelligence for Image Analysis and Understanding; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2017; Chapter 16; pp. 369–390. [Google Scholar] [CrossRef]
  9. Giannoglou, V.G.; Stavrakoudis, D.G.; Theocharis, J.B. IVUS-based characterization of atherosclerotic plaques using feature selection and SVM classification. In Proceedings of the 2012 IEEE 12th International Conference on Bioinformatics & Bioengineering (BIBE), Larnaca, Cyprus, 11–13 November 2012; pp. 715–720. [Google Scholar] [CrossRef]
  10. Chen, X.; Fu, Y.; Lin, J.; Ji, Y.; Fang, Y.; Wu, J. Coronary Artery Disease Detection by Machine Learning with Coronary Bifurcation Features. Appl. Sci. 2020, 10, 7656. [Google Scholar] [CrossRef]
  11. Antczak, K.; Liberadzki, Ł. Stenosis Detection with Deep Convolutional Neural Networks. MATEC Web Conf. 2018, 210, 04001. [Google Scholar] [CrossRef]
  12. Garcea, F.; Serra, A.; Lamberti, F.; Morra, L. Data augmentation for medical imaging: A systematic literature review. Comput. Biol. Med. 2023, 152, 106391. [Google Scholar] [CrossRef]
  13. Bas, H.; Hugo, J.K.; Gilhuijs, K.G.; Viergever, M.A. Explainable artificial intelligence (XAI) in deep learning-based medical image analysis. Med. Image Anal. 2022, 79, 102470. [Google Scholar] [CrossRef]
  14. Li, Y.; Wang, S.; Tian, Q.; Ding, X. A survey of recent advances in visual feature detection. Neurocomputing 2015, 149, 736–751. [Google Scholar] [CrossRef]
  15. Tessmann, M.; Vega-Higuera, F.; Fritz, D.; Scheuering, M.; Greiner, G. Multi-scale feature extraction for learning-based classification of coronary artery stenosis. In Proceedings of the Medical Imaging 2009: Computer-Aided Diagnosis, Lake Buena Vista, FL, USA, 10–12 February 2009; Karssemeijer, N., Giger, M.L., Eds.; International Society for Optics and Photonics, SPIE: Bellingham, WA, USA, 2009; Volume 7260, p. 726002. [Google Scholar] [CrossRef]
  16. Fazlali, H.R.; Karimi, N.; Soroushmehr, S.M.R.; Sinha, S.; Samavi, S.; Nallamothu, B.; Najarian, K. Vessel region detection in coronary X-ray angiograms. In Proceedings of the 2015 IEEE International Conference on Image Processing (ICIP), Quebec City, QC, Canada, 27–30 September 2015; pp. 1493–1497. [Google Scholar] [CrossRef]
  17. Acharya, U.R.; Sree, S.V.; Krishnan, M.M.R.; Molinari, F.; Saba, L.; Ho, S.Y.S.; Ahuja, A.T.; Ho, S.C.; Nicolaides, A.; Suri, J.S. Atherosclerotic Risk Stratification Strategy for Carotid Arteries Using Texture-Based Features. Ultrasound Med. Biol. 2022, 38, 899–915. [Google Scholar] [CrossRef] [PubMed]
  18. Ricciardi, C.; Valente, A.S.; Edmunds, K.; Cantoni, V.; Green, R.; Fiorillo, A.; Picone, I.; Santini, S.; Cesarelli, M. Linear discriminant analysis and principal component analysis to predict coronary artery disease. Health Inform. J. 2020, 26, 2181–2192. [Google Scholar] [CrossRef] [PubMed]
  19. Barburiceanu, S.; Terebes, R.; Meza, S. 3D Texture Feature Extraction and Classification Using GLCM and LBP-Based Descriptors. Appl. Sci. 2021, 11, 2332. [Google Scholar] [CrossRef]
  20. Murphy, L.M. Linear feature detection and enhancement in noisy images via the Radon transform. Pattern Recognit. Lett. 1986, 4, 279–284. [Google Scholar] [CrossRef]
  21. Cruz-Aceves, I.; Oloumi, F.; Rangayyan, R.M.; Aviña-Cervantes, J.G.; Hernandez-Aguirre, A. Automatic segmentation of coronary arteries using Gabor filters and thresholding based on multiobjective optimization. Biomed. Signal Process. Control 2016, 25, 76–85. [Google Scholar] [CrossRef]
  22. Otsu, N. A Threshold Selection Method from Gray-Level Histograms. IEEE Trans. Syst. Man Cybern. 1979, 9, 62–66. Available online: https://cw.fel.cvut.cz/b201/_media/courses/a6m33bio/otsu.pdf (accessed on 8 January 2024). [CrossRef]
  23. Gil-Rios, M.A.; Chalopin, C.; Cruz-Aceves, I.; Lopez-Hernandez, J.M.; Hernandez-Gonzalez, M.A.; Solorio-Meza, S.E. Automatic Classification of Coronary Stenosis Using Feature Selection and a Hybrid Evolutionary Algorithm. Axioms 2023, 12, 462. [Google Scholar] [CrossRef]
  24. Frangi, A.; Nielsen, W.; Vincken, K.; Viergever, M. Multiscale vessel enhancement filtering. In Proceedings of the Medical Image Computing and Computer-Assisted Intervention—MICCAI’98: First International Conference, Cambridge, MA, USA, 11–13 October 1998; Volume 1496, pp. 130–137. [Google Scholar] [CrossRef]
  25. Salem, N.M.; Salem, S.A.; Nandi, A.K. Segmentation of retinal blood vessels based on analysis of the hessian matrix and Clustering Algorithm. In Proceedings of the 2007 15th European Signal Processing Conference, Poznan, Poland, 3–7 September 2007; pp. 428–432. [Google Scholar]
  26. Eiho, S.; Qian, Y. Detection of coronary artery tree using morphological operator. In Proceedings of the Computers in Cardiology 1997, Lund, Sweden, 7–10 September 1997; pp. 525–528. [Google Scholar] [CrossRef]
  27. Qian, Y.; Eiho, S.; Sugimoto, N.; Fujita, M. Automatic extraction of coronary artery tree on coronary angiograms by morphological operators. In Proceedings of the Computers in Cardiology 1998. Vol. 25 (Cat. No.98CH36292), Cleveland, OH, USA, 13–16 September 1998; pp. 765–768. [Google Scholar] [CrossRef]
  28. Kang, W.; Kang, W.; Li, Y.; Wang, Q. The segmentation method of degree-based fusion algorithm for coronary angiograms. In Proceedings of the 2013 2nd International Conference on Measurement, Information and Control, Harbin, China, 16–18 August 2013; Volume 01, pp. 696–699. [Google Scholar] [CrossRef]
  29. Ricci, E.; Perfetti, R. Retinal Blood Vessel Segmentation Using Line Operators and Support Vector Classification. IEEE Trans. Med Imaging 2007, 26, 1357–1365. [Google Scholar] [CrossRef]
  30. Cruz-Aceves, I.; Cervantes-Sanchez, F.; Avila-Garcia, M.S. A Novel Multiscale Gaussian-Matched Filter Using Neural Networks for the Segmentation of X-Ray Coronary Angiograms. J. Healthc. Eng. 2018, 2018, 5812059. [Google Scholar] [CrossRef] [PubMed]
  31. Gabor, D. Theory of communication. Part 1: The analysis of information. J. Inst. Electr.-Eng. Part Iii Radio Commun. Eng. 1946, 93, 429–441. [Google Scholar] [CrossRef]
  32. Rangayyan, R.M.; Oloumi, F.; Oloumi, F.; Eshghzadeh-Zanjani, P.; Ayres, F.J. Detection of Blood Vessels in the Retina Using Gabor Filters. In Proceedings of the 2007 Canadian Conference on Electrical and Computer Engineering, Vancouver, BC, Canada, 22–26 April 2007; pp. 717–720. [Google Scholar] [CrossRef]
  33. Rangayyan, R.M.; Ayres, F.J.; Oloumi, F.; Oloumi, F.; Eshghzadeh-Zanjani, P. Detection of blood vessels in the retina with multiscale Gabor filters. J. Electron. Imaging 2008, 17, 023018. [Google Scholar] [CrossRef]
  34. Kirkpatrick, S.; Gelatt, C.D.; Vecchi, M.P. Optimization by simulated annealing. Science 1983, 220, 671–680. [Google Scholar] [CrossRef]
  35. Valdez-Peña, S.I.; Hernández, A.; Botello, S. A Boltzmann based estimation of distribution algorithm. Inf. Sci. 2013, 236, 126–137. [Google Scholar] [CrossRef]
  36. Dang, D.C.; Lehre, P.K.; Nguyen, P.T.H. Level-Based Analysis of the Univariate Marginal Distribution Algorithm. Algorithmica 2017, 81, 668–702. [Google Scholar] [CrossRef]
  37. Fix, E.; Hodges, J.J.L. Discriminatory Analysis—Nonparametric Discrimination: Small Sample Performance; Technical Report; Air University, USAF School of Aviation Medecine: San Antonio, TX, USA, 1952; pp. 13–37. [Google Scholar]
  38. Cover, T.; Hart, P. Nearest neighbor pattern classification. IEEE Trans. Inf. Theory 1967, 13, 21–27. [Google Scholar] [CrossRef]
  39. Cortes, C.; Vapnik, V. Support-Vector Networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar] [CrossRef]
  40. Noble, W.S. What is a support vector machine? Nat. Biotechnol. 2006, 24, 1565–1567. [Google Scholar] [CrossRef] [PubMed]
  41. Hari, B.N.; Salankar, S.S.; Bora, V.R. MRI brain cancer classification using Support Vector Machine. In Proceedings of the 2014 IEEE Students’ Conference on Electrical, Electronics and Computer Science, Bhopal, India, 1–2 March 2014; pp. 1–6. [Google Scholar] [CrossRef]
  42. Haralick, R.; Shanmugam, K.; Dinstein, I. Textural Features for Image Classification. IEEE Trans. Syst. Man Cybern. 1973, SMC-3, 610–621. [Google Scholar] [CrossRef]
  43. Welikala, R.; Fraz, M.; Dehmeshki, J.; Hoppe, A.; Tah, V.; Mann, S.; Williamson, T.; Barman, S. Genetic algorithm based feature selection combined with dual classification for the automated detection of proliferative diabetic retinopathy. Comput. Med. Imaging Graph. 2015, 43, 64–77. [Google Scholar] [CrossRef] [PubMed]
  44. Szegedy, C.; Liu, W.; Jia, Y.; Sermanet, P.; Reed, S.; Anguelov, D.; Erhan, D.; Vanhoucke, V.; Rabinovich, A. Going Deeper with Convolutions. arXiv 2014, arXiv:1409.4842. [Google Scholar]
  45. Harouni, A.; Karargyris, A.; Negahdar, M.; Beymer, D.; Syeda-Mahmood, T. Universal multi-modal deep network for classification and segmentation of medical images. In Proceedings of the 2018 IEEE 15th International Symposium on Biomedical Imaging (ISBI 2018), Washington, DC, USA, 4–7 April 2018; pp. 872–876. [Google Scholar] [CrossRef]
  46. Gil-Rios, M.A.; Guryev, I.V.; Cruz-Aceves, I.; Avina-Cervantes, J.G.; Hernandez-Gonzalez, M.A.; Solorio-Meza, S.E.; Lopez-Hernandez, J.M. Automatic Feature Selection for Stenosis Detection in X-ray Coronary Angiograms. Mathematics 2021, 9, 2471. [Google Scholar] [CrossRef]
  47. Hajhosseiny, R.; Bustin, A.; Munoz, C.; Rashid, I.; Cruz, G.; Manning, W.J.; Prieto, C.; Botnar, R.M. Coronary Magnetic Resonance Angiography: Technical Innovations Leading Us to the Promised Land? JACC Cardiovasc. Imaging 2020, 13, 2653–2672. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Comparison of five main diseases around the world: coronary heart disease (C.H.D.), stroke, chronic inflammatory pulmonary disease (C.I.P.D.), lower respiratory infections (L.R.I.) and lung cancer.
Figure 1. Comparison of five main diseases around the world: coronary heart disease (C.H.D.), stroke, chronic inflammatory pulmonary disease (C.I.P.D.), lower respiratory infections (L.R.I.) and lung cancer.
Diagnostics 14 00268 g001
Figure 2. Coronary angiograms with their respective coronary stenosis regions labeled by the specialist.
Figure 2. Coronary angiograms with their respective coronary stenosis regions labeled by the specialist.
Diagnostics 14 00268 g002
Figure 3. Proposed feature extraction, selection, and classification methods to classify coronary stenosis.
Figure 3. Proposed feature extraction, selection, and classification methods to classify coronary stenosis.
Diagnostics 14 00268 g003
Figure 4. Flowchart of the proposed hybrid method to perform automatic feature selection.
Figure 4. Flowchart of the proposed hybrid method to perform automatic feature selection.
Diagnostics 14 00268 g004
Figure 5. Patch samples from a database of 608 images. Rows (1) and (2) correspond to positive stenosis cases. Rows (3) and (4) correspond to negative stenosis cases. Column (a) corresponds to the original patch image. Columns (bi) correspond to the different vessel enhancement method responses, as follows: Frangi, Salem, simple-scale Gabor, multi-scale Gabor, multi-scale linear, multi-scale matched filter, single-scale matched filter, and top-hat operator.
Figure 5. Patch samples from a database of 608 images. Rows (1) and (2) correspond to positive stenosis cases. Rows (3) and (4) correspond to negative stenosis cases. Column (a) corresponds to the original patch image. Columns (bi) correspond to the different vessel enhancement method responses, as follows: Frangi, Salem, simple-scale Gabor, multi-scale Gabor, multi-scale linear, multi-scale matched filter, single-scale matched filter, and top-hat operator.
Diagnostics 14 00268 g005
Figure 6. Patch samples from the second database of 2788 images. Rows (1) and (2) correspond to positive stenosis cases. Rows (3) and (4) correspond to negative stenosis cases. Column (a) corresponds to the original patch image. Columns (bi) correspond to the different enhancement method responses as follows: Frangi, Salem, simple-scale Gabor, multi-scale Gabor, multi-scale linear, multi-scale matched filter, single-scale matched filter, and top-hat operator.
Figure 6. Patch samples from the second database of 2788 images. Rows (1) and (2) correspond to positive stenosis cases. Rows (3) and (4) correspond to negative stenosis cases. Column (a) corresponds to the original patch image. Columns (bi) correspond to the different enhancement method responses as follows: Frangi, Salem, simple-scale Gabor, multi-scale Gabor, multi-scale linear, multi-scale matched filter, single-scale matched filter, and top-hat operator.
Diagnostics 14 00268 g006
Figure 7. Performance evaluation chart considering the FDR and the training accuracy for selected metaheuristics: UMDA, BUMDA, GA, SA, and the hybrid-evolutionary method.
Figure 7. Performance evaluation chart considering the FDR and the training accuracy for selected metaheuristics: UMDA, BUMDA, GA, SA, and the hybrid-evolutionary method.
Diagnostics 14 00268 g007
Figure 8. Frequencies for each selected feature in the best solution achieved by the hybrid-evolutionary method contrasted with the frequency for that same feature in the best solution achieved by the other methods. Selected features are as follows: Frangi filter—mean intensity (F1), Frangi filter—mean std. dev. of segments length in all arterial sections (F2), linear multi-scale filter—gradient mean (F3), top-hat—gradient coefficient of variation (F4).
Figure 8. Frequencies for each selected feature in the best solution achieved by the hybrid-evolutionary method contrasted with the frequency for that same feature in the best solution achieved by the other methods. Selected features are as follows: Frangi filter—mean intensity (F1), Frangi filter—mean std. dev. of segments length in all arterial sections (F2), linear multi-scale filter—gradient mean (F3), top-hat—gradient coefficient of variation (F4).
Diagnostics 14 00268 g008
Figure 9. Selection probability for all 473 features based on the frequency in which each of them was selected in all trials considering all contrasted strategies.
Figure 9. Selection probability for all 473 features based on the frequency in which each of them was selected in all trials considering all contrasted strategies.
Diagnostics 14 00268 g009
Table 1. Summary of the different feature types extracted from the original images and their corresponding vessel-enhancement responses.
Table 1. Summary of the different feature types extracted from the original images and their corresponding vessel-enhancement responses.
Feature TypeRequire the Original ImageRequire Vessel-EnhancementQuantity
IntensityYesYes45
TextureYesNo28
ShapeYesYes400
Total Extracted Features:473
Table 2. Parameter settings and statistical analysis of the feature decreasing rate (FDR) for the feature selection performance of distinct search metaheuristics of 30 independent trials.
Table 2. Parameter settings and statistical analysis of the feature decreasing rate (FDR) for the feature selection performance of distinct search metaheuristics of 30 independent trials.
MethodPop. SizeMax. Gens.Min. FDRMax. FDRMedian FDRMean FDRStd. Dev. FDR
UMDA10010000.480.980.920.830.17
BUMDA10010000.570.500.570.530.02
GA10010000.940.910.940.930.01
SA 1-10,0000.870.970.920.920.03
Hybrid Evolutionary 2100100/5000 10.920.990.970.970.03
1 Simulated annealing ( T m a x = 1 , T m i n = 0 , T s t e p = 0.0001 ). 2 BUMDA was configured with 100 maximum number of generations. Subsequently, the inner SA strategy was configured with ( T m a x = 1 , T m i n = 0 , T s t e p = 0.005 ) . In total, the hybrid-evolutionary strategy iterated 5 × 10 5 times.
Table 3. Description of the test results for each compared strategy using the KNN and the SVM classifiers, including the number of selected features (NSF) with their corresponding feature decreasing rate (FDR), the achieved classification accuracy, and Jaccard coefficient (JC), using the first image database containing 100 balanced instances.
Table 3. Description of the test results for each compared strategy using the KNN and the SVM classifiers, including the number of selected features (NSF) with their corresponding feature decreasing rate (FDR), the achieved classification accuracy, and Jaccard coefficient (JC), using the first image database containing 100 balanced instances.
MethodNSFFDRClassifierAccuracyJC
GLNet [44]0.850.76
UNet [45]0.850.75
CNN-16C [11]0.860.76
UMDA [46]100.98KNN0.810.75
SVM0.800.67
BUMDA2050.57KNN0.810.70
SVM0.820.72
GA290.94KNN0.790.65
SVM0.790.66
SA160.97KNN0.810.65
SVM0.850.75
Hybrid-Evolutionary40.99KNN0.870.76
SVM0.860.75
Table 4. Description of the testing results for each compared strategy using the KNN and the SVM classifiers, including the number of selected features (NSF) with their corresponding feature decreasing rate (FDR), the achieved classification accuracy and Jaccard coefficient, using a dataset formed from the Antczak image database, containing 960 balanced instances.
Table 4. Description of the testing results for each compared strategy using the KNN and the SVM classifiers, including the number of selected features (NSF) with their corresponding feature decreasing rate (FDR), the achieved classification accuracy and Jaccard coefficient, using a dataset formed from the Antczak image database, containing 960 balanced instances.
MethodNSFFDRClassifierAccuracyJC
GLNet [44]0.720.63
UNet [45]0.760.72
CNN-16C [11]0.860.74
UMDA [46]100.98KNN0.850.75
SVM0.830.71
BUMDA2050.57KNN0.860.77
SVM0.740.56
GA290.94KNN0.850.75
SVM0.870.78
SA160.97KNN0.850.75
SVM0.850.75
Hybrid Evolutionary40.99KNN0.840.74
SVM0.890.80
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

Gil-Rios, M.-A.; Cruz-Aceves, I.; Hernandez-Aguirre, A.; Moya-Albor, E.; Brieva, J.; Hernandez-Gonzalez, M.-A.; Solorio-Meza, S.-E. High-Dimensional Feature Selection for Automatic Classification of Coronary Stenosis Using an Evolutionary Algorithm. Diagnostics 2024, 14, 268. https://doi.org/10.3390/diagnostics14030268

AMA Style

Gil-Rios M-A, Cruz-Aceves I, Hernandez-Aguirre A, Moya-Albor E, Brieva J, Hernandez-Gonzalez M-A, Solorio-Meza S-E. High-Dimensional Feature Selection for Automatic Classification of Coronary Stenosis Using an Evolutionary Algorithm. Diagnostics. 2024; 14(3):268. https://doi.org/10.3390/diagnostics14030268

Chicago/Turabian Style

Gil-Rios, Miguel-Angel, Ivan Cruz-Aceves, Arturo Hernandez-Aguirre, Ernesto Moya-Albor, Jorge Brieva, Martha-Alicia Hernandez-Gonzalez, and Sergio-Eduardo Solorio-Meza. 2024. "High-Dimensional Feature Selection for Automatic Classification of Coronary Stenosis Using an Evolutionary Algorithm" Diagnostics 14, no. 3: 268. https://doi.org/10.3390/diagnostics14030268

APA Style

Gil-Rios, M. -A., Cruz-Aceves, I., Hernandez-Aguirre, A., Moya-Albor, E., Brieva, J., Hernandez-Gonzalez, M. -A., & Solorio-Meza, S. -E. (2024). High-Dimensional Feature Selection for Automatic Classification of Coronary Stenosis Using an Evolutionary Algorithm. Diagnostics, 14(3), 268. https://doi.org/10.3390/diagnostics14030268

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