Next Article in Journal
Multi Set-Point Explicit Model Predictive Control for Nonlinear Process Systems
Next Article in Special Issue
An Improved Hybrid Aquila Optimizer and Harris Hawks Algorithm for Solving Industrial Engineering Optimization Problems
Previous Article in Journal
Multivariable Robust Regulation of Alkalinities in Continuous Anaerobic Digestion Processes: Experimental Validation
Previous Article in Special Issue
Improved NSGA-III with Second-Order Difference Random Strategy for Dynamic Multi-Objective Optimization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Evolutionary Arithmetic Optimization Algorithm for Multilevel Thresholding Segmentation of COVID-19 CT Images

1
Faculty of Computer Sciences and Informatics, Amman Arab University, Amman 11953, Jordan
2
School of Computer Sciences, Universiti Sains Malaysia, Pulau Pinang 11800, Malaysia
3
Division of Engineering, New York University Abu Dhabi, Saadiyat Island, Abu Dhabi 129188, United Arab Emirates
4
Department of Civil and Urban Engineering, Tandon School of Engineering, New York University, Brooklyn, NY 11201, USA
5
Faculty of Engineering and Information Technology, University of Technology Sydney, Ultimo, NSW 2007, Australia
*
Author to whom correspondence should be addressed.
Processes 2021, 9(7), 1155; https://doi.org/10.3390/pr9071155
Submission received: 5 May 2021 / Revised: 16 June 2021 / Accepted: 18 June 2021 / Published: 2 July 2021
(This article belongs to the Special Issue Evolutionary Process for Engineering Optimization)

Abstract

:
One of the most crucial aspects of image segmentation is multilevel thresholding. However, multilevel thresholding becomes increasingly more computationally complex as the number of thresholds grows. In order to address this defect, this paper proposes a new multilevel thresholding approach based on the Evolutionary Arithmetic Optimization Algorithm (AOA). The arithmetic operators in science were the inspiration for AOA. DAOA is the proposed approach, which employs the Differential Evolution technique to enhance the AOA local research. The proposed algorithm is applied to the multilevel thresholding problem, using Kapur’s measure between class variance functions. The suggested DAOA is used to evaluate images, using eight standard test images from two different groups: nature and CT COVID-19 images. Peak signal-to-noise ratio (PSNR) and structural similarity index test (SSIM) are standard evaluation measures used to determine the accuracy of segmented images. The proposed DAOA method’s efficiency is evaluated and compared to other multilevel thresholding methods. The findings are presented with a number of different threshold values (i.e., 2, 3, 4, 5, and 6). According to the experimental results, the proposed DAOA process is better and produces higher-quality solutions than other comparative approaches. Moreover, it achieved better-segmented images, PSNR, and SSIM values. In addition, the proposed DAOA is ranked the first method in all test cases.

1. Introduction

One of the most often used image segmentation techniques is multilevel thresholding. It is divided into two types: bi-level and multilevel [1,2]. Multilevel thresholding is used to separate complex images, which can generate several thresholds, such as tri-level or quad-level thresholds, which break pixels into several identical parts depending on size. Bi-level thresholding divides the image into two levels, while multilevel thresholding divides the image into two classes [3,4]. When there are only two primary gray levels in an image, bi-level thresholding yields acceptable results; however, when it is expanded to multilevel thresholding, the main drawback is the time-consuming computation [5]. Bi-level thresholding cannot precisely find the optimum threshold, due to the slight variation between the target and the context of a complex image [6,7].
Medical imaging, machine vision, and satellite photography all use image segmentation [8,9,10]. The primary aim of image segmentation is to divide an image into relevant regions for a specific mission. The process of finding and isolating points of interest from the rest of the scene is known as the segmentation of pattern recognition systems [11,12]. Following image segmentation, certain features from objects are removed, and then objects are grouped into specific categories or classes, depending on the extracted features. Segmentation is used in medical applications to detect organs, such as the brain, heart, lungs, and liver, in CT or MR images [13,14]. It is also used to tell the difference between abnormal tissue, such as a tumor, and healthy tissue. Image segmentation techniques, such as image thresholding, edge detection, area expanding, stochastic models, Artificial Neural Network (ANN), and clustering techniques, have all been used, depending on the application [15,16].
Tsallis, Kapur, and Otsu procedures are the most widely used thresholding strategies [17,18]. The Otsu method maximizes the between-class variance function to find optimum thresholds, while the Kapur method maximizes the posterior entropy of the segmented groups. Due to exhaustive search, Tsallis and Otsu’s computational complexity grows exponentially as the number of thresholds increases [19]. Many researchers have worked on image segmentation over the years. Image segmentation has been tackled using a variety of approaches and algorithms [20]. Examples of the used optimization algorithms are the Bat Algorithm (BA) [21], Firefly Algorithm (FA) [22], Genetic Algorithm (GA) [23], Gray Wolf Optimizer (GWO) [24,25], Dragonfly Algorithm (DA) [26], Moth–Flame Optimization Algorithm (MFO) [27], Marine Predators Algorithm (MPA) [28], Arithmetic Optimization Algorithm (AOA) [29], Aquila Optimizer (AO) [30], Krill Herd Optimizer (KHO) [31], Harris Hawks Optimizer (HHO) [32], Red Fox Optimization Algorithm (RFOA) [33], Artificial Bee Colony Algorithm (ABC) [34], and Artificial Ecosystem-based Optimization [35]. Many other optimizers can be found in [36,37].
The paper [38] used Kapur and Otsu’s approaches to adjust the latest Elephant Herding Optimization Algorithm for multilevel thresholding. Its performance was compared to four other swarm intelligence algorithms, using regular benchmark images. The Elephant Herding Optimization Algorithm outperformed and proved more stable than other methods in the literature. Sahlol et al. in [39] introduced an improved hybrid method for COVID-19 images by merging the strengths of convolution neural networks (CNNs) to remove features and the MPA feature selection algorithm to choose the most important features. The proposed method exceeds several CNNs and other well-known methods on COVID-19 images.
The multi-verse optimizer (MVO), based on the multi-verse theorem, is a new algorithm for solving real-world multi-parameter optimization problems. A novel parallel multi-verse optimizer (PMVO) with a coordination approach is proposed in [40]. For each defined iteration, the parallel process is used to randomly split the original solutions into multiple groups and exchange the various groups’ details. This significantly improves individual MVO algorithm cooperation and reduces the shortcomings of the original MVO algorithm, such as premature convergence, search stagnation, and easy trapping into the local optimal search. The PMVO algorithm was compared to methods under the CEC2013 test suite to validate the proposed scheme’s efficiency. The experimental findings show that the PMVO outperforms the other algorithms under consideration. In addition, using minimum cross entropy thresholding, PMVO is used to solve complex multilevel image segmentation problems. In comparison with different related algorithms, the proposed PMVO algorithm seems to achieve better quality image segmentation.
For image segmentation, a modified artificial bee colony optimizer (MABC) is proposed [41], which balances the tradeoff between the search process by using a pool of optimal foraging strategies. MABC’s main goal is to improve artificial bee foraging behaviors by integrating local search with detailed learning, using a multi-dimensional PSO-based equation. With detailed learning, the bees combine global best solution knowledge into the solution quest equation to increase exploration. Simultaneously, local search allows the bees to thoroughly exploit across the promising field, providing a good combination of exploration and exploitation. The proposed algorithm’s feasibility was shown by the experimental findings comparing the MABC to several popular EA and SI algorithms on a series of benchmarks. The experimental findings verify the suggested algorithm’s efficacy.
For solving the image segmentation problem, a novel multilevel thresholding algorithm based on a meta-heuristic Krill Herd Optimization (KHO) algorithm is proposed in [42]. The optimal threshold values are calculated, using the Krill Herd Optimization technique to maximize Kapur’s or Otsu’s objective function. The suggested method reduces the amount of time it takes to calculate the best multilevel thresholds. Various benchmark images are used to illustrate the applicability and numerical performance of the Krill Herd Optimization-based multilevel thresholding. To demonstrate the superior performance of the proposed method, a detailed comparison with other current bio-inspired techniques based on multilevel thresholding techniques, such as Bacterial Foraging (BF), Particle Swarm Optimization (PSO), Genetic Algorithm (GA), and Moth-Flame Optimization (MFO), was performed. The results confirmed that the proposed method achieved better results than other methods.
This paper presents a modified version of the Manta Ray Foraging Optimizer (MRFO) algorithm to deal with global optimization and multilevel image segmentation problems [43]. MRFO is a meta-heuristic technique that simulates the behaviors of manta rays to find food. The performance of the MRFO is improved by using fractional-order (FO) calculus during the exploitation phase. In this experiment, a variant of natural images is used to assess FO-MRFO. According to different performance measures, the FO-MRFO outperformed the compared algorithms in global optimization and image segmentation.
The concept “optimization” refers to the process of identifying the best solutions to a problem while keeping those constraints in mind [44,45]. The used optimization in solving the image segmentation problem is the method of finding the best threshold values for a given image. Swarm intelligence (SI) algorithms are used widely for multilevel thresholding problems to determine the optimal threshold values, using various objective functions to solve the problems of the computational inefficiency of traditional thresholding techniques. The primary motivation behind this paper is to find the optimal threshold values for image segmentation problems. At the same time, to address the weakness of the original AOA, it suffers from the local optimal problem and premature coverage in some cases. In this paper, an improved version of the Arithmetic Optimization Algorithm (AOA) by using Differential Evolution, called DAOA, is proposed. The proposed method uses Differential Evolution to tackle the conventional Arithmetic Optimization Algorithm’s weaknesses, such as being trapped in local optima and fast convergence. Thus, Differential Evolution is used to enhance the performance of the Arithmetic Optimization Algorithm. The proposed DAOA assists by using eight standard test images from different groups: two-color images, two gray images, two normal CT COVID-19 images, and two confirmed COVID-19 CT images. Peak signal-to-noise ratio (PSNR), structural similarity index test (SSIM), and fitness function (Kapur’s) are used to determine the accuracy of segmented images. The proposed DAOA method’s efficiency is evaluated and compared to other multilevel thresholding methods. The findings are presented with a number of different threshold values (i.e., 2, 3, 4, 5, and 6). According to the experimental results, the proposed DAOA process is better and produces higher-quality solutions than other approaches. The encouraging findings suggest that using the DAOA-based thresholding strategy has potential and is helpful.
The rest of this paper is organized as follows. Section 2 presents the procedure of the proposed DAOA method. Section 3 presents the definitions and procedures of the image segmentation problem. The experiments and results are given in Section 4. Finally, in Section 5, the conclusions and potential future work directions are given.

2. The Proposed Method

In this section, we present the conventional Arithmetic Optimization Algorithm (AOA), Differential Evolution (DE), and the proposed Evolutionary Arithmetic Optimization Algorithm (DAOA).

2.1. Arithmetic Optimization Algorithm (AOA)

In this section, we describe the exploration and exploitation phases of the original AOA [29], which is motivated by the main operators in math science (i.e., multiplication (M), division (D), subtraction (S), and addition (A)). The main search methods of the AOA are presented in Figure 1, which are illustrated in the following subsections.
The AOA should choose the search process before beginning its work (i.e., exploration or exploitation). So, in the following search steps, the Math Optimizer Accelerated ( M O A ) function is a coefficient determined by Equation (1).
M O A ( C _ I t e r ) = M i n + C _ I t e r × M a x M i n M _ I t e r
where M O A (C_Iter) means the value at the t th iteration of MOA function, determined by Equation (1). C _ I t e r is the current iteration: [1 ……… M _ I t e r ]. M i n and M a x are the accelerated function values (minimum and maximum), respectively.

2.1.1. Exploration Phase

