key: cord-0881282-z2r2ki9n authors: Abualigah, Laith; Al-Okbi, Nada Khalil; Elaziz, Mohamed Abd; Houssein, Essam H. title: Boosting Marine Predators Algorithm by Salp Swarm Algorithm for Multilevel Thresholding Image Segmentation date: 2022-03-03 journal: Multimed Tools Appl DOI: 10.1007/s11042-022-12001-3 sha: 88b03bc7aaf74322856a4e3d4454b7e9f7b3e58c doc_id: 881282 cord_uid: z2r2ki9n Pixel rating is considered one of the commonly used critical factors in digital image processing that depends on intensity. It is used to determine the optimal image segmentation threshold. In recent years, the optimum threshold has been selected with great interest due to its many applications. Several methods have been used to find the optimum threshold, including the Otsu and Kapur methods. These methods are appropriate and easy to implement to define a single or bi-level threshold. However, when they are extended to multiple levels, they will cause some problems, such as long time-consuming, the high computational cost, and the needed improvement in their accuracy. To avoid these problems and determine the optimal multilevel image segmentation threshold, we proposed a hybrid Marine Predators Algorithm (MPA) with Salp Swarm Algorithm (SSA) to determine the optimal multilevel threshold image segmentation MPASSA. The obtained solutions of the proposed method are represented using the image histogram. Several standard evaluation measures, such as (the fitness function, time consumer, Peak Signal-to-Noise Ratio, Structural Similarity Index, etc.…) are employed to evaluate the proposed segmentation method’s effectiveness. Several benchmark images are used to validate the proposed algorithm’s performance (MPASSA). The results showed that the proposed MPASSA got better results than other well-known optimization algorithms published in the literature. Image segmentation is an important and challenging process [31] . It is involved in various fields, including pattern recognition, digital image processing, and artificial intelligence [32] . It is considered one of the most critical image processing steps that separate the image components into different parts by merging the same pixel groups. It is also used in the process of extracting object-related features from an image [37] . This is not an easy task; usually, the images contain some unwanted background noise. Various techniques are used in the segmentation process to subtract objects from images and separate or split them from backgrounds [21] . Image segmentation is also an essential step in computer vision. Also, it represents a necessary technique in image analysis, especially for medical image analysis. The images are divided into two types: color and gray ideas. Each of them has different methods of segmentation [46] . Image segmentation methods are classified into layer-based segmentation methods and block-based segmentation methods, as given in Fig. 1 . Among the most popular segmentation techniques used are Threshold, Histogram, Edge detection, Watershed Transformation, Clustering, and Region-based methods. As for separating the objects from the background, the Threshold method is used. The segmentation method, which uses threshold, is characterized by its simplicity and accuracy compared to other methods [23] . So, it has received the attention of many researchers in this field. The threshold method is widely used for dividing pixels in an image into different classes [49, 50] . There are two methods to determine image thresholds: the first is bi-level thresholding, and the second is multilevel thresholding. The first method divides the image into two categories based on the threshold value. It groups pixels densest from a threshold value into one class and pixels less dense than a threshold value into another level. They are often used to separate the foreground and background of the images. The process of dividing into two categories is insufficient, especially when the image is more complex and contains several objects with similar gray levels [45] . Therefore, it is necessary to extend the threshold levels from bi-level to multilevel [44] . The multilevel threshold (MTH) splits an image into several classes belonging to multiple objects Fig. 1 Methods of image segmentation in the image. In this method, it is possible to reveal more information and more items from the segmented image. Additionally, several histogram methods have been developed, such as Otsu's method and Kapoor's method. The basis for their work is to maximize the variance between classes and entropy to achieve homogeneity between categories. One of the problems caused by these methods is the complexity of the mathematical operations, primarily as the number of threshold levels increases to obtain better threshold values. Therefore, this problem's solution lies in combining these Metaheuristic algorithms. Several algorithms are used by researchers, such as particle swarm optimization, social-spider optimization [40] , and others [13] . In this case, conflict problems arise to achieve a specific objective; this target may not be suitable for image types. This is because each goal has a particular kind of image that makes it powerful. In the Otsu-based method, the best value of the threshold is maximizing the variance between classes. So, the choice of optimal thresholds in multilevel was a challenge over decades. Which affects segmentation accuracy. To address this weakness, this paper proposes an efficient multilevel threshold method depending on an enhanced Marine Predators Algorithm for image segmentation. The Salp Swarm Algorithm is used as a local search to improve the basic MPA's performance in this proposed method. The proposed algorithm, MPASSA, is a hybrid optimization algorithm for MTH that beats on the shortcomings of individual optimization algorithms using the force of both MPA and SSA. The candidate solutions for the proposed method are represented using an image histogram. The hybrid technique combines the characteristics of two different methods (MPA and SSA). Therefore, the proposed hybrid algorithm (MPASSA) averts getting stuck on a local optimum and has a high ability to find the optimal solution for the image. Several standard evaluation metrics, such as the function of fitness, signal-to-noise ratio, structural similarity index, consumer time, etc., are used to evaluate the proposed segmentation method's effectiveness. The MPASSA is assessed by using two experiments that have a set of images. In the first experiment set, two gray images are used. These images were widely applied in different studies to assess various segmentation algorithms. At the second experiment set, three-color images are used. Also, these images were involved in other studies to evaluate various segmentation algorithms. The results of MPASSA are compared with several optimization algorithms, such as WOA, PSO, AOA, SSA, MPA, and it shows significant performance improvement. The proposed method MPASSA has proved its effectiveness in all test cases. Both experiment sets' performance indicates that the MPASSA is an effective segmentation method, and it can be used in various segmentation applications. The results reveal that the proposed MPASSA is better than other known optimization algorithms published in the literature. The main contributions of this study are summarized as follows. 1) This research proposes a new image segmentation method using a multilevel threshold based on using optimization method. 2) The SSA operators are used to enhance the exploitation ability of the MPA, called MPASSA. 3) This research tests the proposed method's performance in two experiment sets using two popular grayscale images and three-color popular images. 4) The proposed method is compared with several state-of-the-art methods. The design of this paper is given as follows. Section 2 contains a literature review of the past and current image segmentation by multilevel thresholding algorithms. Section 3 presents the research methodology and techniques used in Image segmentation by multilevel thresholding algorithms. Section 4 illustrates the implemented proposed algorithms, the evolution measurements that have been used to test and evaluate the proposed algorithms, and research discussion and results. Finally, Section 5 presents the conclusion of this research, followed by suggestions and motivations for future work in the field of Hybrid Algorithms. One of the most critical image analysis problems and processing is the pixel segment into its various categories. The goal of segmentation is to identify elements and isolate them from the background and distinguish between pixels to improve contrast [41] . Several common methods can be used to solve this problem [5, 8, 12, 20, 26, 43, 53] . The most common segmentation techniques that researchers still use nowadays are Threshold, Histogram, Edge detection, Watershed Transformation, clustering, and Region-based methods [1] . There are two types of thresholds; threshold one value is called bi-level thresholds. If more than one threshold value is used, it is called a multilevel point, the second type. In the first type, pixels are divided into two categories. Pixels with an intensity level higher than the threshold value are classified as objects, while the remaining pixels are classified as a part of the background [38] . However, the second one used a multilevel threshold to divide the grayscale images into sections where it became necessary to extend the bi-level threshold levels to the multilevel. When the image is more complex and contains several objects with similar gray levels, it is also used with pictures with colored backgrounds. The multilevel threshold mechanism, which has more than one threshold, divides the image into several classes belonging to multiple objects. Therefore, it results in multiple items with one background. In this way, more information and things can be revealed from the segmented image [2] . Among the most critical methods used with threshold segmentation are the Otsu and Kapur methods [30] . The Otsu method increases the contrast between image classes. In contrast, the Kapur method maximizes entropy as a measure of homogeneity between categories. Both approaches have proven to be one of the most widely used forms of image processing. They are also considered accurate and effective alternatives to pixilation in the bi-level threshold technique. It can be used in a multilevel threshold, but its accuracy decreases as the thresholds' number increases, leading to increased complexity. Many other approaches adopt entropy, such as Tallies entropy, Fuzz,y entropy, the Shannon entropy, and Renyi's entropy. The same problem of increasing computational complexity appears at a multilevel threshold [38] . Therefore, to solve such issues, Meta-Heuristic algorithms provide highly accurate results in most cases. The Meta-Heuristic algorithms are essential tools for solving optimization problems' intricacies with high accuracy [4] . According to the official data provided by those countries, a new forecasting method is proposed in [19] to forecast the number of people infected with Covid-19 in some countries. The proposal offers an improvement to the ANFIS model to predict the number of injured people, depending on another optimizer, the Marine Predator algorithm. The algorithm is used to improve ANFIS parameters, as well as enhance prediction performance. The results of the prediction performance evaluation of the proposed method MPA-ANFIS were compared with several other ways. The result was that the proposed approach significantly outperforms almost all techniques and measures of performance. Among the performance measures that were compared are, mean absolute error (MAE), Root Mean Squared Relative Error (RMSRE), Mean Absolute Percentage Error (MAPE), and others [9, 14] . In [25] , a new variant of MPA is proposed. The proposed method's performance is tested on several problems, including the real-world engineering design and CEC-2017 tests. The algorithm was compared to other optimization algorithms: the first is GA and PSO, which is one of the algorithms that have been well studied, the second is GSA, CS, and SSA, as the algorithms developed in the recent period, third is CMA-ES, SHADE, and LSHADE -cnEpSin, are the algorithms that are as optimizers for performance. They are also IEEE CEC competition winners. The results were that the Marine Predator algorithm ranked second as the best performing method, and it showed competitive results compared to LSHADE-cnEpSin. The Marine Predator algorithm is one of the winning algorithms in the CEC 2017 competition. The statistical analysis presented in that paper shown that the Marine Predator algorithm can work as a high-performance optimizer. It also proved that the Marine Predator algorithm is significantly superior to the SSA, GSA, CMA-ES, PSO, CS, and GA algorithms. The study also showed the similarity of performance statistically for each Marine Predator, SHADE, and LSHADE-cnEpSin. In [34] , a novel hybrid algorithm is presented. The Proposed method ensembles Salp Swarm and Multi-objective Salp Swarm in solving the optimization problems for both single and multiple objectives. The Salp's primary behavior is a crowd of strains while on the move to search for food [6] . The researcher tested both algorithms on several mathematical optimization tasks to observe, monitor, and confirm their behavior in finding the best solution. The results showed that the Salp swarm algorithm could optimize the initial random solutions and are close to being the optimum. The results also showed that the multi-objective Salp Swarm algorithm could converge Pareto optimal solutions. The study also showed that both proposed algorithms could solve complex and computationally costly engineering design problems (such as designing marine propellers). The research also indicates that the proposed algorithms have many advantages that distinguish them in solving many real-world situations. In [3] , an effective multilevel threshold (MTH) method is proposed for image segmentation, including medical image segmentation including COVID-19 CT images. Suggested is the Marine Predatory Animal (MPA) algorithm. MPA is a new SI method. According to the researchers' study, the proposal provides the first application for use in image segmentation. The (MFO) algorithm was used in the MPA optimization process. The new proposal has been called MPAMFO. The request was evaluated with various images, as it included crosssectional images of Coronavirus . The results showed stable and good performance in all tests. Also, several comparisons were applied to prove the MPAMFO proposal's superiority over many algorithms, including PSO, GWO, and CS, in terms of SSIM, fitness value, and PSNR. The proposed algorithm has proven to be highly effective. Therefore, it is possible to improve it and apply it in various improvement processes, including data clustering, machine job scheduling, time series prediction, cloud computing, etc. A new proposal combining the SSA algorithm with the PSO algorithm is given in [28] . This proposal is to enhance SSA's ability for exploration and exploitation by using PSO characteristics to improve the quality of the SSA for searching for results. Thus, the rate of convergence increases. The proposed algorithm was evaluated with two experiments; first, it was tested on 15 standard groups. The second was applied to determine the optimal subset of the features from among ten UCI groups to increase classification accuracy. The results from the two experiments were compared to several algorithms, including SSA. The work SSA with the PSO algorithm gives better results in performance measures, including chosen feature ratio, processing unit time, accuracy, and evaluating the fitness function. This means that hybridization of the SSA gives better results than it does separately. In [47] , a novel objective function and a new application of the MPA approach are presented to properly extract the nine parameters of the PV module's TDPV model. This study aims to obtain an accurate PV model for any commercial PV panel, which plays an essential role in the grid-connected PV power systems' simulation studies. In this study, the optimization problem's main objective is to minimize a function representing the current error. The MPA technology was successfully utilized to reduce the objective function, obtaining the PV module's TDPV model's nine parameters. Various comparisons were exhibited to check the efficacy of the offered TDPV model using the MPA technology. The proposed algorithm was successfully employed to optimally design the parameters of two marketable KC200GT and MSX-60 PV modules. The optimal parameters, realized using the MPA-TDPV model, are coherent with those achieved using other algorithms. The MPA approach has recorded lower optimal fitness values of 1.245e-14 and 7.458e-13 for KC200GT and MSX-60 PV modules. Furthermore, the simulation outcomes of the MPA-TDPV model concur with the measured data for these well-known PV panels under several environmental situations. The MPA-TDPV model's ACE indicates a lower value than other PV models for the marketable PV panels. This points out the proposed approach's superiority, efficacy, and robustness for achieving a precise TDPV model-based PV panel. With the help of the MPA approach, accurate modeling of any marketable PV panel can be realized. Moreover, the MPA technology can be further extended to solve various optimization problems in power system applications, energy storage devices, and smart grids. In [52] , an innovative objective function with a robust and reliable optimization algorithm named the marine predators' algorithm (MPA) is proposed to provide the optimal pattern structure for three dimensions of PV arrays nine × 9, 16 × 16, and 25 × 25. The MPA is tested with several shade patterns and compared with manta ray foraging optimization (MRFO), Harris hawk optimizer (HHO), and particle swarm optimizer (PSO), as well as the total-cross, tied (TCT) connection. Several quality and statistical measures are computed, such as mismatch power loss, fill factor, percentage power loss, and Wilcoxon signed-rank test to assess the performance of the proposed approach. The research showed that the I-V and P-V characteristics were used to investigate the proposed MPA's applicability compared with the other counterparts. Moreover, the mean execution time has been evaluated. The results reveal that MPA enhanced the PV array power by the percentage of 28.6%, 2.7%, and 5.7% in cases of 9 × 9, 16 × 16, and 25 × 25 PV arrays, respectively, and a uni-peak PV characterizes is achieved as well with lowest execution time and highest consistency in the results across the number of independent runs. Therefore, the authors recommend MPA as an efficient and applicable algorithm for PV reconfiguration systems at any dimension of PV array structures. In [29] , a node localization algorithm has been proposed based on Salp Swarm Algorithm (SSA), which handled the node localization problem as an optimization problem. The proposed algorithm has been implemented and validated in various WSN deployments using other target nodes and anchor nodes. Moreover, the proposed algorithm has been evaluated compared to four well-known optimization algorithms, namely PSO, BOA, FA, and GWO, in terms of localization accuracy, computing time, and several localized nodes. +e obtained simulation results have proved the proposed algorithm's superiority compared to the other localization algorithms regarding the various performance metrics. The proposed approach can be hybridized with different algorithms to reduce the localization error in future work. In [27] , an improved version of the salp swarm algorithm is introduced for predicting chemical compound activities. A set of assessment indicators are used to evaluate and compared with different algorithms, including particle swarm optimization (PSO), Grasshopper Optimization Algorithm (GOA), Grey Wolf Optimizer (GWO), Sine Cosine Algorithm (SCA), Whale Optimization Algorithm (WOA) using three initialization method, and a superior accuracy was obtained with our proposed approach. Also, compared with other algorithms that used the same data, this research's system has a higher performance using fewer features. The previous algorithms (GOA, GWO, PSO, SSA, SCA, and WOA) are compared. Three different methods were used to initialize the various optimization algorithms to ensure the other optimizers' capability to converge from various initial positions, namely mixed initialization, short initialization, and extensive initialization. In [10] , a new hybrid Arithmetic Optimization Algorithm (AOA) and Deferential Evolution (DE) algorithm (DAOA) is proposed for multi-thresholding image segmentation. The proposed algorithm employs the AOA algorithm to optimize the threshold and then uses this thresholding value to partition the images through DE. So, the DAOA integrates AOA global optimization and DE fast convergence. The experiment results were compared against four algorithms. The DAOA achieved better results than other methods; also, the DAOA provided a faster convergence with relatively lower CPU time. In the future, the DAOA can be applied to other complex image segmentation problems. We concluded that the given methods in the literature are employed to solve the image segmentation problems by finding the optimal number of threshold values. However, finding these values is a complicated problem and needs an efficient method to determine them [48, 17] . Several optimization techniques proved their ability to solve this proposed and demonstrated that an efficient method is required. Thus, we proposed an effective method to address this problem and finding the optimal threshold values. Before presenting the proposed hybrid algorithm (MPASS), we need to define the problem and understand the mathematical model for the main algorithms: Marine Predators Algorithm and Salp Swarm Algorithm. In this section, we present the formulation of the problem, MPA and SSA, and their main items. In this section, the problem formulation of multilevel thresholding is illustrated, which gives a mathematical definition through considering a gray level image I. The image I is tested to be segmented consisting K + 1 classes. Each segmentation process requires division pixels of the image I into subregions or classes as in Eq. (1). This can be done by determining K thresholds {t 1 , t 2 , . ., t K } [22] . In given Equations, G K represents K-th class of image I, I (i, j) is the value of gray level to the pixel (i, j), t K where (k = 1,…, K) that define k-th threshold value, and L represents the gray levels of I, all levels within the range (0, 1…L -1). Multilevel thresholding aimed to locate the optimal threshold values that divide I into several groups; Optimization methods can determine these values. One of the standard methods and well-known optimization function is Otsu's method [39] . It is used in bi-level segmentation problems as well as in multilevel. In multilevel thresholding, Otsu's technique locates the optimum threshold values of I by maximizing as in Eq. (2) . It can be defined as the following: where η 1 pointed to mean in tensity of image I, with t 0 = 0 and t K+1 = L. h(i) is frequency and P j is probability of the ith gray level. N is the all pixels in the evaluated I. This part will explain the Marine Predators Optimization Algorithm (MPA) development process as an efficient and straightforward meta-heuristic optimization method [19] . The basis of the algorithm's work is population. The first solution is distributed randomly on the search area in a uniform manner, as in Eq. (6): where X min and X max are the minimum and maximum bound for variables, the rand is a random vector within range (0, 1). Based on the theory of survival of the fittest, an array can be constructed to include the best solution for the best predator. The best predator is the one that has the most talent for foraging. This matrix is called the elite. The array is represented mathematically as in Eq. (7). Elite ¼ X I In Eq. (8), the Xi,j represents the jth dimension of (ith ) prey. It is noted that the optimization process is entirely linked directly and mainly to these two matrices. The process of optimizing the algorithm and moving it from one phase to another depends heavily on the predator's speed relative to the prey. Therefore, the process of improving the Marine Predators algorithm can be divided into three primary stages, depending on the different speed ratios for both the prey and the predator and the life of each of them. The steps are as follows: 1-The movement of the prey relative to the predator is faster with a high rate. 2-The movement of both predator and prey at equal or close speed. 3-The movement of the predator to the prey is faster, with a low rate. Phase No. 1 The exploration stage is where the prey moves at high speed to discover new search areas. At this stage, the predator should remain and monitor the prey's movements. This stage's performance is in the first third of the total iterations of development (i.e., 1/3 * Max_iter) as in Eq. (9). Where CF is an adaptive variable for controlling the step volume of a predator. R B represents random numbers and is represented by a vector. These numbers are regular in [0,1] and represent Brownian motion. P = 0.5 represents a vector and constant number for random numbers. The variable Iter represents the current iteration, while Max_iter represents the maximum iteration. Phase No. 2 At this point, both the predator and the prey are moving at the same speed. This part occurs in the middle stage of improvement in Eq. (10) . Depending on the rule, the unit pace ratio (v ≈ 1), if prey moves in Lévy, the best-organized predator is Brownian motion. This can be considered a transitional stage between exploration and exploitation. Both exploitation and exploration are essential in this section. Where R L refers to the random numbers, where these numbers are following the Lévy movement in distribution, the first half represents exploitation (represented by predators), Eq. (11) is applied to it. While the prey that represents the second half of the population assumes: Max_iter is the overall number of generations. The R B and Elite simulate the predator's movement in a Brownian motion, while the prey changes its position depending on the predator movement in Brownian fashion. This phase is the last of the optimization process. They represented when the predator's movement is faster than a prey. The predator exploits the prey that is detected and attack it very quickly to get it. So, the exploitability is often high. The best motion for the predator in a low pace ratio (v = 0.1) is Lévy. This phase is executed on the last third of the iteration numbers (Iter > 2/3 Max_iter), the predator follows Lévy. The predator will update its position in Eq. (12) . where, Step size is added to the Elite site for the predator's movement to update the prey site. The doubling of the R L and Elite also simulates the predator's movement within Lévy's strategy. Another point that causes a behavioral change in marine predators is environmental issues such as the eddy formation or Fish Aggregating Devices (FADs) effects. As mentioned in [28] , Sharks spend over 80% of their time in the vicinity of FADs. The others are taking the long jumps in different dimensions to find other prey. The effects of FADs can be mathematically expressed in Eq. (14) . Where FADs = 0.2 represent the probability of affected FADs on the process of the optimization. U is a vector of binary with matrices that have contained zero and one. This is done by generating an arbitrary vector in [0,1] and update the array to zero if it's less than 0.2, and one if it's more. r is the uniform arbitrary number [0,1]. Both (X min , X max ) are the vectors consist of (lower, upper) bounds of the dimensions. r 1 and r 2 represent indexes of random prey. Figure 2 illustrates the second phase of optimization, in which the predator adopts the Brownian strategy better in searching for its prey within the field, as indicated in blue. When optimization is in its final phase, algorithms need a high ability for exploitation. The third phase represents the last phase of optimization when the predator changes its behavior from Brown's strategy to Levi so that the search process will be more efficient in a particular area. At the same time, the Convergence Factor (CF) has an excellent role for predators. This limits the search in several parts of the specific area for exploitation. It also avoids the effort in the search that extends from long steps as a result of using Levi's strategy for non-promising areas in the field. The flowchart of the MPA algorithm is represented as shown in Algorithm 1. The Salp algorithm is one of the meta-heuristic algorithms that has been successfully used in solving many optimization problems in various fields [34] . Mathematically, as we know, that the swarm is divided into a leader and followers. The leader guides his followers in their movements. The X represents a swarm of n of the Salps as Eq. (15), represented by a twodimensional matrix. The fitness of each Salp is calculated to determine the best among them (meaning, defines the leader for the swarm). The leader positions will be updated by using Eq. (16). . . . Wherexi1represents the position for the leader (first) Salp in the dimension ith, and y i is the food site in the ith dimension. lb i and ub i is the lower and upper bound at the ith dimension, respectively, and the coefficient r 1 is calculated by Eq. (17) . r 2 and r 3 are random numbers between [0,1]. Where L is the upper iterations and l is the running iteration. Note that the coefficient r1 is very significant in SSA because the balances of exploration and exploitation depend on it during the entire search process. As for the followers, Eq. (18) shows the updated positions. Where j >= 2, xi1pointed to the site of the nth Salp in the ith dimension, 0 represents an initial speed, t is the time, and 0 ¼ xÀx0 t . The time in optimization points to the iteration. So, the difference between iterations is equal to 1. Considering the proposition that o ¼ 0 , the following equation (in Eq. (19)) is employed for this issue. Where j >= 2. In some cases, the Salp leaving outside of the search space, how we can bring them back to the search space in Eq. (20) . Exploration and exploitation, diversification and intensification, global search and local searches, few pairs of words are common in optimization algorithms. At least the optimization algorithms include one of those pairs. Here we will talk about exploration and exploitation. Exploration aims to discover new areas in the search area. This makes the search process not stop within one level, leading to discovering the optimal solution level. The exploitation aims to reach a better solution within the suitable and discovered solutions. A coefficient responsible for the balance between exploration and exploitation, it is called the coefficient c 1 . It is calculated from the equation Eq. (2). The Pseudo-code of the SSA is given in Algorithm 2. The research introduces a hybrid algorithm used to segment images. The proposed hybrid algorithm determines the optimal multilevel threshold values that maximize Otsu's objective function. The proposed algorithm is called MPASSA and is based on the MPA and SSA algorithms and uses Otsu's method as an objective function. The proposed hybrid algorithm depends on boosting the performance of MPA via using the SSA. The MPA is used to reduce the search region by determining the best solution; then, SSA optimizes each agent less than the limited base. Therefore, the ABCSCA algorithm starts by computing the histogram of an input image and then generates a random population of pop_size solutions (threshold values). Then, the MPA updates these populations using its levels representing by three faces of motion, then Eddy Formation and FADs' Effect. The optimal solution is then determined from the population based on Otsu's method (the best solutions from the MPA algorithm). The SSA begins determining the minimum threshold value by using the MPA's output (worst solution) and optimizing the solutions of the population via the strategy discussed before. The optimal global solution (the SSA algorithm's output) is determined, and all the previous steps are repeated until the stopping conditions are met. The final stage of the MPASSA is illustrated in Fig. 3 . It can be observed that the proposed method introduces an arithmetic complex which is O (t (nd Cof * n)), where n represents a number of factors, t is the iteration, d is the problem dimension, and Cof is the cost of function. In this part, the performance of the MPASSA is evaluated by two experiments (gray and colored). It is compared with five algorithms: original Marine Predators Algorithm MPA, original Salp Swarm Algorithm (SSA), Whale Optimization Algorithm (WOA), Particle Swarm Optimization (PSO), and Arithmetic Optimization Algorithm (AOA). Besides, two types of images are used to test the proposed method (2 grays and 3 color images). The parameters settings of the employed methods are takes from the original methods. The supplied results is expressed as the best fitness function, PSNR, and SSIM. Furthermore, for all of the acquired findings, the Friedman ranking test is used to demonstrate how significantly the suggested approach outperforms the comparable methods for all of the examined issues. The methods are developed and implemented in MATLAB (2015a) on a machine with the following specifications: CPU 2.3 GHz (an Intel Core i7 platform), 16 GB RAM, 2400 MHz, and DDR4. The comparison of all algorithms used in the research has the same cessation conditions (maximum iterations set to 100), with a total run of 30 per algorithm and the population's size (25) . For performance evaluation of all experiments on the test images are done by the number of thresholds: 2, 3, 4, and 5 as in [18] . The selection for such thresholds was to illustrate the performance of the proposed algorithm (MPASSA) compared with the traditional segmentation algorithms based on the swarm. All algorithms are programmed in "Matlab2017" and implemented on a "Windows 10 -64bit" environment on a computer with Intel Corei7 8th Gen (1.99 GHz) processor" 1T" memory and SSD (256). In order to evaluate the fitness of the segmented image, four measures are used. The measures that applied to the proposed algorithms, as follows: 1. The execution time. 2. The fitness value is calculated to evaluate and assess each solution based on its current position, as mention in chapter four. 3. The Peak Signal-to-Noise Ratio (PSNR) measure [51] : used to measure the variance between the reference image and segmented image, and it depends on the value of intensity in the image, and this refers to the fitness of the reconstructed image. The PSNR is defined as: RMSE refer to the root-mean-squared error, detect as: (Seg) image, defined as: µI and µSeg refer to mean intensity of (I) and (Seg) respectively; while σI and σSeg refer to standard deviation of (I) and (Seg) respectively; σI, Seg represent a variance of (I) and (Seg). c1 and c2 are constants, that c 1 = 6.5025 and c 2 = 58.52252 [36] . The maximum value of SSIM refers to better performance. The Friedman test is the non-parametric alternative to the one-way ANOVA with repeated measures. It is used to test for differences between groups when the dependent variable being measured is ordinal. It can also be used for continuous data that has violated the assumptions necessary to run the one-way ANOVA with repeated measures [7] . Table 1 show the parameters settings for each algorithm that are used in the following experiments. In addition, the general parameters are set as follows. This paper tested all algorithms on the five images: two gray images and three-color images. These images are from the dataset of Berkeley University [33] , namely, lena, baboom for gray images as shown in Fig. 4(a, b) , and lena, baboom, and peppers for color images as shown in Fig. 5(a, b, c) . With size (512 × 512) for all. This section discusses the proposed hybrid MPASSA compared to other popular algorithms in the segmentation domain. The comparisons are conducted using several measures (i.e., fitness value, PSNR, SSIM, and time). Two grayscale of images and three-color images with four numbers of thresholds (2, 3, 4, and 5) are applied and compared the results with five comparative algorithms (i.e., MPA [25] , SSA [34] , WOA [1] , PSO [15] , and AOA [11] ) The results of the comparison are given in terms of PSNR, SSIM, CPU time, Fitness values, in addition to the Friedman ranking test for the comparative methods using the PSNR and SSIM Mean values for each image. These are standard performance evaluation measures in image segmentation, especially for evaluating multilevel thresholds methods. The discussion depends on each criterion are given as follows. In terms of PSNR measure, the testing results of the threshold two as in Tables 6 and 7 , the MPASSA got the best results in all images (gray and color) followed by MPA at some times which got the best results, followed by WOA, PSO, SSA, and AOA. The results of threshold 2 indicate that the MPASSA obtained the best values in all images. The MPASSA outperformed all other algorithms in thresholds 3, 4, and 5. It achieved the best results in all images, followed by the WOA, MPA, PSO, and SSA. From the results, we can derive that MPASSA is the best algorithm according to the PSNR measure, followed by MPA and WOA, and these algorithms outperformed the others when the threshold is greater than 4. In Tables 8 and 9 , the comparative algorithms' results are presented in terms of SSIM for gray and color images, respectively. According to the results obtained by applying two thresholds, the proposed MPASSA got the best results in all images (i.e., gray and color), followed by SSA, which got the best results in 3 images out of 5, While the values of each AOA, MPA, WOA and PSO algorithms are comparative according to the tested images. At the threshold level 3, the proposed MPASSA method obtained the best results in all images, followed by MPA, WOA, SSA, and AOA. Also, at the threshold levels 4 and 5, the proposed method obtained the best results in all images. Simultaneously, each WOA, SSA, AOA, PSO, and MPA are similar according to the SSIM values. From these results, we can conclude that MPASSA is the best algorithm among all comparative algorithms. The CPU time is considered for all algorithms and listed in Tables 10 and 11 to show all algorithms' performance in terms of exhaustion time. To compare the time complexity of the proposed MPASSA method against other algorithms, Table 10 presents the CPU time for the gray level, and Table 11 presents the CPU time for the color images. These tables illustrated that the proposed method consuming average time between the other algorithms. It was ranked fourth, better than MPA and PSO algorithms, and worse than SSA, AOA, and WOA. While SSA was the lowers consuming of time execution followed by AOA and then WOA at all levels with all images (gray and color). The Friedman ranking test results using PSNR and SSIM for all the tested methods using different tested images (gray and color) are given in Tables 12, 13, 14, 15, 16, 17, 18, 19, 20 For the gray images, in Table 12 of PSNR value, the MPASSA got the first ranking followed by WOA, it got the second-ranking, and SSA has the third ranking. Simultaneously, the AOA and PSO ranked fourth, and MPA got the last ranking. These results proved the ability of the proposed MPASSA in solving the image segmentation problem. In Table 13 of the PSNR value, the MPASSA and SSA got the first ranking followed by AOA; it got the second-ranking, MPA has the third ranking. In contrast, WOA got the fourth ranking, and PSO has the fifth the last ranking in the table. This result illustrated that the proposed MPASSA has a promising outcome compared to other methods. Tables 14, 15 and 16 shows the ranking results based on using PSNR values for color images. In Table 14 , the MPASSA got the first ranking, followed by each of WOA and SSA, which have got the second-ranking, MPA has the third-ranking. In contrast, PSO has the fourth-ranking, and AOA got the last ranking in the table. In Table 15 , the MPASSA got the first ranking followed by PSO; it got the second-ranking, SSA has the third-ranking, while Table 21 The results of the Friedman ranking test for the comparative methods using the SSIM values for image peppers SSIM K MPASSA WOA SSA AOA PSO MPA peppers 2 1 5 2 6 4 3 3 1 6 2 3 5 4 4 1 2 5 6 4 3 5 2 6 1 3 4 5 SUM 5 19 10 18 17 Fig. 9 Segmented lena gray images using Otsu's (a2-f2) represent threshold2, (a3-f3) represent threshold3, (a4-f4) represent threshold4, (a5-f5) represent threshold5 WOA got the fourth-ranking, followed by AOA has the fifth ranking, and MPA at the last ranking in the table. In Table 16 , the MPASSA got the first ranking, followed by SSA. It got the second-ranking, MPA has the third-ranking, while AOA got the fourth-ranking, followed by PSO have the fifth ranking and WOA at the last ranking in the table. The given results revealed that the modified version of MPA using SSA got better results than other comparative methods and avoided the basic version of MPA's weaknesses. For the gray images in Table 17 of SSIM value, the MPASSA got the first ranking followed by PSO; it got the second-ranking, each of WOA and SSA has the third-ranking. In contrast, MPA has the fourth-ranking, and AOA got the last ranking. In Table 18 of SSIM value, the MPASSA got the first ranking followed by SSA; it got the second-ranking, MPA has the thirdranking. In contrast, AOA got the fourth-ranking, and WOA has the fifth ranking and PAO at the last ranking. This result illustrated that the proposed MPASSA has a better performance compared to other methods. In Table 19 , the MPASSA got the first ranking for the color images, followed by SSA, which got the second-ranking with WOA. At the same time, MPA has the third-ranking, and PSO got the fourth-ranking, and at the end ranking, AOA has the last ranking. In Table 20 , the MPASSA got the first ranking followed by PSO; it got the second-ranking, SSA has the third- Fig. 10 Segmented baboom gray images using Otsu's (a2-f2) represent threshold2, (a3-f3) represent threshold3, (a4-f4) represent threshold4, (a5-f5) represent threshold5 ranking, while WOA got the fourth-ranking, followed by AOA has the fifth ranking and MPA at the last ranking. In Table 21 , the MPASSA got the first ranking, followed by SSA; it got the second-ranking, MPA has the third-ranking, while PSO got the fourth-ranking, followed by AOA has the fifth ranking, and WOA at the last ranking. The given results showed that the proposed version of MPA using SSA got promising results compared to other comparative methods and avoided the basic version of MPA. Figures 9, 10, 11 , 12, and 13 show the segmented test images (gray and color) obtained by the comparative methods using various threshold values (k= 2, 3, 4, and 5). The proposed MPASSA got the best-segmented images compared to other comparative methods. Figure 14 shows the gray and color images histogram distribution estimated by the proposed hybrid algorithm MPASSA at k equal to 2, 3, 4, and 5. The proposed MPASSA got the optimum results in most cases. This because each of the images is considered as a varied optimization problem. Moreover, because of the randomness of the swarm methods, the results can change in some cases. Also, depending on No-Free Lunch Theorem NFL that illustrates, it is hard to find one algorithm suitable for many optimization problems. As well as, each image has a varied gray or colored level histogram. Also, image segmentation becomes a complex operation due to Fig. 11 Segmented lena color images using Otsu's (a2-f2) represent threshold2, (a3-f3) represent threshold3, (a4-f4) represent threshold4, (a5-f5) represent threshold5 the histogram's multimodality. The segmented image is constructed by determining each class a grey or the colored level value specified from the class members' mean. So, for each image, the proposed MPASSA performed better in terms of the PSNR, SSIM, and fitness values at each level. The MPASSA may fail to segment the image in some cases and at different levels. For these reasons, we used statistical comparisons to determine the best algorithm for several images. Multilevel thresholding is one of the methods used in image segmentation. It is considered a preprocessing phase in many applications. In this research, the problem of defining the optimal/ best threshold for image segmentation in the case of multilevel thresholding. This problem has been considered as an 'optimization problem'. Otsu's function has been applied as a fitness function. So, a novel hybrid algorithm in image segmentation of multilevel threshold was proposed. It draws from the properties of both Marine Predators and Salp Swarm algorithms. The proposed algorithm has been used to solve this problem. Each one aims to define the optimal threshold value, which increases Otsu's function. The proposed Fig. 12 Segmented baboom color images using Otsu's (a2-f2) represent threshold2, (a3-f3) represent thresh-old3, (a4-f4) represent threshold4, (a5-f5) represent threshold5 novel algorithm results have been compared with MPA, SSA, WOA, AOA, and PSA algorithms. The performance of algorithms has been evaluated based on the following measures: fitness values, Peak Signal-to-Noise Ratio (PSNR), Structural Similarity Index (SSIM), CPU-time, and Friedman ranking test. We used five benchmark images for gray and color images in all experiments. The study concluded that the proposed hybrid MPASSA method's performance is better than that of the MPA and SSA algorithms and other algorithms. The leading cause for this supremacy is that hybrid technology combines the characteristics of different methods. Therefore, the proposed hybrid algorithm averts getting stuck on a local optimum. All of this indicates that hybrid techniques have a high ability to find the optimal solution for the image segmentation problem. The proposed method can solve various problems and applications as of image processing, such as visualization, computer vision, computer-aided diagnostics, image classification, etc. In all of these applications, the finesses of images are a critical problem in image segmentation. Low resolutions, high noises, or lower contrast, can cause this problem. So, the preprocessing of the images is an essential point to improve the quality of that status and multi-objective cases to gain optimal segmentation results. Other recent optimizers, such as Arithmetic Optimization Algorithm (AOA), can solve the given problem with new modifications in the future. Fig. 13 Segmented peppers color images using Otsu's (a2-f2) represent threshold2, (a3-f3) represent thresh-old3, (a4-f4) represent threshold4, (a5-f5) represent threshold5 Authors contribution Laith Abualigah: Conceptualization, methodology, writing, review, editing, and supervision. Nada Khalil Al-Okbi: Experimental, draft writing and validation. Mohamed Abd Elaziz: Experimental, Writing, review, editing, and supervision. Essam H. Houssein: Writing, review, and editing.Data availability -All the used data in this Human and animal rights This article does not contain any studies with human participants performed by any of the authors. No animal studies were carried out by the authors for this article Gy and Color images histogram at K=2, 3, 4, 5 obtained by MPASSA References Whale optimization algorithm and moth-flame optimization for multilevel thresholding image segmentation Multi-level thresholding-based grey scale image segmentation using multiobjective multi-verse optimizer An improved Marine Predators algorithm with fuzzy entropy for multi-level thresholding: Real world example of COVID-19 CT image segmentation Hyper-heuristic method for multilevel thresholding image segmentation Advanced optimization technique for scheduling IoT tasks in cloud-fog computing environments Salp swarm algorithm: a comprehensive survey Nature-inspired optimization algorithms for text document clustering-a comprehensive analysis Aquila Optimizer: A novel meta-heuristic optimization Algorithm Group search optimizer: a nature-inspired meta-heuristic optimization algorithm with its results, variants, and applications A novel evolutionary arithmetic optimization algorithm for multilevel thresholding segmentation of COVID-19 CT images The arithmetic optimization algorithm Amended hybrid multi-verse optimizer with genetic algorithm for solving task scheduling problem in cloud computing A comprehensive survey of the Grasshopper optimization algorithm: results, variants, and applications Advances in sine cosine algorithm: a comprehensive survey Unsupervised text feature selection technique based on hybrid particle swarm optimization algorithm with genetic operators for the text clustering A new feature selection method to improve the document clustering using particle swarm optimization algorithm A novel instance segmentation algorithm based on improved deep learning algorithm for multi-object images A study on particle swarm optimization and artificial bee colony algorithms for multilevel thresholding Marine predators algorithm for forecasting confirmed cases of COVID-19 in Italy, USA, Iran and Korea Solving capacitated vehicle routing problem using cooperative firefly algorithm Medical image segmentation algorithm based on multilayer boundary perceptionself attention deep learning model Cuckoo search algorithm and wind driven optimization based study of satellite image segmentation for multilevel thresholding using Kapur's entropy Monarch butterfly optimization algorithm for computed tomography image segmentation Marine predators algorithm for optimal allocation of active and reactive power resources in distribution networks Marine Predators Algorithm: A nature-inspired metaheuristic Development and application of slime mould algorithm for optimal economic emission dispatch Swarming behaviour of salps algorithm for predicting chemical compound activities Improved salp swarm algorithm based on particle swarm optimization for feature selection Abualigah L (2022) Dwarf Mongoose Optimization Algorithm A new method for gray-level picture thresholding using the entropy of the histogram. Computer vision, graphics, and image processing Various image segmentation techniques: a review Image segmentation via multi dimensional color transform and consensus based region merging A database of human segmented natural images and its application to evaluating segmentation algorithms and measuring ecological statistics Salp Swarm Algorithm: A bio-inspired optimizer for engineering design problems The whale optimization algorithm A hybrid differential evolution for optimal multilevel image thresholding A general framework for complex network-based image segmentation Multilevel thresholding by fuzzy type II sets using evolutionary algorithms A threshold selection method from gray-level histograms Social spiders optimization and flower pollination algorithm for multilevel image thresholding: a performance study A survey: image segmentation techniques. In: Harmony Search and Nature Inspired Optimization Algorithms A new arithmetic optimization algorithm for solving real-world multiobjective CEC-2021 constrained optimization problems: diversity analysis and validations A novel deep learning-based feature selection model for improving the static analysis of vulnerability detection Multi-level image thresholding using Otsu and chaotic bat algorithm A generalized Masi entropy based efficient multilevel thresholding method for color image segmentation A study of image segmentation algorithms for different types of images Marine predators algorithm for parameters identification of triple-diode photovoltaic models A novel deep learning pipeline architecture based on CNN to detect Covid-19 in chest X-ray images A joint convolutional neural networks and context transfer for street scenes labeling Weakly supervised adversarial domain adaptation for semantic segmentation in urban scenes Multilevel minimum cross entropy threshold selection based on particle swarm optimization A robust strategy based on marine predators algorithm for large scale photovoltaic array reconfiguration to mitigate the partial shading effect on the performance of PV system Hybrid Aquila optimizer with arithmetic optimization algorithm for global optimization tasks