key: cord-0955496-ry9f19gn
authors: Houssein, Essam H.; Emam, Marwa M.; Ali, Abdelmgeid A.
title: Improved manta ray foraging optimization for multi-level thresholding using COVID-19 CT images
date: 2021-07-07
journal: Neural Comput Appl
DOI: 10.1007/s00521-021-06273-3
sha: 94dc0a1d6bcfdd0fd05b655b86de3661a9ea3fa3
doc_id: 955496
cord_uid: ry9f19gn

Coronavirus disease 2019 (COVID-19) is pervasive worldwide, posing a high risk to people’s safety and health. Many algorithms were developed to identify COVID-19. One way of identifying COVID-19 is by computed tomography (CT) images. Some segmentation methods are proposed to extract regions of interest from COVID-19 CT images to improve the classification. In this paper, an efficient version of the recent manta ray foraging optimization (MRFO) algorithm is proposed based on the oppositionbased learning called the MRFO-OBL algorithm. The original MRFO algorithm can stagnate in local optima and requires further exploration with adequate exploitation. Thus, to improve the population variety in the search space, we applied Opposition-based learning (OBL) in the MRFO’s initialization step. MRFO-OBL algorithm can solve the image segmentation problem using multilevel thresholding. The proposed MRFO-OBL is evaluated using Otsu’s method over the COVID-19 CT images and compared with six meta-heuristic algorithms: sine-cosine algorithm, moth flame optimization, equilibrium optimization, whale optimization algorithm, slap swarm algorithm, and original MRFO algorithm. MRFO-OBL obtained useful and accurate results in quality, consistency, and evaluation matrices, such as peak signal-to-noise ratio and structural similarity index. Eventually, MRFO-OBL obtained more robustness for the segmentation than all other algorithms compared. The experimental results demonstrate that the proposed method outperforms the original MRFO and the other compared algorithms under Otsu’s method for all the used metrics.

The coronavirus disease 2019 (COVID-19) is one of the newest viruses [1, 2] . It had arisen from a virus named SARSCoV-2. The first cases of COVID-19 were recorded in Wuhan, China, in December 2019 [3] . COVID-19 is a deadly pandemic that the world is currently facing. COVID-19 has impacted approximately 11 million globally, including over 500,000 people in Europe and the USA [4] . COVID-19 has become a worldwide pandemic because of an exponential growth rate and poorly understood transmission mechanism [5, 6] . Cough, fever, sore throat, headache, exhaustion, and discomfort in the muscles are the previously reported symptoms of COVID-19 [1] . This virus also contributes to pneumonia, an infection that inflames the lungs' air sacs. COVID-19 influences countries' health, and its effects (e.g., economic and psychological) are crucial. Early detection may be mirrored in early therapy. Accordingly, COVID-19 is an epidemic with a wide range of challenges to tackle. Based on the knowledge mentioned above, there is a need for tools to detect this fatal disease on time. One of the critical problems is using multiple tests to diagnose COVID-19. The real-time polymerase chain reaction is the most commonly used test. It is invasive and time-consuming, and the diagnosis is false-negative. Chest computed tomography (CT) is another test that plays an essential role in COVID-19 diagnosis [7] .

CT offers a guide to pathophysiology, which might shed light on the diagnosis and development of specific disease stages. It evolves as a useful diagnostic method for the medical practice of lung disease associated with COVID-19 [5] . Preliminary studies suggested that chest CT has a high detection sensitivity for COVID-19 lung pathology. Several groups have demonstrated the capacity for diagnosis based on CAD systems, with accuracy as high as 95% [7] . Medical imaging has been used recently for multiple diagnoses of the disease. Medical imaging methods may be used as a critical diagnostic to guide a possible disease.

Image segmentation is considered a critical process in image processing. It plays a significant role in image research in different applications such as computer vision, pattern recognition, and medical images [8, 9] . The classification system may fail if the segmentation result is incorrect. Segmentation aimed to split the image into several homogeneous regions or segments containing similar features, such as texture, color, brightness, contrast, form, and size, regarding a specific thresholding value(s) [10] . Several methods (e.g., edge-based, regionalbased, threshold-based, and feature-based clustering) have been suggested to solve the current problems and improve research quality in the literature.

Thresholding is extremely important in the image segmentation domain. It is categorized into two types: bi-level and multilevel thresholding. The simplest thresholding technique is bi-level thresholding segmentation that divides an image into two classes by searching a single thresholding value. Multilevel thresholding divides the image into several distinct parts by establishing several threshold values [11] . Multilevel thresholding is the most common image segmentation method that attempts to group the pixels that share characteristics into a finite number of pixels and segment an image into more than two regions based on pixel intensity [12] .

Most image segmentation based on multilevel thresholding methods uses the image histogram as an input to determine the thresholds by maximizing or minimizing the objective functions [13] , for example, fuzzy entropy [14] , Kapur's entropy [15] , Tsallis entropy [16] , and Otsu's method [17] .

However, multilevel thresholding segmentation methods have some limitations in computational costs. The computational time will increase for multilevel thresholding problems with increasing thresholding values as they search for the best thresholding values exhaustively to optimize the objective function. There are two types of approaches for finding the appropriate threshold values in multilevel thresholding segmentation: parametric and nonparametric. A parametric approach assumes that each class of the image can identify by using probability density distributions. These classes are all used to represent the pixels in an image. This approach takes a long time to compute. However, the nonparametric methods use discriminated rules to separate the pixels into homogeneous regions. The threshold values are determined by a statistical criterion [18] . Otsu's method [17] presented the maximization of the variance between classes to obtain the optimal threshold values. In this paper, we focused on a parametric approach, in which we measured the specific statistical parameters for each class in the image. To achieve the desired threshold values, we used Otsu's method by maximizing between variances. Otsu suggested a common approach that can determine thresholds to optimize the class variance of the foreground and background intensity levels. It is one of the most suitable threshold selection methods for real-world images [19] . However, there were challenges in these approaches when achieving the desired threshold values for multilevel thresholding. The computational time is also expensive, which grows exponentially with the number of threshold values. The process of searching for optimum thresholding values for multilevel thresholding segmentation can be considered a constrained optimization problem. The optimization algorithms, especially the meta-heuristic methods, are commonly used in related research to solve these problems [20] .

