Next Article in Journal
An Evaluation System for COVID-19 Vaccine Transportation Quality Based on Fuzzy Analytic Hierarchy Process
Previous Article in Journal
Effluent Quality-Aware Event-Triggered Model Predictive Control for Wastewater Treatment Plants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Levy Flight and Chaos Theory-Based Gravitational Search Algorithm for Image Segmentation

National Institute of Technology Warangal, Warangal 506004, Telangana, India
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(18), 3913; https://doi.org/10.3390/math11183913
Submission received: 22 June 2023 / Revised: 26 July 2023 / Accepted: 1 August 2023 / Published: 14 September 2023
(This article belongs to the Section Mathematics and Computer Science)

Abstract

:
Image segmentation is one of the pivotal steps in image processing due to its enormous application potential in medical image analysis, data mining, and pattern recognition. In fact, image segmentation is the process of splitting an image into multiple parts in order to provide detailed information on different aspects of the image. Traditional image segmentation techniques suffer from local minima and premature convergence issues when exploring complex search spaces. Additionally, these techniques also take considerable runtime to find the optimal pixels as the threshold levels are increased. Therefore, in order to overcome the computational overhead and convergence problems of the multilevel thresholding process, a robust optimizer, namely the Levy flight and Chaos theory-based Gravitational Search Algorithm (LCGSA), is employed to perform the segmentation of the COVID-19 chest CT scan images. In LCGSA, exploration is carried out by Levy flight, while chaotic maps guarantee the exploitation of the search space. Meanwhile, Kapur’s entropy method is utilized for segmenting the image into various regions based on the pixel intensity values. To investigate the segmentation performance of ten chaotic versions of LCGSA, firstly, several benchmark images from the USC-SIPI database are considered for the numerical analysis. Secondly, the applicability of LCGSA for solving real-world image processing problems is examined by using various COVID-19 chest CT scan imaging datasets from the Kaggle database. Further, an ablation study is carried out on different chest CT scan images by considering ground truth images. Moreover, various qualitative and quantitative metrics are used for the performance evaluation. The overall analysis of the experimental results indicated the efficient performance of LCGSA over other peer algorithms in terms of taking less computational time and providing optimal values for image quality metrics.

Graphical Abstract

1. Introduction

The process of optimization involves locating the ideal solution within a challenging search space. Numerous engineering issues in various disciplines can be categorized as optimization challenges. For instance, computer engineers use image processing techniques to produce the best pixel images, civil engineers design structures that are both affordable and durable, mechanical engineers design machines that have high component life or low production costs, and chemical engineers design the processing facility that ensures the best production rate, while electrical engineers create communication networks that provide the shortest time possible for communication from one node to another, and industrial engineers create systems that minimize the overall turnaround time. Therefore, the best solution for any problem remains a research area that has remained popular for years. Additionally, using algorithms created in this discipline has proven to be a successful way to solve engineering challenges, and consequently, new novel algorithms have been introduced as a result of efforts to optimize more complicated situations.
Meta-heuristic algorithms, which are utilized to address optimization problems, are frequently motivated by simple notions such as physical facts, animal behavior, and evolutionary patterns. In comparison to traditional optimization techniques, they can find the optimum in a wide variety of problem types because of their simplicity, non-derivative nature, flexibility, and avoidance of local optima. However, there is no meta-heuristic method that gives the optimum result for every problem type [1]. Some metaheuristics may work best for one type of problem but may be inadequate for another. The fact that the most successful algorithm is different for each problem type allows for the creation of new meta-heuristic algorithms or the improvement of existing ones. In the literature, it is possible to come across studies on improved versions of optimization algorithms [2], hybrid algorithms created by combining conventional ones [3,4], or newly created algorithms [5,6].
The heuristic algorithms can generally be examined under two basic headings: single solution-based and population-based. A candidate solution is used to begin the single solution-based procedure. Iterations are used to generate this solution [7]. On the other hand, in population-based meta-heuristics, the search begins with a set of solutions (population) and is developed through iterations [8,9,10,11,12,13,14,15,16,17]. Moreover, population-based meta-heuristics have a full comprehension of the solution area. Furthermore, they have the benefit of being able to avoid local optima more easily than single solution-based techniques. Normally, population-based metaheuristics can be grouped as Swarm Intelligence [18,19,20,21], Physics-based [22,23,24], Mathematics-based [5,25], Evolution-based, etc. algorithms [26,27,28].
It has been seen that there are few vaccinations and medicines produced to prevent or treat COVID-19 [29,30,31,32,33,34,35]. Yet, usually, a new virus strand is added to the ever-growing list of COVID-19 variants [36]. The vaccine’s efficiency for this modified virus strand is uncertain [37]. Accurate and early diagnosis of COVID-19 cases is crucial in the early phases of medical treatment and prevention, as infected patients who are not correctly identified continue to disseminate the virus to healthy individuals, allowing the pandemic to spread at an uncontrollable rate. As a result of these factors, physicians are in severe need of an early diagnosis method for COVID-19 [38]. Reverse Transcription Polymerase Chain Reaction (RT-PCR) is a widespread but flawed method of detecting COVID-19 infections. It has a relatively low positive rate and may fail to detect the virus. In addition, the procedure is costly and time-consuming [39]. Medical imaging techniques (X-ray, CT) have an accelerated detection rate and can show patients affected by the virus within hours of infection. Therefore, medical imaging techniques are much more reliable than RT-PCR in the early diagnosis of COVID-19 [40,41]. Automatic recognition of COVID-19 disease utilizing chest X-ray and CT-scan images will aid in reducing the pandemic’s influence on human society. However, these scans have many similar imaging features that make distinguishing COVID-19 from other kinds of pneumonia difficult. Moreover, because of the uneven forms, varying sizes, and indistinct borders between normal and infected tissues, reliably segmenting COVID-19-infected lesions on CT scans remains a difficult process [38]. This challenge may be overcome by employing a variety of image-processing approaches to differentiate such characteristics in terms of similarities and differences.
Image segmentation is one of the best techniques that can be used in medical imaging because it provides detailed information that can be vital to the success of any medical procedure. This is especially true when it comes to MRI and CT scans. The area of concern can be determined during the examination using image segmentation, which physicians can use to diagnose the illness [42,43,44]. In fact, image segmentation is the process of splitting an image into multiple parts in order to provide detailed information on different aspects of the image [45]. This is performed by separating a larger image into smaller, more manageable sections and then identifying the specific characteristics in each section. The particular section represents a different feature or class. Image segmentation algorithms are used in a wide range of areas, such as computer vision, machine learning, pattern recognition, medical imaging, and many more [46,47,48]. They provide an invaluable tool for simplifying complex tasks and enhancing efficiency. As such, they make it possible for researchers and students alike to perform analysis on large amounts of data while offering much greater accuracy than manual rectifications would ever provide. In computer vision and machine learning, image segmentation is a technique that extracts characteristics from images and uses them for data analysis. Medical image segmentation is often applied to MRI and CT scans to identify the part of the body that needs treatment or attention [49,50]. As a result, healthcare providers can focus on treating the area where there is a problem. It can also be used as a way of organizing information since the binary images created can be chosen and classified dynamically by sorting them based on various criteria like storage space, robustness against noise, or speed of computation [51].
Image segmentation may be accomplished using a number of approaches, such as edge detection, thresholding, etc. [52]. Edge detection is the process of identifying which parts of an image are edges and which areas contain mostly noise or flat surfaces [53]. Its typical task is usually to detect edges in an image without relying on complicated mathematical methods like morphology, which generally require multiple measurements from different points in the image. Moreover, edge detection can be used to identify objects in an image, find the edges of text, distinguish between foreground and background, or extract a contour map [54].
Thresholding, on the other hand, is a set of techniques used for image segmentation when an image is represented by pixels. Thresholding is extensively utilized because of its simplicity and robustness [55]. This technique can be used in image enhancement or filtering as well as geometric segmentation [56]. It usually defines different regions in an image with a simple threshold value, allowing us to separate objects that may be too similar. Unlike other segmentation methods, it does not require any form of complex data processing or image manipulation beforehand. There are two types of thresholding processes: bi-level thresholding and multi-level thresholding. These approaches refer to how many levels the algorithm uses to divide the image into segments, rather than how it calculates the level between segments [57].
Bi-level thresholding is similar to binary thresholding in that after a given intensity value is assigned, all pixels below that level are set to that single intensity value. It creates a binary image with transparent pixels, or white ones, and opaque pixels, or black ones. This means that any pixel that has both “black” and “white” layers is either completely black or completely white. In contrast, multilevel thresholding assigns multiple intensities in a hierarchical manner. It is a more complex form of bi-level thresholding where multiple levels are set throughout the image rather than just one level, as in binary images [58,59,60]. Otsu’s method and Kapur’s scheme are the two well-known bi-level thresholding techniques [61]. These approaches show consistent results while still having more flexibility than other approaches. Both methods use a binary decision rule to decide whether a pixel is above or below a given threshold [62]. But these two techniques differ in their treatment of edge pixels, which makes them a perfect case study of how different approaches can yield distinct results. Otsu’s method maximizes class variance, whereas Kapur’s scheme maximizes histogram entropy [63]. While they may perform well in simple structures or with low threshold numbers, as the number of thresholds increases, so does the processing cost. In other words, computational overhead is the main problem with traditional thresholding approaches. Heuristic algorithms, on the other hand, with their simplicity and high convergence speed, can minimize computing costs while increasing decision accuracy. They can be successfully applied to graph theory, optimization, computer vision, and machine learning [10,25,64,65,66,67].
The main motivation of our research is to improve the segmentation capability of standard image segmentation techniques like Kapur’s entropy method. It has been seen that traditional image segmentation techniques suffer from local minima and premature convergence issues while exploring the complex pixel search space. Moreover, these techniques also take considerable time to find the optimal pixels as the threshold levels are increased. Therefore, in order to overcome the computational overhead and convergence problems of the segmentation process, we have employed an efficient hybrid optimizer, namely LCGSA. The LCGSA is able to provide widely expected segmentations at a faster speed and with less computing cost. In LCGSA, the parameters are adjusted by utilizing both Levy flight and Chaos theory in order to improve segmentation results. The efficiency of LCGSA is benchmarked using numerous state-of-the-art techniques. The method presented is fast, precise, and reliable, and it can be employed in a complicated background environment. Furthermore, we employed the LCGSA approach to chest CT scans in order to rapidly and efficiently evaluate the severity of COVID-19 disease. It is expected that by using this novel technique, many more people will be able to receive treatment on time, as the patients most likely to suffer from COVID-19 pneumonia will be diagnosed before clinical symptoms appear.
The following are the primary contributions of this paper:
  • A novel hybrid image segmentation technique, namely LCGSA, is developed to overcome the inadequacies of traditional segmentation approaches and provide predicted segmented output at a faster speed and reduced processing cost.
  • To enhance segmentation results, the algorithm parameters are updated using Levy’s flight and Chaos theory.
  • The algorithm incorporates the Levy flight to enhance exploration capabilities and obtain a suitable balance between the exploration and exploitation stages.
  • Chaos theory prevents the algorithm from getting trapped in local optima and, hence, increases the chances of locating feasible regions of the search space.
  • The proposed LCGSA approach is applied to two benchmark images from the USC-SIPI database.
  • Moreover, LCGSA is also applied to three chest CT scan images in order to quickly and efficiently assess the severity of COVID-19 disease.
  • An ablation study is carried out on COVID-19 images and infection masks to further authenticate the optimal performance of LCGSA.
  • LCGSA’s performance is evaluated and compared with 12 state-of-the-art heuristic algorithms.
The other parts of this paper are structured as follows: Section 2 deals with the literature survey of heuristic approaches for multilevel thresholding. Section 3 covers the methodology related to GSA, chaotic maps, and Levy flight. Section 4 explains the LCGSA and its application in image segmentation. Subsequently, the experimental results of the benchmark and CT scan images are discussed in Section 5. Moreover, Section 6 presents the ablation study in which COVID-19 and ground truth images are analyzed for segmentation purposes. Section 7 illustrates the overall analysis of the experimental results. Lastly, the conclusion and future scope of this study are presented in Section 8.

2. Literature Survey

In the recent period when optimization studies have increased, heuristic algorithms have been widely employed for image segmentation method development. Segmentation involves breaking down an image into clusters of meaningful, non-overlapping, homogeneous parts. Kandhway et al. [68] carried out the segmentation of standard images by using a recently proposed heuristic technique, namely the water cycle algorithm. Meanwhile, in order to find the optimal pixels in the complex pixel search space, they employed the Masi and Tsallis methods. The experimental results confirmed the efficient performance of the Masi entropy-based water cycle algorithm. Similarly, Reptile Search Optimizer (RSO), Arithmetic Optimization Algorithm (AOA), and Aquilia Optimizer (AO) have also been employed for the multilevel thresholding of the benchmark and medical images. The authors utilized the K-means algorithm, Kapur’s entropy, and Otsu’s variance schemes to find the best pixels in the problem space. The simulation results showed the optimal performance of AOA, AO, and RSO over other competitive algorithms [69,70,71]. Furthermore, Su et al. [72] suggested a multilayer thresholding image segmentation approach that was based on horizontal and vertical search processes by utilizing an updated Ant Bee Colony (CCABC) algorithm to increase the efficacy of the traditional ABC method. They utilized the above approach to segment COVID-19 X-ray images. The novel technique’s performance was compared with fifteen different algorithms over benchmark functions. They claimed that the proposed strategy produced higher-quality results.
Chakraborty et al. [73] utilized Kapur’s entropy-based fitness function to segment six benchmark images and three different COVID-19 chest X-ray images in a multilevel thresholding strategy to assess the efficacy of the modified WOA (mWOAPR). Based on the findings, they claimed that the suggested method outperforms several metaheuristics such as WOA, Heap-Based Optimizer (HBO), Hunger Games Search (HGS), SMA, and some variant algorithms of WOA.
CLACO, a new ant colony optimization method developed by Liu et al. [43], was created by merging the Cauchy mutation with the greedy Levy mutation. They used the algorithm to segment COVID-19 X-ray images, utilizing Kapur’s entropy as a fitness function, and found that it performed better than other approaches. In addition, they stated that CLACO outperformed other algorithms considering 30 benchmark functions in terms of search capability and convergence speed.
Singh et al. [38] presented FFQOAK, a novel image segmentation approach built on the K-means clustering method and the Fast Forward Quantum Optimization Algorithm (FFQOA). They aimed to segment CT scan images of the chest in order to properly detect infected regions. They evaluated the suggested strategy on multiple chest CT scan images of COVID-19 patients using multiple comparative image segmentation techniques and concluded that the FFQOAK technique outperformed the others based on various performance assessment parameters.
Zhang et al. [74] designed GBSFSSSA, a new segmentation approach that combines SSA with Gaussian barebone and stochastic fractal searches. They used it for the image segmentation of COVID-19 CT scans and compared it with other techniques on several benchmark problems. They stated that GBSFSSSA is a more trustworthy and efficient approach than other methods after assessing the outcomes using three distinct metrics: PSNR, SSIM, and FSIM.
Houssein et al. [75] presented an enhanced form of the Equilibrium Optimizer that blends conventional operators with dimension learning hunting (I-EO). They evaluated the technique against a set of benchmark functions. They contended that the findings validate the suggested algorithm’s resilience when compared to other optimization approaches. In addition, they employed I-EO for the segmentation of a set of COVID-19 CT images using multi-level thresholding. It was stated that the suggested technique is effective for image segmentation.
Zhao et al. [76] aimed to develop an automated method for segmenting lung CT images by assessing and comprehending the tissue properties of the segmented regions and investigating clinically interpretable information. The aim of this study was to help radiologists diagnose COVID-19 disease. For this purpose, they proposed a new method (SP-V-Net) that integrates a three-dimensional V-Net with Shape Priorities and detects COVID-19 using interpretable characteristics obtained from segmentation findings. They noted that the suggested method performs well and facilitates the automated identification of COVID-19 disease on chest CT images.
Munusamy et al. [39] developed the FractalCovNet model, which employs fractal blocks and U-Net to segment and classify chest CT scan images. They evaluated segmentation results with models such as U-Net, DenseUNet, and other models, as well as classification results from ResNet5-, Xception, and other models. They stated that, when compared to previous techniques, the suggested model can reliably predict COVID-19 infection with high accuracy values.
Jin et al. [77] developed a self-correction model based on field adaptation (DASC-Net) for detecting COVID-19 infection from CT images. DASC-Net is comprised of a new domain adaption model for dealing with domain changes and a self-correcting learning mechanism for improving segmentation outcomes. They claimed that DASC-Net outperformed other coronavirus infection segmentation algorithms in extended trials on three COVID-19 CT datasets.
Nama [78] suggested a new Quasi-Reflected SMA (QRSMA) that merges the SMA with a quasi-reflection-based learning system. QRSMA’s performance was evaluated by comparing it with other algorithms using various benchmarking functions, and it was stated that the results can considerably increase QRSMA’s convergence speed and solution accuracy. Furthermore, the method was applied to COVID-19 X-ray images. The simulation outcomes revealed that QRSMA is a more efficient multilevel thresholding technique than other recent methods. In Table 1, a summary of the related works dealing with heuristic image segmentation approaches is provided. It is evident that Kapur’s entropy scheme is a widely used objective function for image segmentation. It is because HAs treat image segmentation as an optimization problem in which an objective function is necessary for checking the quality of the searcher agents. The survey also depicts that heuristic segmentation methods have high exploitation capability and require less computational overhead to find the best pixels from the segmented image. Moreover, USC-SIPI image database benchmarks like Lena, Cameraman, Aeroplane, Hunter, etc. are often used by researchers for performance evaluation. Furthermore, it can also be seen that MSE, PSNR, FSIM, and SSIM are commonly used performance metrics for multi-level thresholding.

3. Methodology

The LCGSA is a hybrid strategy that has powerful exploration and exploitation capabilities. Actually, LCGSA is a combination of standard GSA, Levy flight distribution, and ten chaotic maps. In this section, the mathematical foundation of the LCGSA technique is laid out.

3.1. Gravitational Search Algorithm (GSA)

Each optimization algorithm draws inspiration from natural, physical, anthropological, or chemical processes. There are also a lot of heuristic techniques that are inspired by nature, such as PSO, GWO, ACO, etc. The gravitational search algorithm (GSA) is one of the heuristic techniques that draws inspiration from physics. In the GSA optimization process, mass initialization is the first stage since searcher agents adopt the shape of masses. It is based on Newton’s rule of global gravity and motion, which states that “the gravitational force between two masses is proportional to their product and inversely proportional to the square of the distance between them”. If we consider a system of N masses, then the position of the ith mass in the search space is given by Equation (1).
X i = x 1,1 x 1 , i x 1 , n 1 x 1 , n x 2,1 x 2 , i x 2 , n x N 1,1 x N 1 , i x N 1 , n x N , 1 x N , i x N , n 1 x N , n
At a specific time ‘t’, we define the force, F i j d acting on mass ‘i’ from mass ‘j’, as follows:
F i j d ( t ) = G ( t )   M p i ( t ) M a j ( t ) R i j t +   ( x i d t + x j d t )
The active and passive attractive masses in Equation (2) are M𝑎j(t) and M𝑝i(t). Actually, M𝑎j(t) denotes the attracting force acting on a point mass. Likewise, M𝑝i(t) indicates the attractive force exerted by a point mass in a gravitational field. In addition, 𝑅ij(𝑡) and ∈ represent Euclidean distance and a small constant, respectively.
Having the correct balance between diversification and intensification stages is vital in GSA. Consequently, ‘G’, a gravitational constant, assists in locating the solution space’s feasible regions. Furthermore, it helps in creating consistency in the solutions during the optimization operation. It is expressed by Equation (3).
G ( t ) = G ( t 0 )   e α C I   M I
where G(t0) and G(t), respectively, are the initial and final values of G, and α signifies a minor coefficient. Additionally, CI and MI represent the current iteration and the maximum number of iterations, respectively. In addition, it is important to calculate the masses existing in the solution area. Aside from the active and passive masses, there is also an inertia mass that quantifies particle resistance to external forces. It is obvious that the bigger the mass, the greater the gravitational attraction. The value of gravitational mass (Mi) is calculated as shown in Equation (6), when the values of active (mai), passive (mpi), and inertial (mii) masses are equal as indicated in Equation (4).
Mai = Mpi = Mii = Mi i = 1, 2, 3,…, N
m i ( t ) = f i t i ( t )     w o r s t ( t )   b e s t ( t )     w o r s t ( t )
M i ( t ) = m i ( t ) y = 1 m m i ( t )
The fitness function represented by 𝑚i(𝑡) in Equation (5) measures the quality of the masses in a particular environment, whereas best (t) and worst (t) are parameters of f i t i ( t ) that determine whether the problem is minimization or maximization one. In Equation (6), ‘m’ represents point mass candidate solutions in the search space. Likewise, Newtonian mechanics utilizes Equation (7) to determine the total gravitational force that will be used to find the masses that can attract other masses toward themselves.
F i d ( t )   = j = 1 , j i m γ y F i j d ( t )  
In Equation (7), γy is a random variable. It is also apparent that, according to Equation (2), heavier masses will have a stronger gravitational field. Furthermore, in the whole search domain, feasible neighborhoods will be discovered. As a consequence, the quality of the solutions is maintained by employing the kbest (cardinality constraint) strategy, as demonstrated in Equation (8).
F i d ( t )   = j = k b e s t , j i m γ j F i j d ( t )  
When a physical system undergoes acceleration, it necessarily produces a force. Masses in the solution space constantly exert a force on each other, which causes acceleration and directs solutions toward feasible regions. In Equation (9), F i d t   is the force generated by the masses on each other.
a i d ( t )   = F i d ( t )   m i i ( t )
In contrast, inertial mass is represented by M i i ( t ) in GSA. Each point mass has a position and velocity. Meanwhile, at the end of the iteration process, just one mass with a significant gravitational field remains. Because of this, it is critical to calculate the velocity v i d ( t ) and position x i d ( t )   in order to find an optimal solution, as shown in Equations (10) and (11).
v i d ( t + 1 )   = γ j v i d ( t )   + a i d ( t )  
x i d ( t + 1 ) = x i d ( t ) + v i d ( t + 1 )  