The exploration operators of AOA are modeled in Equation (2). The exploration phase uses the D or M operators conditioned by r 1 > MOA. The D operator is prepared by r 2 < 0.5, or, otherwise, by the M operator. r 2 is a random number. The position updating process is determined as follows.
x i , j ( C _ I t e r + 1 ) = b e s t ( x j ) ÷ M O P × ( ( U B j L B j ) × μ + L B j ) , r 2 < 0.5 b e s t ( x j ) × M O P × ( ( U B j L B j ) × μ + L B j ) , o t h e r w i s e
where x i (C_Iter + 1 ) is the i th next solution, x i , j (C_Iter) is the j th location of the i th solution, and b e s t ( x j ) is the j th location in the best solution. μ is a control value (0.5) to tune the exploration search.
M O P ( C _ I t e r ) = 1 C _ I t e r 1 / α M _ I t e r 1 / α
where M O P (C_Iter) denotes the coefficient value at the t th iteration. α is a control value (5) to tune the exploration search.

2.1.2. Exploitation Phase

The exploitation searching phase uses the S and A operators conditioned by the MOA function value. Subtraction (S) and addition (A) search strategies are represented in Equation (4).
x i , j ( C _ I t e r + 1 ) = b e s t ( x j ) M O P × ( ( U B j L B j ) × μ + L B j ) , r 3 < 0.5 b e s t ( x j ) + M O P × ( ( U B j L B j ) × μ + L B j ) , o t h e r w i s e
The intuitive and detailed process of AOA is shown in Figure 2.

2.2. Differential Evolution (DE)

In [46], Storn and Price introduced the DE as the first version to solve multiple optimization problems in 1997. DE stands out for its versatility, quick execution time, rapid acceleration pattern, and fast and accurate local operators [47,48]. In DE, the optimization process begins with a random selection of solutions for finding the majority of the points in the search space (initialization phase). The solutions can then be improved, using a series of operators called mutation and crossover. The new solution can be accepted if it has a higher objective value. For the current solution X i , the mathematical model of the mutation operator Z i t can be applied as follows:
Z i , j = X D r a n d 1 + F × ( X D r 2 X D r 3 ) ,
where r 1 , r 2 , and r 3 are random numbers, F is the mutation balancing factor, and F is greater than 0.
For the crossover operator, Equation (6) represents the new solution V i , which is produced using the mutated operator through the crossover Z i . The crossover is considered a mixture process among vectors Z i and X D i .
V i , j = Z i , j i f r a n d C r X D i , j o t h e r w i s e
C r is the crossover probability.
The DE algorithm improves its selected solutions according to the objective function values, where the generated V i , C _ I t e r is replaced with the current one if it obtained a better fitness value, which is as follows.
X D i , j = V i , j i f f ( V i , j ) < f ( X D i , j ) X D i , j o t h e r w i s e

2.3. The Proposed DAOA

In this section, the procedure of the proposed Evolutionary Arithmetic Optimization Algorithm (DAOA) is presented as follows.

2.3.1. Initialization Phase

When using the AOA, the optimization procedure begins with a number of random solutions (X) as designated in matrix (8). The best solution is taken in each iteration as the best-obtained solution.
X = x 1 , 1 x 1 , j x 1 , n 1 x 1 , n x 2 , 1 x 2 , j x 2 , n x N 1 , 1 x N 1 , j x N 1 , n x N , 1 x N , j x N , n 1 x N , n

2.3.2. Phases of the Proposed DAOA

In this section, the main details and procedures of the proposed Evolutionary Arithmetic Optimization Algorithm (DAOA) are given.
The DAOA is introduced mainly to develop the original AOA’s convergence ability, the quality of solutions, and the ability to avoid the local optima problem. Thus, the DE technique is introduced into the conventional AOA to form DAOA. This proposed DAOA method is introduced to perform the exploration search by the AOA and exploitation search by the DE. This also makes an excellent balance between the search strategies and guarantees that the proposed method averts the local optima.
Figure 3 depicts the proposed DAOA approach in this section. The DAOA procedure begins with (1) determining the values of the used algorithms’ parameters, (2) generating candidate solutions, (3) calculating fitness functions, (4) selecting the best solution, (5) if a given condition is true, the AOA is executed to update the solutions; otherwise, the DE is executed to update the solutions, and (6) then another condition is given to stop or continue the optimization process. Figure 3 shows the flowchart for the proposed DAOA. The pseudo-code of the DAOA algorithm is given in Algorithm 1.
Algorithm 1 Pseudo-code of the DAOA algorithm
1:
Initialize the Arithmetic Optimization Algorithm parameters α , μ .
2:
Initialize the solutions’ positions randomly. (Solutions: i = 1, …, N.).
3:
Calculate the Fitness values.
4:
while ( C _ I t e r < M _ I t e r ) do
5:
 Find the best solution (determined best so far).
6:
 Update the MOA value using Equation (1).
7:
 Update the MOP value using Equation (3).
8:
 Calculate the Fitness Function ( F F ) for the given solutions.
9:
for (i = 1 to S o l u t i o n s ) do
10:
  if rand < 0.5 then
11:
   Generate a random values between [0, 1] ( r 1 , r 2 , and r 3 )
12:
   if  r 1 > MOA then
13:
    if  r 2 > 0.5 then
14:
     Update the i th solutions’ positions using the first rule in Equation (2).
15:
    else
16:
     Update the i th solutions’ positions using the second rule in Equation (2).
17:
    end if
18:
   else
19:
    if  r 3 > 0.5 then
20:
     Update the i th solutions’ positions using the first rule in Equation (4).
21:
    else
22:
     Update the i th solutions’ positions using the second rule in Equation (4).
23:
    end if
24:
   end if
25:
  else
26:
   if rand < 0.5 then
27:
    Update the i th solutions’ positions using Mutation operator as given in Equation (5).
28:
   else
29:
    Update the i th solutions’ positions using Crossover operator as given in Equation (6).
30:
   end if
31:
  end if
32:
end for
33:
 C_Iter = C_Iter + 1
34:
end while
35:
Return the best solution (x).

3. Definitions of the Multilevel Thresholding Image Segmentation Problems

In this section, we describe the main problem of multilevel thresholding. Let us suppose that I is a gray or color image that needs to be processed, where K + 1 presents the classes that need to be produced. For segmenting the given image (I) into K + 1 classes, the k thresholds’ values are required to progress in the image segmentation procedure; { t k , k = 1 . . . . . . . . . . K } , and this can be expressed as follows [1,7,49].
C 0 = { I i , j 0 I i , j t 1 1 } , C 1 = { I i , j t 1 I i , j t 2 1 } , C K = { I i , j t K I i , j L 1 }
where L indicates the highest gray levels and C K indicates the kth class of the image I. The t k is the k-th threshold, with I i , j being the gray level at the ( i , j ) th pixel. Furthermore, in Equation (10), multilevel thresholding is identified as a maximization optimization problem that needs to find the optimal threshold values.
K multilevel threshold values can be presented as follows.
t 1 , , t 2 , , , t K , = arg max t 1 , , t K F i t ( t 1 , , t K )

3.1. Fitness Function (Kapur’s Entropy)

For the purpose of thresholding, consider a digital image I with N pixels and L gray levels. Via thresholds, these L number of gray levels are divided into classes: Class1, Class2, …, Classk [1].
In this proposed DAOA, Kapur’s entropy is utilized for achieving optimum threshold values. Measurement of the bi-level thresholds needs the optimization process’s objective function, as shown in Equation (11).
F i t ( t 1 , , t K ) = k = 1 , K H i
H k = i = 0 L 1 p i × μ k ( i ) P k × l n ( p i × μ k ( i ) P k ) ,
P k = i = 0 L 1 p i × μ k ( i )
μ 1 ( l ) = 1 l a 1 l c 1 a 1 c 1 a 1 l c 1 0 l > c 1 μ K ( l ) = 1 l a K 1 l a K c K a K a K 1 < l c K 1 0 l > c K 1
where p i is the probability distribution, h ( i ) is the numbers of pixels for the used gray level L, and N p is the total numbers of pixels of the image I. p i presents the probability value for the distribution, determined as p i = h ( i ) / N p ( 0 < i < L 1 ). h ( i ) and N k are the numbers of pixels for the used gray level L and total pixel of the image I. a 1 , c 1 , . , a k 1 , c k 1 are the used fuzzy parameters, and 0 a 1 c 1 a K 1 c K 1 .
Then, t 1 = a 1 + c 1 2 , t 2 = a 2 + c 2 2 , . . . , t K 1 = a K 1 + c K 1 2 . The best fitness function obtained is the highest value.

3.2. Performance Measures

We assess the proposed DAOA method performance, using three performance measures: the fitness function value, the Structural Similarity Index (SSIM), and the Peak Signal-to-Noise Ratio (PSNR) [50,51]. The following equations compute SSIM and PSNR:
S S I M ( I , I S ) = ( 2 μ I μ I S + c 1 ) ( 2 σ I , I S + c 2 ) ( μ I 2 + μ I S 2 + c 1 ) ( σ I 2 + σ I S 2 + c 2 )
where μ I S ( σ I S ) and μ I ( σ I ) are the images’ mean intensity of I S and I, respectively, where σ I , I S is the governance of I and I S , and c 1 and c 2 coefficient values are equal to 6.5025 and 58.52252, respectively [1].
P S N R = 20 l o g 10 ( 255 R M S E ) , R M S E = i = 1 N r j = 1 N c ( I i , j I S i , j ) 2 N r × N c
where the R M S E is the root-mean-squared error of each pixel, and M × N depicts the image’s size. I i , j is the gray pixel value of the initial image, and I s i j is the gray value of the pixel in the obtained segmented image.

4. Experiments and Results

4.1. Benchmark Images

In this section, the benchmark image data sets are presented in Figure 4 and Figure 5. Two image types were used in this paper’s experiments, taken from nature (as seen in Figure 4) and medical CT images (as seen in Figure 5). We chose eight images: two-color images from nature (i.e., Test 1 and Test 2), two gray images from nature (i.e., Test 3 and Test 4), two COVID-19 CT images (i.e., Test 5 and Test 6), and two normal COVID-19 CT images (i.e., Test 7 and Test 8). These benchmarks were taken from the Berkeley Segmentation Data Set: Images and BIMCV-COVID19 [52].

4.2. Comparative Algorithms and Parameter Setting

The proposed DAOA is analyzed and compared with six recently well-known algorithms, including Aquila Optimizer (AO) [30], Whale Optimization Algorithm (WOA) [53], Salp Swarm Algorithm (SSA) [54], Arithmetic Optimization Algorithm (AOA) [29], Particle Swarm Optimization (PSO) [55], Marine Predators Algorithm (MPA) [56], and Differential Evolution (DE) [57].
These algorithms’ parameters are set in the same way as they were in their original papers. The values of different parameter settings used in the tested algorithms are shown in Table 1. These sensitive parameters can be tuned for further investigation to show the effect of each parameter on the performance of the tested methods. The algorithms are executed by using the MATLAB 2015a software. These algorithms are run on an Intel Core i7 1.80 GHz 2.30 GHz processor with 16 GB RAM. The number of solutions used is twenty-five. For a systematic comparison, the maximum number of iterations is set to one hundred. Each competitor algorithm generates thirty independent runs.

4.3. Performance Evaluation