Optimization is the process of determining the best possible solution(s) for a particular problem. Trying to find the fastest road to get to a destination or coordinate work activities to minimize the idle time between each task are examples of optimization problems. Surely, everyone uses optimization daily. Optimization intends to reach a satisfying or optimal solution concerning different objectives for a specific problem. Optimization is received a huge increased concentration in the last decade since it has ubiquitous nature with all issues [21] . It appears everywhere and almost in every application like Computer Science, engineering problems [22, 23] , and Finance, feature selection [24] , fuzzy control systems [25] , chemoinformatics, bioinformatics [26] , engineering optimization problems [27] , improving the collective construction systems [28] , and image segmentation problems [29] . So many fields of science need optimization, different optimization methods have been developed. Optimization algorithms can be defined as processes where possible solutions are evaluated to choose the most suitable solution for the problem at hand. The quality of the optimized solution depends on the chosen algorithm to solve the problem. As a result, one needs to select the most appropriate algorithm to solve the problem at hand [30] . An optimization algorithm may start and perform the optimization process by single or multiple random solutions. In the first case, the optimization process begins with a single random solution, and it is iteratively improved over the iterations. In the second case, a set of solutions is created and enhanced during optimization. Many factors have an impact on the quality of the optimization algorithm, such as the number of features, diversity of training sets, the number of samples, and representativeness [31] . Optimization algorithms can be classified into two categories: deterministic and stochastic methods. Deterministic algorithms may be a good choice when the gradient is accessible. This variety of methods is based on predetermined calculations, implying that the optimization process is replicable. A reproducible process means that the optimization path and the final solution will be identical every time the algorithm is executed with the same initial states. In contrast, stochastic methods have random features that generate different optimization pathways, leading to other optimized solutions, even if the initial conditions are identical at each run. The main advantage of stochastic algorithms is that they do not need gradient information. Gradient information is hard to access due to the problems becoming more and more complex. Researchers have been developing stochastic algorithms in which no gradient information is needed. The development of more advanced stochastic methods has produced the category of metaheuristics algorithms [30] .

Meta-heuristics algorithms have gained enormous interest in multilevel image thresholding [32] . For the past years, researchers have widely demonstrated meta-heuristic algorithms' ability to solve many kinds of difficult optimization problems in engineering, communications, transportation, industry, and social sciences because these algorithms have high performance and are easy to implement [33] . Researchers have proposed meta-heuristic algorithms to solve the multilevel thresholding problems by maximizing some specific criteria for research. Some of the most widespread meta-heuristic algorithms include a genetic algorithm (GA) [34] , cuckoo search (CS) [35] , grey wolf optimization (GWO) [36] , social engineering optimizer (SEO) [37] , differential evolution (DE) [38] , colony optimization (ACO) [39] , sine cosine algorithm (SCA) [40] , Harris hawks optimization (HHO) [41] , and moth flame optimization (MFO) [42] . Zhao et al. [43] introduced manta ray foraging optimization (MRFO) as one of the most recent additions. MRFO showed more high-quality results than several classical and the latest counterparts.

Most optimization algorithms recorded exemplary achievements in multilevel thresholding, as all can get the optimal threshold result; however, they meet particular challenges because of the trapping in local regions, early conversion, and missing global searchability. These lay grounds for researchers to propose modified and hybrid versions and improve techniques. For the optimization problems, evaluating a candidate and its opposite solution simultaneously can increase the convergence rate toward a globally optimal solution. Opposition-based learning (OBL) [44] is one of the healthiest approaches to enhance the search performance of the meta-heuristic algorithms [45] . In [46] , OBL was implemented to increase the efficiency of meta-heuristic optimization methods. OBL has regraded the suggested solutions built through a stochastic repetition system and their ''opposite solutions'' located in the search space's opposite regions. It has been combined with several bioinspired optimizations to provide smaller predicted distances to the globally optimal relative to randomly experimented solution pairs [47] , such as cuckoo optimization algorithm [48] , shuffled complex evolution algorithm [49] , and fireworks algorithm [50] . In [51] , the authors applied the OBL to enhance the Equilibrium Optimizer adopted for image segmentation using Otsu's function to obtain the best threshold values for image segmentation.

A new meta-heuristic swarm optimization algorithm called MRFO algorithm has been applied in structural optimization [43] . Despite its efficacy in search mechanisms, certain areas can be improved to prove its efficiency on challenging optimization problems, like image thresholding, because the algorithm may miss some critical search regions. We integrated the MRFO algorithm with OBL (MRFO-OBL) to generate solutions from potential regions and explore search space more rigorously to address this issue. More importantly, we improved its local searchability with the help of solutions generated from around the promising regions, avoiding traps in local optima. Consequently, the proposed MRFO variant comes with a trade-off equilibrium among exploration and exploitation. We then applied the MRFO-OBL to a multilevel image thresholding problem.

The following are the contributions of this paper:

• Boosted MRFO using OBL to solve image segmentation problems has been proposed. • We applied the proposed segmentation algorithm to COVID-19 CT images. • Different segmentation levels are used to assess the stability and the efficiency of the proposed method. • Several well-known meta-heuristics are compared with the efficiency of the proposed method. • Validate segmentation efficiency according to peak signal-to-noise ratio (PSNR) and structural similarity index (SSIM) values. • The proposed method can be generalized to different medical imaging diagnoses and applied to different benchmark images.