3.2. Levy Flight and Chaos Theory-Based Gravitational Search Algorithm (LCGSA)

This section provides a solid groundwork for the suggested LCGSA technique by explaining Levy flight and Chaos theory concepts. In LCGSA, Levy flight is utilized to guarantee global exploration of the search space and maintain a proper balance between exploration and exploitation phases. Moreover, chaotic maps help in the convergence of candidate solutions towards the global optimum.

3.2.1. Levy Flight

Exploration is crucial for handling complicated and large-scale optimization problems, according to a number of investigations [79,80]. In actuality, optimization and premature convergence problems are brought on by a decrease in the diversity of candidate solutions. Additionally, GSA encounters difficulties with exploration while attempting to solve multidimensional and challenging multimodal problems. Therefore, Levy flight distribution has been combined with the gravitational constant of GSA to solve diversity difficulties in HAs and ease the exploration problem of GSA [81,82]. In general, numerous studies have shown sufficient experimental support for the use of Levy flight to address diversity issues in HAs [81,82].
A random walking strategy called Levy flight [83] is based on the probabilistic distribution of position changes that take place during the motions of living objects. The magnitude of the step made by the moving element varies in dynamical systems. The Levy distribution, which is based on the Fourier transform, as indicated in Equation (12), determines this variability.
F ( k ) = e x p [ α | k | β ] ,   0 < β     2
where α is a scaling coefficient, β is the Levy index, and k is the characteristic equation variable.
I f   β = 2 ,   t h e n   F ( k ) = e x p [ α k 2 ]  
whose inverse Fourier is consistent with a Gaussian distribution.
O r   i f   β = 1 ,   t h e n   F ( k ) = e x p [ α | k | ]
which is consistent with a Cauchy distribution [84].
Brownian random walks are less effective than Levy flights [85]. Levy-like flights or movements have been observed in a variety of animals, including insects, monkeys, and more. In addition, there are other physical processes that, under the right circumstances, exhibit Levy-flight behavior, such as the spreading of fluorescent molecules. The reason for this is that Levy flights can maximize the effectiveness of resource discovery in unstable environments [86,87].
Levy flight is used in the proposed LCGSA technique to offer stability between diversification and intensification, which solves the problem of local minima. Additionally, the infinite variation of the Levy distribution ensures that the GSA’s difficulties with sensitive initialization and slipping into local minima will be resolved. In fact, Levy flight with big step sizes enhances the feasible solution diversity, while with small step sizes, there is a strong convergence of the solutions to the optimal neighborhood.

3.2.2. Chaos Theory

Most metaheuristic methods make significant use of long-period random number sequences. The likelihood of the algorithm becoming trapped in local optima may rise by amassing a collection of randomly produced numbers in a particular region or by producing the same values. The numbers generated should not be identical and must have a spread spectrum in order to address these limitations [88].
Chaos-based methods are based on a class of functions called chaotic maps. Discretized-time systems with chaotic behavior are known as chaotic maps. It has been demonstrated that chaotic maps generate numbers that are unpredictable and non-periodic in nature. Values derived from chaotic visualizations are used instead of random variables, which are frequently chosen in heuristic algorithms that incorporate them [89,90].
A chaotic sequence is represented by the total of the chaotic variables employed during a particular iteration. The use of chaotic sequences demonstrates the versatility of the algorithm by allowing it to break free of local minima while looking for the global minimum. By doing so, HAs can circumvent inaccessible portions of the search space [91]. Therefore, it is anticipated that creating viable solutions for optimization issues using metaheuristic algorithms with chaotic maps will be quicker and more efficient when the initial random number string is determined [92,93,94,95,96,97,98,99,100,101]. Ten chaotic maps were used in this effort to improve GSA’s performance and fix its optimization issues. The chaotic maps employed in this work are shown in Figure 1, among the numerous described in the literature.

4. Image Segmentation Using LCGSA

In this part, a modified version of the classic GSA that is based on two very intriguing mathematical methods, Levy flight and Chaos theory, is proposed for the image segmentation challenge. The GSA suffers from the disadvantage of skipping real solutions during the optimization process, as well as from delayed convergence and entrapment in local minima issues. Two methods are used to address the aforementioned problems.
The Levy flight distribution is used in the first strategy to address the standard GSA’s diversity problem. Actually, the Levy distribution’s unlimited variance and adjustable step size aid in solving the local optimum problem. In other words, it broadens the diversity of the search process. Equation (15) illustrates the mathematical computation of the Levy flight [83] using Mantegna’s technique.
L e v y   ( u ,   v ) = c u     v 1 β σ u
where c′ is a multiplicative constant having a value of 0.01; u and v are normal distributions; and β is a Levy index with a value of 1.5.
In the second strategy, ten different chaotic maps are employed to overcome slow convergence and the local searching issues of standard GSA. The chaotic maps create huge changes in the output when the initial conditions of the maps are modified. This helps search agents move out of the local minima traps. Moreover, chaotic normalization aids in the proper balance between exploration and exploitation. It is mathematically calculated as shown in Equation (16).
C i n o r m ( t ) = C i t a ( D c ) ( b a ) + c  
In Equation (18), (a, b) is the range of the chaotic map; i represents a chaotic index value from 1 to 10 because ten chaotic maps were considered; and (c, D) is the chaotic normalized interval where c has a value of zero, while D is calculated using Equation (17).
D = M I t   M I ( M a x M i n )  
Here, MI and t represent the maximum number of iterations and the current iteration, respectively. In addition, adaptive intervals are indicated by Max and Min with values of 20 and 1 × 10−10, respectively. Therefore, the chaotic normalization equations of ten chaotic maps can be written as follows:
C 1 n o r m ( t ) = C 1 t a ( D c ) ( b a ) + c ( C h e b y s h e v )  
so on to…
C 10 n o r m ( t ) = C 10 t a ( D c ) ( b a ) + c ( T e n t )
According to Equation (3), the gravitational constant (G) is the primary parameter in standard GSA that governs the gravitational field’s intensity. The right balance between the phases of exploration and exploitation depends on this parameter. In reality, the value of G declines exponentially throughout the preliminary iteration phase, ensuring solution diversification. Additionally, the value of G varies gradually through the final iterations, encouraging the adoption of potential solutions in the direction of the overall optimum. Therefore, G is selected since it is a key governing parameter of the conventional GSA that enables the exploration and exploitation stages to move more quickly. Levy flight and chaotic sequences are integrated with the GSA gravitational constant in the proposed LCGSA. Consequently, the sum of Equations (3), (15) and (16) yields the Levy–Chaotic gravitational constant (GLC(t)).
G L C ( t ) = L e v y   ( u , v ) + C i n o r m ( t ) + G ( t 0 ) e α t   M I
On the basis of the chaotic normalization of the search space, ten chaotic versions of LCGSA are also illustrated below:
G L C 1 ( t ) = L e v y   ( u , v ) + C 1 n o r m ( t ) + G ( t 0 ) e α t   M I   ( C h e b y s h e v V e r s i o n )
G L C 2 ( t ) = L e v y   ( u , v ) + C 2 n o r m ( t ) + G ( t 0 ) e α t   M I ( C i r c l e V e r s i o n )
so on to….
G L C 10 ( t ) = L e v y   ( u , v ) + C 10 n o r m ( t ) + G ( t 0 ) e α t   M I ( T e n t V e r s i o n )
There is no doubt that the Levy–Chaotic Gravitational constant, GLC(t), possesses the intriguing traits of Levy randomness, chaotic stochasticity, and heuristic adaptive learning capacity. In general, GLC(t) possesses all the necessary traits for resolving standard GSA’s entrapment in local minima, intensification, and diversification problems. Figure 2 shows the flowchart of LCGSA algorithm while Figure 3 displays the multilevel image thresholding scheme, which is based on the LCGSA.

5. Experimental Results and Discussion

The simulation analysis of image segmentation based on ten chaotic versions of LCGSA, that is, LCGSA1 to LCGSA10, has been benchmarked using two standard images and three COVID-19 chest CT scan images. The benchmark images, namely Airport and Boat, were taken from the USC-SIPI database, while COVID-19 chest images, namely CT1, CT2, and CT3, were acquired from the Kaggle database. The images have a symmetrical pixel distribution, that is, a pixel range of 0–255. The histogram representation of the standard test images and CT scan images is shown in Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8, respectively. Both traditional and robust HAs were employed for the empirical analysis, such as PSO, GSA, PSOGSA, SCA, SSA, DE, BBO, CPSOGSA, MFO, ABC, GWO, and SMA. The initial values of the competitive algorithms were acquired from the base papers of the algorithms.
The simulation’s results were documented with a constant population size of 20 and a maximum iteration count of 300. Ten trials of the algorithms were completed before recording the experimental data. Additionally, the participating HAs’ optimization process comes to an end when they produce identical results for 10% of the total iterations. The source codes will be available online at https://github.com/SAJADAHMAD1 (accessed on 5 June 2023). Moreover, the simulation analysis was performed using the 3.40 GHz i7 Intel core CPU and the R2013a MATLAB version.
The quality of the output segmented image was evaluated using a variety of performance measures, including Peak Signal-to-Noise Ratio (PSNR), Structural Similarity Index Measure (SSIM), and Feature Similarity Index Measure (FSIM). In reality, PSNR examines the segmented image’s dependability by taking into account the threshold values in subsequent iterations. Equation (18) illustrates how PSNR is expressed mathematically.
PSNR = 10   l o g 10 255 2 M S E
MSE = 1 R C i = 1 R j = 1 C ( I   i , j   O   i , j ) 2
The matrix’s rows and columns are represented by R and C, respectively, in Equation (25). I and O stand for the standard input image and segmented output image, respectively. SSIM, another image quality metric, measures how uniform and similar the segmented image is to the input standard image. Given that it indicates the presence of high-quality pixels in the segmented image, the value of SSIM should be high. Equation (26) illustrates how SSIM is expressed mathematically.
SSIM   ( x ,   y ) = ( 2 μ x μ y + c 1 ) ( 2 σ x y + c 2 ) μ x 2 + μ y 2 + c 1     ( σ x 2 + σ y 2 + c 2 )
Such that μ x   a n d   μ y are the mean intensities of input and segmented images; σ x 2   a n d   σ y 2 represents standard deviation; σ x y is the covariance; and <c1, c2> are the small constants.
FSIM is another image assessment metric that deals with examining the quality of the segmented image by assessing the pixels in the local neighborhood of the image. Meanwhile, FSIM inspects the values of Phase Congruency (PC) and Gradient Magnitude (GM), which are essential for locating the best pixel values to determine the accuracy of the output image. The mathematical formulation of FSIM is shown in Equation (27).
FSIM   ( x ) = x Ω S L ( x ) P C m ( x ) x Ω P C m ( x )
Whereas the domain of the pixel search space is represented by   Ω , and S L (x) is the similitude of the image. If P C 1 and P C 2 are the symmetrical phases of two mathematical functions, then PC is given by Equation (28).
P C m   ( x ) = max   ( P C 1   ( x ) ,   P C 2 ( x ) )
Likewise, for the similitude content of the segmented image, S L (x) is formulated in Equation (29).
S L ( x ) = [ S P C ( x ) ] α [ S G ( x ) ] β
Such that
S P C ( x ) = 2 P C 1 ( x ) P C 2 ( x ) + T 1 P C 1 2   x   P C 2 2   x + T 1
S G ( x ) = 2 G 1 ( x ) G 2 ( x ) + T 2 G 1 2   x   G 2 2   x + T 2
where S G   and S P C are the similitude of Gradient and Phase Congruency while constants are denoted by α , β , T 1 ,   a n d   T 2 . As the image segmentation task is a maximization problem, the higher value of FSIM shows the capability of locating the best pixels in the solution space.

5.1. Experimental Analysis of Benchmark Images

The LCGSA is a hybrid heuristic technique that is used to find the optimal solution in a non-linear search space. It is clear that LCGSA has robust exploration and exploitation operators to handle complex problem spaces. Therefore, in order to authenticate the problem-solving capability of LCGSA, it was applied to two benchmark images, namely the Airport and the Boat. The LCGSA will be mainly tested and benchmarked for its capability of finding efficient values for the image thresholds and providing optimal values for various image quality assessment measures like PSNR, SSIM, and FSIM. Moreover, it will be interesting to see how LCGSA will deal with the computational overhead problem as the values of the image thresholds increase.

5.1.1. Simulation Results of the Airport Image

The simulation outcomes of competitive algorithms for the Airport image are presented in Table 2 and Table 3. It can be clearly seen that all ten versions of LCGSA have better values for the image thresholds at k = 2, 4, 6, 8, and 10. The minimum values of standard deviation (STD) and MSE for LCGSA show the stability of the optimization process and the high quality of the pixels. It can also be observed that LCGSA has acquired large values for the mean as the threshold levels are increased. It indicates that LCGSA was successful in obtaining better pixels in successive generations. Moreover, LCGSA versions, particularly LCGSA2, LCGSA4, and LCGSA8, have large values for PSNR, FSIM, and SSIM, which shows superior quality of the output image and efficient segmentation capability. Furthermore, LCGSA takes less runtime (in seconds) to find the best pixels in the search space. It can also be seen that SCA, BBO, and ABC provide appreciable values for the image quality metrics. In contrast, MFO, ABC, SSA, SMA, and GWO provide sub-optimal results for the statistical measures and image quality metrics. As far as best values for Kapur’s objective function are concerned, it can be seen that SCA, DE, and BBO provide better values, while GSA, PSOGSA, PSO, CPSOGSA, MFO, ABC, SSA, SMA, and GWO provide smaller values, indicating optimization problems while handling uneven optimization landscapes.
The segmented images and histogram curves of LCGSA are shown in Figure 9. It is clear that as the number of thresholds increases, the clarity and contrast of the output segmented image also increase. Further, the convergence curves and box plot graphs are depicted in Figure 10 and Figure 11, respectively. The convergence curves show that LCGSA takes less computational time to find the best pixels in the problem space. Moreover, LCGSA exhibits appreciable exploitation capability as its convergence speed is better than all other peer algorithms. Besides, the convergence curves of GWO and SMA are at the bottom, indicating lower values for the objective function and the presence of outliers in the segmentation output. It also indicates that both SMA and GWO have serious issues while handling the complex pixel problem spaces and are unable to locate the optimal pixels during the optimization process. Likewise, SSA, MFO, ABC, PSOGSA, GWO, and SMA have small values for Kapur’s fitness function, indicating issues in the optimization capability. The box plots also validated the optimal performance of LCGSA, as it has the highest symmetrical fitness values in the range of 50 for the image pixels. It is worthwhile to mention that GWO, SMA, SSA, GSA, MFO, PSOGSA, ABC, PSO, and CPSOGSA give sub-optimal values for the average and inter-quartile range, highlighting difficulties in the exploration.

5.1.2. Simulation Results of the Boat Image

The simulation outcomes for the Boat image are recorded in Table 4 and Table 5. It can be observed that LCGSA versions show smaller values for standard deviation and mean square error, implying coherence in the image pixel values. It can also be noted that standard GSA has minimum values for image quality metrics like PSNR, SSIM, and FSIM, indicating difficulties in handling non-linear problem spaces. Further, GSA has <15.39, 21.06, 26.35, 32.39, 39.77> values for Kapur’s objective function, which are nowhere near the values provided by LCGSA, such as LCGSA1 <18.10, 24.24, 35.75, 44.86, 52.55>. Moreover, when we closely look at the mean values of the LCGSA and peer algorithms, it is clear that the LCGSA has better mean values, and they are also close to Kapur’s objective function value. It shows that LCGSA has lower presence of outliers and noise in the segmented output. Furthermore, LCGSA takes less runtime to locate the optimal pixels, while ABC, BBO, DE, GSA, and SSA take more computational time, indicating optimization issues.
The segmented output and histogram curves of LCGSA at k = 10 are shown in Figure 12. It can be seen that histogram values are scattered across the whole frequency spectrum of the image. It means that LCGSA has optimal values for the thresholds, which help in the efficient segmentation of the image. The improvement in the segmented output can be seen as the number of thresholds increases. It shows that LCGSA is obtaining optimal pixels as the complexity of the segmentation increases. Moreover, convergence curves, as shown in Figure 13, depict that LCGSA has higher local optimization power while GWO, SMA, SSA, BBO, and MFO have a slower convergence rate. It is clear that GWO, SMA, and SSA are facing optimization problems while countering the uneven search spaces and local minima terrains. The box plots, as shown in Figure 14, convey that LCGSA has consistency and proximity in the objective function values. It is because LCGSA has the maximum value for Kapur’s objective function, while SMA, GWO, ABC, SSA, and MFO have smaller values, indicating sub-optimal performance. Moreover, SMA and SSA contain whiskers and outliers in the output, implying that values are away from the central mean, and it also indicates that these techniques have serious exploration problems.

5.2. COVID-19 Case Study: Experimental Analysis of COVID-19 CT Scan Images

COVID-19 is a respiratory viral disease that severely affects the normal functioning of the human lungs. Researchers have applied various heuristic approaches in order to segment the CT images of COVID-19 patients [43,73,102]. Normally, axial non-enhanced chest CT images, also called lung windows, are used for segmentation purposes. Therefore, we have considered three CT images of COVID-19-affected patients from the Kaggle database, namely CT1, CT2, and CT3. The LCGSA versions were applied to chest CT images in order to properly identify the consolidation and ground glass opacity areas of the lung window images. Moreover, we have also used colormap images at different threshold levels to clearly show the segmentation areas in the output image.

5.2.1. Simulation Results of the CT1 Image

The experimental results of the CT1 image are recorded in Table 6 and Table 7. The algorithms that have provided optimal values for the thresholds include LCGSA versions, BBO, DE, ABC, and SCA. It can also be seen that LCGSA, BBO, and DE have suitable values for PSNR, SSIM, and FSIM. Meanwhile, GWO, SMA, PSO, and SSA have large values for the standard deviation and MSE, indicating the presence of outliers in the segmented output. However, ABC, BBO, DE, and SSA take substantial computational time to find the optimal pixels in the problem space. Furthermore, it is obvious that LCGSA versions take less CPU time to segment the image. As far as mean values are concerned, LCGSA4, LCGSA8, and LCGSA9 provide large mean values, indicating these versions have the best values for Kapur’s objective function. When we closely look at the mean and best objective function values of the LCGSA versions, it is clear that the mean values at k = 10 (LCGSA4 (47.70), LCGSA8 (49.19), and LCGSA9 (51.07)) are close to the best objective function values at the same threshold level (LCGSA4 (53.14), LCGSA8 (53.71), and LCGSA9 (53.77)). It shows that LCGSA versions have fewer outliers present in the segmented output because objective function outcomes have fewer value differences across consecutive generations. On the other hand, MFO, SSA, and CPSOGSA have large differences in the mean and Kapur’s objective function values, indicating an improper balance between the exploration and exploitation stages during the optimization process.
The multi-level thresholding output of the CTI image is depicted in the form of segmented images, colormap images, and histograms, as shown in Figure 15. In the grayscale segmented output image, the hazy white part in both lungs is clearly visible, and it becomes clearer as the number of thresholds increases. It shows that LCGSA was successful in finding feasible pixels in the complex search space environment. Moreover, it also implies that LCGSA has reduced the impact of infeasible pixels or noise on the segmented output. The obscured consolidated parts can also be seen in the colormap images, which show the presence of COVID-19 disease in the lungs. The reason behind the optimal performance of LCGSA is its intelligent hybrid framework, in which Chaos theory handles local exploitation issues while variable step size and infinite variance of Levy flight help in the global exploration of the search space. In simpler terms, LCGSA has a robust optimization combination that helps it counter difficult optimization terrains like uneven pixel search space with proper ease and without getting stuck in the local minima regions. Similarly, the local exploitation capability of the peer algorithms is benchmarked through convergence curves, as shown in Figure 16. It can be observed that LCGSA and SCA have fast convergence rates, while SSA, BBO, CPSOGSA, GWO, SMA, and PSO take more runtime to converge towards the optimal regions of the pixel search space, indicating local minima issues. Furthermore, the box plot analysis is depicted in Figure 17, and it portrays the best performance of LCGSA because it has efficient values for Kapur’s objective function. It is also clear that the boxplots of GWO, SMA, CPSOGSA, and SSA are at the bottom, implying local minima and exploration problems.