A comparison of the proposed DAOA for multilevel thresholding segmentation, using eight different images, is presented in this section. The following tables show the max, mean, min, and standard deviation of each test case’s PSNR and SSIM. Moreover, the summation, mean rank, and final ranking are given, using the Friedman ranking test to prove the proposed method’s significant improvement [58,59].
The PSNR and SSIM results of Test 1 are given in Table 2 and Table 3. It is clear that the proposed DAOA obtained excellent results in almost all the test cases in terms of PSNR. For threshold 2, the proposed DAOA obtained the best results, and it ranked first when compared to all other comparative methods, followed by AOA, SSA, PSO, WOA, MPA, AO, and finally, DE. In addition, for threshold 6, the proposed method obtained promising results compared to other methods. DAOA obtained the first rank, followed by WOA, PSO, SSA, AOA, AO, SSA, MPA, and DE. Overall, we can see that the proposed method obtained the first ranking, followed by AOA, PSO, SSA, WOA, AO, MPA, and DE. The obtained results in this table prove the ability of the proposed DAOA to solve the given problems effectively.
For threshold 2 in Table 3, the proposed DAOA obtained the best results, and it ranked as the first method, compared to all other comparative methods, followed by PSO, SSA, DE, AOA, WOA, MPA, and finally, AO. In addition, for threshold 3, the proposed method obtained promising results, compared to other methods. DAOA obtained the first ranking, followed by PSO, SSA, DE, AOA, MPA, WOA, and AO. Overall, we can see that the proposed DAOA method obtained the first ranking, followed by PSO, AOA, SSA, DE, WOA, AO, and MPA. The obtained results in this table prove the ability of the proposed DAOA to solve the given problems effectively.
The PSNR and SSIM results of Test 2 are given in Table 4 and Table 5. The proposed DAOA achieved excellent results in almost all the test cases in terms of PSNR. For threshold 4, the proposed DAOA obtained the best results, and it ranked as the first method, compared to all other comparative methods, followed by AO, MPA, AOA, SSA, DE, PSO, and finally, WOA. For threshold 5, the proposed method obtained promising results, compared to other methods. DAOA obtained the first rank, followed by AO, DE, PSO, WOA, AOA, SSA, and MPA. Overall, we can see that the proposed method obtained the first ranking, followed by AO, DE, SSA, PSO, MPA, WOA, and AOA. The achieved results in this table demonstrate the ability of the proposed DAOA to solve the given problems efficiently.
For threshold 4 in Table 5, the proposed DAOA obtained the best results, and it ranked as the first method, compared to all other comparative methods, followed by AO, AOA, MPA, SSA, DE, PSO, and finally, WOA. For threshold 3, the proposed method obtained promising results, compared to other methods. DAOA obtained the first ranking, followed by WOA, DE, AO, PSO, SSA, MPA, and AOA. Overall, we can see that the proposed DAOA method obtained the first ranking, followed by AO, MPA, SSA, DE, PSO, WOA, and AOA. The obtained results in this table confirm the performance of the proposed DAOA and its ability to solve the given problems efficiently.
The PSNR and SSIM results of Test 3 are given in Table 6 and Table 7. The proposed DAOA obtained new, promising results in almost all the test cases in terms of PSNR. For threshold 5, the proposed DAOA obtained the best results, and it ranked as the first method, compared to all other comparative methods, followed by SSA, AO, PSO, AOA, MPA, and finally, DE. For threshold 6, the proposed method obtained promising results compared to other methods. DAOA obtained the first rank, followed by AO, PSO, DE, WOA, AOA, SSA, and MPA. Overall we can see that the proposed method obtained the first ranking, followed by WOA, DE, AO, AOA, MPA, PSO, and SSA. The achieved results in this table demonstrate the ability of the proposed DAOA to solve the given problems efficiently. As well, it is clear the proposed DAOA has this ability at different threshold levels.
For threshold 2 in Table 7, the proposed DAOA obtained the best results, and it ranked as the first method, compared to all other comparative methods, followed by MPA, WOA, DE, AO, SSA, AOA, and finally, PSO. For threshold 5, the proposed method obtained promising results compared to other methods. DAOA obtained the first ranking, followed by PSO, AOA, MPA, WOA, AO, PSO, and DE. Overall, we can see that the proposed DAOA method obtained the first ranking, followed by MPA, AO, SSA, AOA, WOA, PSO, and DE. The obtained results in this table confirm the performance of the proposed DAOA and its ability to solve the given problems efficiently. The following results prove and support that the proposed algorithm’s ability to solve such problems is strong and that it is capable of finding robust solutions in this field.
The PSNR and SSIM results of Test 4 are given in Table 8 and Table 9. The proposed DAOA obtained new promising results in almost all the test cases in terms of PSNR, as shown in Table 8. The proposed DAOA obtained the best results for two threshold values (i.e., 3 and 4 levels). For threshold 3, the proposed DAOA obtained the best results, and it ranked as the first method, compared to all other comparative methods, followed by PSO, SSA, WOA, AOA, MPA, DE, and finally, AO. For threshold 4, the proposed method obtained promising results, compared to other methods. DAOA obtained the first rank, followed by WOA, SSA, AOA, MPA, AO, DE, and AOA. Overall, we can see that the proposed method obtained the first ranking, followed by WOA, PSO, SSA, MPA, AO, AOA, and DE. The achieved results in this table demonstrate the ability of the proposed DAOA to solve the given problems efficiently. As well, it is clear the proposed DAOA has this ability at different threshold levels.
Table 9 shows that the proposed DAOA method obtained better results in almost all the test cases in terms of SSIM for Test 4. The proposed DAOA obtained the best results for two threshold values (i.e., 2 and 3 levels). For threshold 2, the proposed DAOA obtained the best results, and it ranked as the first method, compared to all other comparative methods, followed by DE, AOA, WOA, MPA, PSO, SSA, and finally, AOA. For threshold 3, the proposed method obtained promising results, compared to other methods. DAOA obtained the first ranking, followed by AOA, AO, PSO, WOA, MPA, SSA, and DE. Overall, we can see that the proposed DAOA method obtained the first ranking, followed by PSO, AO, WOA, AOA, MPA, SSA, and DE. The obtained results in this table confirm the performance of the proposed DAOA to solve the given problems efficiently. The following results prove and support that the proposed algorithm’s ability to solve such problems is strong and that it is capable of finding robust solutions in this field.
In Table 10 and Table 11, the PSNR and SSIM results of Test 5 are shown. As shown in Table 10, the proposed DAOA yielded new promising PSNR results in almost all test cases. For two threshold values, the proposed DAOA gave the best results (i.e., 2 and 5 levels). For threshold 2, the proposed DAOA produced the best results, placing it first among all other comparative methods, ahead of SSA, AOA, MPA, DE, AO, WOA, and PSO. In addition, when compared to other methods, the proposed method produced positive results for threshold 5. DAOA came first, followed by AOA, MPA, AO, WOA, PSO, and DE. Overall, we can see that DAOA came first, followed by SSA, DE, AO, AOA, PSO, MPA, and WOA. The obtained results in this table demonstrate the proposed DAOA’s ability to solve the given problems efficiently. Furthermore, it is evident that the proposed DAOA has the potential to operate at various threshold levels.
In terms of SSIM for Test 5, Table 11 shows that the proposed DAOA system obtained better results in almost all test cases. For two threshold values, the proposed DAOA produced the best results (i.e., 2 and 5 levels). For threshold 2, the proposed DAOA received the best results, placing it first among SSA, AO, MPA, AOA, DE, PSO, and WOA. In addition, when compared to other methods, the proposed method produced positive results for threshold 5. The first-place winner was DAOA, followed by AO, WOA, SSA, AOA, DE, MPA, and PSO. Overall, we can see that the proposed DAOA method obtained the first ranking, followed by AO, SSA, WOA, PSO, MPA, DE, and AOA. The obtained results in this table confirm the performance of the proposed DAOA and its ability to solve the given problems efficiently. The following results prove and support that the proposed algorithm’s ability to solve such problems is strong and that it is capable of finding robust solutions in this field.
In Table 12 and Table 13, the PSNR and SSIM results of Test 6 are shown. As shown in Table 12, the proposed DAOA yielded new promising PSNR results in nearly all test cases. For two threshold values, the proposed DAOA gave the best results (i.e., 5 and 6 levels). For threshold 5, the DE produced the best results, placing it first among all other comparative approaches, ahead of DAOA, SSA, PSO, AOA, WOA, MPA, and AO. In addition, when compared to other methods, the proposed method produced positive results for threshold 6. DAOAO came first, followed by SSA, WOA, AOA, MPA, DE, AO, and PSO. Overall, we can see that DAOA came first, followed by AOA, DE, PSO, SSA, WOA, AO, and MPA. The obtained results in this table demonstrate the proposed DAOA’s ability to solve the given problems efficiently. Furthermore, it is evident that the proposed DAOA has the potential to operate at various threshold levels.
In terms of SSIM for Test 6, Table 13 shows that the proposed DAOA method obtained better results in almost all test cases. For one threshold value, the proposed DAOA produced the best results (i.e., three levels). For threshold 3, the proposed DAOA received the best results, placing it first, followed by AOA, AO, PSO, MPA, SSA, WOA, and DE. Overall, we can see that the proposed DAOA method obtained the first ranking, followed by AOA, SSA, MPA, AO, WOA, PSO, and DE. The obtained results in this table confirm the performance of the proposed DAOA to solve the given problems efficiently. The following results prove and support that the proposed algorithm’s ability to solve such problems is strong and that it is capable of finding robust solutions in this field.
The PSNR and SSIM results of Test 7 are given in Table 14 and Table 15. The proposed DAOA obtained new promising results in almost all the test cases in terms of PSNR, as shown in Table 14. The proposed DAOA obtained the best results for three threshold values (i.e., 3, 4, and 6 levels). For threshold 3, the proposed DAOA obtained the best results, and it ranked as the first method, compared to all other comparative methods, followed by PSO, WOA, MPA, PSO, SSA, AOA, and finally, DE. For threshold 6, the proposed method obtained promising results, compared to other methods. DAOA obtained the first rank, followed by MPA, AO, AOA, SSA, PSO, DE, and WOA. Overall, we can see that the proposed method obtained the first ranking, followed by MPA, AO, AOA, SSA, PSO, WOA, and DE. The achieved results in this table demonstrate the ability of the proposed DAOA to solve the given problems efficiently. Furthermore, it is obvious that the proposed DAOA has the potential to operate at various threshold levels.
Table 15 shows that the proposed DAOA method obtained better results in almost all the test cases in terms of SSIM for Test 7. The proposed DAOA obtained the best results for two threshold values (i.e., 2 and 4 levels). For threshold 2, the proposed DAOA obtained the best results, and it ranked as the first method, compared to all other comparative methods, followed by MPA, PSO, AOA, AO, WOA, SSA, and finally, DE. For threshold 4, the proposed method obtained promising results compared to other methods. DAOA obtained the first ranking, followed by PSO, SSA, AO, AOA, WOA, DE, and MPA. Overall, we can see that the proposed DAOA method obtained the first ranking, followed by PSO, AOA, AO, MPA, WOA, SSA, and DE. The obtained results in this table confirm the performance of the proposed DAOA to solve the given problems efficiently. The presented results demonstrate and declare the proposed algorithm’s ability to solve such problems and find reliable solutions in this area.
The PSNR and SSIM results of Test 8 are given in Table 16 and Table 17. The proposed DAOA obtained got new promising results in almost all the test cases in terms of PSNR, as shown in Table 16. The proposed DAOA obtained the best results for four threshold values (i.e., 2, 3, 4, and 5 levels). For threshold 2, the proposed DAOA obtained the best results, and it ranked as the first method, compared to all other comparative methods, followed by DE, AOA, PSO, SSA, AO, MPA, and finally, WOA. For threshold 5, the proposed method obtained promising results, compared to other methods. DAOA obtained the first rank, followed by AO, DE, PSO, WOA, AOA, MPA, and SSA. Overall, we can see that the proposed method obtained the first ranking, followed by AO, AOA, PSO, DE, WOA, MPA, and SSA. The achieved results in this table demonstrate the ability of the proposed DAOA to solve the given problems efficiently. Furthermore, it is obvious that the proposed DAOA has the potential to operate at various threshold levels.
Table 17 shows that the proposed DAOA method obtained better results in almost all the test cases in terms of SSIM for Test 8. The proposed DAOA obtained the best results for two threshold values (i.e., 3 and 5 levels). For threshold 3, the proposed DAOA obtained the best results, and it ranked as the first method, compared to all other comparative methods, followed by AOA, SSA, DE, WOA, AO, PSO, and finally, MPA. For threshold 5, the proposed method obtained promising results compared to other methods. DAOA obtained the first ranking, followed by AO, DE, SSA, PSO, AOA, WOA, and MPA. Overall, we can see that the proposed DAOA method obtained the first ranking, followed by AOA, AO, WOA, SSA, PSO, DE, and MPA. The obtained results in this table confirm the performance of the proposed DAOA and its ability to solve the given problems efficiently. The presented results demonstrate and declare the proposed algorithm’s ability to solve such problems and find reliable solutions in this area. We added a Wilcoxon sign test to show the significant improvement for test case number 8 as shown in Table 16 and Table 17. It is clear that the proposed method is more effective and better than the other methods.
The segmentation results (segmented images) of the proposed DAOA and the other comparative methods for Test 8 are shown in Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10. Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10 show the segmented images for all the tested methods, when the threshold values are 2, 3, 4, 5, and 6, respectively. According to these figures, we can recognize that the proposed DAOA showed good segmented images for various images (CT COVID-19 medical images) under different thresholds. Additionally, these figures prove that the segmented images are better in terms of quality when the threshold value is higher.
The thresholds are shown in Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15, applied over the selected images. In Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15, the histogram images are given with the best threshold values obtained by the comparative methods for Test 8, where the threshold values are taken (i.e., 2, 3, 4, 5, and 6). The X and Y axes present the threshold values and Kapur measure values, respectively. It is feasible to recognize that the histogram classes are uniformly created, even in complex situations from such images. This means that the proposed method has an excellent ability to find always the same threshold values. The complexity is different from case to case because of the various peaks displayed in the pixels’ distribution, which could create multiple classes or even carefully obtain the selection of the optimal thresholds.
Figure 16 shows the convergence curves of the proposed DAOA and its comparative optimization algorithms on eight tested images (i.e., Test 1 to Test 8); it can be seen that the proposed DAOA performs better than all involved other optimization methods in Test 8 when the threshold value is 6. For almost all the test images, the excellent optimized performance with accelerated convergence and more reliable accuracy achieved by the proposed DAOA can be seen as being remarkably smoothing behavior in the convergence curve. Moreover, we recognize that the curves of the proposed method always converge smoothly, reflecting the proposed DAOA’s ability to avoid the common problem (local optima). In the end, the proposed DAOA reached the best solutions almost in all the tested cases, compared to the other comparative methods, as clearly shown in Figure 16.