The rest of the paper is structured as follows: Sect. 2 explains some related works. Section 3 confers an overview of the used methods. Section 4 introduces the proposed method. The experimental results derived from the proposed method to address multilevel thresholding are discussed in Sect. 5. Section 6 concludes the proposed method and discusses future work.

The image segmentation problem plays a crucial role in image processing and computer vision, which are present in many fields, such as medical diagnosis, object recognition, satellite image processing, and more. There are more methods used to solve the problem of image segmentation. Multilevel thresholding segmentation is considered the best one for solving the image segmentation problem. Due to the time complexity problem with the increased threshold levels, traditional techniques failed to solve the image segmentation problem. Thus, the meta-heuristic algorithms have been more used to solve these problems and overcome the time complexity. This section introduces an overview of the recent works for overcoming the image segmentation problem of medical imaging and standard imaging and shows some recent works of segmentation of COVID-19 medical imaging. Table 1 summarizes some optimization algorithms used for image segmentation. Moreover, several researchers work on medical imaging segmentation using optimization algorithms; for example, in [52] , the authors proposed a method MRI image segmentation using the LASHED optimization algorithm. The statistically verified results demonstrate that the suggested approach improves consistency and segmentation quality. Furthermore, [53, 54] also using the optimization algorithms for medical image segmentation. The researchers in [55] proposed a multilevel thresholding method for medical image segmentation based on a partitioned and cooperative quantum-behaved PSO. The proposed method was tested with four stomach CT images and compared with two modified PSO algorithms. In [56] , the authors presented an evolutionary grey gradient algorithm-based multilevel thresholding technique for brain MRI image segmentation (EGGA). To improve the fitness function, they used an adaptive swallow swarm optimization (ASSO) technique. They evaluated the ASSO using twenty-five MR scans, which revealed that it performed better than the original SSO. While in [57] , the authors presented an enhanced method of the FPA algorithm to tackle the problem of medical image segmentation. They applied Otsu's method as an objective function and tested the proposed method using eight CT images. This method outperformed several meta-heuristics algorithms (the original FPA, PSO, GA, and DE). In [58] , the authors proposed a segmentation method for MRI brain images using the artificial bee colony algorithm (ABC). They enhanced the ABC algorithm by using the mean bestguided approach. They test the proposed method on 20 MRI images.

Also, the multilevel thresholding has been used for COVID-19 medical imaging segmentation. Therefore, in this paper, some methods for COVID-19 imaging segmentation are presented. In [7] , the authors proposed a segmentation method for CT-COVID-19 images based on marine predators algorithm (MPA) and MFO algorithm using the fuzzy entropy fitness function for multilevel thresholding segmentation. They enhanced the MPA algorithm using the MFO algorithm and used 13 CT images to test its performance. In [59] , the researchers presented a multilevel thresholding segmentation method for COVID-19 chest X-ray images by hybridized two metaheuristics algorithms (Slime mould algorithm and WOA algorithm). They applied the Kapur method as an objective function and compared it with five optimization algorithms. Also, in [60] , the authors applied the MPA algorithm to diagnose the COVID-19 X-ray images and obtain a good performance. While the authors in [61] proposed a model to predict the COVID-19 using the machine learning techniques with the nature-inspired algorithms. The proposed method was compared against more existing methods and tested on the same datasets.

As conferred from the previous researches mentioned above, image segmentation of medical images, especially in COVID-19 images, is critical to detect the infection with COVID-19.

This section discusses the materials and methods needed for developing the proposed method. The basic MRFO algorithm and its structure are clarified, and some fundamental concepts of the OBL strategy and the objective function are described.

This subsection explains the standard method's mathematical model used in image thresholding; namely, Otsu's method [17] . The segmentation method is based on the image histogram [75] , which carries the pixel's distribution over the image. This method takes the image histogram as the input and then finds the optimal threshold values to its criteria. Otsu's method is one of the segmentation methods used to achieve the image's best threshold values by maximizing the variance between classes and determining the optimum values for the thresholds that divide the image into many classes. This method identifies the gray image intensity levels L v and determines the probability distribution in Eq. (1) [8] .

where i l is the intensity level specified in ð0 i l L v À 1Þ, PN represents the total number of pixels, h i denotes the intensity number i l in the image histogram. The histogram Year Algorithm Technique [62] 2017 MFO-WOA Both MFO and WOA algorithms were used for multilevel thresholding segmentation. The proposed method used Otsu's as the fitness function and tested both WOA and MFO using several images [63] 2020 ABC-SCA In this method, a hybrid of the ABC algorithm and the SCA algorithm was proposed for multilevel thresholding image segmentation. The SCA is employed as a local search for the ABC to boost its performance. This model obtains good performances compared to several existing meta-heuristics methods [10] 2020 HHO The HHO algorithm is used for image segmentation and applied the minimum crossentropy as a fitness function. The performance of the algorithm has been tested in standard images and digital mammograms. The proposed method is verified based on other comparable optimizers and two machine learning algorithms (K-means and the Fuzzy IterAg) [64] 2015 FFO The FFO algorithm has been proposed to maximize Otsu's variance to solve timeconsuming and low-accuracy problems in multilevel thresholding image segmentation [65] 2020 EO The EO algorithm was used to find the optimal threshold value for a grayscale image and applied the Kapur entropy as a fitness function. The performance of this algorithm is compared with seven other algorithms [66] 2016 CS This paper introduced the comparative performance study of different objective functions using cuckoo search and other optimization algorithms to solve the color image segmentation problem using Otsu or Kapur's method [67] 2018 ABC This method presented an Otsu segmentation method based on the ABC algorithm [68] 2020 PSO This technique was used to segment the color images [32] 2019 WOA-GWO-PSO This method used three meta-heuristics algorithms for multilevel thresholding image segmentation to maximize the Otsu method. It tested on 20 benchmark test images using six different thresholds [69] 2018 Firefly algorithm (FA) This is a technique for multilevel color image thresholding used the fuzzy entropy as a fitness function and enhanced the FA algorithm by Levy flight (LF) strategy [70] 2020 PSO This paper proposed a non-revisiting quantum-behaved PSO (NrQPSO) algorithm to find the optimal multilevel thresholds for gray-level images using Kapur's entropy as an objective function [71] 2020 Teaching learning based optimization algorithm (TLBO)