5.2.2. Simulation Results of the CT2 Image

The CT2 is another CT scan image of the COVID-19 patient that we have considered for benchmarking the segmentation ability of the heuristic algorithms. The experimental analysis is shown in Table 8 and Table 9. As far as best threshold values are concerned, the LCGSA1, LCGSA2, LCGSA10, DE, and SCA have efficient pixel values. At the same time, PSO, PSOGSA, and CPSOGSA depicted large standard deviation values, implying pixels are distant from the central optima region. Meanwhile, SMA, GWO, PSO, and PSOGSA also show substantial outcomes for MSE, indicating the presence of outliers and noise in the segmented output. It can also be noted that LCGSA6 and LCGSA7 perform better than other LCGSA versions because they have better values for PSNR, SSIM, and FSIM. Further, the mean and Kapur’s objective function values of LCGSA versions are close to each other, indicating symmetry in the segmented output and less deviation from the mean objective value. On the other hand, SMA, GWO, PSO, PSOGSA, and CPSOGSA are again showing suboptimal results for the multilevel thresholding problem. In fact, when we closely look at the simulation results of these techniques, we find that they have large values for the STD and MSE but small values for the mean, Kapur’s objective function, PSNR, SSIM, and FSIM.
The colormap images shown in Figure 18 clearly show the existence of ground glass opacity and fibrous bands in the lower part of the lungs. It can also be seen that optimal pixel values cover most of the frequency spectrum of the CT2 image at respective threshold values. It indicates the exploration capacity of LCGSA and its competence in handling local minima regions. Likewise, Figure 19 indicates the efficient performance of LCGSA versions in locating feasible pixels in the complex search space. When we closely look at the convergence curves, it is clear that the SMA and GWO curves are at the bottom, implying premature convergence problems. On the other hand, the LCGSA convergence curve is at the top, portraying optimal fitness values at consecutive iterations. It is also evident that SCA, DE, and PSO show appreciable convergence performance. Furthermore, box plots in Figure 20 also authenticate the outstanding performance of LCGSA because it has a higher concentration of optimal fitness values around the mean. Moreover, the LCGSA has fewer outliers, and its mean fitness values are better than other competitive algorithms. The boxplots also convey that SMA and GWO have very small values for the fitness function, implying serious optimization issues and less segmentation power. In contrast, SCA and DE again showed optimal performance because their fitness values are very close to LCGSA, indicating competitive segmentation capability and symmetry in the segmented output.

5.2.3. Simulation Results of the CT3 Image

The experimental analysis of the CT3 image is presented in Table 10 and Table 11. It is clear that LCGSA, DE, and SCA are the best-performing algorithms, as they have efficient values for the pixels. It is worthwhile to note that ABC, PSO, PSOGSA, and CPSOGSA gave competitive values for the image quality metrics. On the other hand, SSA, MFO, SMA, and GWO again provide suboptimal values for the thresholds. The reason behind the poor performance of the above techniques is that they lack diversity in the optimal threshold values. It shows that they got stuck in the local minima in the pixel search space, which resulted in premature convergence. Moreover, when we further investigate their evaluation metrics values, it is clear that they have large values for STD and MSE and small values for the mean and Kapur’s fitness function. In contrast, the LCGSA versions provided very high values for the mean and fitness functions indicating optimal performance and segmentation proficiency. Meanwhile, the LCGSA versions take very little runtime to reach the global optimum. It indicates that LCGSA was successful in resolving the computational overhead problem of the traditional Kapur’s method while handling high threshold levels. At the same time, when we look at the run times of peer algorithms, GSA, PSOGSA, SSA, ABC, and SCA take more time to find the feasible pixels in the search space. Moreover, DE and BBO take less runtime than the above-mentioned heuristic techniques, indicating appropriate optimization capability.
The air bronchogram is clearly visible in the segmented output images provided by LCGSA, as shown in Figure 21. Evidently, the colormap images also depict the heavy consolidation density in both lungs. The distortion of the lung architecture is obvious in the case of the left lung because it is excessively affected by the Coronavirus infection. Moreover, in the colormap image, the yellowish color shows air-filled bronchi with fluid or pus, while the light blue color represents fibrosis, which is opacification of the alveoli, and the dark blue color indicates healthy parts of the lungs. In Figure 22, the convergence rate of LCGSA is faster than other competitive approaches. At the same time, the convergence curves of SMA, GWO, and SSA are at the bottom, indicating premature convergence issues in handling difficult search spaces. However, the convergence curves of DE, SCA, and MFO are close to the LCGSA curve. It indicates these techniques have also obtained the best values for Kapur’s objective function in consecutive generations. Similarly, Figure 23 conveys that LCGSA has the highest segmentation capability. It is because LCGSA has a maximum value for the fitness function, and its simulation values are concentrated around the mean. The box plot analysis also shows that GWO, SMA, MFO, SSA, and ABC have lower values for Kapur’s objective function, while DE, SCA, and BBO have average values for the fitness function.

5.3. Statistical Analysis of the Results

The minimum value of standard deviation and mean square error does not imply that an algorithm is more efficient than others [103]. In fact, statistical tests should be performed on the simulation results to find the optimal competitive algorithm. Therefore, a pairwise non-parametric signed Wilcoxon rank-sum test [104] was performed at a 5% significance level to statistically validate the simulation results between LCGSA and other peer algorithms. The reason behind selecting a Wilcoxon rank-sum test is that it uses the median as a statistical measure. Moreover, in the Wilcoxon rank-sum test, the distribution of the dataset is not considered.
The null hypothesis (H0) states that the LCGSA does not provide the best pixel values for the image, whereas the alternate hypothesis (H1) includes that the LCGSA provides efficient pixel values for the image quality metrics. The p-values of all the competitive algorithms are computed through MATLAB simulation analysis. If an algorithm provides a p-value less than 0.05, then the null hypothesis is rejected and H1 is accepted. Moreover, if the p-value is equal to one, it simply means the performance of the peer algorithm is consistent with the best-performing algorithm. Table 12 shows the statistical analysis of the five grayscale images in which LCGSA values are statistically compared with the other twelve competitive algorithms. It can be clearly seen that LCGSA has better values for the image quality measures than its competitors because p-values are less than 0.05, indicating the null hypothesis is rejected. Moreover, it implies that LCGSA was successful in providing optimal values to the pixels at all the threshold levels (k = 2, 4, 6, 8, and 10) in the complex search space. In short, the results of the Wilcoxon statistical test specify that LCGSA is an efficient optimizer compared to other peer algorithms.

6. Ablation Study