5. Conclusions and Future Works

The most crucial aspect of image segmentation is multilevel thresholding. However, multilevel thresholding displays require increasingly more computational complexity as the number of thresholds grows. In order to address this weakness, this paper proposes a new multilevel thresholding approach based on using an improved optimization-based evolutionary method.
The Arithmetic Optimization Algorithm (AOA) is a recently proposed optimization technique to solve different complex optimization problems. An enhanced version of the Arithmetic Optimization Algorithm is proposed in this paper to solve multilevel thresholding image segmentation problems. The proposed method combines the conventional Arithmetic Optimization Algorithm with the Differential Evolution technique, called DAOA. The main aim of the proposed DAOA is to improve the local search of the Arithmetic Optimization Algorithm and to establish an equilibrium among the search methods (exploration and exploitation).
The proposed DAOA method was applied to the multilevel thresholding problem, using Kapur’s measure between class variance functions as a fitness function. The proposed DAOA evaluated eight standard test images from two different groups: nature images and CT medical images (i.e., COVID-19). The Peak Signal-to-Noise Ratio (PSNR) and Structural Similarity Index Test (SSIM) were used to determine the segmented images’ accuracy. The proposed DAOA method’s efficiency was evaluated and compared to other multilevel thresholding methods, including the Aquila Optimizer (AO), Whale Optimization Algorithm (WOA), Salp Swarm Algorithm (SSA), Arithmetic Optimization Algorithm (AOA), Particle Swarm Optimization (PSO), and Marine Predators Algorithm (MPA). The findings were presented, using a number of different threshold values (i.e., 2, 3, 4, 5, and 6). According to the experimental results, the proposed DAOA produced higher quality solutions than the other approaches. It achieved better results in almost all the tested cases, compared to other methods.
For future work, other fitness functions, evaluation measures, and benchmark images can be used. The conventional Arithmetic Optimization Algorithm can be improved, using other different optimization operations to enhance its performance further. As well, the proposed DAOA method can be used to solve other problems, such as text clustering, feature selection, photovoltaic parameter estimations, task scheduling in fog and cloud computing, appliances management in smart homes, advanced benchmark functions, text classification, text summarization, data clustering, engineering design problems, industrial problems, image construction, short-term wind speed forecasting, fuel cell modeling, damage identification, the prediction of the software vulnerability, knapsack problems, and others.

Author Contributions