In this paper, LebTLBO was applied on ten standard test images and used the Otsu and Kapur's entropy objective functions for image segmentation and compared with the MTEMO, GA, PSO, and BF algorithms for both Otsu and Kapur's entropy methods.

The results demonstrated that the LebTLBO outperforms the compared algorithms [72] 2020 DE This paper proposed a beta differential evolution (BDE)-based fast color image multilevel thresholding method using two objective functions (Kapur's and Tsallis entropy). The efficiency of the proposed method is examined over existing multilevel thresholding methods such as artificial bee colony, particle swarm optimization, wind-driven optimization, and differential evolution [73] 2021 Black Widow optimization algorithm (BWO)

The BWO algorithm used to find the best threshold value to solve the multilevel thresholding image segmentation and used both Otsu and Kapur methods as objective functions [74] 2019 Elephant herding optimization (EHO)

The authors enhance the EHO algorithm by the OBL and dynamic Cauchy mutation (DCM) to solve the multilevel image thresholding problem for image segmentation by maximizes two objective functions: Kapur's entropy and between-class variance Neural Computing and Applications has standardized in Ph i . Depending on the probability distribution or thresholding value (th), the classes for bilevel segmentation have determined as:

where x 0 ðthÞ and x 1 ðthÞ are additive probability distributions for C 1 and C 2 , respectively, as shown in Eq. (3).

It is mandatory to find the average intensity levels l 0 and l 1 using Eq. (4). When these values have been calculated, the Otsu-based between-class r 2 B is defined in Eq. (5).

In Eq. (5), r 1 and r 2 are the C 1 variances specified as follows:

where l T ¼ x 0 l 0 þ x 1 l 1 and x 0 þ x 1 ¼ 1 based on r 1 and r 2 . Equation (7) defines the objective function:

where r 2 B ðthÞ is the Otsu's variance for a specified th value. The objective function F otsu ðthÞ in Eq. (7) is updated for multiple thresholds as:

where TH=½th 1 ; th 2 ; . . .; th qÀ1 is a vector containing multiple thresholds and L denotes maximum gray level. The variances are calculated by Eq. (9).

where i defines a particular class and x i and l j are the probability of occurrence and the mean of a level. For multilevel thresholding, such values are obtained as follows:

for mean values:

OBL is a useful search strategy in bypassing stagnancy in candidate solutions [44] . Tizhoosh [46] proposed the OBL's basic idea, which improves the search mechanism's exploitation worth. Generally, when the initial solutions are near the optimal location in meta-heuristic algorithms, convergence happens instantly; otherwise, late convergence is expected. Here, the OBL strategy obtained better solutions by considering opposite search regions closer to the global optimum. The OBL works by searching both directions in the search space. One of the original solutions is used in these two ways, whereas the opposite solution defines the other direction. The OBL then takes the most suitable solution from all solutions [76] .

• Opposition number The idea of opposite numbers can be described by explaining OBL. An opposition-based number can be identified as follows: Consider z 0 is a real number on an interval: z 0 2 ½q; w. The opposite number of Z 0 is defined by Eq. (12) [76]

Moreover, the opposite number in D-dimensional space is determined by Eqs. (13) and (14) z ¼ z 1 ; z 2 ; z 3 ; . . .. . .; z D ð13Þ

The values of all items in z are calculated using Eq. (15)

• Opposition-based optimization In the optimization process, the opposite point z 0 is substituted by the corresponding solution z 0 as regards the fitness function. If f ðz 0 Þ is better than f ðz 0 Þ, then z 0 does not alter, oppositely, z 0 = z 0 . So, the solutions are updated as regards the best value of z and z [77].

Zhao et al. [43] designed the MRFO algorithm, a bioinspired optimizer. MRFO assumes the actions of manta rays in catching the prey. MRFO utilizes three foraging procedures of manta rays: chain foraging, cyclone foraging, and somersault foraging. Chain foraging mimics the process of necessary food searching. Foraging manta rays systematically grasp up to capture the disappeared or undetected prey in the chain by the last manta ray. Cyclone foraging occurs when the collection amount of the prey is significant. The head is paired with the manta ray's tail, making a spiral to create an edge in a cyclone's eye. In somersault foraging, manta rays perform backward rotation and circle movements throughout as the prey planktons move them into their open lips [78] . Like various metaheuristic algorithms, MRFO's initialization start is defined randomly to generate random positions for a set of agents X, followed by getting the most desirable agent X best with the best fitness value. The agents are updated according to the aforementioned three strategies. The following are the mathematical models of these strategies.

1. Chain foraging In this process, the agent X k is updated at iteration k using Eq. (16)

where x t k is the k th position in time k, r is a random vector in [0, 1], a is the weight coefficient, and x t best is the high-concentration plankton [43] . 2. Cyclone foraging In this process, the agent X k is updated following Eq. (18)

where b is the weight coefficient. Furthermore, based on a random position provided in the search space, the cyclone foraging step agents can change their position to improve the MRFO's exploration. Thus, the current agent location is updated using the following equation:

where Low and Up are the lower and upper limits of the search space. 3. Somersault foraging The mathematical formulation used to update the agent in this phase is defined in the following equation:

where s is the somersault factor equal to 2, and r 1 and r 2 are the random numbers in [0, 1].

To conclude, the MRFO algorithm is initialized by generating a random population within specified permissible boundaries. The position-update procedure depends on the individual manta ray at the front of the recent one and the estimated pivot position. Changing from exploration to exploitation phase depends on the (itr/maxitr) ratio value. The exploitation phase is determined when ðitr=maxitr
Þ, in which the current most suitable position is recognized as a pivot position. The algorithm moves to the exploration phase when ðitr=maxitr [ rÞ. Furthermore, the algorithm can shift between chain foraging and cyclone foraging based on a randomly generated number. Somersault foraging then takes steps to update the individuals' current status through the current best approach. These three distinct foraging mechanisms are conducted interchangeably to simultaneously achieve the optimization problem's optimal global solution and meet the predefined end criterion.

This section discusses in detail the implementation of the proposed MRFO-OBL algorithm for segmenting the CT chest for COVID-19. The MRFO-OBL model is illustrated in detail in the flowchart presented in Fig. 1 . In the proposed method, the OBL is used to boost the convergence of the MRFO toward the global solution. OBL is a local search technique aimed to avoid the random population's disadvantages and increase the algorithm's convergence by enhancing the variety of its solutions. Consequently, the smaller steps enable the complete quest to scan the promising area rigorously. The proposed MRFO-OBL multilevel thresholding image segmentation method has explained with more details in the following phases:

1. Initialization phase The first step in the proposed method is to read the input image from the COVID-19 dataset and store it as the gray-scale image. It then obtains the input image histogram and calculates the probability distribution following Eq. (1). The next step initializes the MRFO-OBL parameters, such as maximum iterations it max , population size N, and problem dimensions D. Like many other meta-heuristic algorithms, MRFO-OBL starts with the first population's random initialization x 0 and saving results. The OBL strategy can be applied to the MRFO in the initialization phase to improve the search process as follows:

where Opp i is a vector maintaining solution generated by applying OBL and lb j and ub j are the lower and upper bounds of the j t h component of a vector X. 

It is important to check the pixel classification's accuracy for multilevel segmentation. The metrics used in this paper to measure the quality of the segmented images include the peak signal-to-noise ratio (PSNR) and the structure similarity index (SSIM). These metrics are defined briefly as: PSNR It is applied to validate the similarities among the original and segmented images using the Root Mean Square Error (RMSE) [79] of each pixel. PSNR and RMSE are represented as:

where I og is the original image, I segg is the segmented image, while r ro and c co are the maximum numbers of rows and columns of the images, respectively. SSIM It is another measure of the similarity among the original and segmented image [80] . It is identifiable as follows:

where l I and l segg are the mean intensities of the original image I og and the segmented image I segg , r I and r segg are the standard deviations of I og and I seg . r I og ;I segg is the covariance of I og and I segg , and c 1 and c 2 are constants. Moreover, to validate the quality of threshold values, the fitness value is used. To statistically illustrate the significance of the proposed method, the Wilcoxon rank-sum test has been utilized to determine the difference between MRFO-OBL and all other methods used in the comparison [81] . It is a nonparametric test to compare the outcomes of each pair of methods [82] . This test is based on two assumptions: (i) null assumption, which designates no difference between the ranks of a pair of algorithms extracted from the results, and (ii) alternative assumption, which affirms a discrepancy between the results' ranks produced by a pair of algorithms. It is based on a 5% significant level. According to the P values, if ðP [ 0:05Þ, then the null assumption is true. If ðP:05Þ, then the alternative assumption is true. 

Return the best solution, which contains the best thresholds, and apply them to the CT images.

In this section, the experiments have conducted to verify the proposed MRFO-OBL method's performance and compare their performance with some of the state-of-theart algorithms and the original MRFO algorithm to tackle the image segmentation problem. The experiments were conducted on a set of chest CT COVID-19 images, namely CT-image1, CT-image2, CT-image3, CT-image4, CT-im-age5, CT-image6, CT-image7, CT-image8, CT-image9, and CT-image10 as shown in Fig. 3 . The proposed method was evaluated in terms of the fitness function, also to assess the quality of the segmented image, we used a set of performance metrics, including the PSNR and SSIM. The reminder of this section is organized as follows: Sect. 5.1 shows the dataset used in the experiments. Section 5.2 explains the environmental setup. Section 5.3 describes the experimental settings. Section 5.4 shows Otsu's fitness, PSNR, SSIM, and Wilcoxon rank-sum test results.

In this paper, computed tomography (CT) images for COVID-19 from [83, 84] datasets used to test the proposed method. The CT COVID-19 dataset has 349 CT images containing clinical findings of COVID-19 from 216 patients. The COVID-19 images were collected from patients with ages ranging from 40 to 84 of both genders. Figure 2 shows case study of samples of the CT COVID-19 dataset and metadata for each sample. The proposed method evaluates ten test images from this dataset for different patients to test the proposed method's performance. The test images are labeled as CT-image1, CT-image2, ..., CT-image10. Figure 3 displays the selected test images and their corresponding histograms.

The efficiency of the proposed method is compared to six meta-heuristic algorithms which are MFO [42] , whale optimization algorithm (WOA) [85] , SCA [40] , slap swarm algorithm (SSA) [86] , equilibrium optimization (EO) [87] , and the original MRFO [43] . The COVID-19 CT images are segmented with 7, 8, 9, and 10 thresholds. Also, the proposed method's assessment is achieved according to fitness, PSNR, and SSIM. Both algorithms have been coded and run on MATLAB R2016b on a computer running Windows 10 (64 bit), Intel Core I7, and 8 GB RAM. Table 2 presents the default parameters values of each algorithm that are applied in the experiments. As mentioned in [88] , default parameter values is a fair parametrization. Moreover, employing default values reduce comparison bias risks as no algorithm could be advantaged with a better parametrization. To ensure fair benchmarking comparison, the proposed and counterparts algorithms are evaluated over 30 independent runs with a maximum number of 500 iterations for each test image, and the population size is set to 30.

This subsection displays and discusses the proposed method results. The Otsu's method in Eq. (8) has been used as an objective function. Tables 3 and 4 show the segmented COVID-19 CT images obtained from the MRFO-OBL with different numbers of thresholds nTh = 7, 8, 9, 10 outlined over histograms. Table 5 shows the gained threshold values at level 7 outlined over the CT-image3 histogram and the segmented images at the same level acquired by the compared algorithms. Table 6 shows the gained threshold values at level 9 outlined over the CT-image10 histogram and the segmented images at the same level acquired by the compared algorithms.

Tables 7, 8, and 9 include the mean results of fitness, PSNR, and SSIM evaluation matrices, respectively. A higher value means a more reliable and more effective algorithm. In summary, the following observation from the experiments are worth mentioning:

• In terms of fitness values Table 7 shows the fitness values for each level in the comparison of MRFO-OBL with all other algorithms. The best values are presented in bold. Higher fitness values were achieved using the MRFO-OBL and MRFO algorithms, but one gave better results than the other for a set of test images. For CT-image6, CT-image9, and CT-image10 images, the MRFO-OBL algorithm values are the best for all or at least three levels. Moreover, for CT-image1, CT-image4, CT-image5, and CT-image8, MRFO-OBL has higher values for only two levels, while the CT-image3 image obtained only one higher value for level 9. In MRFO, the CT-image7 image acquired higher values for each level, and the CT-image2 image acquired higher values for three levels (level 8, 9, and 10). Also, for CT-image4, CT-image5, and CT-image8, the MRFO algorithm has higher fitness values at only two levels comparing with other algorithms, while for images CT-image1, CT-image3, and CT-image9 has obtained higher values at only one level. The best values were acquired for level 7 and level 8 at only three test images in the SSA algorithm. On the other hand, the EO, WOA, SCA, and MFO algorithms have not produced any best value in either image. In general, the MRFO-OBL and MRFO outperform all algorithms. • In terms of PSNR values Table 8 shows the algorithms' PSNR values. It provides an estimate of the similarity between the segmented image and the original, as mentioned above, where a higher value represents better quality segmentation. Higher PSNR values were generally achieved using the MRFO-OBL and EO algorithms, but one gave improved results than the other for a set of test images. In MRFO-OBL, the PSNR obtained for the CT-image2 image gives higher values for all levels than the compared algorithms, while the CT-image3 and CT-image4 images show the best results compared to the other algorithms for all levels except level 10. Also, in CT-image1 and CT-image10, it has higher values for only two levels (level 9, level 10), and (level 8, level 9), respectively. For CT-image7 and CT-image9, the MRFO-OBL has only one higher value at level 10 and level 9, respectively. In EO, the images CT-image6, CT-image7, and CT-image9 acquired higher values for three levels, while the images CT-image1 and CT-image8 have higher values for only two levels (level 7, level 8) and (level 8, level 10), respectively. In CT-image3, CT-image4, and CT-image10, the EO has higher values than other comparable algorithms at only one level. In WOA, the CT-image5 image acquired for all levels than the compared algorithms, while in CT-image10 has only one higher value at level 10. The best PSNR values were achieved in the SSA algorithm for two levels only in the CT-image8 image. In the MRFO algorithm, the higher PSNR value is achieved only one at level 8 for image CT-image6. The MFO and SCA algorithms do not give any high PSNR results in any of the images. As a result, the highest PSNR values are generated through the MRFO-OBL and EO algorithms, displaying the most tested images' greatest efficiency, indicating a great affiliation between the original and the segmented image. • In terms of SSIM values According to the SSIM values in Table 9 , it was perceived that the obtained SSIM values in the MRFO-OBL are better than those all algorithms compared, implying significance in most images, particularly in the CT-image1, CT-image2, CT-image3, CT-image4, CT-image6, CT-image7, and CT-image10. In the EO algorithm, the SSIM higher values were obtained for image CT-image9, while for image CT-image5, the SSIM higher values were achieved by the WOA algorithm. The MFO algorithm has only one high SSIM value for CT-image4 image at level 10. The MRFO and SCA algorithms do not have any higher SSIM values at any test image.

The following points can be observed from this analysis: Table 10 shows the Wilcoxon rank-sum test results for the best values applied on PSNR between MRFO-OBL and each counterpart. As mentioned above, according to the P values, if ðP [ 0:05Þ, then the null assumption is true. If ðP:05Þ, then the alternative assumption is true. To clarify the discussion of these values, the symbols ?? and -are involved. The þþ symbol means a significant Image nTh = 7 nTh = 8 nTh = 9 nTh = 10 Particularly, the major outcomes of the results reported in Table 10 are summarized as follows: the proposed MRFO-OBL method has a better quality for segmentation.

With the spread of COVID-19 worldwide since December 2019, the entire world has moved to find techniques to help distinguish infected persons from normal ones. After many efforts and confirmation by the medical specialists, CT images could significantly identify whether the suspected patients have been infected. Various segmentation methods are applied to extract regions of interest from CT images essential to improve the classification methods. Medical image segmentation considers a crucial step for many medical applications that need to be done correctly for practical image analysis. One of the most primary and essential techniques for image segmentation is thresholding. In this paper, finding the optimum thresholding values in multilevel thresholding image segmentation was considered an optimization problem where Otsu's method has been used as an objective function. So, a new enhanced version of the MRFO algorithm has been introduced to solve this problem. This algorithm aims to determine the best threshold values that maximize Otsu's function. This method depends on the OBL strategy to improve MRFO's ability to reach the optimum threshold value. The proposed method has been used for COVID-19 CT image segmentation. The performance evaluation of MRFO-OBL is evaluated using 10 COVID-19 CT images with threshold numbers nTh from 7 to 10 and compared to six metaheuristic algorithms: MFO, WOA, SCA, SSA, EO, and the original MRFO. The proposed method's performance has been assessed based on three measures, the best fitness values, PSNR, and SSIM metrics. The results showed that the proposed method gets good results compared to the other competed algorithms. In the future, we will extend the proposed method to test in other applications such as feature selection and color image segmentation and increase the number of thresholds to attain more reliable results. We will also hybridize the proposed method with different meta-heuristic optimization algorithms to improve the segmentation results and applied more activation functions such as Kapur entropy and Fuzzy entropy. Table 5 Segmented image and its histograms acquired by the compared algorithms over CT-image3 Table 6 Segmented image and its histograms acquired by the compared algorithms over CT-image10 

Conflict of interest The authors declare that there is no conflict of interest.

Ethical standard This article does not contain any studies with human participants or animals performed by any of the authors.

CT-image10

836E-02 -8 6.308E-01 -9

CT-image2 7 5.741E-01 -1.650E-03 ?? 2.993E-11 ?? 5.640E

972E-03 ?? 7.585E-09 ?? 3.020E-11 ?? 2.277E-05 ?? 1.679E-03 ?? 8.760E-02 -10 1.292E-05 ?? 1.404E-04 ?? 3.018E-11 ?? 1.428E-08 ?? 6.520E

CT-image3 7 6.794E-03 ?? 8.026E-04 ?? 3.012E-11 ?? 3.767E-02 ?? 7.900E-01 -6.350E-01 -8 1.202E-01 -7.031E-05 ?? 2.977E-11 ?? 1.551E-05 ?? 6.721E

059E-02 ?? 1.044E-09 ?? 3.018E-11 ?? 2.952E-05 ?? 7.007E-01 -2.736E-03 ?? 10 5.827E-03 ?? 7.187E-02 -3.016E-11 ?? 1.805E-05 ?? 5.554E-02 -7.534E-02

CT-image4 7 6.149E-01 -1.132E-04 ?? 2.403E-09 ?? 1.426E-04 ?? 6.650E-02 -3.135E-02 ?? 8 1.205E-01 -9

709E-05 ?? 8.110E-09 ?? 3.020E-11 ?? 9.188E-06 ?? 7.897E-05 ?? 2.704E-02 ?? 10 2

469E-02 -6

CT-image6 7 2.064E-02 ?? 2.342E-05 ?? 6.503E-07 ?? 7.935E-06 ?? 9.646E-01 -2.911E-05

318E-03 ?? 9.472E-02 -3.508E-07 ?? 3.083E-02 ?? 4.732E-01 -2.156E-01 -9 1.958E-01 -1.698E-04 ?? 9.243E-09 ?? 1.070E-01 -1

CT-image7 7 5.150E-01 -1.134E-02 ?? 5.070E-03 ?? 3.535E-02 ?? 3.946E-01 -5.417E-02 -8 9.882E-01 -7.945E-03 ?? 2.274E-05 ?? 3.631E-01 -5.249E-01 -6.194E-01 -9 7.974E-02 -1.122E-02 ?? 1.056E-03 ?? 1.909E-02 ?? 9.234E-01 -1.215E-03 ?? 10 1.857E-03 ?? 1.669E-01 -5

CT-image8 7 9.528E-01 -2.573E-03 ?? 2.911E-09 ?? 1.750E-02 ?? 3.152E-03 ?? 1.274E-03

3.310E-03 ?? 6.244E-02 -1.409E-01 -9 5.541E-01 -8.994E-01 -2.958E-11 ?? 7.103E-06 ?? 9.166E-03 ?? 2.889E-01 -10 1.659E-02 ?? 1.409E-02 ?? 2.999E-11 ?? 8.072E-01 -1.554E-03 ?? 5.486E-01

CT-image9 7 7.429E-01 -2.839E-02 ?? 5.191E-10 ?? 9.215E-02 -3.057E-02 ?? 7.297E-01 -8 7.124E-02 -2.136E-07 ?? 8.338E-08 ?? 1.131E-04 ?? 2.035E-01 -1.129E-03 ?? 9 1.785E-01 -2.051E-03 ?? 1.359E-07 ?? 2.951E-05 ?? 6.458E-02 -6

CT-image10 7 1.082E-01 -8.764E-08 ?? 2.571E-05 ?? 2.278E-06 ?? 4.460E

005E-02 ?? 7.677E-03 ?? 1.976E-06 ?? 1.063E-03 ?? 3.891E-02 ?? 3.125E-01 -9 1.051E-02 ?? 6.202E-01 -1.871E-07 ?? 4.840E-04 ?? 7.726E

Detection of coronavirus (covid-19) associated pneumonia based on generative adversarial networks and a fine-tuned deep transfer learning model using chest X-ray dataset

A review on covid-19 forecasting models

Covid-19 X-ray images classification based on enhanced fractional-order cuckoo search optimizer using heavytailed distributions

Perspectives on the definition of data visualization: a mapping study and discussion on coronavirus (covid-19) dataset

Artificial intelligence for the detection of covid-19 pneumonia on chest CT using multinational datasets

Krishnamurthi R (2021) Dynamic models and control techniques for drone delivery of medications and other healthcare items in covid-19 hotspots

An improved marine predators algorithm with fuzzy entropy for multilevel thresholding: real world example of covid-19 CT image segmentation

A novel black widow optimization algorithm for multilevel thresholding image segmentation

Efficient solution of otsu multilevel image thresholding: a comparative study

An efficient Harris Hawks-inspired image segmentation method

An efficient krill herd algorithm for color image multilevel thresholding segmentation problem

A local fuzzy thresholding methodology for multiregion image segmentation

Image thresholding segmentation based on a novel beta differential evolution approach

Fuzzy entropy and conditioning

A new method for graylevel picture thresholding using the entropy of the histogram

Moment-preserving thresolding: a new approach

A threshold selection method from gray-level histograms

Image segmentation by minimum cross entropy using evolutionary methods

A survey of thresholding techniques

Toward an optimal kernel extreme learning machine using a chaotic moth-flame optimization strategy with applications in medical diagnoses

An enhanced opposition-based salp swarm algorithm for global optimization and engineering problems

Henry gas solubility optimization: a novel physicsbased algorithm

Archimedes optimization algorithm: a new metaheuristic algorithm for solving optimization problems

Enhanced Harris Hawks optimization with genetic operators for selection chemical descriptors and compounds activities

Szedlak-Stinean A-I (2021) Slime mould algorithm-based tuning of cost-effective fuzzy controllers for servo systems

A modified henry gas solubility optimization for solving motif discovery problem

Lévy flight distribution: a new metaheuristic algorithm for solving engineering optimization problems

A hybrid swarm algorithm for collective construction of 3d structures

Hybrid sine cosine artificial bee colony algorithm for global optimization and image segmentation

Hybrid slime mould algorithm with adaptive guided differential evolution algorithm for combinatorial and global optimization problems

A BA-based algorithm for parameter optimization of support vector machine

A new heuristic for multilevel thresholding of images

Introduction to optimization and metaheuristic methods. In: Recent metaheuristics algorithms for parameter identification

Genetic algorithms

Cuckoo search via lévy flights

Grey wolf optimizer

The social engineering optimizer (SEO)

Differential evolution-a simple and efficient heuristic for global optimization over continuous spaces

Ant colony optimization

SCA: a sine cosine algorithm for solving optimization problems

Harris hawks optimization: algorithm and applications

Moth-flame optimization algorithm: a novel nature-inspired heuristic paradigm

Manta Ray foraging optimization: an effective bio-inspired optimizer for engineering applications

Local search in combinatorial optimization

A survey and classification of opposition-based metaheuristics

Opposition-based learning: a new scheme for machine intelligence

How to speed up optimization? Opposite-center learning and its application to differential evolution

A cuckoo optimization algorithm using elite opposition-based learning and chaotic disturbance

A shuffled complex evolution algorithm with opposition-based learning for a permutation flow shop scheduling problem

Opposition-based adaptive fireworks algorithm

Oppositionbased Laplacian equilibrium optimizer with application in image segmentation using multilevel thresholding

Improving the segmentation of magnetic resonance brain images using the lshade optimization algorithm

Paip 2019: liver cancer segmentation challenge

A novel reformed histogram equalization based medical image contrast enhancement using krill herd optimization

Partitioned-cooperative quantum-behaved particle swarm optimization based on multilevel thresholding applied to medical image segmentation

An evolutionary gray gradient algorithm for multilevel thresholding of brain MR images using soft computing techniques

A hybrid flower pollination algorithm based modified randomized location for multi-threshold medical image segmentation

An improved artificial bee colony algorithm based on mean best-guided approach for continuous optimization problems and real brain MRI images segmentation

Hsma\_woa: a hybrid novel slime mould algorithm with whale optimization algorithm for tackling the image segmentation problem of chest X-ray images

Covid-19 image classification using deep features and fractional-order marine predators algorithm

Covid-19 cases prediction by using hybrid machine learning and beetle antennae search approach

Whale optimization algorithm and moth-flame optimization for multilevel thresholding image segmentation

Improved artificial bee colony using sine-cosine algorithm for multi-level thresholding image segmentation

A method of twodimensional otsu image threshold segmentation based on improved firefly algorithm

A novel equilibrium optimization algorithm for multi-thresholding image segmentation problems

A novel color image multilevel thresholding based segmentation using nature inspired optimization algorithms

A multilevel thresholding image segmentation based on an improved artificial bee colony algorithm

A multimodal particle swarm optimization-based approach for image segmentation

A new technique for multilevel color image thresholding based on modified fuzzy entropy and Lévy flight firefly algorithm

A non-revisiting quantum-behaved particle swarm optimization based multilevel thresholding for image segmentation

A multilevel thresholding algorithm using LebTLBO for image segmentation

A novel beta differential evolution algorithm-based fast multilevel thresholding for color image segmentation

A novel black widow optimization algorithm for multilevel thresholding image segmentation

Neural Computing and Applications

Oppositional elephant herding optimization with dynamic Cauchy mutation for multilevel image thresholding

An analysis of histogram-based thresholding algorithms

Improved salp swarm algorithm based on opposition based learning and novel local search algorithm for feature selection

An improved oppositionbased sine cosine algorithm for global optimization

A novel chaotic Manta-Ray foraging optimization algorithm for thermo-economic design optimization of an air-fin cooler

A study on particle swarm optimization and artificial bee colony algorithms for multilevel thresholding

Image quality assessment: from error visibility to structural similarity

Individual comparisons by ranking methods

Gene selection using Wilcoxon rank sum test and support vector machine for cancer classification

Covid-CT-dataset: a CT scan dataset about covid-19

Covid-19 image data collection: prospective predictions are the future

The whale optimization algorithm

Salp swarm algorithm: a bio-inspired optimizer for engineering design problems

Equilibrium optimizer: a novel optimization algorithm

Parameter tuning or default values? An empirical investigation in search-based software engineering