The ablation study was carried out to further benchmark the performance of LCGSA in solving segmentation problems. In this study, we have considered grayscale and ground truth images from the famous COVID-19 CT scan lesion segmentation dataset from Kaggle (https://www.kaggle.com/datasets/maedemaftouni/covid19-ct-scan-lesion-segmentation-dataset, accessed on 5 June 2023). All the images have 512 × 512 pixel dimensions, and the search space consists of 0–255 pixels. In Figure 24, the grayscale images, ground truth images, and their corresponding histograms are presented. The LCGSA was applied to all three images, and the simulation results were compared with twelve state-of-the-art algorithms. The simulation setup and system specifications are the same as those mentioned in Section 3.

6.1. Performance Metrics

In this study, we have considered five commonly used performance metrics in semantic segmentation for evaluation purposes. Our goal is to score the similarity between the predicted (prediction) and annotated segmentation (ground truth). The five evaluation metrics are Pixel Accuracy (Rand Index), Precision, Recall, Dice Coefficient (Dice Score or F1-Score), and Jaccard Index (Intersection over Union (IoU)). All presented metrics are based on the computation of a confusion matrix for a binary segmentation mask, which contains the number of true positive (TP), false positive (FP), true negative (TN), and false negative (FN) predictions. The value ranges of all presented metrics span from zero (worst) to one (best).
Pixel Accuracy: The accuracy score, also known as the Rand index, is the number of correct predictions, consisting of correct positive and negative predictions, divided by the total number of predictions, as shown in Equation (32).
Accuracy = T P + T N T P + T N + F N + F P  
Precision: The precision score is the number of true positive results divided by the number of all positive results, as shown in Equation (33).
Precision = T P T P + F P  
Recall: It is also known as Sensitivity or true positive rate, is the number of true positive results divided by the number of all samples that should have been identified as positive. It is mathematically calculated as shown in Equation (34).
Recall = T P T P + F N  
Dice Coefficient: It is also called the F1-measure or F-score. It is one of the most widely used scores for performance measurement in computer vision and in MIS (Medical Image Segmentation). The dice coefficient is calculated from the precision and recall of a prediction. It scores the overlap between predicted segmentation and ground truth. It also penalizes false positives, which is a common factor in highly class-imbalanced datasets like MIS. It is a harmonic mean of precision and recall. In other words, it is calculated by 2 times intersection divided by the total number of pixels in both images, as shown in Equation (35).
Dice Coefficient = 2 T P 2 T P + F P + F N  
Intersection over Union (IoU): It is also referred to as the Jaccard Index. It is essentially a method to quantify the percent overlap between the target mask and our segmented output. Quite simply, the IoU metric measures the number of pixels common between the target and prediction masks divided by the total number of pixels present across both masks, as shown in Equation (36).
IoU = T P T P + F P + F N  
This metric ranges from 0 to 1 (0–100%), with 0 signifying no overlap and 1 signifying perfectly overlapping segmentation.

6.2. Quantitative and Qualitative Analysis of the Results

The experimental results of CT-g1, CT-g2, and CT-g3 images are shown in Table 13, Table 14 and Table 15, respectively. It can be clearly seen that the LCGSA has better values for pixel accuracy, dice coefficient, and Jaccard index. In Table 13, the LCGSA versions have dice scores of <0.97, 0.96, 0.97, 0.90, 0.98, 0.90, 0.98, 0.98, 0.97, 0.98>, which shows that LCGSA-based segmentation has more similarity with the ground truth image. As far as IoU values are concerned, LCGSA has obtained <0.95, 0.93, 0.95, 0.83, 0.96, 0.82, 0.97, 0.97, 0.94, 0.96> values that are close to 1, indicating optimal overlap between the segmented image and the ground truth mask. It can also be seen that BBO, DE, and SCA also provide better values for the dice coefficient, Jaccard index, and pixel accuracy.
Similarly, in Table 14, the results of the CT-g2 image show the optimal performance of LCGSA versions in the segmentation of the grayscale image considering the annotated infection mask. When we compare the results of standard GSA and LCGSA values, it can be observed that LCGSA has obtained better results for accuracy (GSA(0.96), LCGSA(0.98)), dice score (GSA(0.96), LCGSA(0.98)), and Jaccard index (GSA(0.93), LCGSA(0.97)). The optimal values of the dice score and Jacquard index indicated that LCGSA has the capability of finding feasible regions in the complex pixel search space. Moreover, the LCGSA runtime values clearly indicate a proper balance between the exploration and exploitation stages.
Table 15 presents the results of LCGSA and other peer algorithms for the CT-g3 image. It can be observed that PSO, PSOGSA, CPSOGSA, BBO, and SCA provide optimal values for the image quality metrics. On the other hand, GSA, DE, MFO, ABC, GWO, and SMA have sub-optimal values for pixel accuracy, dice coefficient, and Jaccard index, indicating a higher false positive rate and issues in handling complex search spaces. Moreover, most of the LCGSA versions have shown feasible results for segmentation. However, it is also obvious that LCGSA2, LCGSA6, LCGSA8, and LCGSA10 have infeasible results, indicating optimization issues while handling complex segmentation problem spaces.
Figure 25 shows the colormap segmentation output of LCGSA at threshold level ten. It can be seen that segmented output has a total resemblance to the ground truth. It indicates that LCGSA was successful in extracting the optimal pixels from the original image. Moreover, it also shows that LCGSA has obtained optimal values for the Jaccard index and dice coefficient because the segmented image has maximal overlap with the ground truth image. Furthermore, Figure 26 shows the box plot analysis of the simulation results of CT-g1, CT-g2, and CT-g3 images. It is obvious that LCGSA has obtained large values for Kapur’s objective function, while other peer algorithms have obtained smaller values. It can also be seen that DE, SCA, and BBO provide feasible values for the objective function.

6.3. Statistical Analysis of the Results

The non-parametric statistical test, namely the signed Wilcoxon rank-sum test, was employed to verify the simulation results statistically at the 5% significance level. We have carried out the pairwise Wilcoxon rank-sum test between LCGSA and other peer algorithms separately. It can be seen in Table 16 that the p-values of all the peer algorithms are less than 0.05. In simpler terms, it implies that the null hypothesis is rejected while the alternate hypothesis is accepted. The null hypothesis is that LCGSA has issues obtaining better results for the image quality measures, and the alternate hypothesis is that LCGSA has the potential to obtain optimal pixels in the complex pixel search space. The p-values are clearly indicating that LCGSA values are statistically superior to other competitive algorithms in consecutive generations.

7. Overall Analysis of Simulation Results

The experimental analysis of benchmark images and CT scan images clearly showed that LCGSA is a superior and robust image segmentation technique. It is because LCGSA provided optimal values for image pixels at various threshold levels. It was quite interesting to see that all LCGSA versions took very little computational overhead to segment the benchmark images and locate the feasible pixels in the complex search space environment. At the same time, it can be noted that LCGSA10 struggled as far as runtime is concerned. Meanwhile, the CPU time results of LCGSA versions were better than those of other heuristic algorithms. Further, it is also evident that LCGSA versions have better values for PSNR, SSIM, and FSIM, indicating segmentation competence. We have found that in all the benchmark image results, the threshold values of LCGSA3 were infeasible and comparatively sub-optimal as compared to other LCGSA versions. Additionally, LCGSA6, LCGSA5, and LCGSA7 were the best optimizers in locating the optimal pixels and providing efficient values for Kapur’s objective function.
Similarly, when we closely look at the simulation outcomes of competitive heuristic algorithms for benchmark images, we can observe that DE, SCA, and BBO provided suitable values for the pixels. Meanwhile, DE performed better because its results were statistically superior to those of other heuristic approaches, indicating optimization proficiency and the potential for handling complex solution spaces. Moreover, DE and SCA took less runtime to segment the benchmark images. Also, we have noticed that SMA, GWO, SSA, MFO, PSO, and PSOGSA results show the presence of infeasible fitness values portraying noise, attenuation, and outliers. It also conveys convergence issues and difficulty in overcoming stagnation in local minima. Moreover, it was also surprising to see ABC, GSA, BBO, DE, and SSA taking more CPU time to reach the global optimum. Furthermore, the performance of GWO, SMA, PSO, PSOGSA, CPSOGSA, and ABC was ordinary because they provided appropriate values for PSNR, SSIM, and FSIM, but their results for other measures were unsatisfactory.
In the ablation study, the experiments were conducted on three COVID-19 chest CT scan images by considering infection masks. We have found that LCGSA versions are providing optimal results for dice coefficient, Jaccard index, and pixel accuracy metrics. It shows LCGSA has the potential to detect the abnormal portions in the CT images with high accuracy. Moreover, LCGSA again took less computational time to find the best pixels in the search space. It also indicates a proper balance between exploration and exploitation in LCGSA during the optimization process. Furthermore, the Wilcoxon rank-sum test verified the statistical superiority of the LCGSA simulation results over other peer algorithms. On the other hand, we have observed that SMA, GWO, CPSOGSA, and SSA provide suboptimal results for the image quality metrics. Meanwhile, DE and SCA provided better results for the dice coefficient, Kapur’s objective function, Jaccard index, and pixel accuracy.
Now, keeping in view the simulation results of both LCGSA and heuristic algorithms, we can unanimously declare LCGSA a robust optimizer for multilevel thresholding because its outcomes were superior and more effective than those of other peer algorithms. Moreover, the LCGSA versions took very little computational time to find the best pixels. They had statistically excellent results, and quite interestingly, the performance got better with the increase in the threshold levels. Lastly, the LCGSA versions were successful in segmenting the CT scan images of COVID-19 patients and were also proficient in finding the areas where the lungs were affected by the virus. Furthermore, the ablation study further authenticated the optimal performance of LCGSA in segmentation as it efficiently segmented the chest CT scan images while considering the gold standard (ground truth) of the original images.
Similarly, our study has also faced some challenges. Firstly, we have used a simple experimental setup consisting of 8 GB RAM and an i7 processor. However, we can obtain better simulation results if we use a system with advanced configurations like 128 GB RAM and GPUs for the simulation analysis. Secondly, in the experimental analysis, we found that certain LCGSA versions, such as LCGSA3 and LCGSA10, are facing premature convergence issues and taking more computational time to find the optimal pixels in the images. It shows that there is certainly a room for further improvement in these versions in specific areas like improving local exploitation capability to avoid convergence issues. Moreover, we have observed that other peer algorithms are also obtaining better results for image quality metrics such as PSNR, SSIM, and FSIM as the threshold values are increased. It shows that LCGSA needs more exploration capability in order to beat its competitor algorithms by a wide margin in the comparative analysis. Therefore, we will be using other learning strategies like opposition-based learning and the Gbest strategy to further enhance the segmentation and optimization capabilities of LCGSA.

8. Conclusions and Future Scope

An efficient hybrid strategy, namely the Levy flight and Chaos theory-based Gravitational Search Algorithm (LCGSA), was applied for the multilevel thresholding of the grayscale images. Kapur’s entropy scheme was combined with LCGSA in order to find the best pixels in the search space. Two benchmark images, namely Airport and Boat, and three chest CT scan images, such as CT1, CT2, and CT3, of COVID-19 patients were employed for the empirical analysis. Various performance metrics like Dice Coefficient, Pixel Accuracy, SSIM, FSIM, PSNR, CPU Time, Jaccard Index, and so on were utilized for the performance evaluation. The simulation results of LCGSA were compared with 12 state-of-the-art peer algorithms. Moreover, an ablation study consisting of various chest CT scan images and infection masks, was carried out to further evaluate the segmentation performance of LCGSA. The simulation results clearly revealed that LCGSA provides better values for the image pixels at various threshold intensity levels. Further, LCGSA requires less computational overhead to converge toward a feasible neighborhood. Moreover, the LCGSA has optimal values for pixel accuracy, PSNR, SSIM, and FSIM, indicating symmetry, quality, and consistency in the segmented output images. Furthermore, the LCGSA versions were proficient in locating the consolidated and abnormal patchy areas in the CT scan images, which can aid doctors in adequately diagnosing COVID-19-symptomatic patients. Meanwhile, it conveys the applicability of LCGSA to solving real-world image processing problems. In addition, we noticed that DE, SCA, and BBO also provided suitable values for the image thresholds.
Lastly, if we think about the future prospectus of LCGSA in image processing, we can see lots of possibilities, such as it will be interesting to apply LCGSA to segment the color RGB benchmark images in place of grayscale images. Similarly, the LCGSA can be utilized for the segmentation of crack images in the civil engineering field. Moreover, X-rays and other medical images, such as brain tumor images, heart images, lupus nephritis images, and so on, can be considered for evaluating the efficiency of LCGSA in medical imaging. In fact, it will be fascinating to employ Otsu’s variance scheme and the Rényi entropy method as fitness functions for LCGSA-based image segmentation. Furthermore, opposition learning-based LCGSA can be proposed and applied for the multilevel thresholding task. In addition, a version of the LCGSA for solving multi-objective problems is included in our future work. Moreover, a binary version of LCGSA for solving the feature selection problems will be given in the future. Finally, we can also explore the potential of LCGSA in deep learning. In fact, LCGSA can be used to optimize the hyper-parameters, such as learning rate, number of epochs, selection of an activation function, batch size, number of hidden layers, and so on of the deep neural architectures, such as CNN, U-Net, and LSTM.

Author Contributions

Conceptualization, S.A.R.; Software, S.A.R.; Validation, S.A.R.; Formal analysis, S.A.R.; Investigation, S.A.R.; Resources, S.D.; Writing—original draft, S.A.R.; Writing—review & editing, S.D.; Visualization, S.A.R.; Supervision, S.D.; Project administration, S.D.; Funding acquisition, S.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

MATLAB source codes will be publicly available on the GitHub platform (https//:github.com/SAJADAHMAD1, accessed on 5 June 2023) and on MathWork’s webpage (https://in.mathworks.com/matlabcentral/profile/authors/6240015-sajad-ahmad-rather, accessed on 5 June 2023).

Conflicts of Interest

The authors clearly state that there are no conflict of interest, whether financial or professional, regarding the publication of this work.

Nomenclature

ISImage Segmentation
MTMultilevel Thresholding
HAHeuristic Algorithm
PSNRPeak Signal-to-Noise Ratio
STDStandard Deviation
SSIMStructural Similarity Index Measure
FSIMFeature Similarity Index Measure
MSEMean Square Error
BVBest Value
PSOParticle Swarm Optimization
CPSOGSAConstriction Coefficient-based PSO and GSA
GSAGravitational Search Algorithm
SSASalp Swarm Optimizer
BBOBiogeography-Based Optimizer
DEDifferentıal Evolution
SCASine–Cosine Algorithm
MFOMoth Flame Optimizer
ABCArtificial Bee Colony Algorithm
GWOGray Wolf Optimizer
SMASlime Mould Algorithm
MISMedical Image Segmentation

References

  1. Wolpert, D.H.; Macready, W.G. No free lunch theorems for optimization. IEEE Trans. Evol. Comput. 1997, 1, 67–82. [Google Scholar] [CrossRef]
  2. Alinaghian, M.; Tirkolaee, E.B.; Dezaki, Z.K.; Hejazi, S.R.; Ding, W. An augmented Tabu search algorithm for the green inventory-routing problem with time windows. Swarm Evol. Comput. 2021, 60, 100802. [Google Scholar] [CrossRef]
  3. Tirkolaee, E.B.; Mardani, A.; Dashtian, Z.; Soltani, M.; Weber, G.-W. A novel hybrid method using fuzzy decision making and multi-objective programming for sustainable-reliable supplier selection in two-echelon supply chain design. J. Clean. Prod. 2020, 250, 119517. [Google Scholar] [CrossRef]
  4. Bansal, P.; Kumar, S.; Pasrija, S.; Singh, S. A hybrid grasshopper and new cat swarm optimization algorithm for feature selection and optimization of multi-layer perceptron. Soft Comput. 2020, 24, 15463–15489. [Google Scholar] [CrossRef]
  5. Mirjalili, S. SCA: A Sine Cosine Algorithm for solving optimization problems. Knowl.-Based Syst. 2016, 96, 120–133. [Google Scholar] [CrossRef]
  6. 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]
  7. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey Wolf Optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef]
  8. Erdal, F.; Doğan, E.; Saka, M.P. Optimum design of cellular beams using harmony search and particle swarm optimizers. J. Constr. Steel Res. 2011, 67, 237–247. [Google Scholar] [CrossRef]
  9. Abualigah, L.; Yousri, D.; Elaziz, M.A.; 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]
  10. Abualigah, L.; Diabat, A.; Mirjalili, S.; Abd Elaziz, M.; Gandomi, A.H. The arithmetic optimization algo-rithm. Comput. Methods Appl. Mech. Eng. 2021, 376, 113609. [Google Scholar] [CrossRef]
  11. Ezugwu, A.E.; Agushaka, J.O.; Abualigah, L.; Mirjalili, S.; Gandomi, A.H. Prairie Dog Optimization Algorithm. Neural Comput. Appl. 2022, 34, 20017–20065. [Google Scholar] [CrossRef]
  12. Agushaka, J.O.; Ezugwu, A.E.; Abualigah, L. Dwarf Mongoose Optimization Algorithm. Comput. Methods Appl. Mech. Eng. 2022, 391, 114570. [Google Scholar] [CrossRef]
  13. Oyelade, O.N.; Ezugwu, A.E.-S.; Mohamed, T.I.A.; Abualigah, L. Ebola Optimization Search Algorithm: A New Nature-Inspired Metaheuristic Optimization Algorithm. IEEE Access 2022, 10, 16150–16177. [Google Scholar] [CrossRef]
  14. Abualigah, L.; Abd Elaziz, M.; Sumari, P.; Geem, Z.W.; Gandomi, A.H. Reptile Search Algorithm (RSA): A na-ture-inspired meta-heuristic optimizer. Expert Syst. Appl. 2022, 191, 116158. [Google Scholar] [CrossRef]
  15. Santos, M.R.; Costa, B.S.J.; Bezerra, C.G.; Andonovski, G.; Guedes, L.A. An evolving approach for fault diagnosis of dynamic systems. Expert Syst. Appl. 2022, 189, 115983. [Google Scholar] [CrossRef]
  16. Precup, R.-E.; David, R.-C.; Roman, R.-C.; Szedlak-Stinean, A.-I.; Petriu, E.M. Optimal tuning of interval type-2 fuzzy controllers for nonlinear servo systems using Slime Mould Algorithm. Int. J. Syst. Sci. 2021, 1–16. [Google Scholar] [CrossRef]
  17. Bojan-Dragos, C.-A.; Precup, R.-E.; Preitl, S.; Roman, R.-C.; Hedrea, E.-L.; Szedlak-Stinean, A.-I. GWO-Based Optimal Tuning of Type-1 and Type-2 Fuzzy Controllers for Electromagnetic Actuated Clutch Systems. IFAC-Pap. 2021, 54, 189–194. [Google Scholar] [CrossRef]
  18. Khurma, R.A.; Aljarah, I.; Sharieh, A.; Mirjalili, S. EvoloPy-FS: An Open-Source Nature-Inspired Optimization Framework in Python for Feature Selection. Evol. Mach. Learn. Tech. Algorithms Appl. 2020, 131–173. [Google Scholar] [CrossRef]
  19. Karaboga, D. An Idea Based on Honeybee Swarm for Numerical Optimization; Erciyes University: Kayseri, Türkiye, 2005; Volume 200, pp. 1–10. [Google Scholar]
  20. Mirjalili, S. Moth-flame optimization algorithm: A novel nature-inspired heuristic paradigm. Knowl. Based Syst. 2015, 89, 228–249. [Google Scholar] [CrossRef]
  21. Kennedy, J.; Eberhart, R. Particle Swarm Optimization. IEEE Int. Conf. Neural Netw. 1995, 4, 1942–1948. [Google Scholar]
  22. Pereira, J.L.J.; Francisco, M.B.; Diniz, C.A.; Oliver, G.A.; Cunha, S.S.; Gomes, G.F. Lichtenberg algorithm: A novel hybrid physics-based meta-heuristic for global optimization. Expert Syst. Appl. 2021, 170, 114522. [Google Scholar] [CrossRef]
  23. Erol, O.K.; Eksin, I. A new optimization method: Big Bang–Big Crunch. Adv. Eng. Softw. 2006, 37, 106–111. [Google Scholar] [CrossRef]
  24. Rashedi, E.; Nezamabadi-Pour, H.; Saryazdi, S. GSA: A Gravitational Search Algorithm. Inf. Sci. 2009, 179, 2232–2248. [Google Scholar] [CrossRef]
  25. Khalilpourazari, S.; Naderi, B.; Khalilpourazary, S. Multi-Objective Stochastic Fractal Search: A powerful algorithm for solving complex multi-objective optimization problems. Soft Comput. 2020, 24, 3037–3066. [Google Scholar] [CrossRef]
  26. Holland, J.H. Genetic Algorithms. Sci. Am. 1992, 267, 66–73. [Google Scholar] [CrossRef]
  27. Ma, H.; Simon, D. Biogeography-based optimization with blended migration for constrained optimization problems. In Proceedings of the 12th Annual Conference on Genetic and Evolutionary Computation, Portland, OR, USA, 7–11 July 2010; pp. 417–418. [Google Scholar] [CrossRef]
  28. 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]
  29. Wu, F.; Zhao, S.; Yu, B.; Chen, Y.-M.; Wang, W.; Song, Z.-G.; Hu, Y.; Tao, Z.-W.; Tian, J.-H.; Pei, Y.-Y.; et al. A new coronavirus associated with human respiratory disease in China. Nature 2020, 579, 265–269. [Google Scholar] [CrossRef]
  30. Khalilpourazari, S.; Doulabi, H.H.; Çiftçioğlu, A.; Weber, G.-W. Gradient-based grey wolf optimizer with Gaussian walk: Application in modelling and prediction of the COVID-19 pandemic. Expert Syst. Appl. 2021, 177, 114920. [Google Scholar] [CrossRef]
  31. Huang, C.; Wang, Y.; Li, X.; Ren, L.; Zhao, J.; Hu, Y.; Zhang, L.; Fan, G.; Xu, J.; Gu, X.; et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet 2020, 395, 497–506. [Google Scholar] [CrossRef]
  32. Chan, J.F.; Yuan, S.; Kok, K.H.; To, K.K.; Chu, H.; Yang, J.; Xing, F.; Liu, J.; Yip, C.C.; Poon, R.W.; et al. A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: A study of a family cluster. Lancet 2020, 395, 514–523. [Google Scholar] [CrossRef]
  33. World Health Organization. Laboratory Testing for Coronavirus Disease 2019 (COVID-19) in Suspected Human Cases; WHO: Geneva, Switzerland, 2020; pp. 1–7. [Google Scholar]
  34. Kumar, S.; Sharma, P.P.; Shankar, U.; Kumar, D.; Joshi, S.K.; Pena, L.; Durvasula, R.; Kumar, A.; Kempaiah, P.; Poonam Rathi, B. Discovery of New Hydroxyethylamine Analogs against 3CLpro Protein Target of SARS-CoV-2: Molecular Docking, Molecular Dynamics Simulation, and Structure–Activity Relationship Studies. J. Chem. Inf. Model. 2020, 60, 5754–5770. [Google Scholar] [CrossRef]
  35. Le, T.T.; Andreadakis, Z.; Kumar, A.; Román, R.G.; Tollefsen, S.; Saville, M.; Mayhew, S. The COVID-19 vaccine development landscape. Nat. Rev. Drug Discov. 2020, 19, 305–306. [Google Scholar] [CrossRef]
  36. V’kovski, P.; Kratzel, A.; Steiner, S.; Stalder, H.; Thiel, V. Coronavirus biology and replication: Implications for SARS-CoV-2. Nat. Rev. Microbiol. 2021, 19, 155–170. [Google Scholar] [CrossRef]
  37. Toyoshima, Y.; Nemoto, K.; Matsumoto, S.; Nakamura, Y.; Kiyotani, K. SARS-CoV-2 genomic variations associated with mortality rate of COVID-19. J. Hum. Genet. 2020, 65, 1075–1082. [Google Scholar] [CrossRef] [PubMed]
  38. Singh, P.; Bose, S.S. A quantum-clustering optimization method for COVID-19 CT scan image segmentation. Expert Syst. Appl. 2021, 185, 115637. [Google Scholar] [CrossRef] [PubMed]
  39. Munusamy, H.; Muthukumar, K.J.; Gnanaprakasam, S.; Shanmugakani, T.R.; Sekar, A. FractalCovNet architecture for COVID-19 Chest X-ray image Classification and CT-scan image Segmentation. Biocybern. Biomed. Eng. 2021, 41, 1025–1038. [Google Scholar] [CrossRef] [PubMed]
  40. Feng, H.; Liu, Y.; Lv, M.; Zhong, J. A case report of COVID-19 with false negative RT-PCR test: Necessity of chest CT. Jpn. J. Radiol. 2020, 38, 409–410. [Google Scholar] [CrossRef] [PubMed]
  41. Ai, T.; Yang, Z.; Hou, H.; Zhan, C.; Chen, C.; Lv, W.; Tao, Q.; Sun, Z.; Xia, L. Correlation of Chest CT and RT-PCR Testing for Coronavirus Disease 2019 (COVID-19) in China: A Report of 1014 Cases. Radiology 2020, 296, E32–E40. [Google Scholar] [CrossRef]
  42. Sarkar, A.; Vandenhirtz, J.; Nagy, J.; Bacsa, D.; Riley, M. Identification of Images of COVID-19 from Chest X-rays Using Deep Learning: Comparing COGNEX VisionPro Deep Learning 1.0TM Software with Open Source Convolutional Neural Networks. SN Comput. Sci. 2021, 2, 130. [Google Scholar] [CrossRef]
  43. Liu, L.; Zhao, D.; Yu, F.; Heidari, A.A.; Li, C.; Ouyang, J.; Chen, H.; Mafarja, M.; Turabieh, H.; Pan, J. Ant colony optimization with Cauchy and greedy Levy mutations for multilevel COVID 19 X-ray image segmentation. Comput. Biol. Med. 2021, 136, 104609. [Google Scholar] [CrossRef]
  44. Wang, G.; Liu, X.; Li, C.; Xu, Z.; Ruan, J.; Zhu, H.; Meng, T.; Li, K.; Huang, N.; Zhang, S. A Noise-Robust Framework for Automatic Segmentation of COVID-19 Pneumonia Lesions from CT Images. IEEE Trans. Med Imaging 2020, 39, 2653–2663. [Google Scholar] [CrossRef] [PubMed]
  45. Luo, S.; Li, Y.; Gao, P.; Wang, Y.; Serikawa, S. Meta-seg: A survey of meta-learning for image segmentation. Pattern Recognit. 2022, 126, 108586. [Google Scholar] [CrossRef]
  46. Wang, X.; Li, Z.; Huang, Y.; Jiao, Y. Multimodal medical image segmentation using multi-scale context-aware network. Neurocomputing 2021, 486, 135–146. [Google Scholar] [CrossRef]
  47. Oskouei, A.G.; Hashemzadeh, M.; Asheghi, B.; Balafar, M.A. CGFFCM: Cluster-weight and Group-local Feature-weight learning in Fuzzy C-Means clustering algorithm for color image segmentation. Appl. Soft Comput. 2021, 113, 108005. [Google Scholar] [CrossRef]
  48. Civit-Masot, J.; Luna-Perejón, F.; Corral, J.M.R.; Domínguez-Morales, M.; Morgado-Estévez, A.; Civit, A. A study on the use of Edge TPUs for eye fundus image segmentation. Eng. Appl. Artif. Intell. 2021, 104, 104384. [Google Scholar] [CrossRef]
  49. Fournel, J.; Bartoli, A.; Bendahan, D.; Guye, M.; Bernard, M.; Rauseo, E.; Khanji, M.Y.; Petersen, S.E.; Jacquier, A.; Ghattas, B. Medical image segmentation automatic quality control: A multi-dimensional approach. Med Image Anal. 2021, 74, 102213. [Google Scholar] [CrossRef] [PubMed]
  50. Cui, X.; Chang, S.; Li, C.; Kong, B.; Tian, L.; Wang, H.; Huang, P.; Yang, M.; Wu, Y.; Li, Z. DEAttack: A differential evolution based attack method for the robustness evaluation of medical image segmentation. Neurocomputing 2021, 465, 38–52. [Google Scholar] [CrossRef]
  51. Shu, X.; Yang, Y.; Wu, B. A neighbor level set framework minimized with the split Bregman method for medical image segmentation. Signal Process. 2021, 189, 108293. [Google Scholar] [CrossRef]
  52. Cai, Y.; Mi, S.; Yan, J.; Peng, H.; Luo, X.; Yang, Q.; Wang, J. An unsupervised segmentation method based on dynamic threshold neural P systems for color images. Inf. Sci. 2022, 587, 473–484. [Google Scholar] [CrossRef]
  53. Chen, X.; Huang, H.; Heidari, A.A.; Sun, C.; Lv, Y.; Gui, W.; Liang, G.; Gu, Z.; Chen, H.; Li, C.; et al. An efficient multilevel thresholding image segmentation method based on the slime mould algorithm with bee foraging mechanism: A real case with lupus nephritis images. Comput. Biol. Med. 2021, 142, 105179. [Google Scholar] [CrossRef]
  54. Dai, M.; Baylou, P.; Humbert, L.; Najim, M. Image segmentation by a dynamic thresholding using edge detection based on cascaded uniform filters. Signal Process. 1996, 52, 49–63. [Google Scholar] [CrossRef]
  55. Chakraborty, S.; Mali, K. Biomedical image segmentation using fuzzy multilevel soft thresholding system coupled modified cuckoo search. Biomed. Signal Process. Control. 2022, 72, 103324. [Google Scholar] [CrossRef]
  56. Wu, T.; Shao, J.; Gu, X.; Ng, M.K.; Zeng, T. Two-stage image segmentation based on nonconvex ℓ2−ℓp approximation and thresholding. Appl. Math. Comput. 2021, 403, 126168. [Google Scholar] [CrossRef]
  57. Kalyani, R.; Sathya, P.; Sakthivel, V. Trading strategies for image segmentation using multilevel thresholding aided with minimum cross entropy. Eng. Sci. Technol. Int. J. 2020, 23, 1327–1341. [Google Scholar] [CrossRef]
  58. Abdel-Basset, M.; Mohamed, R.; AbdelAziz, N.M.; Abouhawwash, M. HWOA: A hybrid whale optimization algorithm with a novel local minima avoidance method for multi-level thresholding color image segmentation. Expert Syst. Appl. 2022, 190, 116145. [Google Scholar] [CrossRef]
  59. Houssein, E.H.; Hussain, K.; Abualigah, L.; Elaziz, M.A.; Alomoush, W.; Dhiman, G.; Djenouri, Y.; Cuevas, E. An improved opposition-based marine predators algorithm for global optimization and multilevel thresholding image segmentation. Knowl.-Based Syst. 2021, 229, 107348. [Google Scholar] [CrossRef]
  60. Zhao, D.; Liu, L.; Yu, F.; Heidari, A.A.; Wang, M.; Oliva, D.; Muhammad, K.; Chen, H. Ant colony optimization with horizontal and vertical crossover search: Fundamental visions for multi-threshold image segmentation. Expert Syst. Appl. 2021, 167, 114122. [Google Scholar] [CrossRef]
  61. Cao, X.; Li, T.; Li, H.; Xia, S.; Ren, F.; Sun, Y.; Xu, X. A Robust Parameter-Free Thresholding Method for Image Segmentation. IEEE Access 2019, 7, 3448–3458. [Google Scholar] [CrossRef]
  62. 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]
  63. Kotte, S.; Kumar, P.R.; Injeti, S.K. An efficient approach for optimal multilevel thresholding selection for gray scale images based on improved differential search algorithm. Ain Shams Eng. J. 2018, 9, 1043–1067. [Google Scholar] [CrossRef]
  64. Khalilpourazari, S.; Pasandideh, S.H.R. Modeling and optimization of multi-item multi-constrained EOQ model for growing items. Knowl.-Based Syst. 2019, 164, 150–162. [Google Scholar] [CrossRef]
  65. Rather, S.A.; Bala, P.S. Swarm-based chaotic gravitational search algorithm for solving mechanical engineering design problems. World J. Eng. 2020, 17, 97–114. [Google Scholar] [CrossRef]
  66. Rather, S.A.; Bala, P.S. Analysis of Gravitation-Based Optimization Algorithms for Clustering and Classification. In Handbook of Research on Big Data Clustering and Machine Learning; IGI Global: Hershey, PA, USA, 2020; pp. 74–99. [Google Scholar] [CrossRef]
  67. Rather, S.A.; Bala, P.S. Hybridization of Constriction Coefficient Based Particle Swarm Optimization and Gravitational Search Algorithm for Function Optimization. In Proceedings of the International Conference on Advances in Electronics, Electrical & Computational Intelligence (ICAEEC), Prayagraj, India, 31 May–1 June 2019. [Google Scholar]
  68. Kandhway, P.; Bhandari, A.K. A Water Cycle Algorithm-Based Multilevel Thresholding System for Color Image Segmentation Using Masi Entropy. Circuits Syst. Signal Process. 2019, 38, 3058–3106. [Google Scholar] [CrossRef]
  69. Jamazi, C.; Manita, G.; Chhabra, A.; Manita, H.; Korbaa, O. Mutated Aquila Optimizer for assisting brain tumor segmentation. Biomed. Signal Process. Control. 2023, 105089. [Google Scholar] [CrossRef]
  70. Abualigah, L.; Diabat, A.; Sumari, P.; Gandomi, A. A Novel Evolutionary Arithmetic Optimization Algorithm for Multilevel Thresholding Segmentation of COVID-19 CT Images. Processes 2021, 9, 1155. [Google Scholar] [CrossRef]
  71. Abualigah, L.; Habash, M.; Hanandeh, E.S.; Hussein, A.M.; Al Shinwan, M.; Abu Zitar, R.; Jia, H. Improved Reptile Search Algorithm by Salp Swarm Algorithm for Medical Image Segmentation. J. Bionic Eng. 2023, 20, 1766–1790. [Google Scholar] [CrossRef]
  72. Su, H.; Zhao, D.; Yu, F.; Heidari, A.A.; Zhang, Y.; Chen, H.; Li, C.; Pan, J.; Quan, S. Horizontal and vertical search artificial bee colony for image segmentation of COVID-19 X-ray images. Comput. Biol. Med. 2022, 142, 105181. [Google Scholar] [CrossRef] [PubMed]
  73. Chakraborty, S.; Saha, A.K.; Nama, S.; Debnath, S. COVID-19 X-ray image segmentation by modified whale optimization algorithm with population reduction. Comput. Biol. Med. 2021, 139, 104984. [Google Scholar] [CrossRef] [PubMed]
  74. Zhang, Q.; Wang, Z.; Heidari, A.A.; Gui, W.; Shao, Q.; Chen, H.; Zaguia, A.; Turabieh, H.; Chen, M. Gaussian Barebone Salp Swarm Algorithm with Stochastic Fractal Search for medical image segmentation: A COVID-19 case study. Comput. Biol. Med. 2021, 139, 104941. [Google Scholar] [CrossRef] [PubMed]
  75. Houssein, E.H.; Helmy, B.E.; Oliva, D.; Jangir, P.; Premkumar, M.; Elngar, A.A.; Shaban, H. An efficient mul-ti-thresholding based COVID-19 CT images segmentation approach using an improved equilibrium optimizer. Biomed. Signal Process. Control 2022, 73, 103401. [Google Scholar] [CrossRef]
  76. Zhao, C.; Xu, Y.; He, Z.; Tang, J.; Zhang, Y.; Han, J.; Shi, Y.; Zhou, W. Lung segmentation and automatic detection of COVID-19 using radiomic features from chest CT images. Pattern Recognit. 2021, 119, 108071. [Google Scholar] [CrossRef] [PubMed]
  77. Jin, Q.; Cui, H.; Sun, C.; Meng, Z.; Wei, L.; Su, R. Domain adaptation based self-correction model for COVID-19 infection segmentation in CT images. Expert Syst. Appl. 2021, 176, 114848. [Google Scholar] [CrossRef] [PubMed]
  78. Nama, S. A novel improved SMA with quasi reflection operator: Performance analysis, application to the image seg-mentation problem of COVID-19 chest X-ray images. Appl. Soft Comput. 2022, 118, 108483. [Google Scholar] [CrossRef]
  79. Dimitrov, D.; Abdo, H. Tight independent set neighborhood union condition for fractional critical deleted graphs and ID deleted graphs. Discret. Contin. Dyn. Syst.-S 2018, 12, 711. [Google Scholar] [CrossRef]
  80. Gao, W.; Guirao, J.L.; Basavanagoud, B.; Wu, J. Partial multi-dividing ontology learning algorithm. Inf. Sci. 2018, 467, 35–58. [Google Scholar] [CrossRef]
  81. Jensi, R.; Jiji, G.W. An enhanced particle swarm optimization with levy flight for global optimization. Appl. Soft Comput. 2016, 43, 248–261. [Google Scholar] [CrossRef]
  82. Li, Y.; Li, X.; Liu, J.; Ruan, X. An improved bat algorithm based on lévy flights and adjustment factors. Symmetry 2019, 11, 925. [Google Scholar] [CrossRef]
  83. Mandelbrot, B.B. The Fractal Geometry of Nature; WH Freeman: New York, NY, USA, 1982. [Google Scholar]
  84. Yang, X.-S. Chapter 3-Random Walks and Optimization. In Nature-Inspired Optimization Algorithms; Yang, X.-S., Ed.; Elsevier: Oxford, UK, 2014; pp. 45–65. [Google Scholar]
  85. Gutowski, M. Levy Flights as an underlying mechanism for global optimization algorithms. arXiv 2001, arXiv:0106003. [Google Scholar]
  86. Pavlyukevich, I. Lévy flights, non-local search and simulated annealing. J. Comput. Phys. 2007, 226, 1830–1844. [Google Scholar] [CrossRef]
  87. Ramos-Fernández, G.; Mateos, J.L.; Miramontes, O.; Cocho, G.; Larralde, H.; Ayala-Orozco, B. Levy walk patterns in the foraging movements of spider monkeys (Ateles geoffroyi). Behav. Ecol. Sociobiol. 2004, 55, 223–230. [Google Scholar]
  88. Mirjalili, S.; Gandomi, A.H. Chaotic gravitational constants for the gravitational search algorithm. Appl. Soft Comput. J. 2017, 53, 407–419. [Google Scholar] [CrossRef]
  89. Alatas, B. Chaotic bee colony algorithms for global numerical optimization. Expert Syst. Appl. 2010, 37, 5682–5687. [Google Scholar] [CrossRef]
  90. Li, C.; Zhou, J.; Xiao, J.; Xiao, H. Parameters identification of chaotic system by chaotic gravitational search algorithm. Chaos Solitons Fractals 2012, 45, 539–547. [Google Scholar] [CrossRef]
  91. Mingjun, J.; Huanwen, T. Application of chaos in simulated annealing. Chaos Solitons Fractals 2004, 21, 933–941. [Google Scholar] [CrossRef]
  92. Gandomi, A.H.; Yang, X.-S.; Talatahari, S.; Alavi, A.H. Firefly algorithm with chaos. Commun. Nonlinear Sci. Numer. Simul. 2013, 18, 89–98. [Google Scholar] [CrossRef]
  93. Wang, N.; Liu, L.; Liu, L. Genetic algorithm in chaos. OR Trans. 2001, 5, 1–10. [Google Scholar]
  94. Li-Jiang, Y.; Tian-Lun, C. Application of Chaos in Genetic Algorithms. Commun. Theor. Phys. 2002, 38, 168–172. [Google Scholar] [CrossRef]
  95. Jothiprakash, V.; Arunkumar, R. Optimization of Hydropower Reservoir Using Evolutionary Algorithms Coupled with Chaos. Water Resour. Manag. 2013, 27, 1963–1979. [Google Scholar] [CrossRef]
  96. Zhenyu, G.; Bo, C.; Min, Y.; Binggang, C. Self-Adaptive Chaos Differential Evolution. In Advances in Natural Computation: Second International Conference ICNC, Xi’an, China, 24–28 September 2006; Springer: Berlin/Heidelberg, German, 2006; pp. 972–975. [Google Scholar]
  97. Saremi, S.; Mirjalili, S.M.; Mirjalili, S. Chaotic Krill Herd Optimization Algorithm. Procedia Technol. 2014, 12, 180–185. [Google Scholar] [CrossRef]
  98. Wang, G.-G.; Guo, L.; Gandomi, A.H.; Hao, G.-S.; Wang, H. Chaotic Krill Herd algorithm. Inf. Sci. 2014, 274, 17–34. [Google Scholar] [CrossRef]
  99. Peitgen, H.; Jurgens, H.; Saupes, D. Chaos and Fractals; Springer: New York, NY, USA, 1992. [Google Scholar]
  100. Li, Y.; Deng, S.; Xiao, D. A novel Hash algorithm construction based on chaotic neural network. Neural Comput. Appl. 2011, 20, 133–141. [Google Scholar] [CrossRef]
  101. Ott, E. Chaos in Dynamical Systems; Cambridge University Press: Cambridge, UK, 2002. [Google Scholar]
  102. Zhao, S.; Wang, P.; Heidari, A.A.; Chen, H.; Turabieh, H.; Mafarja, M.; Li, C. Multilevel threshold image seg-mentation with diffusion association slime mould algorithm and Renyi’s entropy for chronic obstructive pulmonary disease. Comput. Biol. Med. 2021, 134, 104427. [Google Scholar] [CrossRef] [PubMed]
  103. 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]
  104. Wilcoxon, F. Individual Comparisons by Ranking Methods. Biom. Bull. 1945, 1, 80. [Google Scholar] [CrossRef]