L.A.: data curation, formal analysis, investigation, methodology, resources, software, supervision, validation, visualization, writing—original draft, and writing—review and editing. A.D.: formal analysis, investigation, supervision, writing—original draft, and writing—review and editing. P.S.: investigation, supervision, and writing—review and editing. A.H.G.: investigation, supervision, and writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Abd El Aziz, M.; Ewees, A.A.; Hassanien, A.E. Whale optimization algorithm and moth-flame optimization for multilevel thresholding image segmentation. Expert Syst. Appl. 2017, 83, 242–256. [Google Scholar] [CrossRef]
  2. Shubham, S.; Bhandari, A.K. A generalized Masi entropy based efficient multilevel thresholding method for color image segmentation. Multimed. Tools Appl. 2019, 78, 17197–17238. [Google Scholar] [CrossRef]
  3. Abd Elaziz, M.; Lu, S. Many-objectives multilevel thresholding image segmentation using knee evolutionary algorithm. Expert Syst. Appl. 2019, 125, 305–316. [Google Scholar] [CrossRef]
  4. Pare, S.; Bhandari, A.K.; Kumar, A.; Singh, G.K. An optimal color image multilevel thresholding technique using grey-level co-occurrence matrix. Expert Syst. Appl. 2017, 87, 335–362. [Google Scholar] [CrossRef]
  5. Bao, X.; Jia, H.; Lang, C. A novel hybrid harris hawks optimization for color image multilevel thresholding segmentation. IEEE Access 2019, 7, 76529–76546. [Google Scholar] [CrossRef]
  6. Abd Elaziz, M.; Ewees, A.A.; Oliva, D. Hyper-heuristic method for multilevel thresholding image segmentation. Expert Syst. Appl. 2020, 146, 113201. [Google Scholar] [CrossRef]
  7. Houssein, E.H.; Helmy, B.E.d.; Oliva, D.; Elngar, A.A.; Shaban, H. A novel Black Widow Optimization algorithm for multilevel thresholding image segmentation. Expert Syst. Appl. 2021, 167, 114159. [Google Scholar] [CrossRef]
  8. Gill, H.S.; Khehra, B.S.; Singh, A.; Kaur, L. Teaching-learning-based optimization algorithm to minimize cross entropy for Selecting multilevel threshold values. Egypt. Inform. J. 2019, 20, 11–25. [Google Scholar] [CrossRef]
  9. Tan, Z.; Zhang, D. A fuzzy adaptive gravitational search algorithm for two-dimensional multilevel thresholding image segmentation. J. Ambient. Intell. Humaniz. Comput. 2020, 11, 4983–4994. [Google Scholar] [CrossRef]
  10. Yousri, D.; Abd Elaziz, M.; Abualigah, L.; Oliva, D.; Al-Qaness, M.A.; Ewees, A.A. COVID-19 X-ray images classification based on enhanced fractional-order cuckoo search optimizer using heavy-tailed distributions. Appl. Soft Comput. 2021, 101, 107052. [Google Scholar] [CrossRef] [PubMed]
  11. Srikanth, R.; Bikshalu, K. Multilevel thresholding image segmentation based on energy curve with harmony Search Algorithm. Ain Shams Eng. J. 2021, 12, 1–20. [Google Scholar] [CrossRef]
  12. Duan, L.; Yang, S.; Zhang, D. Multilevel thresholding using an improved cuckoo search algorithm for image segmentation. J. Supercomput. 2021, 77, 6734–6753. [Google Scholar] [CrossRef]
  13. Jia, H.; Lang, C.; Oliva, D.; Song, W.; Peng, X. Hybrid grasshopper optimization algorithm and differential evolution for multilevel satellite image segmentation. Remote Sens. 2019, 11, 1134. [Google Scholar] [CrossRef] [Green Version]
  14. Khairuzzaman, A.K.M.; Chaudhury, S. Multilevel thresholding using grey wolf optimizer for image segmentation. Expert Syst. Appl. 2017, 86, 64–76. [Google Scholar] [CrossRef]
  15. He, L.; Huang, S. Modified firefly algorithm based multilevel thresholding for color image segmentation. Neurocomputing 2017, 240, 152–174. [Google Scholar] [CrossRef]
  16. Li, Y.; Bai, X.; Jiao, L.; Xue, Y. Partitioned-cooperative quantum-behaved particle swarm optimization based on multilevel thresholding applied to medical image segmentation. Appl. Soft Comput. 2017, 56, 345–356. [Google Scholar] [CrossRef]
  17. Manic, K.S.; Priya, R.K.; Rajinikanth, V. Image multithresholding based on Kapur/Tsallis entropy and firefly algorithm. Indian J. Sci. Technol. 2016, 9, 89949. [Google Scholar] [CrossRef]
  18. Bhandari, A.K.; Kumar, A.; Singh, G.K. Modified artificial bee colony based computationally efficient multilevel thresholding for satellite image segmentation using Kapur’s, Otsu and Tsallis functions. Expert Syst. Appl. 2015, 42, 1573–1601. [Google Scholar] [CrossRef]
  19. Liang, H.; Jia, H.; Xing, Z.; Ma, J.; Peng, X. Modified grasshopper algorithm-based multilevel thresholding for color image segmentation. IEEE Access 2019, 7, 11258–11295. [Google Scholar] [CrossRef]
  20. Abualigah, L. Group search optimizer: A nature-inspired meta-heuristic optimization algorithm with its results, variants, and applications. Neural Comput. Appl. 2020, 33, 2949–2972. [Google Scholar] [CrossRef]
  21. Alsalibi, B.; Abualigah, L.; Khader, A.T. A novel bat algorithm with dynamic membrane structure for optimization problems. Appl. Intell. 2021, 51, 1992–2017. [Google Scholar] [CrossRef]
  22. Ewees, A.A.; Abualigah, L.; Yousri, D.; Algamal, Z.Y.; Al-qaness, M.A.; Ibrahim, R.A.; Abd Elaziz, M. Improved Slime Mould Algorithm based on Firefly Algorithm for feature selection: A case study on QSAR model. Eng. Comput. 2021, 1–15. [Google Scholar] [CrossRef]
  23. Şahin, C.B.; Dinler, Ö.B.; Abualigah, L. Prediction of software vulnerability based deep symbiotic genetic algorithms: Phenotyping of dominant-features. Appl. Intell. 2021, 1–17. [Google Scholar] [CrossRef]
  24. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey wolf optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef] [Green Version]
  25. Safaldin, M.; Otair, M.; Abualigah, L. Improved binary gray wolf optimizer and SVM for intrusion detection system in wireless sensor networks. J. Ambient. Intell. Humaniz. Comput. 2021, 12, 1559–1576. [Google Scholar] [CrossRef]
  26. Alshinwan, M.; Abualigah, L.; Shehab, M.; Abd Elaziz, M.; Khasawneh, A.M.; Alabool, H.; Al Hamad, H. Dragonfly algorithm: A comprehensive survey of its results, variants, and applications. Multimed. Tools Appl. 2021, 80, 14979–15016. [Google Scholar] [CrossRef]
  27. Shehab, M.; Abualigah, L.; Al Hamad, H.; Alabool, H.; Alshinwan, M.; Khasawneh, A.M. Moth–flame optimization algorithm: Variants and applications. Neural Comput. Appl. 2020, 32, 9859–9884. [Google Scholar] [CrossRef]
  28. Al-Qaness, M.A.; Ewees, A.A.; Fan, H.; Abualigah, L.; Abd Elaziz, M. Marine predators algorithm for forecasting confirmed cases of COVID-19 in Italy, USA, Iran and Korea. Int. J. Environ. Res. Public Health 2020, 17, 3520. [Google Scholar] [CrossRef]
  29. Abualigah, L.; Diabat, A.; Mirjalili, S.; Abd Elaziz, M.; Gandomi, A.H. The arithmetic optimization algorithm. Comput. Methods Appl. Mech. Eng. 2021, 376, 113609. [Google Scholar] [CrossRef]
  30. Abualigah, L.; Yousri, D.; Abd Elaziz, M.; Ewees, A.A.; Al-qaness, M.A.; Gandomi, A.H. Aquila Optimizer: A novel meta-heuristic optimization algorithm. Comput. Ind. Eng. 2021, 157, 107250. [Google Scholar] [CrossRef]
  31. Gandomi, A.H.; Alavi, A.H. Krill herd: A new bio-inspired optimization algorithm. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 4831–4845. [Google Scholar] [CrossRef]
  32. Jouhari, H.; Lei, D.; Al-qaness, M.A.; Elaziz, M.A.; Damaševičius, R.; Korytkowski, M.; Ewees, A.A. Modified Harris Hawks Optimizer for Solving Machine Scheduling Problems. Symmetry 2020, 12, 1460. [Google Scholar] [CrossRef]
  33. Połap, D.; Woźniak, M. Red fox optimization algorithm. Expert Syst. Appl. 2021, 166, 114107. [Google Scholar] [CrossRef]
  34. Mernik, M.; Liu, S.H.; Karaboga, D.; Črepinšek, M. On clarifying misconceptions when comparing variants of the artificial bee colony algorithm by offering a new implementation. Inf. Sci. 2015, 291, 115–127. [Google Scholar] [CrossRef]
  35. Sahlol, A.T.; Abd Elaziz, M.; Tariq Jamal, A.; Damaševičius, R.; Farouk Hassan, O. A novel method for detection of tuberculosis in chest radiographs using artificial ecosystem-based optimisation of deep neural network features. Symmetry 2020, 12, 1146. [Google Scholar] [CrossRef]
  36. Abualigah, L.; Diabat, A. Advances in sine cosine algorithm: A comprehensive survey. Artif. Intell. Rev. 2021, 54, 2567–2608. [Google Scholar] [CrossRef]
  37. Abualigah, L.; Diabat, A. A comprehensive survey of the Grasshopper optimization algorithm: Results, variants, and applications. Neural Comput. Appl. 2020, 32, 15533–15556. [Google Scholar] [CrossRef]
  38. Tuba, E.; Alihodzic, A.; Tuba, M. Multilevel image thresholding using elephant herding optimization algorithm. In Proceedings of the 2017 14th International Conference on Engineering of Modern Electric Systems (EMES), Oradea, Romania, 1–2 June 2017; pp. 240–243. [Google Scholar]
  39. Sahlol, A.T.; Yousri, D.; Ewees, A.A.; Al-Qaness, M.A.; Damasevicius, R.; Abd Elaziz, M. COVID-19 image classification using deep features and fractional-order marine predators algorithm. Sci. Rep. 2020, 10, 1–15. [Google Scholar] [CrossRef] [PubMed]
  40. Wang, X.; Pan, J.S.; Chu, S.C. A parallel multi-verse optimizer for application in multilevel image segmentation. IEEE Access 2020, 8, 32018–32030. [Google Scholar] [CrossRef]
  41. Gao, Y.; Li, X.; Dong, M.; Li, H.P. An enhanced artificial bee colony optimizer and its application to multi-level threshold image segmentation. J. Cent. South Univ. 2018, 25, 107–120. [Google Scholar] [CrossRef]
  42. Resma, K.B.; Nair, M.S. Multilevel thresholding for image segmentation using Krill Herd Optimization algorithm. J. King Saud. Univ.-Comput. Inf. Sci. 2018, 33, 528–541. [Google Scholar]
  43. Abd Elaziz, M.; Yousri, D.; Al-qaness, M.A.; AbdelAty, A.M.; Radwan, A.G.; Ewees, A.A. A Grunwald–Letnikov based Manta ray foraging optimizer for global optimization and image segmentation. Eng. Appl. Artif. Intell. 2021, 98, 104105. [Google Scholar] [CrossRef]
  44. Sörensen, K. Metaheuristics—The metaphor exposed. Int. Trans. Oper. Res. 2015, 22, 3–18. [Google Scholar] [CrossRef]
  45. García-Martínez, C.; Gutiérrez, P.D.; Molina, D.; Lozano, M.; Herrera, F. Since CEC 2005 competition on real-parameter optimisation: A decade of research, progress and comparative analysis’s weakness. Soft Comput. 2017, 21, 5573–5583. [Google Scholar] [CrossRef]
  46. Storn, R.; Price, K. Differential evolution—A simple and efficient heuristic for global optimization over continuous spaces. J. Glob. Optim. 1997, 11, 341–359. [Google Scholar] [CrossRef]
  47. Tanabe, R.; Fukunaga, A. Success-history based parameter adaptation for differential evolution. In Proceedings of the 2013 IEEE Congress on Evolutionary Computation, Cancun, Mexico, 20–23 June 2013; pp. 71–78. [Google Scholar]
  48. Črepinšek, M.; Liu, S.H.; Mernik, M. Replication and comparison of computational experiments in applied evolutionary computing: Common pitfalls and guidelines to avoid them. Appl. Soft Comput. 2014, 19, 161–170. [Google Scholar] [CrossRef]
  49. Maitra, M.; Chatterjee, A. A hybrid cooperative–comprehensive learning based PSO algorithm for image segmentation using multilevel thresholding. Expert Syst. Appl. 2008, 34, 1341–1350. [Google Scholar] [CrossRef]
  50. Yin, P.Y. Multilevel minimum cross entropy threshold selection based on particle swarm optimization. Appl. Math. Comput. 2007, 184, 503–513. [Google Scholar] [CrossRef]
  51. Zhou, W. Image quality assessment: From error measurement to structural similarity. IEEE Trans. Image Process. 2004, 13, 600–613. [Google Scholar]
  52. Sumari, P.; Syed, S.J.; Abualigah, L. A Novel Deep Learning Pipeline Architecture based on CNN to Detect Covid-19 in Chest X-ray Images. Turk. J. Comput. Math. Educ. (TURCOMAT) 2021, 12, 2001–2011. [Google Scholar]
  53. Mirjalili, S.; Lewis, A. The whale optimization algorithm. Adv. Eng. Softw. 2016, 95, 51–67. [Google Scholar] [CrossRef]
  54. Mirjalili, S.; Gandomi, A.H.; Mirjalili, S.Z.; Saremi, S.; Faris, H.; Mirjalili, S.M. Salp Swarm Algorithm: A bio-inspired optimizer for engineering design problems. Adv. Eng. Softw. 2017, 114, 163–191. [Google Scholar] [CrossRef]
  55. Eberhart, R.; Kennedy, J. A new optimizer using particle swarm theory. In Proceedings of the MHS’95. Proceedings of the Sixth International Symposium on Micro Machine and Human Science, Nagoya, Japan, 4–6 October 1995; pp. 39–43. [Google Scholar]
  56. Faramarzi, A.; Heidarinejad, M.; Mirjalili, S.; Gandomi, A.H. Marine Predators Algorithm: A nature-inspired metaheuristic. Expert Syst. Appl. 2020, 152, 113377. [Google Scholar] [CrossRef]
  57. Price, K.; Storn, R.M.; Lampinen, J.A. Differential Evolution: A Practical Approach to Global Optimization; Science & Business Media: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  58. Derrac, J.; García, S.; Molina, D.; Herrera, F. A practical tutorial on the use of nonparametric statistical tests as a methodology for comparing evolutionary and swarm intelligence algorithms. Swarm Evol. Comput. 2011, 1, 3–18. [Google Scholar] [CrossRef]
  59. Abualigah, L.M.Q. Feature Selection and Enhanced Krill Herd Algorithm for Text Document Clustering; Springer: Berlin/Heidelberg, Germany, 2019. [Google Scholar]
Figure 1. The search phases of the Arithmetic Optimization Algorithm.
Figure 1. The search phases of the Arithmetic Optimization Algorithm.
Processes 09 01155 g001
Figure 2. Flowchart of the conventional AOA.
Figure 2. Flowchart of the conventional AOA.
Processes 09 01155 g002
Figure 3. The flowchart of the proposed DAOA.
Figure 3. The flowchart of the proposed DAOA.
Processes 09 01155 g003
Figure 4. The nature benchmark images that have been used.
Figure 4. The nature benchmark images that have been used.
Processes 09 01155 g004
Figure 5. The CT benchmark images that were used.
Figure 5. The CT benchmark images that were used.
Processes 09 01155 g005
Figure 6. The segmented image (Test 8) by the comparative methods when the threshold value is 2.
Figure 6. The segmented image (Test 8) by the comparative methods when the threshold value is 2.
Processes 09 01155 g006
Figure 7. The segmented image (Test 8) by the comparative methods when the threshold value is 3.
Figure 7. The segmented image (Test 8) by the comparative methods when the threshold value is 3.
Processes 09 01155 g007
Figure 8. The segmented image (Test 8) by the comparative methods when the threshold value is 4.
Figure 8. The segmented image (Test 8) by the comparative methods when the threshold value is 4.
Processes 09 01155 g008
Figure 9. The segmented image (Test 8) by the comparative methods when the threshold value is 5.
Figure 9. The segmented image (Test 8) by the comparative methods when the threshold value is 5.
Processes 09 01155 g009
Figure 10. The segmented image (Test 8) by the comparative methods when the threshold value is 6.
Figure 10. The segmented image (Test 8) by the comparative methods when the threshold value is 6.
Processes 09 01155 g010
Figure 11. The histogram image (Test 8) by the comparative methods when the threshold value is 2.
Figure 11. The histogram image (Test 8) by the comparative methods when the threshold value is 2.
Processes 09 01155 g011
Figure 12. The histogram image (Test 8) by the comparative methods when the threshold value is 3.
Figure 12. The histogram image (Test 8) by the comparative methods when the threshold value is 3.
Processes 09 01155 g012
Figure 13. The histogram image (Test 8) by the comparative methods when the threshold value is 4.
Figure 13. The histogram image (Test 8) by the comparative methods when the threshold value is 4.
Processes 09 01155 g013
Figure 14. The histogram image (Test 8) by the comparative methods when the threshold value is 5.
Figure 14. The histogram image (Test 8) by the comparative methods when the threshold value is 5.
Processes 09 01155 g014
Figure 15. The histogram image (Test 8) by the comparative methods when the threshold value is 6.
Figure 15. The histogram image (Test 8) by the comparative methods when the threshold value is 6.
Processes 09 01155 g015
Figure 16. The convergence behavior of the comparative methods in solving Test 8 when the threshold value is 6.
Figure 16. The convergence behavior of the comparative methods in solving Test 8 when the threshold value is 6.
Processes 09 01155 g016
Table 1. Parameter settings.
Table 1. Parameter settings.
No.AlgorithmReferenceParameterValue
1AO[30] α 0.1
/ δ 0.1
2WOA[53] α Decreased from 2 to 0
b2
3SSA[54] v 0 0
4AOA[29] α 5
μ 0.5
5PSO[55]TopologyFully connected
Cognitive and social constant(C1, C2) 2, 2
Inertia weightLinear reduction values [0.9 0.1]
Velocity limit10% of dimension range
6MPA[56] γ γ > 1
P0.0
7DE[57]Co0.5
Mu0.5
Table 2. The PSNR results of the test case 1.
Table 2. The PSNR results of the test case 1.
ThresholdMetricComparative Methods
AOWOASSAAOAPSOMPADEDAOA
2Max12.7447913.2036813.7252113.2753913.7008512.0261513.5155114.55374
Mean11.4368111.9774612.5872812.6157912.5133911.7940711.2585212.76203
Min10.6069710.1622210.9846812.0218411.8595211.6000010.3336511.86358
STD1.146311.604011.428130.629351.030130.215600.658851.55167
Ranking75324681
3Max16.4401114.2753516.0814016.8901314.1881315.1321014.4144015.73708
Mean15.1986313.2860514.4351014.6845613.4928614.3789512.2251414.61015
Min14.2850711.3666512.1875611.7352812.5331712.9104111.0221513.64585
STD1.114311.662512.015342.656680.858581.271940.566981.05506
Ranking17426583
4Max15.1765617.9569117.4723517.1643517.3283316.3005715.6585417.94836
Mean14.0018315.2204116.8726315.5912116.5222515.7283814.2548416.04748
Min13.3003712.5535415.9018014.3754416.0946915.1186113.9555815.09213
STD1.023722.702360.848491.428390.698520.591880.474471.64622
Ranking86152473
5Max16.7262216.4271016.2411017.9031216.3742016.3025616.3225418.67014
Mean15.5495316.0209615.5479116.9225915.8413815.2457115.2254115.86760
Min14.4954315.6180714.8844215.6238515.3476314.5795514.0255414.01993
STD1.120430.404520.678831.172480.514400.925570.655582.46778
Ranking52614783
6Max19.4358220.6194219.5234420.4318719.9683818.8674417.9510120.03906
Mean18.3878118.7539117.8551218.2343918.7172816.9271616.2587019.23425
Min16.3761317.0704014.7826114.8895617.4471314.5785515.3365217.83410
STD1.742671.781492.664142.943911.260732.173401.254121.21708
Ranking42653781
Summation2522201519293911
Mean Rank54.4433.85.87.82.2
Final Ranking65423781
Table 3. The SSIM results of the test case 1.
Table 3. The SSIM results of the test case 1.
ThresholdMetricComparative Methods
AOWOASSAAOAPSOMPADEDAOA
2Max0.2697170.3621460.3807570.2646720.4545870.2205160.3744540.385465
Mean0.1735750.2232270.2560720.2319060.2577140.19740.235550.277138
Min0.1166570.0197210.1845880.2049930.1333670.16060.0189850.122664
STD0.0837310.1801180.1083670.0302670.1724550.0322170.154150.137343
Ranking86352741
3Max0.4173740.6750720.6738590.5803740.6313890.5606410.645440.631938
Mean0.3536850.441140.5558180.5117860.5802120.5042520.514450.588044
Min0.3058380.2176710.4784850.4179650.5545520.4389450.485540.53806
STD0.057430.228880.1038540.0840940.044320.0613370.222520.047235
Ranking87352641
4Max0.4173740.6750720.6319380.5803740.6313890.5606410.588870.673859
Mean0.3536850.441140.5880440.5117860.5802120.5042520.544140.555818
Min0.3058380.2176710.538060.4179650.5545520.4389450.5011410.478485
STD0.057430.228880.0472350.0840940.044320.0613370.088850.103854
Ranking87152643
5Max0.5774960.5005920.5354510.6860320.6850140.4555190.552410.606479
Mean0.4778990.4612720.4728330.6046130.4707270.4064420.435250.483806
Min0.4019350.3900550.3540470.5451240.2904760.3174590.401250.387526
STD0.0901350.0617870.1029220.0729690.1994590.0771980.454520.111837
Ranking36415872
6Max0.7162010.8269430.7511830.8023340.7683440.767270.598580.790498
Mean0.6340750.6742390.5743380.6415460.6790480.5751390.565550.736605
Min0.5413540.5227760.2792560.3940270.6088910.3881210.525550.669011
STD0.0879040.1520870.2572230.2175320.0814310.1896260.044140.061891
Ranking53742681
Summation322918201333278
Mean Rank6.45.83.642.66.65.41.6
Final Ranking76342851
Table 4. The PSNR results of the test case 2.
Table 4. The PSNR results of the test case 2.
ThresholdMetricComparative Methods
AOWOASSAAOAPSOMPADEDAOA
2Max13.7749112.6388114.529710.5123312.6412413.9921413.2514513.27729
Mean11.9828111.581712.1130510.2569211.5477513.5918912.222112.59842
Min10.7361210.4049510.159989.82336910.9514913.2473112.0121112.09539
STD1.5911211.1217282.2214430.3774460.9482810.3755220.252120.610259
Ranking56487132
3Max15.8638816.4411315.8297514.4764414.8727514.3549315.3252116.8866
Mean14.9466415.470315.4983212.835213.8841113.2251214.1419114.508
Min13.1577614.6411214.962311.645713.2492511.8884413.9547812.02299
STD1.5493810.9083280.4685191.4684450.8676451.246192.251412.433551
Ranking32186754
4Max16.5380113.2008517.1271816.2801217.2772816.8410316.547416.96062
Mean15.6968512.2993315.2883915.5305514.8973415.6625115.2514516.84254
Min14.1484211.101313.3489714.7700912.9299314.7832214.2511416.71792
STD1.3426541.0807251.8911160.7550742.2028421.061040.254960.121484
Ranking28547361
5Max17.9479916.7761216.5875917.7246817.6607815.724618.2564120.50293
Mean17.333216.1448615.6930815.8153116.9731615.1332917.1495417.45356
Min16.8560115.4730814.7414712.3377516.2347714.2904616.2541515.26742
STD0.558840.6524610.9243853.0164860.7143590.7494262.333652.722404
Ranking25764831
6Max18.4164116.7333317.6332317.8641716.9889617.6920117.5484520.23421
Mean18.3513516.2680115.5621216.3365516.4010515.760916.3665219.55858
Min18.2381415.4654313.5708514.1392715.275312.9092914.9585418.29391
STD0.0984060.6979992.0323671.9506560.9752542.520731.369451.096094
Ranking26853741
Summation142725312726219
Mean Rank2.805.405.006.205.405.204.201.80
Final Ranking26486531
Table 5. The SSIM results of the test case 2.
Table 5. The SSIM results of the test case 2.
ThresholdMetricComparative Methods
AOWOASSAAOAPSOMPADEDAOA
2Max0.4062680.3727860.4883130.0987760.3810660.4166080.356520.388003
Mean0.2471710.2289510.264020.0625890.2606230.4061560.325410.299891
Min0.0555040.1190710.0903680.0283410.1561910.3972450.323360.153948
STD0.1776360.1302210.2037480.0352570.1132890.0097730.454850.127294
Ranking67485123
3Max0.5425020.5919260.622080.4691520.5562730.4296740.552410.577128
Mean0.485680.5157040.3858010.2803540.4077580.3222740.512540.543097
Min0.3822380.4633660.1169140.1330110.3199450.1165970.465240.479899
STD0.0897290.0675260.2541570.1718620.129330.178180.514250.054785
Ranking42685731
4Max0.6435060.3723270.6494260.5690080.6167090.6302610.536520.609209
Mean0.5406820.2514430.5156420.5362430.4697590.5336840.514140.591323
Min0.3439840.0501870.3731230.4848710.3454580.4839080.465850.565499
STD0.1704050.1754660.1383590.0450490.1370360.0836510.258540.022911
Ranking28537461
5Max0.5682730.5351960.5632170.6133460.6119660.8101830.584750.647781
Mean0.530320.5278670.5297250.4880840.5787730.652710.545410.618783
Min0.4551180.5230470.4670380.2395540.5563870.5078660.512450.604212
STD0.0651280.0064510.0543310.2152350.0293230.1515530.254140.025113
Ranking57683142
6Max0.7289850.5896240.7917030.6926760.550240.646450.425410.666054
Mean0.6880550.5433630.7510290.585070.5050910.450920.425450.490001
Min0.6581360.4980460.6893680.44740.4271970.1832860.401210.306101
STD0.0366860.0457970.0542990.1253710.0677420.2398530.154240.180105
Ranking24135786
Summation1928223025202313
Mean Rank3.85.64.46544.62.6
Final Ranking27486351
Table 6. The PSNR results of the test case 3.
Table 6. The PSNR results of the test case 3.
ThresholdMetricComparative Methods
AOWOASSAAOAPSOMPADEDAOA
2Max16.4385115.5232916.1717813.3402610.7041713.8015611.2545417.02237
Mean11.7733713.6269811.425810.909438.94024311.9769712.5456213.0293
Min6.91949910.526218.8299518.3621727.3223410.2536311.6585610.1937
STD4.7623112.7075364.1161752.4910851.6956361.7761322.665253.558437
Ranking51678432
3Max16.3221814.3762919.880917.0816718.0169516.3207817.5454818.7061
Mean15.5293313.160918.0833215.1535414.4081212.5975815.3652514.62128
Min14.2958111.7374817.0678312.6157811.105788.57687713.5254111.41963
STD1.0826791.3316491.5611182.2945093.4657683.8805143.254143.722657
Ranking27146835
4Max17.0970217.5533317.684719.7312719.4528119.1968118.5695820.754
Mean16.2391714.1666516.124918.2865115.0284115.1775116.5256518.02446
Min14.7375510.4324313.8853916.014189.87408512.6280415.9684113.89724
STD1.304843.5731461.9887671.9919444.8309013.5224252.597163.635777
Ranking48517632
5Max20.4729520.0870220.2156119.4276517.860820.664618.4937120.41022
Mean17.970217.7664718.0440716.9488817.7417816.3960616.4674318.58232
Min15.0663414.6209715.0998614.7863217.5716713.974415.4554716.23457
STD2.7255342.8248632.6439522.3367680.1511863.7078172.654782.135811
Ranking34265871
6Max21.1634121.3918521.7095520.9359323.0905819.1969420.1215421.98094
Mean19.1237419.8505816.8450417.5693616.887616.9517718.1541420.41186
Min16.520718.8202712.8951613.9138912.5185815.3371116.3698719.1292
STD2.3720711.3598064.4778113.5199215.5194552.0056751.648561.447282
Ranking32857641
Summation1722222333322011
Mean Rank3.404.404.404.606.606.404.002.20
Final Ranking24468731
Table 7. The SSIM results of the test case 3.
Table 7. The SSIM results of the test case 3.
ThresholdMetricComparative Methods
AOWOASSAAOAPSOMPADEDAOA
2Max0.8107970.7842230.7377010.8079840.7298050.7849640.745480.810173
Mean0.7013460.7265860.6792750.6525620.6525030.7483630.712540.777743
Min0.510970.6143910.6426330.5124590.5448590.7074090.695840.738525
STD0.1654860.0971760.0511420.1483570.0961330.038960.025140.036303
Ranking53678241
3Max0.8587150.849550.8790720.8699360.8464740.8109390.8014540.862403
Mean0.8291130.821170.8465870.8107680.8035860.7148670.741250.826664
Min0.7949890.7990880.8119910.75670.722530.6563830.702150.798218
STD0.0321030.0258140.0335910.056790.0702370.0838550.021930.032708
Ranking24156873
4Max0.8357990.7862650.8891580.896390.8316340.8774770.812560.842558
Mean0.8165230.7652470.8351750.8637570.8023870.8569420.765850.820005
Min0.7801920.7485030.7828240.8074970.7519950.8240780.713690.776304
STD0.0314830.0192410.0531860.048930.0438280.0287550.214540.037853
Ranking58316274
5Max0.8698990.8628510.8644410.8953980.8890660.866250.856450.901137
Mean0.8468350.8583030.8210750.8631230.8721820.8625270.810210.874055
Min0.8300220.8504940.7882810.8364860.8547540.8560390.756450.837004
STD0.020660.0067930.0391650.0298580.0171620.0056390.0211140.033208
Ranking65732481
6Max0.8970950.8824730.9107290.9202570.8984520.8847570.841450.890175
Mean0.8939140.8688490.8690010.8504370.8257760.871710.795680.874336
Min0.8921610.8489290.8453060.7876370.7588870.8483420.765820.855048
STD0.0027590.0176360.0362480.0665880.0699620.0202830.0294470.017816
Ranking15467382
Summation1925212229193411
Mean Rank3.854.24.45.83.86.82.2
Final Ranking26457281
Table 8. The PSNR results of the test case 4.
Table 8. The PSNR results of the test case 4.
ThresholdMetricComparative Methods
AOWOASSAAOAPSOMPADEDAOA
2Max0.8107970.7842230.7377010.8079840.7298050.7849640.745480.810173
Mean0.7013460.7265860.6792750.6525620.6525030.7483630.712540.777743
Min0.510970.6143910.6426330.5124590.5448590.7074090.695840.738525
STD0.1654860.0971760.0511420.1483570.0961330.038960.025140.036303
Ranking53678241
3Max0.8587150.849550.8790720.8699360.8464740.8109390.8014540.862403
Mean0.8291130.821170.8465870.8107680.8035860.7148670.741250.826664
Min0.7949890.7990880.8119910.75670.722530.6563830.702150.798218
STD0.0321030.0258140.0335910.056790.0702370.0838550.021930.032708
Ranking24156873
4Max0.8357990.7862650.8891580.896390.8316340.8774770.812560.842558
Mean0.8165230.7652470.8351750.8637570.8023870.8569420.765850.820005
Min0.7801920.7485030.7828240.8074970.7519950.8240780.713690.776304
STD0.0314830.0192410.0531860.048930.0438280.0287550.214540.037853
Ranking58316274
5Max0.8698990.8628510.8644410.8953980.8890660.866250.856450.901137
Mean0.8468350.8583030.8210750.8631230.8721820.8625270.810210.874055
Min0.8300220.8504940.7882810.8364860.8547540.8560390.756450.837004
STD0.020660.0067930.0391650.0298580.0171620.0056390.0211140.033208
Ranking65732481
6Max0.8970950.8824730.9107290.9202570.8984520.8847570.841450.890175
Mean0.8939140.8688490.8690010.8504370.8257760.871710.795680.874336
Min0.8921610.8489290.8453060.7876370.7588870.8483420.765820.855048
STD0.0027590.0176360.0362480.0665880.0699620.0202830.0294470.017816
Ranking15467382
Summation1925212229193411
Mean Rank3.854.24.45.83.86.82.2
Final Ranking26457281
Table 9. The SSIM results of the test case 4.
Table 9. The SSIM results of the test case 4.
ThresholdMetricComparative Methods
AOWOASSAAOAPSOMPADEDAOA
2Max0.5131760.4467390.4693470.3037350.4680330.4726750.4658550.493285
Mean0.4205370.4116450.3798390.2682410.3955910.4079720.4321650.450537
Min0.2973930.348530.3126740.2387320.3539130.2858250.4154410.390067
STD0.1110780.0547730.0806910.0329120.0629740.1058450.0651350.053842
Ranking34786521
3Max0.6399690.5188220.5705670.6428740.5502950.5169350.5416850.638931
Mean0.5392640.4967940.4760320.5590450.5223480.4802530.4516840.565094
Min0.4836930.4535090.4063940.4426660.5071160.4217830.4015130.427189
STD0.0873690.0374880.0848710.1039970.0242350.0511810.1651520.119529
Ranking35724681
4Max0.6355840.5209550.6474370.6290880.5833090.5675890.5454380.59171
Mean0.5588470.4950810.6243640.526970.5421970.5409230.4841530.566602
Min0.4960350.4754770.6123180.4555410.496950.4941020.3515350.530519
STD0.0708090.0233780.0199890.0907530.0433280.0406780.913510.032038
Ranking37164582
5Max0.6254590.7226890.6260880.6787750.6958460.6824120.6468490.752727
Mean0.5741590.655050.6089050.6237350.6816370.5934520.5684350.654905
Min0.5446520.5763030.5890480.5633860.6614640.4689450.5154640.571449
STD0.0445940.0738230.0186640.0578780.017950.1110840.513540.091489
Ranking72541683
6Max0.6786990.7610360.6550620.6938750.749220.7216030.7115430.765842
Mean0.6555020.6525310.5583920.5772770.7355730.6277110.6584350.721667
Min0.6133770.5185190.4606510.5011790.7273220.5039220.6155340.656329
STD0.0365430.1232550.097210.1025340.0119050.1118780.1531120.057742
Ranking45871632
Summation202328271628299
Mean Rank44.65.65.43.25.65.81.8
Final Ranking34652681
Table 10. The PSNR results of the test case 5.
Table 10. The PSNR results of the test case 5.
ThresholdMetricComparative Methods
AOWOASSAAOAPSOMPADEDAOA
2Max15.9398514.485415.9986117.1470515.097716.0583116.5574917.53916
Mean14.1061313.0007215.29414.7914912.7287714.4608714.3598915.82038
Min12.3938711.2729814.7230711.840111.5144713.6266311.975213.42614
STD1.7761071.6199430.6481912.7031782.0517711.3838742.2968712.138085
Ranking67238451
3Max19.4098118.8611517.4826616.5574917.9527917.6913919.4466319.93627
Mean18.7227917.2821316.5211314.3598916.2785717.194119.2687718.13099
Min18.2915815.9513715.9220911.975214.0719416.5155119.0160215.35808
STD0.601411.4706940.8410662.2968711.9944560.6085470.2248572.437657
Ranking24687513
4Max18.1338619.122318.9086117.4574519.104318.2340718.8373518.38066
Mean16.9473916.4032917.9492915.5089318.117616.6717117.2830417.01377
Min15.9849814.8129116.6831113.3773516.420614.9310315.5197615.35199
STD1.0918192.3660251.144042.0462041.4761271.6587241.4941981.535719
Ranking57281634
5Max20.2139820.0240421.2721420.1824518.9149620.688620.2554621.56663
Mean18.9910318.6700420.5372720.0468318.2684519.6935618.0534720.9167
Min17.3681117.8649919.608919.8529617.7980418.9790317.5135420.14592
STD1.4644891.1795710.8483320.1722970.5789070.8886340.154340.718025
Ranking56237481
6Max21.0947621.5713721.343422.1873323.8982419.4466321.5499922.99988
Mean20.4181120.4252920.4782721.6081322.1023319.2687720.2598721.55413
Min19.4980519.0729419.9081321.1634721.0231519.0160219.6485619.98576
STD0.8257161.2619260.7617470.5250191.5658310.2248570.166551.510799
Ranking65421873
Summation2429162424272412
Mean Rank4.805.803.204.804.805.404.802.40
Final Ranking38233731
Table 11. The SSIM results of the test case 5.
Table 11. The SSIM results of the test case 5.
ThresholdMetricComparative Methods
AOWOASSAAOAPSOMPADEDAOA
2Max0.6705590.6631850.6701910.6567070.6389990.6413280.6254540.652777
Mean0.6158080.5772090.6271770.6027370.5882120.6078570.589440.638984
Min0.5544930.4976570.5700640.5375080.5475860.5753530.5235650.622532
STD0.0583110.0829510.0515310.0603920.0465460.0329980.0513510.015297
Ranking38257461
3Max0.7289180.6771490.7148390.6636480.6679380.7170540.6468410.704409
Mean0.7191250.6609350.6665360.596340.6292530.6647450.5518440.671047
Min0.7126560.6293150.62190.5105370.5760460.6228260.4935450.639151
STD0.0086260.0273870.0465780.0782130.0476360.0479650.0503150.032654
Ranking15376482
4Max0.7484350.708220.6874220.7297960.7068170.7292980.7154340.720751
Mean0.6767130.6680070.6755210.6616290.6819850.6748470.6984340.678177
Min0.6335310.6279580.6604190.5900720.6363360.6014960.6513540.614625
STD0.0625430.0401310.0137830.0699240.0395840.0659640.051340.056087
Ranking47582613
5Max0.7603160.7445950.7602280.7222690.7010880.7132940.7154690.754593
Mean0.7352560.7210880.7075330.7065240.6863570.6944030.6946850.740905
Min0.7194290.6937680.6807230.6805770.6659940.6652440.6451350.720782
STD0.0219520.0256270.0456370.0226410.0182120.0256180.1003510.0178
Ranking23458761
6Max0.7786060.8023380.7870560.7592680.7751170.7548720.7143540.759996
Mean0.7456320.7678510.7281880.7434220.7605610.7150350.694560.757846
Min0.7171570.7289730.6852190.7304380.7377180.6823790.6467580.755781
STD0.0309710.0368790.0527470.0146270.020030.0367760.143530.002109
Ranking41652783
Summation1424203025282910
Mean Rank2.84.84655.65.82
Final Ranking24385671
Table 12. The PSNR results of the test case 6.
Table 12. The PSNR results of the test case 6.
ThresholdMetricComparative Methods
AOWOASSAAOAPSOMPADEDAOA
2Max14.2900314.9966613.2307614.258613.6727613.8869913.2566913.62076
Mean12.1312413.567912.0079513.5345212.6933412.5635412.7665812.88169
Min11.0230910.8662510.7115913.0077711.9001810.7597411.8157711.91395
STD1.869792.3409951.2611970.6483420.9008471.6180150.7275730.876084
Ranking71825643
3Max17.1360913.2626515.8629716.7032916.3993716.7727615.5643516.70329
Mean14.8029912.7042114.2781814.7378115.5462815.6824514.5643514.73781
Min11.8852511.8172113.1934913.3031914.9515613.8115213.5453113.30319
STD2.6737920.7767161.403251.7611110.7576961.6276520.355311.761111
Ranking38742164
4Max19.7209416.6540518.5366818.3635218.1821916.817116.1615317.29736
Mean17.4221515.637216.418217.4020316.1703114.6640415.6153316.81609
Min14.4788314.3073214.0852916.5222714.9760413.3174514.6584415.89593
STD2.6798281.2042742.2334270.9233421.7524831.8840580.5131530.797164
Ranking16425873
5Max19.9597118.4127719.4627220.7747519.7192718.0980120.1561521.7643
Mean16.2165317.1675318.1213317.2953317.3021116.2344919.5665218.92589
Min12.1534816.0753317.4447414.3399615.6084115.0877118.9143416.95589
STD3.9129321.1762111.1616913.2492462.1488011.6281140.443452.519088
Ranking86354712
6Max19.7033818.9966619.1575219.9779617.437220.5473120.4861821.16495
Mean17.5559418.5168718.9564618.3894116.7303617.6837317.6716419.65434
Min14.9232317.685218.6203817.4598916.3364115.1449414.5616917.72153
STD2.4267380.7230840.2929231.3823530.6134852.7157893.513551.760107
Ranking73248561
Summation2624241724272413
Mean Rank5.204.804.803.404.805.404.802.60
Final Ranking73323831
Table 13. The SSIM results of the test case 6.
Table 13. The SSIM results of the test case 6.
ThresholdMetricComparative Methods
AOWOASSAAOAPSOMPADEDAOA
2Max0.5851710.5758480.5073340.5842040.5547580.5832740.5454580.566448
Mean0.5297890.5111520.474660.5581410.5101730.5170280.532510.551596
Min0.4514290.4215180.4404680.5286480.4685010.4183390.5105120.540869
STD0.0697690.080130.0334590.0279360.0432020.0871230.0516510.013279
Ranking46817532
3Max0.6336020.5664450.5939760.6431070.6246650.5873160.5845470.620738
Mean0.5732180.530210.5394260.5868450.5612260.5543220.525480.6091
Min0.5316030.4761570.4457840.5572980.514240.5007380.495220.590855
STD0.0535270.0477090.0814640.0487450.0570220.0468170.154790.015999
Ranking37624581
4Max0.6258980.5997590.6886340.6470840.5880070.6192830.601470.651858
Mean0.6204550.5848450.6315810.5841620.5612470.5772850.5645490.631427
Min0.6167390.5618250.5760990.5413420.5209710.5524710.5146250.594978
STD0.0048180.0202240.0562840.0556640.0355030.0365710.515560.031643
Ranking34158672
5Max0.6393090.6115990.6703740.7063680.6624260.6219190.614440.688438
Mean0.5803970.5934070.6260180.6590530.6155590.6011870.532550.657403
Min0.4814860.571930.589530.5960540.5525830.5850430.501440.634943
STD0.0861790.0200380.0409920.0568050.0566660.0188610.2545160.027759
Ranking76314582
6Max0.7327610.7211160.7020310.6394470.6694560.7412370.621560.670238
Mean0.6170620.6758430.660370.6146230.6230980.6980530.601560.621129
Min0.5004420.6281450.6154040.5822160.5934190.6324550.5816850.555772
STD0.1161620.0465330.0434080.0293590.040670.0577510.0516170.058938
Ranking62374185
Summation2325211627223412
Mean Rank4.654.23.25.44.46.82.4
Final Ranking56327481
Table 14. The PSNR results of the test case 7.
Table 14. The PSNR results of the test case 7.
ThresholdMetricComparative Methods
AOWOASSAAOAPSOMPADEDAOA
2Max13.2051514.8154113.4664315.4953115.6260315.2009910.5616511.17017
Mean12.4186811.7862113.2745912.5565813.5268813.848878.684689.014158
Min11.329619.31536313.006269.4429639.67579212.682717.516546.574355
STD0.9736972.7922110.2394343.0299613.3396571.2693972.265582.311011
Ranking56342187
3Max16.5669517.5788616.290816.5806716.5362615.8871916.9897118.01761
Mean15.6999415.4668814.5061514.1924114.8415415.1714114.1555616.2308
Min15.153314.2529713.4516110.4479912.2482214.0316910.4416213.04052
STD0.7593281.8358351.5540453.2835712.2808850.9977553.2550332.769505
Ranking23675481
4Max16.9096317.9775217.7810517.4368916.9491819.3431717.4554319.60255
Mean15.0881614.6929216.4521414.221112.888416.9938514.2266917.64468
Min11.4457811.3766614.6807711.048978.25561113.2974611.426516.42001
STD3.1543973.3005441.5967973.1941844.3749213.2401573.1437371.71328
Ranking45378261
5Max19.9565119.1138317.6000720.1124319.9177620.233319.1148220.05669
Mean17.8966615.9625816.5868918.855317.8866718.6673615.9601818.17262
Min16.7525510.4190915.404716.9134315.0505817.3494710.4294416.15425
STD1.787534.8157841.10741.7058892.5314721.4578314.397281.954685
Ranking47615283
6Max21.435420.2960621.0597722.6183520.0363521.4768420.2433821.70695
Mean19.439917.2588619.5810219.8962519.3783818.8426917.3544520.29721
Min17.785114.2829917.6144217.4866918.9935915.2411814.2324919.17387
STD1.8488423.0070071.7737252.5800680.5725433.228423.2342341.290597
Ranking48325671
Summation1929212125153713
Mean Rank3.805.804.204.205.003.007.402.60
Final Ranking37446281
Table 15. The SSIM results of the test case 7.
Table 15. The SSIM results of the test case 7.
ThresholdMetricComparative Methods
AOWOASSAAOAPSOMPADEDAOA
2Max0.732780.7122880.5601220.725050.7333050.7235760.5648650.730313
Mean0.6436060.6177960.4857290.6457120.6760130.6834730.4486420.710998
Min0.5677980.4789260.340290.5190340.5670550.6396780.3486490.691961
STD0.08330.1228480.1259650.1108670.0944020.0420710.1844760.019178
Ranking56743281
3Max0.7733420.7380510.6932110.7299150.7471020.7063630.6948640.760805
Mean0.7706040.7211440.678170.650640.6778950.7002460.6764630.725786
Min0.7682110.6986010.6651590.4955050.5939470.6930750.6664560.672091
STD0.0025820.020320.0141360.1343620.0776340.0067070.0177430.047213
Ranking13586472
4Max0.7421820.748470.7278270.7325080.7822840.7461160.7288070.790683
Mean0.6872770.6646430.7012160.6692050.7392910.5794210.6527130.763686
Min0.6101960.5498270.6590580.6326160.7014940.3949560.629960.719773
STD0.0687310.1028830.0369270.0550440.0406450.1762530.0507620.038362
Ranking46352871
5Max0.790070.8056450.8226520.8192850.7485880.7404070.7915010.798539
Mean0.763360.6997540.7722850.7804820.7145240.7284730.7744410.775382
Min0.7333480.5088250.7362180.7513910.654350.7078170.7324240.760992
STD0.0285050.1656780.0449560.0349730.0522640.017960.0360230.020251
Ranking58417632
6Max0.7674820.8070740.7691070.7846270.8035150.7974390.7626090.800054
Mean0.7482830.776610.7482590.7688660.7911770.7669870.7452710.762369
Min0.7138170.7237760.7109880.7515430.7689490.7429030.7187950.741915
STD0.0299130.0459340.0323530.0165970.0192880.027820.0303720.032675
Ranking62731485
Summation2125262119243311
Mean Rank4.255.24.23.84.86.62.2
Final Ranking36732581
Table 16. The PSNR results of the test case 8.
Table 16. The PSNR results of the test case 8.
ThresholdMetricComparative Methods
AOWOASSAAOAPSOMPADEDAOA
2Max13.3350312.7041813.7258913.8754115.017913.2251413.9346314.2634
Mean11.8700111.3538211.9045512.3002412.078811.7867911.1583613.26591
Min9.6276959.50403510.8572811.29669.7384959.93056410.8536411.58005
STD1.9721111.6574971.5832321.3811422.6901541.6865711.5476911.468158
Ranking57423681
3Max17.2862516.5708113.4927814.9270517.5435215.1341516.0904916.35347
Mean14.7453915.0154412.5053314.4846314.2193214.4329915.0403515.21107
Min9.82240912.0590511.3041714.1886312.354213.6257512.0571313.80324
STD4.2641672.5614711.1098490.3903752.8859870.7597732.5548521.295669
Ranking43857621
4Max17.0062217.036316.1863117.2699815.6695717.8916115.6081918.921
Mean14.9835812.9410414.4914115.4164913.3105714.8713113.3988715.69837
Min13.051138.26819713.0519413.978068.73806911.048778.87934713.14273
STD1.9790874.4124941.5827171.6847543.960573.4912483.7032732.946311
Ranking38527461
5Max19.3365217.7490519.6840819.6673220.930917.5095717.8340719.59432
Mean17.5197216.4675214.9974515.6861217.0327115.3883517.1549717.61533
Min16.5611514.438269.1690888.17969614.1345812.6395315.9855815.74747
STD1.5741891.7775385.3496676.5047493.5067932.4949391.298081.925829
Ranking25864731
6Max22.3520917.9620820.4587519.7350819.5629219.2291119.9104921.32792
Mean19.1804717.2823615.3681617.4752217.4928816.8596615.7510518.93027
Min14.919815.9689512.6047614.8289713.7224312.61538.96319216.39383
STD3.8339641.1376784.4139742.4757713.2705343.6840266.5857672.46997
Ranking15843672
Summation152833192429266
Mean Rank3.005.606.603.804.805.805.201.20
Final Ranking26834751
p-value2.254 × 10 3 3.455 × 10 2 4.254 × 10 2 2.368 × 10 2 6.589 × 10 2 4.554 × 10 2 3.887 × 10 2 NaN
Wilcoxon sign1111111NaN
Table 17. The SSIM results of the test case 8.
Table 17. The SSIM results of the test case 8.
ThresholdMetricComparative Methods
AOWOASSAAOAPSOMPADEDAOA
2Max0.5760230.4548960.606480.6299780.6257020.5673630.6993810.633192
Mean0.4480620.4160270.472570.5428420.4822230.5018810.4443510.525409
Min0.3136770.3436050.3710430.4249050.337310.4003590.3547330.383819
STD0.1312910.0627760.1210130.105950.1442010.0891450.1426740.128078
Ranking68514372
3Max0.625940.6491270.6060390.6210.6737560.4826750.6530150.670905
Mean0.5214450.55970.5779730.591040.5097750.4438290.5704870.644174
Min0.3281920.417910.5337240.5625740.4116280.3945580.522120.615842
STD0.1675460.1241840.0387790.0292410.1429270.0449740.0203420.027566
Ranking65327841
4Max0.6595530.6890730.6510310.6884220.6191710.6544260.6099380.702364
Mean0.5857890.5772150.5726520.6502460.5388350.5696990.4443510.581725
Min0.4573140.4218030.5244730.6021930.4386350.464780.3733310.415275
STD0.1116690.1388560.0684710.0439550.0918930.0964220.1444270.148926
Ranking24517683
5Max0.711790.6531030.6796730.7285780.7390260.6548660.7106140.723238
Mean0.6678340.6080330.6281760.6129950.6259720.5692950.652340.680701
Min0.6453180.5183080.5294370.4052710.5350250.492560.5432640.618906
STD0.0380710.0777050.0855370.1802760.1037820.0815130.0201350.054767
Ranking27465831
6Max0.7498970.7586640.7004910.7399960.7330310.6740320.6212320.726155
Mean0.7001480.7024080.6251820.6755970.6817830.5903690.603270.683228
Min0.6146340.6662460.5675250.5859310.606070.439360.5133410.620501
STD0.0743870.0493770.0682180.080080.0669230.1310320.7807990.05554
Ranking21654873
Summation1825231527332910
Mean Rank3.654.635.46.65.82
Final Ranking35426871
p-value3.856 × 10 2 2.669 × 10 2 2.665 × 10 2 2.814 × 10 2 6.665 × 10 2 3.854 × 10 2 3.225 × 10 2 NaN
Wilcoxon sign1111111NaN
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abualigah, L.; Diabat, A.; Sumari, P.; Gandomi, A.H. A Novel Evolutionary Arithmetic Optimization Algorithm for Multilevel Thresholding Segmentation of COVID-19 CT Images. Processes 2021, 9, 1155. https://doi.org/10.3390/pr9071155

AMA Style

Abualigah L, Diabat A, Sumari P, Gandomi AH. A Novel Evolutionary Arithmetic Optimization Algorithm for Multilevel Thresholding Segmentation of COVID-19 CT Images. Processes. 2021; 9(7):1155. https://doi.org/10.3390/pr9071155

Chicago/Turabian Style

Abualigah, Laith, Ali Diabat, Putra Sumari, and Amir H. Gandomi. 2021. "A Novel Evolutionary Arithmetic Optimization Algorithm for Multilevel Thresholding Segmentation of COVID-19 CT Images" Processes 9, no. 7: 1155. https://doi.org/10.3390/pr9071155

APA Style

Abualigah, L., Diabat, A., Sumari, P., & Gandomi, A. H. (2021). A Novel Evolutionary Arithmetic Optimization Algorithm for Multilevel Thresholding Segmentation of COVID-19 CT Images. Processes, 9(7), 1155. https://doi.org/10.3390/pr9071155

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