Figure 1. Stochastic behavior of chaotic maps.
Figure 1. Stochastic behavior of chaotic maps.
Mathematics 11 03913 g001
Figure 2. Flow chart of the LCGSA algorithm.
Figure 2. Flow chart of the LCGSA algorithm.
Mathematics 11 03913 g002
Figure 3. LCGSA-based image segmentation process.
Figure 3. LCGSA-based image segmentation process.
Mathematics 11 03913 g003
Figure 4. (a) Grayscale background of the Airport image. (b) Histogram of the Airport image.
Figure 4. (a) Grayscale background of the Airport image. (b) Histogram of the Airport image.
Mathematics 11 03913 g004
Figure 5. (a) Grayscale background of the Boat image. (b) Histogram of the Boat image.
Figure 5. (a) Grayscale background of the Boat image. (b) Histogram of the Boat image.
Mathematics 11 03913 g005
Figure 6. (a) Grayscale background of the CT1 image. (b) Histogram of the CT1 image.
Figure 6. (a) Grayscale background of the CT1 image. (b) Histogram of the CT1 image.
Mathematics 11 03913 g006
Figure 7. (a) Grayscale background of the CT2 image. (b) Histogram of the CT2 image.
Figure 7. (a) Grayscale background of the CT2 image. (b) Histogram of the CT2 image.
Mathematics 11 03913 g007
Figure 8. (a) Grayscale background of the CT3 image. (b) Histogram of the CT3 image.
Figure 8. (a) Grayscale background of the CT3 image. (b) Histogram of the CT3 image.
Mathematics 11 03913 g008
Figure 9. LCGSA segmented images, colormap images, and histogram curves for the Airport benchmark image at k = 2, 4, 6, 8, and 10.
Figure 9. LCGSA segmented images, colormap images, and histogram curves for the Airport benchmark image at k = 2, 4, 6, 8, and 10.
Mathematics 11 03913 g009
Figure 10. Convergence curves for the Airport image.
Figure 10. Convergence curves for the Airport image.
Mathematics 11 03913 g010
Figure 11. Box plots for the Airport image.
Figure 11. Box plots for the Airport image.
Mathematics 11 03913 g011
Figure 12. LCGSA segmented images, colormap images, and histogram curves for the Boat benchmark image at k = 2, 4, 6, 8, and 10.
Figure 12. LCGSA segmented images, colormap images, and histogram curves for the Boat benchmark image at k = 2, 4, 6, 8, and 10.
Mathematics 11 03913 g012
Figure 13. Convergence curves for the Boat image.
Figure 13. Convergence curves for the Boat image.
Mathematics 11 03913 g013
Figure 14. Box plots for the Boat image.
Figure 14. Box plots for the Boat image.
Mathematics 11 03913 g014
Figure 15. LCGSA segmented images, colormap images, and histogram curves for the CT1 image at k = 2, 4, 6, 8, and 10.
Figure 15. LCGSA segmented images, colormap images, and histogram curves for the CT1 image at k = 2, 4, 6, 8, and 10.
Mathematics 11 03913 g015
Figure 16. Convergence curves for the CT1 image.
Figure 16. Convergence curves for the CT1 image.
Mathematics 11 03913 g016
Figure 17. Box plots for the CT1 image.
Figure 17. Box plots for the CT1 image.
Mathematics 11 03913 g017
Figure 18. LCGSA segmented images, colormap images, and histogram curves for the CT2 image at k = 2, 4, 6, 8, and 10.
Figure 18. LCGSA segmented images, colormap images, and histogram curves for the CT2 image at k = 2, 4, 6, 8, and 10.
Mathematics 11 03913 g018
Figure 19. Convergence curves for the CT2 image.
Figure 19. Convergence curves for the CT2 image.
Mathematics 11 03913 g019
Figure 20. Box plots for the CT2 image.
Figure 20. Box plots for the CT2 image.
Mathematics 11 03913 g020
Figure 21. LCGSA segmented images, colormap images, and histogram curves for the CT3 image at k = 2, 4, 6, 8, and 10.
Figure 21. LCGSA segmented images, colormap images, and histogram curves for the CT3 image at k = 2, 4, 6, 8, and 10.
Mathematics 11 03913 g021
Figure 22. Convergence curves for the CT3 image.
Figure 22. Convergence curves for the CT3 image.
Mathematics 11 03913 g022
Figure 23. Box plots for the CT3 image.
Figure 23. Box plots for the CT3 image.
Mathematics 11 03913 g023
Figure 24. Grayscale images, ground truth images, and their histograms.
Figure 24. Grayscale images, ground truth images, and their histograms.
Mathematics 11 03913 g024
Figure 25. CT-g1, CT-g2, and CT-g3 images, their ground truths, and LCGSA colormap segmented output.
Figure 25. CT-g1, CT-g2, and CT-g3 images, their ground truths, and LCGSA colormap segmented output.
Mathematics 11 03913 g025
Figure 26. Box plots of the (a) CT-g1 image, (b) CT-g2 image, and (c) CT-g3 image.
Figure 26. Box plots of the (a) CT-g1 image, (b) CT-g2 image, and (c) CT-g3 image.
Mathematics 11 03913 g026aMathematics 11 03913 g026b
Table 1. Meta-heuristic methods for multilevel image thresholding.
Table 1. Meta-heuristic methods for multilevel image thresholding.
ReferenceAlgorithm UsedThresholding
Technique
PerformanceComparative AlgorithmsPerformance Metrics
Abualigah et al., 2023 [71]RSA-SSAOtsu’s variance schemeImproved segmentation of COVID-19 images and reduction in computational overheadAO, WOA, SSA, RSA, MPA, and PSOSSIM, PSNR, Best Fitness values, and statistical tests
Jamazi et al., 2023 [69]AOK-meansImproved brain tumor detectionFuzzy C-means, U-Net, Z-Net, Adaptive K-means, SegNet, and so on PSNR, SSIM, MSE, DSC (Dice Similarity Coefficient), and Sensitivity
Su et al., 2022
[72]
CCABCKapur entropyImproved performance with high threshold valuesABC, SCA, MFO, PSO, SSA, CBA, ACWOA, IWOA, IGWO, and HHOPSNR, SSIM, and FSIM
Nama, 2022
[78]
QRSMAShannon entropyImproved accuracy and convergence speedSMA, MFO, SCA, SHO, SOA, STOA, TSA, and WOAMSE and PSNR
Houssein et al., 2022 [75]I-EOFuzzy entropyIncreased accuracy, PSNR, SSIM, and FSIMAGDE, GWO, MFO, SCA, HHO, and TSAPSNR, SSIM, and FSIM
Abualigah et al., 2021 [70]AOAKapur entropyImproved quality of segmentationAO, WOA, SSA, PSO, MPA, and DE PSNR, SSIM, and Optimal threshold values
Chakraborty et al., 2021 [73]mWOAPRKapur entropyEnhanced performance WOA, HBO, HGS, SMA, and variant algorithms of WOAPSNR and SSIM
Liu et al., 2021
[43]
CLACOKapur entropyImproved performance of search capability and convergence speedGWO, MFO, PSO, ACOR (ant colony optimization (ACO) for continuous domains), SCA, WOA, OBLGWO (boosted GWO), mSCA (modified SCA), and OBSCA (opposition-based SCA)PSNR, SSIM, and FSIM
Singh et al., 2021
[38]
FFQOAKEuclidean distanceImproved MSE, PSNR, and JSCGAK, PSOK, DPSOK, and ACOKMSE, PSNR, Jaccard Similarity Coefficient (JSC), and MSE
Zhang et al., 2021
[74]
GBSFSSSAKapur entropyImproved performance of medical image segmentation, search capability, and convergence speedPSO, SCA, BA, FA, MFO, WOA, and HHOPSNR, SSIM, and FSIM
Zhao et al., 2021
[76]
SP-V-NetSigmoid cross-entropyImproved accuracy, sensitivity, and accelerated convergenceMC-V-Net (multi-channel V-Net) and V-NetOptimal segmentation
Munusamy et al., 2021 [39]FractalCovNetCross-entropyImproved accuracy, precision, and recallU-Net, DenseUNet, Segnet, FCN, ResnetUNet, ResNet5, Xception, Inception-
ResNetV2, and VGG-16
F-measure and Dice Coefficient
Jin et al., 2021
[77]
DASC-NetCross-entropyImproved segmentationU-Net, U2-Net, AdaptSegNet, and ADVENTSensitivity, Specificity, Jaccard, and Dice Coefficient
Kandhway et al., 2019 [68]WCAMasi/Tsallis entropiesConvergence speedBAT, PSO, WDO, MBO, and GOAPSNR, MSE, FSIM, and SSIM
Proposed methodLCGSAKapur entropyTo enhance segmentation and resolve computational issuesGSA, PSO, PSOGSA, CPSOGSA, SCA, SSA, BBO, and so onPSNR, SSIM, MSE, FSIM, BV, STD, and so on
Table 2. Simulation results for the Airport image using classical, hybrid, and recent HAs.
Table 2. Simulation results for the Airport image using classical, hybrid, and recent HAs.
AlgorithmkOptimal ThresholdsMeanSTDMSEPSNRSSIMFSIMBest
Value
Run Time
2121, 13615.940.784898.2611.230.230.7516.764.8409
471, 89, 118, 15817.220.091718.1215.780.600.8820.3113.3771
GSA6197, 206, 147, 221, 176, 20517.820.326441.5810.040.110.6728.7511.7766
8160, 163, 125, 156, 180, 194, 141, 16027.890.525093.2711.060.220.7333.3619.1112
10125, 177, 187, 204, 133, 191, 174, 183, 190, 21035.390.255103.7111.050.220.7437.0319.6661
27, 210.440.415101.2011.050.200.438.414.6734
45, 17, 26, 2910.500.654173.2711.920.250.4315.4511.6513
PSO622, 36, 54, 81, 59, 7623.320.92933.8418.420.710.7523.851.4316
816, 22, 25, 46, 72, 76, 77, 11427.541.11479.4821.320.890.9128.8115.5366
108, 18, 29, 32, 33, 77, 103, 90, 69, 9329.161.18684.6219.770.820.8533.8017.4482
212, 158.350.515881.1410.430.150.429.754.1786
417, 39, 70, 740.370.471303.7916.970.660.7316.5610.9052
PSOGSA614, 23, 31, 34, 36, 5122.720.892231.5114.640.400.5524.759.0639
81, 6, 21, 33, 64, 69, 93, 10619.720.41593.2320.390.850.8731.8812.3254
101, 51, 55, 68, 79, 85, 97, 107, 128, 13430.090.45502.6421.110.840.9034.5615.9034
26, 228.150.394978.2111.160.210.4310.654.5809
47, 7, 34, 5213.750.532183.6514.730.420.5716.7611.5291
CPSOGSA625, 32, 36, 60, 102, 11116.030.24555.3820.680.880.8926.329.1880
87, 10, 32, 36, 57, 60, 61, 6317.890.971663.7415.910.520.6430.6914.0344
102, 14, 22, 29, 31, 48, 64, 101, 109, 11431.470.43414.5121.950.900.9034.3715.9678
2118, 798.120.972371.6714.380.520.8616.936.0754
4146, 94, 44, 22016.490.72654.7319.970.870.9220.5511.2251
BBO663, 194, 59, 33, 215, 14618.290.82969.4718.260.780.8325.1615.9688
882, 24, 255, 105, 198, 198, 134, 3719.810.89495.1621.180.900.9432.6221.9184
10185, 203, 92, 230, 123, 238, 120, 5, 77, 9339.991.111698.7515.820.630.8944.8723.9526
2250, 10612.511.853899.2312.220.320.8016.926.3959
4227, 59, 235, 15020.811.601907.8315.320.600.8123.5111.1228
DE643, 141, 88, 41, 167, 21826.012.18568.2120.580.890.9329.9815.9738
8225, 196, 194, 140, 61, 250, 157, 1730.102.421241.3417.190.750.8736.6422.2046
10117, 40, 188, 219, 137, 39, 41, 21, 74, 734.982.68369.5322.450.930.9544.0523.7699
27, 1439.971.632239.3714.620.550.8717.255.6035
479, 132, 151, 318.861.792218.6914.660.560.8924.4314.3033
SCA6108, 157, 37, 240, 8, 5928.672.56471.5721.390.910.9431.1413.8524
824, 8, 133, 166, 248, 40, 118, 17430.472.811145.6617.540.760.8239.3822.8713
10178, 26, 50, 2, 70, 84, 188, 37, 66, 19440.333.26539.1420.810.840.8547.5421.6741
225, 19.973.974621.2511.480.230.4312.595.1492
497, 22, 131, 9917.411.951585.8016.120.710.8718.2714.1512
SSA6206, 215, 255, 1, 253, 25514.523.197613.139.310.030.4925.6413.2537
8255, 54, 255, 255, 254, 205, 255, 229.623.842285.5414.540.460.6629.2424.3640
10255, 255, 255, 255, 133, 115, 228, 255, 255, 25531.689.634426.8011.660.280.7734.2721.0928
2255, 113.621.467981.009.110.010.4015.723.9914
4255, 240, 1, 110.715.917981.009.110.010.4019.746.3647
MFO61, 2, 255, 1, 255, 119.963.587981.009.110.010.4025.888.5121
8254, 1, 2, 1, 255, 254, 255, 130.812.817818.019.190.020.4126.2610.8618
10254, 255, 255, 141, 1, 255, 225, 255, 1, 25431.235.676075.1810.290.150.7036.0013.6318
2213, 22912.731.947870.289.170.010.4914.3210.2643
4153, 192, 193, 21417.382.096647.109.900.090.6520.2817.4111
ABC6196, 249, 152, 215, 243, 17521.972.246605.169.930.100.6527.3524.8882
8183, 180, 201, 184, 211, 220, 232, 16426.782.336955.699.700.070.6031.7130.8205
10212, 235, 242, 226, 238, 204, 225, 187, 228, 18527.342.867371.759.450.040.5432.5837.8003
233, 4013.231.263022.1213.320.310.4411.263.0042
475, 83, 82, 218.961.832524.0614.100.480.7814.845.5414
GWO63, 4, 0, 0, 0, 76.504.937033.069.650.070.4214.088.0195
80, 32, 9, 1, 25, 0, 0, 112.067.253757.310.270.440.4314.559.9017
100, 1, 2, 0, 0, 0, 1, 5, 3, 00.804.587981.009.110.010.408.1612.9664
2144, 1404.732.336373.5310.080.110.696.663.5162
4200, 201, 180, 1904.902.377692.869.260.020.506.805.8604
SMA61, 1, 1, 2, 5, 104.592.447981.009.110.010.406.778.4221
8197, 191, 189, 188, 190, 178, 197, 1974.792.327621.139.310.020.516.8310.7873
104, 5, 6, 3, 4, 4, 8, 9, 4, 64.772.347498.039.380.040.426.8313.4375
Table 3. Simulation results for the Airport image using ten LCGSA versions.
Table 3. Simulation results for the Airport image using ten LCGSA versions.
AlgorithmkOptimal ThresholdsMeanSTDMSEPSNRSSIMFSIMBest
Value
Run Time
2156, 24212.700.136822.279.790.080.6417.692 × 10−6
492, 160, 239, 2824.031.231236.6817.200.780.8926.162 × 10−6
LCGSA16156, 240, 32, 161, 244, 2631.401.702863.1713.560.490.6534.802 × 10−6
8192, 39, 235, 26, 91, 160, 242, 6136.981.85441.1721.680.900.9243.173 × 10−6
10128, 230, 18, 86, 138, 181, 248, 101, 246, 3045.893.07781.8119.190.840.9352.461 × 10−6
292, 16317.490.573100.4013.210.420.8617.691 × 10−6
4154, 243, 20, 15324.861.323995.8112.110.390.6626.162 × 10−6
LCGSA2695, 161, 243, 37, 206, 5530.411.46459.9421.500.910.9332.262 × 10−6
8143, 232, 28, 148, 245, 30, 213, 3237.651.952542.7714.070.560.7041.463 × 10−6
1092, 159, 235, 29, 85, 142, 188, 241, 49, 19145.072.67367.3722.470.930.9549.661 × 10−6
2130, 22712.230.215577.8110.660.180.7317.672 × 10−6
4142, 232, 160, 6918.720.492099.7314.900.570.8723.821 × 10−6
LCGSA36161, 77, 140, 199, 97, 18623.211.232083.7914.940.550.8832.812 × 10−6
8120, 175, 112, 192, 142, 86, 164, 7928.491.762140.9214.820.540.8838.723 × 10−6
10128, 162, 207, 91, 158, 88, 135, 168, 208, 11441.792.962618.0913.950.470.8644.271 × 10−6
292, 16017.520.503077.7813.240.420.8617.702 × 10−6
4204, 27, 95, 16023.881.411246.8117.170.780.8925.482 × 10−6
LCGSA46148, 245, 32, 156, 239, 3031.222.292779.5013.690.520.6835.202 × 10−6
897, 161, 246, 28, 92, 162, 240, 5540.443.00505.1621.090.890.9342.853 × 10−6
10179, 28, 156, 245, 27, 201, 32, 234, 32, 15643.313.702828.9813.610.490.6552.301 × 10−6
293, 16017.500.703128.8613.170.410.8617.692 × 10−6
491, 160, 238, 3223.801.631073.3817.820.810.9026.162 × 10−6
LCGSA5671, 100, 132, 165, 199, 24229.431.131696.5915.830.610.9035.373 × 10−6
8181, 18, 97, 162, 237, 29, 155, 24435.912.201163.8317.470.790.8943.171 × 10−6
10153, 246, 23, 158, 244, 63, 199, 26, 92, 15947.464.39618.4720.210.880.9351.742 × 10−6
294, 16017.410.923150.1413.140.410.8517.692 × 10−6
497, 161, 249, 2523.991.181417.9616.160.750.8925.482 × 10−6
LCGSA66155, 248, 16, 223, 27, 15731.292.353303.0612.940.450.6532.311 × 10−6
890, 145, 189, 245, 36, 154, 241, 3440.112.73768.1819.270.860.9241.583 × 10−6
1091, 160, 245, 39, 222, 36, 150, 238, 40, 23342.153.22683.6119.780.870.9152.102 × 10−6
2152, 25012.820.606695.829.870.090.6517.692 × 10−6
4227, 30, 237, 1720.210.603918.7312.190.270.4626.151 × 10−6
LCGSA7693, 160, 242, 28, 100, 16334.072.481122.4417.420.790.8935.372 × 10−6
888, 162, 250, 32, 205, 17, 151, 23836.142.07964.2718.280.830.9143.163 × 10−6
1096, 159, 236, 22, 91, 147, 186, 242, 44, 15749.595.04534.6420.850.890.9352.373 × 10−6
2231, 913.540.326569.599.950.100.4417.702 × 10−6
490, 160, 240, 2724.250.731236.7417.200.780.8925.481 × 10−6
LCGSA8691, 162, 237, 19, 154, 24129.741.761663.1315.920.700.8935.372 × 10−6
884, 132, 176, 247, 6, 154, 241, 5739.992.91823.0718.970.800.9242.513 × 10−6
10231, 8, 101, 160, 247, 22, 154, 240, 31, 16048.523.391148.8917.520.790.8850.831 × 10−6
295, 15917.580.533189.7013.090.400.8517.702 × 10−6
495, 160, 246, 2524.240.831426.8616.580.750.8926.162 × 10−6
LCGSA9693, 162, 249, 24, 225, 3030.711.461110.5717.670.800.8935.382 × 10−6
8210, 24, 242, 31, 157, 246, 18, 15137.723.062917.5113.480.500.6743.172 × 10−6
10153, 243, 18, 86, 142, 190, 246, 26, 180, 2546.963.111144.9817.540.800.9152.491 × 10−6
298, 16317.500.643315.1512.920.390.8417.695.6194
4237, 15, 239, 2220.170.824977.7211.160.210.4325.4812.0170
LCGSA10693, 149, 192, 245, 29, 15434.111.811119.7117.630.800.9135.3714.0934
8158, 241, 33, 89, 142, 184, 243, 3939.623.74651.1419.990.870.9242.5320.7450
10153, 243, 22, 150, 242, 34, 205, 9, 175, 3144.402.582615.5113.950.530.6750.9223.7940
Table 4. Simulation results for the Boat image using classical, hybrid, and recent HAs.
Table 4. Simulation results for the Boat image using classical, hybrid, and recent HAs.
AlgorithmkOptimal ThresholdsMeanSTDMSEPSNRSSIMFSIMBest
Value
Run Time
253, 6513.940.246240.4710.170.360.5715.395.3828
4146, 173, 136, 15417.250.193742.7712.390.460.6421.068.6466
GSA679, 150, 111, 160, 149, 8424.210.29821.4118.980.730.8326.3511.7598
8153, 101, 181, 95, 113, 146, 92, 7526.030.09590.2120.420.770.8732.3915.1586
1091, 106, 114, 142, 112, 145, 146, 132, 108, 13032.290.47866.6218.750.730.7939.7717.4477
237, 4512.920.449203.828.490.260.5214.874.2307
424, 34, 66, 6218.810.865985.3610.350.410.5719.467.3875
PSO67, 22, 26, 62, 65, 5224.260.776110.3010.270.400.5626.7110.0751
86, 6, 11, 28, 44, 44, 55, 7325.820.995076.8011.070.450.5930.5212.4567
108, 13, 19, 23, 28, 34, 46, 40, 49, 6932.440.805699.0510.570.430.5837.5914.7040
21, 19.621.8418,744.495.4000.3511.334.2504
44, 11, 13, 3814.850.4510,592.337.880.240.5017.586.6844
PSOGSA628, 44, 45, 60, 71, 8017.590.474385.8011.710.460.6024.389.1286
821, 31, 37, 40, 45, 109, 118, 13830.580.64742.2219.420.800.8033.9711.5332
103, 47, 50, 66, 77, 87, 133, 134, 154, 15632.230.65533.1420.860.810.8638.0314.4190
23, 567.420.567679.369.270.320.5613.583.8311
414, 15, 27, 3518.790.5510,951.407.730.220.4818.086.4918
CPSOGSA616, 22, 71, 94, 97, 9717.210.722793.9713.660.580.6723.569.5377
85, 10, 13, 14, 17, 31, 71, 7319.701.235235.2610.940.460.6031.9511.4005
101, 17, 17, 17, 35, 41, 87, 92, 94, 10933.131.201799.3915.570.660.7236.8814.2903
295, 16912.230.462359.5814.400.590.7517.947.0360
4188, 232, 140, 15221.410.624537.5711.560.410.6223.636.1897
BBO6242, 40, 199, 115, 228, 19425.320.511433.4116.560.700.7926.998.8856
840, 68, 190, 197, 192, 234, 115, 18529.890.551101.7717.700.740.8239.3020.1073
1042, 1, 204, 191, 41, 111, 58, 9, 11, 24328.371.261355.8416.800.740.8141.7025.5141
289, 14116.530.561304.5416.970.660.7717.556.5577
4175, 153, 191, 11118.452.691489.0516.400.660.8024.246.5118
DE6236, 156, 127, 252, 183, 16023.252.992181.0114.740.600.7332.639.0063
816, 144, 113, 12, 240, 19, 190, 25033.261.98945.4818.370.770.8039.4620.5685
1056, 230, 60, 69, 65, 161, 226, 87, 189, 1443.401.411871.3215.400.670.7746.7325.2107
285, 15212.761.421901.6915.330.600.7617.805.5593
4243, 16, 34, 10322.461.862377.2114.370.640.7124.749.8259
SCA664, 104, 216, 9, 163, 20529.262.371205.4917.310.740.8332.3113.4399
8123, 194, 255, 4, 25, 60, 36, 13033.503.05761.5919.310.790.8341.7317.0773
1061, 73, 107, 227, 33, 61, 225, 3, 14, 11541.103.651417.3716.610.710.7747.4119.4771
275, 2259.582.285275.2110.900.420.6214.385.5763
4223, 255, 254, 21623.134.0518,227.455.520.010.4618.349.7626
SSA6143, 1, 1, 117, 65, 7319.405.34781.0219.200.760.8325.0312.9424
8255, 255, 255, 255, 173, 255, 1, 25524.606.3815,122.676.330.090.5123.5018.0670
10112, 88, 76, 164, 137, 7, 146, 93, 1, 825.636.79452.0321.570.850.8734.2318.4091
2245, 25513.202.1018,994.395.3400.4017.353.6887
4255, 13, 255, 321.901.9618,233.185.520.020.4521.576.9507
MFO64, 1, 255, 1, 255, 25521.863.1717,980.765.580.030.4528.289.0449
861, 209, 255, 1, 255, 252, 35, 24921.902.5710,346.167.980.250.5334.4812.9243
101, 255, 255, 1, 77, 255, 1, 255, 211, 121.962.844586.4311.510.480.6741.7214.0383
291, 2113.091.771661.2415.920.620.7317.1310.0483
4164, 127, 160, 16818.151.952281.2814.540.580.7121.3319.0667
ABC6116, 165, 200, 144, 195, 16418.252.281403.2216.650.690.8025.9424.7383
8176, 159, 175, 205, 180, 230, 190, 15418.352.519814.848.210.230.5530.8734.3259
10178, 187, 197, 168, 208, 143, 170, 169, 223, 17518.972.765521.2710.710.370.6135.3739.0609
21, 01.771.5518,488.085.460.010.3711.723.2673
451, 44, 25, 4616.572.858178.599.000.310.5318.805.7354
GWO614, 80, 26, 26, 51, 423.692.084277.5511.810.500.6224.098.0261
819, 60, 40, 46, 70, 20, 32, 7030.022.695443.9410.770.430.5826.6410.3239
100, 0, 1, 1, 1, 1, 0, 0, 0, 13.088.5318,744.495.4000.3632.9014.1903
2185, 1855.362.5016,562.175.930.050.507.183.4958
4126, 120, 123, 1255.302.442881.7113.530.520.657.036.2318
SMA60, 1, 2, 0, 0, 34.852.6819,002.915.34007.188.9006
8179, 178, 177, 179, 179, 178, 179, 1785.522.0415,998.536.090.070.517.1811.6382
10170, 171, 172, 176, 176, 176, 176, 176, 175, 1764.922.6715,734.556.160.070.517.1914.0558
Table 5. Simulation results for the Boat image using ten LCGSA versions.
Table 5. Simulation results for the Boat image using ten LCGSA versions.
AlgorithmkOptimal ThresholdsMeanSTDMSEPSNRSSIMFSIMBest
Value
Run Time
2132, 613.900.343103.4913.210.510.6318.103 × 10−6
4175, 153, 191, 11124.172.691489.0516.400.660.8024.242 × 10−6
LCGSA1663, 124, 237, 6, 116, 23830.021.301212.1417.290.700.7635.752 × 10−6
8106, 233, 9, 114, 239, 9, 69, 12741.344.291010.0518.080.720.7744.862 × 10−6
1060, 113, 180, 239, 164, 4, 58, 97, 128, 18347.633.78440.8621.680.830.8952.552 × 10−6
2112, 24612.870.192737.2513.750.550.6818.131 × 10−6
4119, 245, 6, 11825.201.302325.0314.460.590.6826.871 × 10−6
LCGSA2675, 122, 181, 246, 5, 11834.042.46981.9318.210.740.8235.692 × 10−6
8127, 5, 109, 179, 245, 47, 181, 24635.772.25721.2519.540.780.8442.731 × 10−6
10108, 182, 246, 32, 185, 246, 46, 210, 42, 22341.903.591464.8016.470.720.8051.432 × 10−6
2105, 17617.890.422199.8914.700.610.7518.122 × 10−6
4134, 99, 193, 7817.881.60829.8718.940.740.8324.571 × 10−6
LCGSA36121, 188, 78, 128, 190, 9328.751.63890.1618.630.720.8132.381 × 10−6
8103, 139, 201, 98, 187, 86, 185, 15128.451.28684.5819.770.760.8539.141 × 10−6
10122, 180, 102, 145, 184, 153, 104, 179, 105, 14536.362.91906.7218.550.750.8346.322 × 10−6
2108, 17717.521.642158.9814.780.610.7518.112 × 10−6
471, 125, 245, 824.231.081314.5916.940.710.7725.841 × 10−6
LCGSA46108, 180, 245, 5, 107, 17834.263.221992.5915.130.640.7536.242 × 10−6
8114, 243, 25, 117, 238, 28, 182, 24236.451.981401.1416.660.720.7742.672 × 10−6
10110, 229, 29, 216, 8, 146, 5, 120, 222, 5544.332.89476.0221.350.860.8751.042 × 10−6
269, 12617.570.221395.2316.680.670.7618.121 × 10−6
4113, 236, 5, 10625.131.582311.4714.490.590.6926.871 × 10−6
LCGSA56111, 243, 9, 115, 176, 24130.171.681703.0815.810.670.7535.681 × 10−6
8132, 10, 106, 238, 7, 102, 181, 24436.033.38954.7618.330.750.8143.052 × 10−6
10108, 181, 244, 79, 205, 21, 113, 182, 240, 6246.754.241069.1717.840.750.8253.812 × 10−6
2109, 18117.950.552247.3414.610.600.7418.122 × 10−6
4240, 9, 242, 1121.080.8916,015.896.080.080.4625.841 × 10−6
LCGSA6673, 128, 184, 244, 4, 10634.901.77791.3319.140.760.8435.672 × 10−6
8109, 220, 3, 105, 226, 7, 73, 12541.223.861092.7717.740.710.7843.903 × 10−6
10109, 230, 6, 99, 181, 250, 29, 220, 6, 10149.013.101393.4216.690.720.7950.872 × 10−6
2239, 514.250.3617,709.115.640.040.4418.112 × 10−6
4119, 243, 10, 11625.191.392151.8814.800.610.6826.861 × 10−6
LCGSA76105, 178, 246, 41, 233, 1231.061.881618.1116.040.720.8035.731 × 10−6
8102, 176, 245, 5, 109, 229, 4, 10841.544.301814.5415.540.650.7643.882 × 10−6
10190, 8, 196, 14, 178, 243, 7, 108, 229, 22044.472.811618.3516.040.710.7752.311 × 10−6
2104, 18117.920.792341.6614.430.600.7518.132 × 10−6
468, 123, 180, 24722.570.871055.4917.890.720.8326.851 × 10−6
LCGSA8666, 127, 236, 8, 67, 12334.551.961157.0717.490.710.7736.202 × 10−6
863, 120, 178, 244, 7, 65, 121, 17743.523.62971.3618.250.760.8345.572 × 10−6
10107, 230, 7, 104, 185, 247, 9, 231, 11, 10949.263.841845.1515.470.670.7551.071 × 10−6
2106, 17917.661.452262.7914.580.600.7518.112 × 10−6
4116, 247, 8, 11825.570.802231.5914.640.610.6825.841 × 10−6
LCGSA9667, 125, 249, 9, 226, 1131.401.351212.9417.290.730.7833.912 × 10−6
866, 124, 238, 9, 239, 4, 118, 25036.822.191205.0417.320.710.7645.561 × 10−6
1050, 95, 131, 186, 247, 13, 181, 249, 10, 11949.944.17475.6521.350.830.8752.211 × 10−6
264, 12417.410.811440.9716.540.670.7618.125.6869
4110, 241, 5, 11525.181.882294.5514.520.590.6926.776.4817
LCGSA106108, 178, 242, 8, 247, 631.161.611989.4315.140.660.7532.578.8356
866, 123, 235, 4, 121, 236, 26, 23336.102.821191.9917.360.740.7844.2819.9924
10116, 222, 60, 225, 5, 119, 235, 6, 70, 11947.774.281315.2216.940.700.7753.4824.1814
Table 6. Simulation results for the CT1 image using classical, hybrid, and recent HAs.
Table 6. Simulation results for the CT1 image using classical, hybrid, and recent HAs.
AlgorithmkOptimal ThresholdsMeanSTDMSEPSNRSSIMFSIMBest
Value
Run Time
2132, 18415.800.043289.5412.950.400.6217.194.9900
4109, 143, 92, 12922.120.214190.3611.900.310.5822.259.2151
GSA6184, 142, 212, 154, 173, 17823.960.182980.1513.380.500.7124.4912.9996
8198, 220, 182, 176, 231, 177, 195, 21526.280.544468.8111.620.370.6830.8716.1160
10190, 200, 219, 215, 184, 216, 229, 204, 213, 222 30.450.245511.3610.710.310.6433.7519.3697
211, 195.261.1320,779.204.950.050.2710.764.2840
418, 23, 36, 4113.590.5815,772.396.150.120.3218.577.4031
PSO65, 5, 6, 4, 6, 615.977.5723,917.214.340.020.2721.8211.7522
813, 18, 22, 26, 26, 41, 63, 5620.510.7111,549.667.500.280.5030.2413.0390
101, 15, 22, 31, 38, 54, 45, 53, 76, 8537.411.138435.008.870.400.5738.1216.6038
21, 1811.430.4621,029.844.900.050.2712.353.7307
47, 27, 29, 3613.610.3816,826.145.870.100.3020.126.7921
PSOGSA611, 24, 28, 42, 42, 4916.280.7414,000.126.660.170.4023.4010.3811
87, 12, 20, 30, 34, 44, 62, 6827.960.4010,754.687.810.320.5229.7912.1737
108, 10, 14, 22, 31, 35, 82, 84, 94, 11831.340.794902.7811.220.450.6037.8115.2642
21, 237.830.5619,797.625.160.070.2812.863.7407
44, 33, 45, 5116.320.7513,816.236.720.180.4115.026.8445
CPSOGSA62, 62, 106, 113, 114, 12519.150.684489.0511.600.430.6224.5210.2939
81, 3, 6, 16, 37, 38, 39, 5221.711.5013,631.236.780.190.4332.0712.0088
1012, 15, 21, 23, 49, 64, 75, 80, 85, 10727.510.295791.6110.500.490.6237.6815.3088
2195, 18312.130.285523.0410.700.270.5514.113.7492
452, 162, 222, 19116.631.381023.6218.020.720.8224.4812.3297
BBO688, 174, 128, 208, 151, 2218.180.76934.4018.420.750.8326.7816.4833
8139, 32, 60, 160, 240, 178, 179, 18626.170.68718.0619.560.750.8332.4016.9126
10122, 245, 255, 134, 205, 249, 42, 119, 188, 5223.461.60585.9820.450.780.8645.5723.7955
2130, 5111.952.814054.9212.050.440.6117.043.8522
47, 67, 98, 15224.560.752655.6213.880.480.6624.7511.7624
DE665, 67, 129, 146, 50, 4725.222.402603.8213.970.530.6629.1616.6393
8132, 148, 47, 97, 52, 83, 102, 2636.071.842338.5014.440.570.6938.2817.8548
10188, 177, 140, 244, 106, 95, 73, 176, 114, 6841.222.931154.1517.500.670.7945.6623.8481
228, 9313.791.667783.889.210.350.5417.115.5373
465, 90, 179, 23923.541.821681.0715.870.630.7824.7910.5775
SCA6130, 227, 56, 78, 236, 1429.492.532299.8114.510.670.7731.4715.6037
836, 230, 170, 15, 58, 93, 128, 25133.962.87700.6519.670.820.8838.5718.3722
10195, 65, 165, 64, 155, 33, 125, 252, 57, 18142.943.04523.2220.940.770.8445.3622.6256
21, 115.135.9325,594.164.0400.2515.034.9294
41, 1, 255, 25418.475.8024,740.304.190.010.3220.1610.2434
SSA6134, 255, 255, 255, 1, 25524.524.085597.3210.650.260.5630.8516.2715
81, 1, 21, 47, 1, 2, 1, 1634.969.6514,587.926.490.170.4233.9815.7475
101, 255, 255, 1, 201, 192, 108, 255, 170, 25528.046.882310.8614.490.520.7339.4422.2014
2240, 121.901.3322,230.054.660.050.3815.033.8241
4255, 246, 249, 115.251.8723,488.804.420.030.3622.556.5664
MFO624, 1, 1, 96, 1, 322.113.787574.099.330.340.5432.379.8088
8157, 1, 254, 255, 255, 248, 252, 25529.084.984567.8811.530.300.5727.7611.5291
10188, 186, 273, 216, 142, 223, 142, 100, 6, 10342.022.291733.4715.740.640.8035.4014.3570
2204, 16218.151.733720.1312.420.420.6716.0610.4149
4177, 198, 201, 14817.741.973222.3513.040.440.6519.9417.7692
ABC6234, 154, 250, 213, 219, 23222.592.253548.9112.620.430.6826.6925.6627
898, 152, 134, 137, 148, 170, 108, 11826.992.642850.2313.580.360.6031.3932.8164
10218, 149, 163, 202, 217, 230, 214, 196, 234, 25331.672.743033.1813.310.520.7334.1840.1686
25, 101.583.5625,594.164.0400.2510.553.3954
424, 36, 1, 5514.391.2912,919.727.010.230.4714.856.1411
GWO62, 0, 2, 4, 5, 25.634.9225,027.184.140.010.2619.028.5206
810, 0, 66, 10, 0, 10, 10, 2122.352.4311,330.007.580.290.5218.5110.8290
100, 2, 3, 4, 0, 1, 0, 1, 1, 14.464.1324,191.714.290.010.279.4413.5712
2230, 2325.162.3520,319.375.050.070.407.313.7368
438, 37, 36, 355.072.4416,411.245.970.120.347.436.4121
SMA6160, 159, 159, 160, 158, 1604.562.684888.4511.230.240.507.489.0762
886, 85, 85, 84, 86, 86, 85, 864.922.499568.258.320.230.537.4211.9038
1052, 51, 53, 52, 52, 51, 51, 48, 48, 524.832.5313,602.576.790.220.467.4414.4625
Table 7. Simulation results for the CT1 image using ten LCGSA versions.
Table 7. Simulation results for the CT1 image using ten LCGSA versions.
AlgorithmkOptimal ThresholdsMeanSTDMSEPSNRSSIMFSIMBest
Value
Run Time
2101, 17718.160.522922.3613.470.420.6518.311 × 10−6
498, 172, 254, 2025.241.082012.9515.090.550.7027.021 × 10−6
LCGSA16147, 250, 34, 245, 10, 13332.342.342706.1213.800.490.6435.542 × 10−6
836, 87, 135, 177, 251, 13, 98, 18143.223.45990.4418.170.690.7745.222 × 10−6
10146, 43, 242, 25, 238, 17, 124, 247, 6, 13943.983.052577.9014.010.560.7052.262 × 10−6
297, 17418.150.452952.5013.420.420.6518.301 × 10−6
498, 175, 252, 2125.041.621935.0515.260.570.7126.842 × 10−6
LCGSA26218, 10, 82, 118, 162, 20932.263.081259.2917.120.700.8336.541 × 10−6
8230, 21, 85, 124, 168, 206, 253, 7539.332.54702.3119.650.800.8843.771 × 10−6
10247, 10, 108, 198, 34, 242, 27, 245, 116, 4446.112.851188.7517.370.690.7950.862 × 10−6
2203, 6813.550.524049.6112.050.480.6818.181 × 10−6
4177, 71, 118, 16722.651.531891.4015.360.520.6924.512 × 10−6
LCGSA36169, 110, 191, 78, 133, 18129.131.061663.6215.920.550.7230.092 × 10−6
8117, 187, 94, 175, 156, 93, 184, 8729.901.632033.9615.040.480.6835.902 × 10−6
10132, 85, 124, 167, 66, 96, 124, 145, 177, 11042.203.311581.6616.130.540.6946.352 × 10−6
2135, 25412.930.315615.3410.630.260.5518.311 × 10−6
4251, 15, 251, 1521.601.4020,673.924.970.070.3525.992 × 10−6
LCGSA4699, 174, 253, 38, 245, 1531.272.621369.7116.760.650.7535.371 × 10−6
8218, 10, 115, 249, 32, 250, 12, 15238.862.551464.1216.470.650.7744.161 × 10−6
10212, 42, 246, 62, 172, 251, 44, 251, 32, 24747.703.191011.9218.070.740.8453.143 × 10−6
2102, 17317.871.163052.3313.280.400.6318.301 × 10−6
4130, 251, 25, 13825.411.543740.3512.400.410.5925.962 × 10−6
LCGSA56130, 252, 24, 249, 27, 14531.593.273154.1813.140.440.6136.502 × 10−6
8132, 252, 14, 251, 14, 251, 17, 12539.233.064449.2111.640.380.5944.131 × 10−6
1096, 167, 249, 59, 172, 248, 21, 112, 181, 25338.854.79870.1718.730.710.8051.942 × 10−6
2251, 1614.690.7220,426.785.020.080.3518.311 × 10−6
4130, 252, 42, 16825.181.091644.5115.970.560.6925.982 × 10−6
LCGSA66129, 252, 8, 130, 248, 4331.912.363650.6412.500.490.6635.522 × 10−6
883, 124, 178, 248, 27, 145, 251, 6040.423.14867.2318.740.740.8144.181 × 10−6
10245, 36, 249, 42, 250, 26, 248, 31, 247, 2440.843.5413,970.456.670.200.4153.762 × 10−6
298, 16918.130.513052.1313.280.390.6318.301 × 10−6
496, 176, 253, 1925.141.061947.3815.230.570.7126.862 × 10−6
LCGSA76111, 249, 71, 171, 250, 2232.121.471390.8716.690.660.7736.502 × 10−6
837, 100, 170, 252, 68, 175, 251, 5240.253.521036.6017.970.730.8044.131 × 10−6
10223, 39, 224, 32, 173, 250, 24, 152, 252, 3945.853.561250.9717.080.640.7451.392 × 10−6
296, 17518.170.532889.8213.520.420.6518.301 × 10−6
4102, 172, 254, 3625.360.691555.8316.210.600.7126.011 × 10−6
LCGSA86100, 174, 253, 26, 156, 25129.732.901615.7716.040.570.7036.542 × 10−6
8106, 250, 7, 95, 148, 194, 251, 1240.662.861566.5316.180.640.7844.171 × 10−6
10115, 250, 115, 3, 100, 173, 252, 19, 100, 17349.994.521879.6215.390.560.7053.711 × 10−6
2159, 24712.940.244587.2611.510.300.5818.311 × 10−6
4153, 252, 21, 25320.230.713543.1412.630.380.5725.993 × 10−6
LCGSA96135, 253, 19, 150, 253, 2432.472.363162.7413.130.410.5836.512 × 10−6
8250, 5, 149, 250, 5, 249, 15, 13539.802.623665.1312.480.380.5943.981 × 10−6
10125, 251, 33, 240, 7, 147, 250, 29, 252, 1951.073.122478.0314.180.520.6753.772 × 10−6
297, 17418.160.522962.5513.410.420.6518.313.3164
4148, 253, 19, 25020.150.773779.2612.350.380.5827.0612.2545
LCGSA10692, 144, 194, 253, 71, 17833.943.081335.3716.870.650.7936.5415.0325
835, 100, 166, 251, 19, 120, 249, 3640.803.451555.6616.210.600.7143.8717.9458
10252, 32, 251, 43, 250, 5, 104, 175, 252, 5445.943.141148.6117.520.700.7951.5225.3295
Table 8. Simulation results for the CT2 image using classical, hybrid, and recent HAs.
Table 8. Simulation results for the CT2 image using classical, hybrid, and recent HAs.
AlgorithmkOptimal ThresholdsMeanSTDMSEPSNRSSIMFSIMBest
Value
Run Time
286, 14916.280.023744.0912.390.380.6516.175.1661
4175, 204, 179, 21715.840.153682.3812.460.390.7022.088.8693
GSA6154, 188, 213, 172, 165, 14719.320.103154.1413.140.430.7425.9813.5300
884, 100, 135, 114, 70, 83, 124, 7526.480.303793.9412.330.440.6229.7716.8303
10192, 218, 210, 192, 202, 157, 167, 18, 2, 19232.990.393228.2613.040.450.7433.5620.9540
26, 51.750.8521,156.234.870.030.3211.234.3092
47, 11, 15, 274.140.8816,487.945.950.120.3816.427.2575
PSO616, 20, 19, 42, 76, 8518.291.087405.939.430.470.5821.8412.2000
84, 15, 23, 32, 36, 40, 41, 6625.871.009591.098.310.380.5229.9613.4070
102, 59, 110, 81, 91, 103, 92, 90, 139, 9733.581.422863.7913.560.550.6534.1919.3869
211, 140.480.6719,192.915.290.070.379.103.9591
41, 1, 1, 116.211.7722,712.724.5600.2915.648.2575
PSOGSA616, 30, 42, 51, 58, 6116.360.4610,377.517.960.310.4621.1710.7500
85, 49, 54, 59, 76, 79, 80, 8121.630.997758.789.230.460.5627.8712.8575
101, 5, 23, 37, 79, 99, 124, 129, 130, 13828.990.392776.9413.690.610.6730.7616.8457
21, 313.110.5515,650.866.180.140.378.973.8218
49, 16, 18, 307.411.0815,650.266.180.140.3816.856.7333
CPSOGSA613, 17, 48, 60, 62, 7021.790.769131.618.520.400.5219.2410.9500
841, 65, 78, 80, 79, 112, 113, 16120.670.651283.1217.040.710.7329.7014.0335
1026, 37, 44, 52, 71, 103, 112, 124, 129, 14527.550.362036.9315.040.690.7235.6416.7927
237, 23111.390.6512,557.057.140.220.4215.563.5081
4115, 76, 104, 23818.001.055182.7210.980.410.6523.266.5583
BBO6159, 168, 43, 3, 55, 321.000.921386.6116.710.560.6427.689.2174
8125, 61, 100, 189, 157, 195, 12, 16430.060.55634.3320.100.680.7834.0911.9324
10188, 90, 107, 1, 148, 226, 44, 212, 68, 1531.521.91218.8424.720.880.9336.5714.3068
2217, 18010.802.253992.9612.110.330.6717.103.8954
4183, 111, 137, 22619.251.582749.1713.730.450.7622.946.5492
DE6124, 152, 111, 232, 110, 3625.891.901910.1415.320.620.7230.499.0928
840, 39, 17, 252, 72, 128, 227, 13729.932.452134.6614.830.690.7635.3712.0739
1025, 150, 213, 175, 34, 207, 6, 223, 120, 14337.342.49862.8718.770.690.8041.8114.3418
2231, 113.011.4520,230.265.070.030.4016.856.1295
4138, 246, 35, 10921.641.972838.2513.600.570.6723.2410.5776
SCA6222, 220, 76, 255, 59, 8226.822.235700.8410.570.440.5829.1214.5800
84, 159, 46, 239, 1, 117, 17, 16631.802.621184.7617.390.650.7334.6218.6980
1029, 215, 5, 58, 117, 43, 244, 89, 36, 11540.933.382171.1714.760.670.7342.3322.2553
21, 3911.542.3014,076.236.640.170.3813.905.0792
472, 108, 72, 25511.392.386022.6710.330.390.6216.5210.0753
SSA61, 1, 255, 254, 216, 25524.836.5215,378.186.260.080.4024.2514.2000
8255, 255, 237, 255, 215, 237, 1, 124.046.9814,773.720.430.080.4028.2319.5057
10118, 255, 166, 29, 5, 102, 255, 1, 148, 9228.473.031610.3616.060.600.7137.1122.2910
2253, 412.652.1121,195.234.860.030.3912.363.9982
4255, 254, 254, 20816.012.9710,646.627.850.110.4319.217.1731
MFO63, 6, 108, 84, 255, 25516.543.026445.3110.030.330.6327.338.6488
861, 213, 255, 1, 255, 253, 35, 25032.433.268199.138.990.270.4336.5211.5943
101, 47, 76, 255, 1, 255, 3, 79, 107, 25534.732.386028.6510.320.390.6240.2513.8223
291, 12112.461.745734.2210.540.320.6116.6311.2427
4164, 127, 160, 16817.541.943636.4412.520.310.6020.6518.8235
ABC6102, 151, 182, 132, 176, 13317.902.292778.4813.690.370.6625.4525.7040
8176, 159, 175, 205, 180, 230, 190, 15426.102.563182.4213.100.440.7530.6032.6283
10194, 114, 189, 196, 222, 198, 235, 188, 242, 14630.742.632588.0214.000.480.7836.6739.9652
25, 52.742.2221,667.064.770.020.3110.183.4361
40, 2, 5, 13.432.4621,880.064.730.020.3113.355.4724
GWO61, 0, 1, 2, 3, 00.392.6022,448.304.610.010.2917.368.3194
80, 0, 0, 1, 1, 1, 0, 1, 0, 11.114.6922,979.134.510014.8310.772
101, 3, 0, 1, 0, 1, 0, 1, 1, 64.874.2621,410.644.820.030.3220.4612.8696
21, 14.512.3922,712.724.5600.296.633.8109
438, 37, 36, 384.512.3714,265.636.580.160.376.926.4377
SMA6122, 120, 121, 122, 120, 1214.482.426163.9010.230.270.596.978.7720
8211, 210, 209, 210, 210, 211, 211, 2124.492.5012,457.177.170.090.416.2811.2538
10213, 211, 212, 210, 209, 208, 209, 213, 210, 2114.662.3613,683.726.760.080.406.8913.8051
Table 9. Simulation results for the CT2 image using ten LCGSA versions.
Table 9. Simulation results for the CT2 image using ten LCGSA versions.
AlgorithmkOptimal ThresholdsMeanSTDMSEPSNRSSIMFSIMBest
Value
Run Time
2101, 1817.190.882922.6913.470.380.6717.642 × 10−6
494, 135, 184, 25121.610.952576.3914.020.430.7324.642 × 10−6
LCGSA16249, 28, 100, 188, 90, 18130.301.751098.0717.720.640.7635.091 × 10−6
8101, 181, 252, 23, 103, 177, 248, 7340.392.381092.9817.740.650.7741.641 × 10−6
1097, 190, 76, 201, 15, 107, 187, 79, 199, 7842.092.731247.3317.170.610.7752.721 × 10−6
2100, 17816.571.182992.2713.370.370.6617.641 × 10−6
4101, 182, 252, 8023.511.362300.6814.510.460.7324.561 × 10−6
LCGSA2699, 191, 31, 93, 133, 18432.572.38937.0418.410.650.7433.922 × 10−6
8240, 19, 97, 185, 82, 197, 82, 18336.092.531195.4317.350.640.8141.661 × 10−6
1099, 175, 245, 153, 21, 108, 184, 249, 84, 18145.082.941241.9417.180.620.7751.442 × 10−6
2102, 18811.790.272927.0113.460.380.6917.161 × 10−6
4107, 184, 90, 15521.931.242485.7914.170.420.7025.032 × 10−6
LCGSA36165, 92, 174, 240, 107, 17328.241.532693.8113.820.430.7229.491 × 10−6
8130, 183, 94, 135, 169, 199, 93, 15035.022.972295.9614.520.470.7637.372 × 10−6
1091, 147, 201, 147, 87, 162, 63, 113, 165, 21237.500.91887.3318.640.690.8243.501 × 10−6
2102, 18317.120.962908.9513.490.370.6717.631 × 10−6
480, 111, 150, 18525.271.352176.8314.750.460.7225.762 × 10−6
LCGSA46247, 48, 239, 86, 179, 25330.181.29872.8218.720.700.7833.872 × 10−6
8100, 182, 249, 57, 231, 90, 196, 8036.221.61695.7319.700.710.8341.631 × 10−6
10104, 185, 81, 202, 84, 182, 246, 157, 22, 10340.992.43904.3018.560.720.8551.422 × 10−6
2102, 18317.500.332933.4013.450.370.6717.631 × 10−6
498, 190, 28, 10223.480.911398.1616.670.590.7125.791 × 10−6
LCGSA56101, 187, 28, 98, 193, 3229.871.961159.6917.480.610.7234.111 × 10−6
890, 137, 182, 249, 13, 106, 187, 6639.092.42898.8418.590.660.7941.531 × 10−6
10101, 194, 76, 201, 88, 185, 250, 101, 249, 8642.142.271865.2815.420.530.7850.421 × 10−6
2101, 18017.470.392937.9813.450.370.6717.632 × 10−6
4101, 182, 252, 7923.262.142283.3714.540.470.7325.811 × 10−6
LCGSA6699, 191, 89, 180, 252, 3430.181.79890.5918.630.670.7733.881 × 10−6
8243, 25, 105, 191, 77, 202, 26, 10335.922.20948.1918.360.690.8241.691 × 10−6
1096, 192, 42, 165, 252, 80, 199, 82, 200, 7142.242.18359.5922.570.800.8648.221 × 10−6
2100, 18117.181.362917.9913.480.370.6717.632 × 10−6
4235, 27, 100, 18223.402.001344.5416.840.620.7624.571 × 10−6
LCGSA7698, 180, 249, 96, 179, 25328.871.322815.3713.630.410.7135.232 × 10−6
8100, 188, 43, 231, 82, 199, 79, 19531.141.04593.1520.390.750.8441.441 × 10−6
10232, 33, 103, 189, 22, 104, 182, 250, 173, 8840.343.43805.3319.070.710.8148.541 × 10−6
298, 18217.410.672882.0113.530.380.6717.640
488, 133, 183, 25221.830.812452.7614.230.440.7424.671 × 10−6
LCGSA8698, 184, 253, 25, 103, 18433.882.591458.8916.490.590.7335.211 × 10−6
8247, 53, 235, 23, 98, 138, 33, 25131.332.00686.7919.760.710.7941.681 × 10−6
1094, 140, 183, 251, 83, 200, 28, 100, 192, 8545.202.96775.1019.230.710.8248.641 × 10−6
2104, 18017.460.803006.1913.350.370.6617.642 × 10−6
4104, 181, 253, 3423.921.591221.2417.260.610.7225.830
LCGSA96101, 178, 253, 19, 104, 18034.411.771816.5615.530.540.7135.191 × 10−6
898, 180, 254, 32, 251, 25, 248, 2037.422.661220.0417.260.630.7442.191 × 10−6
1089, 127, 182, 250, 82, 197, 35, 242, 14, 10046.062.50578.0620.510.750.8548.741 × 10−6
2104, 18217.320.982974.3213.390.370.6717.623.4469
4101, 189, 20, 10023.540.75781.9715.620.540.7024.656.1717
LCGSA10697, 191, 14, 102, 183, 25328.691.771785.0515.610.530.7334.098.6448
881, 119, 158, 189, 253, 79, 200, 8237.732.521847.2415.460.540.8041.7811.4366
10248, 91, 181, 248, 86, 183, 248, 19, 101, 18146.723.271489.3916.400.590.7551.4313.9053
Table 10. Simulation results for the CT3 image using classical, hybrid, and recent HAs.
Table 10. Simulation results for the CT3 image using classical, hybrid, and recent HAs.
AlgorithmkOptimal ThresholdsMeanSTDMSEPSNRSSIMFSIMBest
Value
Run Time
2108, 15812.930.063289.2112.950.390.6716.254.7489
4157, 105, 165, 19217.830.042005.0915.100.610.8021.328.4279
GSA6176, 214, 191, 196, 149, 21026.720.272930.6413.460.600.7626.5312.3829
8188, 186, 109, 144, 174, 129, 147, 21024.170.071828.7815.500.690.8230.1221.3364
10142, 155, 108, 136, 166, 119, 115, 135, 54, 15430.630.161518.9116.310.560.7235.1220.8684
214, 148.831.2724,785.004.180.050.3311.384.3324
413, 23, 33, 4516.700.7917,042.085.810.140.3616.817.0285
PSO66, 15, 22, 43, 45, 4321.190.9116,815.095.870.140.3623.9210.3684
82, 12, 44, 52, 69, 61, 65, 9125.741.228351.828.910.380.5830.4915.4691
1018, 32, 66, 71, 76, 126, 93, 122, 109, 10733.281.584161.4611.930.500.6739.3619.0431
224, 315.410.4820,431.315.020.100.349.723.8282
44, 11, 21, 5211.220.5615,305.596.280.180.4115.126.5588
PSOGSA612, 19, 28, 37, 50, 5119.000.2015,496.016.220.160.3821.419.4441
81, 33, 45, 45, 57, 92, 99, 9827.921.357410.029.430.420.6131.5112.6988
1035, 40, 50, 61, 60, 62, 69, 89, 97, 11238.440.585591.0910.650.470.6533.4917.3804
238, 390.400.7618,446.605.470.120.3511.373.9276
41, 1, 1, 113.942.7629,014.483.500.110.3217.597.9724
CPSOGSA67, 12, 30, 31, 43, 7114.651.1411,809.227.400.290.5123.169.5559
82, 17, 39, 90, 95, 114, 112, 12528.180.484487.6211.610.480.6430.0712.0526
104, 21, 35, 37, 45, 60, 89, 92, 95, 11932.440.714953.2011.180.490.6738.5516.8229
2176, 7810.820.462292.9714.520.560.7617.526.6585
428, 160, 145, 20717.260.601744.5015.710.710.7823.5311.2497
BBO6198, 145, 143, 56, 148, 12116.771.28951.9418.340.800.8627.2516.6101
831, 65, 64, 23, 161, 56, 50, 20621.571.07954.0318.330.790.8235.8311.9254
10238, 178, 165, 192, 228, 171, 23, 187, 80, 18933.950.85828.3618.940.820.8938.9015.0964
268, 19512.701.943808.9312.320.530.6917.166.5369
4209, 100, 34, 10121.601.753501.8912.680.640.7624.7111.6915
DE618, 197, 1, 71, 251, 21029.951.483539.3012.640.570.7031.7016.5290
8137, 254, 33, 227, 95, 103, 74, 2734.232.202275.7614.550.720.8339.4011.9878
10236, 134, 1, 111, 200, 63, 73, 166, 17, 1837.612.64441.6421.680.900.9449.2114.7850
2104, 15813.961.553259.6212.990.400.6817.015.3804
4231, 17, 161, 922.642.013194.7913.080.560.7224.209.2408
SCA6174, 6, 40, 222, 1, 4730.232.431990.2115.140.640.7231.5713.5035
859, 141, 212, 88, 161, 60, 186, 1533.852.73431.7521.770.900.9236.2817.8848
10217, 248, 41, 119, 243, 1, 65, 71, 9, 16241.503.22725.2019.520.860.9242.5822.1417
259, 2511.264.9414,106.196.630.220.4610.294.8262
41, 1, 255, 25418.195.9828,355.573.600.010.4118.549.5060
SSA61, 255, 1, 255, 254, 25520.786.4828,355.573.600.010.4124.3314.0087
8255, 1, 255, 254, 255, 254, 1, 19024.367.778543.778.810.2805630.2118.3682
10254, 255, 158, 61, 255, 55, 1, 255, 254, 1728.276.552181.6514.740.560.7334.7222.5992
21, 25511.992.8328,355.573.600.010.4114.865.1625
41, 1, 1, 112.687.0529,014.483.5000.3214.868.8423
MFO639, 1, 255, 1, 1, 25024.162.4117,392.585.720.170.4229.2312.7302
8255, 1, 255, 255, 255, 1, 255, 25526.594.7828,355.573.600.010.4132.8314.2916
101, 4, 92, 255, 1, 255, 3, 98, 15, 25536.922.508293.918.940.360.6740.7513.7680
2108, 17212.981.872742.9413.740.480.7315.9714.0200
4111, 84, 95, 10517.822.076755.959.830.340.6019.7824.5639
ABC6133, 101, 115, 142, 134, 6727.162.203227.6613.040.470.6623.6134.3251
8153, 163, 197, 115, 174, 146, 194, 14922.272.582095.9614.910.590.7731.2842.6460
10189, 110, 189, 190, 220, 198, 235, 188, 240, 14531.502.731962.5615.200.690.8437.1140.1277
20, 25.481.0628,699.383.5500.3210.712.6922
40, 1, 1, 104.972.5225,953.433.980.040.3314.715.6553
GWO60, 0, 1, 1, 0, 110.818.4327,153.853.790.020.3316.477.9962
80, 1, 0, 1, 0, 0, 1, 18.818.8529,014.483.5000.3220.5310.4313
1062, 0, 32, 25, 0, 11, 21, 6, 2, 724.023.5213,302.676.890.240.4725.3012.5223
20, 14.672.5629,331.593.45007.113.6175
420, 19, 21, 204.682.5123,092.374.490.070.337.396.2743
SMA61, 1, 1, 2, 3, 14.672.5428,699.383.5500.327.428.8429
898, 97, 96, 95, 96, 98, 97, 984.512.588697.708.730.290.607.2411.2996
10191, 190, 190, 187, 188, 190, 189, 190, 191, 1914.542.649018.998.570.270.547.3813.9376
Table 11. Simulation results for the CT3 image using ten LCGSA versions.
Table 11. Simulation results for the CT3 image using ten LCGSA versions.
AlgorithmkOptimal ThresholdsMeanSTDMSEPSNRSSIMFSIMBest
Value
Run Time
2253, 2514.650.4621,286.304.840.110.4118.042 × 10−6
4247, 31, 108, 17024.751.431586.7916.120.680.8125.872 × 10−6
LCGSA16223, 23, 157, 254, 29, 16431.892.442145.7714.810.650.7535.933 × 10−6
8246, 21, 104, 168, 253, 42, 253, 2838.891.871327.3916.900.700.8244.411 × 10−6
10140, 253, 40, 234, 47, 249, 25, 102, 171, 25244.152.36894.2218.610.780.8651.581 × 10−6
2252, 3814.370.8917,655.965.660.160.4218.042 × 10−6
488, 127, 170, 21025.851.201512.6816.330.740.8726.522 × 10−6
LCGSA26217, 26, 99, 167, 253, 4631.032.26873.0618.720.830.8933.272 × 10−6
896, 135, 179, 254, 76, 226, 48, 16240.122.87475.6321.350.880.9342.031 × 10−6
1093, 130, 182, 254, 38, 104, 172, 253, 36, 11150.353.94845.4118.860.730.8552.291 × 10−6
2105, 17017.770.452781.5613.680.480.7318.042 × 10−6
4196, 93, 146, 18221.711.911667.1315.910.650.8224.952 × 10−6
LCGSA36118, 169, 223, 86, 130, 16731.612.281621.2316.030.710.8633.492 × 10−6
8124, 177, 97, 165, 86, 128, 165, 20531.412.631301.3116.980.740.8638.561 × 10−6
10103, 133, 165, 205, 151, 98, 135, 171, 211, 8239.973.001169.4117.450.760.8646.381 × 10−6
2253, 1414.620.4824,003.494.320.070.4118.022 × 10−6
4109, 173, 254, 2424.601.421805.8915.560.640.8025.922 × 10−6
LCGSA46102, 173, 253, 36, 151, 24929.882.331176.5617.420.680.8136.012 × 10−6
8224, 40, 243, 20, 159, 252, 29, 24833.442.382000.0015.120.650.7444.051 × 10−6
1096, 171, 251, 33, 108, 216, 42, 160, 251, 4146.544.82674.5019.840.850.9051.011 × 10−6
2103, 17117.681.012718.4613.780.490.7418.042 × 10−6
4250, 28, 253, 3821.780.9917,392.455.720.170.4225.871 × 10−6
LCGSA56139, 250, 14, 146, 253, 3432.241.523154.1913.140.500.6835.111 × 10−6
891, 130, 176, 247, 130, 42, 233, 4335.692.14812.1919.030.830.8944.372 × 10−6
10155, 251, 15, 142, 251, 29, 135, 250, 45, 24443.463.702057.5314.990.580.7251.182 × 10−6
2251, 2714.600.4920,180.095.080.130.4118.053 × 10−6
484, 121, 171, 21225.871.461456.8816.490.750.8826.512 × 10−6
LCGSA66248, 28, 250, 138, 35, 16827.260.821885.1815.370.590.7336.012 × 10−6
8153, 250, 45, 241, 24, 247, 41, 14339.082.912194.8414.710.580.7143.901 × 10−6
10142, 250, 31, 240, 35, 102, 167, 251, 29, 14648.594.101281.3017.050.690.8151.291 × 10−6
2102, 17117.780.862659.4513.880.500.7418.032 × 10−6
4107, 169, 251, 3624.312.311514.5916.320.670.8126.502 × 10−6
LCGSA76103, 171, 252, 43, 163, 25230.151.711280.8617.050.670.8036.022 × 10−6
8223, 31, 248, 38, 238, 31, 130, 25132.992.433099.2513.210.620.7643.682 × 10−6
1084, 126, 172, 249, 38, 98, 148, 190, 249, 4947.335.86413.8021.960.850.9251.161 × 10−6
2166, 25312.950.224681.4511.420.370.6518.042 × 10−6
4253, 31, 159, 25420.020.872976.2913.390.500.6725.911 × 10−6
LCGSA8697, 154, 199, 252, 40, 14434.152.45771.3119.250.840.9035.092 × 10−6
8238, 36, 154, 251, 28, 167, 253, 3339.003.201957.6715.210.600.7243.881 × 10−6
10225, 40, 101, 170, 250, 21, 102, 174, 251, 4148.533.35937.1918.410.810.8851.221 × 10−6
2104, 17217.890.682662.8813.870.500.7418.042 × 10−6
4103, 168, 253, 3624.991.231550.2516.220.660.8126.542 × 10−6
LCGSA96152, 250, 37, 168, 253, 2932.412.182035.4115.040.550.6936.002 × 10−6
8163, 251, 14, 252, 24, 245, 19, 25939.432.812931.4413.460.520.6843.861 × 10−6
10162, 243, 20, 105, 167, 251, 36, 158, 252, 3749.633.911403.4016.650.670.8051.752 × 10−6
2252, 1914.620.5822,704.294.560.090.4118.037.4703
494, 138, 185, 25421.950.731917.9715.300.640.8425.9310.02
LCGSA106123, 250, 27, 108, 167, 25229.661.941683.2615.860.630.7836.0115.4281
8138, 250, 24, 96, 170, 251, 40, 16241.893.691171.3517.440.660.8143.8611.49
10139, 248, 35, 245, 45, 239, 38, 101, 171, 25143.653.09980.0418.210.750.8453.1814.4697
Table 12. Wilcoxon rank-sum test analysis.
Table 12. Wilcoxon rank-sum test analysis.
Grayscale
Image
kLCGSA vs. GSALCGSA vs. PSOLCGSA vs. PSOGSALCGSA vs. CPSOGSALCGSA vs. BBOLCGSA vs. DELCGSA vs. SCALCGSA vs. SSALCGSA vs. MFOLCGSA vs. ABCLCGSA vs. GWOLCGSA vs. SMA
20.01950.00200.00200.00200.01950.00390.00200.00200.00200.01370.00200.0020
40.00390.00200.00200.00200.00200.00200.00190.00200.00200.00200.00200.0020
Airplane60.00200.00200.00200.00200.00200.00200.00200.00200.00200.00200.00200.0020
80.03900.00200.00200.00200.00200.00200.00200.00200.00200.00200.00200.0020
100.00200.00200.00200.00200.00390.00390.02730.00200.00200.00200.00200.0020
20.00980.00390.00200.00200.00590.03710.00220.00390.00200.00200.00200.0020
40.00200.00200.00200.00200.00200.00200.02730.00200.00200.00200.00190.0020
Boat60.00200.00200.00200.00200.00200.00390.01370.00200.00200.00200.00200.0019
80.00200.00200.00200.00200.00200.00200.00390.00200.00200.00200.00230.0020
100.00390.00200.00390.00200.00200.00200.00190.00200.00200.00200.00200.0137
20.00200.00200.00200.00200.00200.00590.00200.00200.00390.02730.00390.0020
40.00200.00200.00200.00200.00390.02730.00200.00200.00200.00200.01270.0020
CT160.00200.00200.00200.00200.00200.00200.00200.00390.00200.00200.00200.0020
80.00390.00200.00200.00200.00200.00200.00200.00200.00200.00200.00200.0020
100.00200.00200.00200.00200.00200.00200.00200.00200.00200.00200.00200.0020
20.03710.00200.00200.00200.00390.01370.00200.00200.00390.02730.00200.0020
40.00390.00200.00200.00200.00200.00200.00390.00200.00200.00390.00200.0020
CT260.00200.00200.00200.00200.00200.00200.00200.00200.00200.00200.00190.0020
80.00200.00200.00200.00200.00200.00200.00200.00200.00200.00200.00200.0019
100.00200.00200.00200.00200.00200.00200.00200.00200.00200.00200.00200.0020
20.03710.00200.00200.00200.00200.00200.00200.00200.00390.00200.00200.0020
40.00200.00200.00200.00200.00200.00590.00200.00200.00200.00200.00200.0020
CT360.00200.00200.01950.00200.00200.00390.00200.00200.00200.00200.00200.0020
80.00200.00200.00200.00200.00200.00200.00200.00200.00200.00200.00200.0020
100.00200.00200.00200.00200.00200.00390.00200.00200.00200.00200.00200.0020
Table 13. Simulation results of the CT-g1 image.
Table 13. Simulation results of the CT-g1 image.
AlgorithmAccuracyPrecisionRecallDice ScoreJaccard IndexRun Time
GSA0.900.8310.900.8316.1872
PSO0.670.5810.730.5814.0780
PSOGSA0.930.8810.930.8811.9169
CPSOGSA0.660.5810.730.5812.7344
BBO0.980.9510.970.9516.5239
DE0.970.9410.950.9417.0990
SCA0.970.9410.970.9416.6560
SSA0.9610.920.960.9217.0693
MFO0.9510.890.940.8916.6198
ABC0.9210.830.910.8347.0522
GWO0.870.7810.880.7816.0203
SMA0.950.9010.950.9017.7029
LCGSA10.9710.950.970.951 × 10−6
LCGSA20.960.9310.960.931 × 10−6
LCGSA30.970.9510.970.951 × 10−6
LCGSA40.9210.830.900.831 × 10−6
LCGSA50.980.9610.980.961 × 10−6
LCGSA60.9210.820.900.822 × 10−6
LCGSA70.980.9710.980.971 × 10−6
LCGSA80.980.9710.980.972 × 10−6
LCGSA90.9710.940.970.941 × 10−6
LCGSA100.9810.960.980.9616.5407
Table 14. Simulation results of the CT-g2 image.
Table 14. Simulation results of the CT-g2 image.
AlgorithmAccuracyPrecisionRecallDice ScoreJaccard IndexRun Time
GSA0.960.9210.960.9312.4194
PSO0.930.8710.930.8711.3754
PSOGSA0.900.8210.900.8210.2045
CPSOGSA0.870.7810.870.7810.6497
BBO0.9010.790.880.7913.3727
DE0.960.9710.950.9913.6886
SCA0.960.9610.970.9815.0199
SSA0.5710.050.090.0513.7934
MFO0.8310.620.760.6214.0207
ABC0.9910.980.990.9839.2531
GWO0.55NaN00013.4122
SMA0.460.4510.620.4518.7029
LCGSA10.960.9310.960.931 × 10−6
LCGSA20.980.9710.980.972 × 10−6
LCGSA30.880.7910.880.791 × 10−6
LCGSA40.8810.740.850.741 × 10−6
LCGSA50.890.8110.890.81‘1 × 10−6
LCGSA60.950.9110.950.911 × 10−6
LCGSA70.960.9310.960.931 × 10−6
LCGSA80.960.9310.960.932 × 10−6
LCGSA90.930.8610.920.862 × 10−6
LCGSA100.910.8410.910.8413.0380
Table 15. Simulation results of the CT-g3 image.
Table 15. Simulation results of the CT-g3 image.
AlgorithmAccuracyPrecisionRecallDice ScoreJaccard IndexRun Time
GSA0.6410.030.050.0313.2421
PSO0.990.9810.990.9812.0374
PSOGSA0.9910.990.990.9910.8915
CPSOGSA0.9810.970.980.9710.8204
BBO0.9710.970.990.9814.1552
DE0.7910.430.600.4313.9826
SCA0.980.9610.980.9613.9312
SSA0.6510.060.110.0614.6093
MFO0.6410.030.060.0315.3080
ABC0.6310.010.020.0140.9452
GWO0.760.6010.750.6013.6832
SMA0.6410.030.060.0314.4526
LCGSA10.980.9710.980.972 × 10−6
LCGSA20.6610.090.170.092 × 10−6
LCGSA30.9810.950.970.951 × 10−6
LCGSA40.990.9710.980.972 × 10−6
LCGSA50.990.9910.990.992 × 10−6
LCGSA60.6610.080.150.081 × 10−6
LCGSA70.990.9910.990.991 × 10−6
LCGSA80.6610.080.150.082 × 10−6
LCGSA90.980.9510.970.952 × 10−6
LCGSA100.6610.080.150.0812.9773
Table 16. Wilcoxon rank-sum test statistical results.
Table 16. Wilcoxon rank-sum test statistical results.
Grayscale ImageLCGSA vs. GSALCGSA vs. PSOLCGSA vs. PSOGSALCGSA vs. CPSOGSALCGSA vs. BBOLCGSA vs. DELCGSA vs. SCALCGSA vs. SSALCGSA vs. MFOLCGSA vs. ABCLCGSA vs. GWOLCGSA vs. SMA
CT-g10.00200.00200.00200.00200.00200.00590.00200.00200.01240.00190.00980.0012
CT-g20.00200.00200.01950.00200.00200.00390.00200.00200.01290.00190.00390.0019
CT-g30.00200.00200.00200.00200.00200.00200.00200.00200.00190.00200.00200.0020
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

Rather, S.A.; Das, S. Levy Flight and Chaos Theory-Based Gravitational Search Algorithm for Image Segmentation. Mathematics 2023, 11, 3913. https://doi.org/10.3390/math11183913

AMA Style

Rather SA, Das S. Levy Flight and Chaos Theory-Based Gravitational Search Algorithm for Image Segmentation. Mathematics. 2023; 11(18):3913. https://doi.org/10.3390/math11183913

Chicago/Turabian Style

Rather, Sajad Ahmad, and Sujit Das. 2023. "Levy Flight and Chaos Theory-Based Gravitational Search Algorithm for Image Segmentation" Mathematics 11, no. 18: 3913. https://doi.org/10.3390/math11183913

APA Style

Rather, S. A., & Das, S. (2023). Levy Flight and Chaos Theory-Based Gravitational Search Algorithm for Image Segmentation. Mathematics, 11(18), 3913. https://doi.org/10.3390/math11183913

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