Next Article in Journal
Generating Red-Edge Images at 3 M Spatial Resolution by Fusing Sentinel-2 and Planet Satellite Products
Next Article in Special Issue
Hybrid Fusion-Based Background Segmentation in Multispectral Polarimetric Imagery
Previous Article in Journal
Gravity Field Model Determination Based on GOCE Satellite Point-Wise Accelerations Estimated from Onboard Carrier Phase Observations
Previous Article in Special Issue
Speckle Noise Reduction Technique for SAR Images Using Statistical Characteristics of Speckle Noise and Discrete Wavelet Transform
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Harris Hawks Optimization with Mutation Mechanism for Satellite Image Segmentation

1
College of Mechanical and Electrical Engineering, Northeast Forestry University, Harbin 150040, China
2
Depto. de Ciencias Computacionales, Universidad de Guadalajara, CUCEI, Av. Revolución 1500, Guadalajara C.P. 44100, Jalisco, Mexico
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(12), 1421; https://0-doi-org.brum.beds.ac.uk/10.3390/rs11121421
Submission received: 2 May 2019 / Revised: 10 June 2019 / Accepted: 12 June 2019 / Published: 14 June 2019
(This article belongs to the Special Issue Image Optimization in Remote Sensing)

Abstract

:
In this paper, a novel satellite image segmentation technique based on dynamic Harris hawks optimization with a mutation mechanism (DHHO/M) is proposed. Compared with the original Harris hawks optimization (HHO), the dynamic control parameter strategy and mutation operator used in DHHO/M can avoid falling into the local optimum and efficiently enhance the search capability. To evaluate the performance of the proposed method, a series of experiments are carried out on various satellite images. Eight advanced thresholding approaches are selected for comparison. Three criteria are adopted to determine the segmentation thresholds, namely Kapur’s entropy, Tsallis entropy, and Otsu between-class variance. Furthermore, four oil pollution images are used to further assess the practicality and feasibility of the proposed method on real engineering problem. The experimental results illustrate that the DHHO/M based thresholding technique is superior to others in the following three aspects: fitness function evaluation, image segmentation effect, and statistical tests.

Graphical Abstract

1. Introduction

Image segmentation is a fundamental and crucial stage in some applications, such as computer vision, pattern recognition, image classification, etc. [1,2,3,4]. Generally speaking, the purpose of the segmentation operation is to partition the grayscale or color image into several non-overlapping regions with unique characteristics. The gray levels of the pixels in the same area are roughly the same, while the pixels from different areas are significantly different. In the last few decades, many segmentation methods have been proposed by scholars, such as clustering, fractal-wavelet modeling, region growing, and thresholding [5,6,7]. Among all available techniques, threshold-based (thresholding) method is extensively used due its simplicity and efficiency [8,9].
Thresholding technique can be mainly divided into two categories: bi-level and multilevel [10]. Bi-level thresholding splits the image into two classes according to the intensity. However, if the image contains much information and multiple objects, the bi-level thresholding is not competent. Therefore, as an extension of bi-level thresholding, multilevel thresholding is proposed to handle the current issues [11]. For example, Díaz–Cortés et al. segment the breast thermograms using multilevel thresholding, which provides a highly reliable clinical decision support and advances the progress of thresholding approach in medical imaging domain [12]. Also, multilevel thresholding is adopted by Bhandari et al. to process the satellite images, which improves the segmentation accuracy and obtains high-quality segmented images [13]. In fact, researchers have introduced numerous approaches based on different criteria to select the segmentation thresholds over the past few decades, such as Kapur’s entropy [14], Tsallis entropy [15], and Otsu between-class variance [16]. Kapur’s entropy measures the homogeneity of segmented regions by maximizing the histogram entropy. Tsallis entropy also determines the thresholds by maximizing the entropy, but there are certain constraints. Otsu’s technique maximizes the between-class variance of the segmented image.
Multilevel thresholding can be considered as a NP-hard problem with high complexity [10]. To be more specific, the computational complexity increases exponentially with the number of thresholds, which greatly affects the efficiency and feasibility. To overcome this drawback, meta-heuristic algorithms are introduced into this domain and combined with the thresholding technique. In 2017, He and Huang proposed a modified firefly algorithm (MFA) based method to select the segmentation thresholds [11]. Kapur’s entropy, Otsu between-class variance, and minimum cross entropy are served as fitness functions. Experimental results reveal that the proposed technique can efficiently reduce the computational complexity and produce segmented images with more details. Khairuzzaman and Chaudhury combined the grey wolf optimizer (GWO) with thresholding for grayscale image segmentation [7]. Kapur’s entropy and Otsu’s method are used to select the segmentation thresholds. It can be found from the results that the proposed approach can accurately determine the segmentation thresholds and significantly reduce the computational complexity. In 2017, Oliva et al. presented an efficient magnetic resonance (MR) image segmentation method based on crow search algorithm (CSA) and minimum cross entropy [17]. The results of the comparison experiment show that the well-delimited regions obtained are easier to distinguish than other techniques [17]. Furthermore, some other algorithms or their modified versions have also been introduced into this domain, such as whale optimization algorithm (WOA) [18], multi-verse optimizer (MVO) [19], grasshopper optimization algorithm (GOA) [20], social spiders optimization (SSO) [21], krill herd optimization (KHO) [22], and cuckoo search (CS) [23] as well as Lévy flight firefly algorithm (LFA) [24], hybrid differential evolution (hjDE) [25], adaptive wind driven optimization (AWDO) [26], etc. These promising results motivate us to apply some other efficient meta-heuristic algorithms to multilevel image thresholding. Besides, as stated in No-Free Lunch Theorem (NFL) [27], there is no algorithm can solve all optimization problems. Thus, the recently proposed and uninvestigated algorithm also has potential.
Harris hawks optimization (HHO) is a high-performance, population-based, gradient-free optimization technique, which was proposed by Heidari et al. in 2019 [28]. Inspired from the cooperative behavior and chasing style of Harris hawks in nature, the mathematical model of HHO was established [28]. Similar to other meta-heuristic algorithms, the exploration and exploitation stages are also included in HHO. Each Harris hawk represents a candidate solution in the search space and the best solution obtained so far is considered as the prey or approximate optimal solution. During the whole process of predation, the Harris hawks patiently wait for the prey, besiege the prey to make it exhausted, and finally attack the prey with a surprise pounce. According to [28], HHO has presented better or occasionally competitive results than well-established techniques on 29 benchmark problems and several real-world engineering tasks. These phenomena illustrate the potential and the superiority of HHO, which also motivates us to apply this powerful algorithm to multilevel image thresholding.
No algorithm is perfect. Each algorithm requires some adjustments to fit the current problem, as does the HHO algorithm. For example, the exploration and exploitation phases of HHO are unbalanced. To be more specific, the absolute value of the parameter E cannot be larger than 1 during the latter half of the iteration (e.g., the last 250 of 500 maximum iterations), that is, the search agent does not perform a global search, although the current region may not be promising. For some complex optimization tasks, there is no guarantee that the population has gathered near the global optimum at the end of the exploration phase, resulting in premature convergence and local optimization. Therefore, two improvement strategies are introduced in this paper to enhance the optimization ability of the original HHO. As for the former, a novel dynamic control parameter strategy is presented to avoid trapping into the local optimum. The essence of the improvement strategy is to add a disturbance term to the update formulation of the escaping energy ( E ) . Sine and cosine functions are used in combination to control where the disturbance peak appears. A Gaussian distribution is also adopted to increase randomness. In the experiment, not only the influence of parameter change on the performance of disturbance term, but also the reason for using the Gaussian distribution instead of the Levy distribution is analyzed. As for the latter, a famous mutation operator, known as DE/best/2 is utilized to improve the global search efficiency and population diversity [29]. Then, the modified version of HHO, known as dynamic Harris hawks optimization with a mutation mechanism (DHHO/M) is applied to multilevel image thresholding.
Satellite imagery contains a wealth of information and is crucial to environmental resource monitoring and evaluation, such as forest cover, wetland resources, climate change, etc. Although remote sensing technology has many advantages (e.g., fast update cycle, fewer interference factors, and saving manpower and material resources), it is still a challenge to extract boundaries, locate objects and separate regions in high-resolution satellite images [30]. To solve the above problems, DHHO/M is combined with thresholding technique for satellite image segmentation. A series of experiment are carried out to evaluate the performance of the proposed method. The experiments can be mainly divided into four parts. For the first part, the impact of both the improvement strategies on performance is investigated and analyzed. For the second and the most important part, the DHHO/M based approach is compared with four advanced thresholding techniques on satellite images. Various performance metrics have been taken into consideration, such as objective function value, standard deviation (Std), peak signal to noise ratio (PSNR), mean squared error (MSE), structural similarity index (SSIM), feature similarity index (FSIM), and convergence property as well as the Wilcoxon’s rank sum test and Friedman test. For the third part, DHHO/M is compared with other thresholding techniques based on different criteria. For the fourth part, four oil pollution images are used to further assess its practicality and feasibility on real engineering problems.
The contributions of this paper have three aspects:
  • Introduce the recently proposed HHO into multilevel image thresholding. To the best of the authors’ knowledge, this attempt has not been made yet.
  • Dynamic control parameter strategy and mutation mechanism are used to improve the search efficiency of the original HHO.
  • Objectively and comprehensively evaluate the performance of the proposed technique.
The structure of this paper is described as follows: in Section 2, the problem statement and definitions of Kapur’s entropy, Tsallis entropy, and Otsu’s method are given. Section 3 briefly reviews the original HHO. In Section 4, the proposed DHHO/M based image thresholding technique is introduced in details. Section 3 presents the experimental results and relevant discussions. Finally, the conclusion and future research direction are represented in the last section.

2. Material and Methods

2.1. Problem Statement

For multilevel image thresholding, the ultimate goal is to select the combination of segmentation thresholds. Assume that the given image is segmented into n + 1 classes { C 0 , C 1 , , C n } by n thresholds { T 1 , T 2 , , T n } , which can be defined as follows [31]:
C 0 = { I ( i , j ) | 0 I ( i , j ) T 1 1 } C 1 = { I ( i , j ) | T 1 I ( i , j ) T 2 1 } C k = { I ( i , j ) | T k I ( i , j ) L 1 }
where I ( i , j ) denotes the gray level of the ( i , j ) th pixel, and L is the number of gray levels in the given image.
Furthermore, the combination of the optimal thresholds is determined by optimizing (maximizing or minimizing) the objective function, which is mathematically expressed as the following equations:
{ T 1 , T 2 , , T n } = arg max T 1 , , T n f ( T 1 , T 2 , , T n ) ,
{ T 1 , T 2 , , T n } = arg min T 1 , , T n f ( T 1 , T 2 , , T n ) ,
where f represents the objective function and { T 1 , T 2 , , T n } is the combination of the optimal thresholds.

2.2. Multilevel Thresholding

2.2.1. Kapur’s Entropy

Kapur’s entropy is a famous thresholding method that selects the thresholds based on the entropy of the segmented classes [14]. The objective function can be defined as:
f K a p u r = i = 0 n H i , H i = j = T i T i + 1 1 p j ω i ln ( p j ω i ) ,
where ω i = j = T i T i + 1 1 p j , p j = h ( j ) / N . ω i indicates the sum of the probabilities of the pixels in the i th class, p j is the proportion of pixels in each gray level to the total, h ( j ) represents the frequency, and N is the number of all pixels. Besides, for multilevel image thresholding, T 0 = 0 and T n + 1 = L .

2.2.2. Tsallis Entropy

Tsallis entropy is extensively used by researchers in the field of image thresholding to choose the optimal combination of thresholds [15]. The objective function is mathematically represented as:
f T s a l l i s = i = 0 n S q c i + ( 1 q ) i = 0 n S q c i , S q c i = 1 j = T i T i + 1 1 ( p j ω i ) q q 1 ,
where q is a constant equal to 4, the value selected is the same as [30,32].

2.2.3. Otsu Between-Class Variance

Otsu’s method was proposed in 1979, which maximizes the between-class variance of segmented classes to determine the thresholds [16]. The objective function is presented as follows:
f O t s u = i = 0 n σ i 2 , σ i 2 = j = T i T i + 1 1 ω i ( μ i μ ¯ ) ,
where μ i = j = T i T i + 1 1 j p j ω i , μ ¯ = j = 0 L 1 j p j . μ ¯ is the mean intensity for whole image.
The mathematical form of above three optimization problems can be expressed as:
Consider x = [ x 1 , x 2 , , x n ] .
Maximize f K a p u r ( x ) = i = 0 n j = x i x i + 1 1 p j ω i ln ( p j ω i )
or
f T s a l l i s ( x ) = i = 0 n 1 j = x i x i + 1 1 ( p j ω i ) q q 1 + ( 1 q ) i = 0 n 1 j = x i x i + 1 1 ( p j ω i ) q q 1
or
f O t s u ( x ) = i = 0 n j = x i x i + 1 1 ω i ( μ i μ ¯ ) .
Subject to x i   must   be   integer , i = 1 , 2 , , n Variable range 0 x i 255 .
As discussed above, the computational complexity of the thresholding techniques will increase exponentially with the number of thresholds. Thus, meta-heuristic algorithms are used to optimize the objective functions in Equations (4)–(6) to improve the efficiency.

2.3. Harris Hawks Optimization

HHO is a novel population-based, gradient-free optimization technique, which was proposed by Heidari et al. in 2019 [28]. HHO simulates the predation, surprise pounce, and attacking behaviors of Harris hawks in nature. Similar to other meta-heuristic algorithms, HHO also includes two optimization stages, namely exploration and exploitation (see Figure 1), which are described in the following subsections.

2.3.1. Exploration Stage

Harris hawks have insightful eyes that help them track and detect prey, however, sometimes the prey is not easy to find. Then, the Harris hawks will perch, wait patiently, and may last for hours. In HHO, above behaviors are modeled as the exploration stage, which can be expressed as:
X i t + 1 = { X r a n d t r 1 | X r a n d t 2 r 2 X ( t ) | q 0.5 ( X rabbit X m t ) r 3 ( L B + r 4 ( U B L B ) ) q < 0.5 ,
where X i t + 1 is the position of i th individual in ( t + 1 ) th iteration, X rabbit denotes the position of the rabbit (prey), and q is a random number in the interval [0,1] that converts the two strategies. r 1 , r 2 , r 3 , and r 4 are also random numbers inside [0,1]. L B and U B are the upper and lower bounds of the given optimization problem respectively. X m t represents the average position of the population that can be calculated as follows:
X m t = 1 N i = 1 N X i t ,
where N is the size of the population, and X i t is the position of i th individual in t th iteration.

2.3.2. Transition from Exploration to Exploitation

The transition between exploration and exploitation is critical to the performance of meta-heuristic algorithms. In HHO, the escaping energy of the rabbit, known as E , is used to convert these two phases. The value of E decreases as the number of iterations increases, which can be mathematically modelled as:
E = 2 E 0 ( 1 t t max ) ,
where E 0 is a random number that changes in the interval [−1,1], t denotes the current iteration, and t max represents the maximum number of iterations.
More specifically, the value of E decreases from [−2,2] to 0 as the iterative process. If | E | 1 , the exploration stage is used to search for the prey; If | E | < 1 , the exploitation stage is adopted to exploit the promising area.

2.3.3. Exploitation Stage

After detecting the prey, the Harris hawks tend to attack it using surprise pounce. The actual predation process is often very complicated, the prey has a chance to escape and the Harris hawks also react differently according to the behavior of the prey. In order to better simulate the real situation, four strategies are used in the exploitation stage. A random number ( r ) is used to determine whether the prey has successfully escaped. The case r < 0.5 indicates successful escape, while r 0.5 shows unsuccessful escape. The escaping energy of the prey ( E ) affects the behavior of the Harris hawks. If | E | 0.5 , the soft besiege occurs; if | E | < 0.5 , the hard besiege happens [28].
A. Soft besiege
Considering r 0.5 and | E | 0.5 , the rabbit has enough energy and keeps trying to escape, even though these acts are in vain and will only make the prey itself exhausted. At the same time, the Harris hawks softly encircle the rabbit and then attack it using surprise pounce. The mathematical model is given as follows:
X i t + 1 = Δ X i t E | J X rabbit X i t |
Δ X i t = X rabbit X i t ,
where Δ X i t indicates the difference between the position vector of the rabbit and the current individual. r 5 is a random number inside [0,1], and J = 2 ( 1 r 5 ) indicates the random jump strength of the rabbit during the escape process.
B. Hard besiege
Considering r 0.5 and | E | < 0.5 , the rabbit becomes so exhausted and does not have enough escaping energy. Meanwhile, the Harris hawks hardly encircle the rabbit to conduct the final pounce. This process can be modeled as:
X i t + 1 = X rabbit E | Δ X i t | .
C. Soft besiege with progressive rapid dives
Considering r < 0.5 and | E | 0.5 , the rabbit has enough energy and managed to escape. In the meantime, a more intelligent soft besiege is performed before the surprise pounce. Note that the position update of the Harris hawks in this situation is a two-step process. If the first step is not an improved move, then the second step will be executed. For the former part, it can be mathematically modelled as:
Y = X rabbit E | J X rabbit X i t | .
For the latter part, the Levy flight is used to simulate the rapid, abrupt, and irregular movement of Harris hawks when chasing the rabbit. The formulation for position update is presented as:
Z = Y + S × Levy ( d i m ) .
where d i m is the dimensions of the optimization problem, S is a random vector of size 1 × d i m , and Levy denotes the Levy distribution that is defined as:
Levy ( x ) = 0.01 × u × σ | v | 1 β , σ = ( Γ ( 1 + β ) × sin ( π β 2 ) Γ ( 1 + β 2 ) × β × 2 ( β 1 2 ) ) 1 β ,
where u and v are random values in the interval [0,1]. β is a constant equal to 1.5.
Then, for a problem to be minimized, the mathematical model of the whole process at this stage is:
X i t + 1 = { Y i f   f ( Y ) < f ( X i t ) Z i f   f ( Z ) < f ( X i t ) ,
where f is the fitness function for the given optimization problem.
D. Hard besiege with progressive rapid dives
Considering r < 0.5 and | E | < 0.5 , the rabbit does not have enough energy but it successfully escapes. In this situation, the Harris hawks perform a hard besiege before the surprise pounce, and they try to bring the whole group closer to the prey, not just the individual. Therefore, the formulation for position update of the Harris hawks is given as follows:
X i t + 1 = { Y i f   f ( Y ) < f ( X i t ) Z i f   f ( Z ) < f ( X i t ) ,
where
Y = X rabbit E | J X rabbit X m t |
Z = Y + S × Levy ( d i m )
Furthermore, the pseudocode of HHO is represented in Algorithm 1.
Algorithm 1 Pseudocode of Harris hawks optimization (HHO) algorithm
Input: The size of population N , maximum number of iterations t max .
Output: The position of the rabbit and the corresponding fitness function value.
  • Initialize the position of the hawks X i and the rabbit X rabbit .
  • Initialize the fitness values of the hawks f i and the rabbit f rabbit .
  • Set the dimensions of the optimization problem   d i m .
  • While (termination condition is not met ( t < t max ) ) do
  •  Check the boundary and evaluate the fitness value of each hawk f i .
  •  Update the location X rabbit and fitness value f rabbit of rabbit if there is a better one.
  • For (each hawk ( i = 1 : N )) do
  •   Update the energy of the rabbit E using Equation (9).
  •   If ( | E | 1 ) then      % Exploration Stage
  •    Update the position using Equation (7).
  •   Else If ( | E | < 1 ) then    % Exploitation Stage
  •    If ( r 0.5   and   | E | 0.5 . ) then        % Soft besiege
  •     Update the position using Equation (10).
  •    Else If ( r 0.5   and   | E | < 0.5 ) then       % Hard besiege
  •     Update the position using Equation (12).
  •    Else If ( r < 0.5   and   | E | 0.5 ) then   % Soft besiege with progressive rapid dives
  •     Update the position using Equation (16).
  •    Else If ( r < 0.5   and   | E | < 0.5 ) then   % Hard besiege with progressive rapid dives
  •     Update the position using Equation (17).
  •    End If
  •   End If
  • End For
  • End While

2.4. Proposed Dynamic Harris Hawks Optimization with Mutation Mechanism

In this section, the proposed DHHO/M algorithm based image segmentation technique is introduced in detail. For multilevel color image segmentation problem, the input to the method was the histogram of each color component, and the output was the segmentation thresholds. This process was implemented by optimizing the objective functions in Equations (4)–(6) using the DHHO/M algorithm. The following subsections will discuss the two improvement strategies respectively, and then give the algorithm steps. Finally, the time complexity of DHHO/M and HHO is analyzed.

2.4.1. Dynamic Control Parameter Strategy

In the original HHO algorithm, the transition of the exploration and exploitation phase is based on the escaping energy ( E ) of the rabbit. More specifically, the absolute value of the parameter E decreased from 2 to 0 throughout the whole iteration. If | E | > 1 , HHO explored the search space to determine the promising area. On the contrary, if | E | < 1 , the exploitation strategy was adopted to enhance the local search efficiency. However, the main drawback of this transition process is that the absolute value of the parameter E cannot be larger than 1 during the latter half of the iteration (e.g., the last 250 of 500 maximum iterations), that is, the search agents do not perform global search during the second half of the iteration, although the current region may not be promising. For some complex optimization tasks, there was no guarantee that the population has gathered near the global optimum at the end of the exploration phase, resulting in premature convergence and local optimal values.
Therefore, a novel dynamic control parameter strategy is introduced with the purpose of jumping out of the local optimum under the premise of ensuring accuracy. The essence of the improvement strategy is to add a disturbance term to the update formulation of E in Equation (9). The disturbance term and modified update formulation of E are presented as follows:
δ = r a n d n ( sin α ( π 2 t t max ) + cos ( π 2 t t max ) 1 )
E = 2 E 0 ( 1 t t max ) + δ ,
where sin and cos denote the sine and cosine functions respectively. r a n d n is a random number subject to Gaussian distribution. t and t max represent the current iteration and the maximum number of iterations respectively. α is a constant that determines where the disturbance peak appears.
Figure 2 illustrates the variation of the disturbance term δ under different values of parameter α . It can be found that the amplitude decreases as the value of the parameter α increases, and the high amplitude disturbance occurs earlier when α = 2 .
In order to more intuitively illustrate the impact of dynamic control parameter strategy on the escaping energy of the rabbit, the variation of the parameter E with the number of iterations is also given in Figure 3. It is worth mentioning that the values of E was not only the key to global and local phases transition, but also affects the convergence of the algorithm because it was part of the position update formulation. To be more specific, a larger value indicates more exploration, while a smaller value indicates more exploitation. As can be observed from the figure, when α = 2 , the absolute value of E in the late iteration can be abruptly greater than 1, but the amplitude of the disturbance term is too large, which seriously affects the convergence performance. In the case of α = 3 , the amplitude of the perturbation term is too small, which cannot enhance the ability of jumping out of the local optimum, but will not affect the performance of the algorithm itself. Considering α = 2.5 , it can be found that the convergence property was not greatly affected, and the global search strategy also had an opportunity to be adopted in the later iteration. As far as the authors are concerned, for optimizing some fitness functions that contain local best, setting α to 3 may reduce the accuracy; setting α to 2 made the results obtained similar to the original algorithm; while setting α to 2.5 can improve the performance of the algorithm. These hypotheses will be verified in the experimental section.
In addition, it was necessary to explain the reason why the Gaussian distribution is used instead of the Levy distribution in the disturbance term. Although the characteristic of usually small steps and occasional large jumps in Levy distribution fits well with the idea of the improvement strategy, sometimes there are some extreme cases as shown in Figure 4, which can seriously affect the performance of the algorithm. More specifically, the disturbance of uncertain position makes the population evolve slowly and even degenerate.
Reducing the step size of Levy flight can theoretically improve this situation, but in fact, it is difficult to determine an appropriate value to meet the desired requirements. In our study, if the step size is larger than 0.05, disturbances with large amplitude will always appear; if the step size is smaller than 0.01, the disturbance is not enough to affect the value of E ; while the value between the two is still not ideal. Therefore, the Gaussian distribution is more suitable for the current problem, and the amplitude generated by it is acceptable compared with the Levy distribution.
Furthermore, according to our experimental statistics, the disturbance peaks of our proposed strategy usually appearred during the 300th to 400th iterations (the maximum number of iterations was 500), especially the 350th to the 400th iterations.

2.4.2. Mutation Mechanism

In order to enhance the global search efficiency of HHO, mutation mechanism was introduced into the exploration stage. The famous DE/best/2 mutation operator was adopted to replace the original location update strategy in Equation (7) [29]. The formulation is given as follows:
X i t + 1 = X best t + F ( X r 1 t X r 2 t ) + F ( X r 3 t X r 4 t ) ,
where x i t + 1 indicates the position of the i th individual at the t th iteration, and F is the scaling factor. r 1 , r 2 , r 3 , and r 4 are different integers selected from the range [ 1 , N ] and distinct from i . N denotes the size of population. Thus, Equation (7) can be rewritten as follows:
X i t + 1 = { X best t + F ( X r 1 t X r 2 t ) + F ( X r 3 t X r 4 t ) q 0.5 ( X rabbit X m t ) r a n d 1 × ( L B + r a n d 2 × ( U B L B ) ) q < 0.5 ,
where r a n d 1 and r a n d 2 are random numbers in the interval [0,1].
The operators of differential evolution (DE) have been widely used as part of modified or hybrid algorithm to improve the optimization capability and search efficiency. This is also the motivation for us to adopt the mutation operator mentioned above. For example, Xiong et al. proposed a hybrid algorithm of whale optimization algorithm (WOA) and DE, known as DE/WOA, for extracting parameters of solar photovoltaic models [33]. DE/best/2 is utilized in this proposed algorithm with the purpose of improving the ability of exploring the search space and locating the region of global optimum [33]. Also, Jadon et al. used DE to modify the onlooker bee phase of artificial bee colony (ABC) algorithm [34]. DE/best/1 is adopted to accelerate the convergence speed and maintain population diversity. Furthermore, DE has also been embedded or merged into some other improved versions of meta-heuristic algorithm, such as moth search algorithm based DE (MSDE) [35], improved whale optimization algorithm (IWOA) [36], self-adaptive gravitational search algorithm and DE (SGSADE) [37], improved dragonfly algorithm (IDA) [38], gbest-guided ABC (GABC) [39], hybrid whale optimization algorithm with DE (WOA-DE) [18], chaotic opposition-based grey wolf optimization algorithm based on DE and disruption operator (COGWO2D) [40], etc.
In fact, there are other types of mutation and crossover operators that have been proposed by researchers to enhance the optimization ability. In 2019, Xu et al. introduced a series of new variants of moth-flame optimization (MFO) by combining MFO with Gaussian mutation (GM), Cauchy mutation (CM), Lévy mutation (LM) or the combination of GM, CM and LM [41]. The experimental results demonstrated that the three strategies can significantly boost exploration and exploitation capabilities of the basic MFO. For the crossover operators, Kita et al. analyzed the unimodal normal distribution crossover (UNDX) by discussing the importance of the distribution and statistics of the offspring yielded by a crossover operator for its evaluation [42]. Akimoto et al. proposed an adaptive real-coded ensemble crossover (AREX) that combines the adaptation of expansion rate technique and the crossover mean descent technique in 2009 [43]. Experimental results showed that the proposed algorithm enlarges the classes of functions that original methods can solve. In 2015, Ariyarit and Kanazaki introduced a modified version of genetic algorithm (GA) with multi-modal distribution crossover (MMDX) [44]. Blend crossover (BLX) [45] and UNDX were used as comparison algorithms. It can be found from the experimental results that MMDX is more competitive and can maintain higher population diversity.
In addition, experiments are performed on eight satellite images (see Figure 5) to determine the optimal parameter settings of the proposed algorithm. The sum of ranks with different combination of α and F is given in Table 1. The lower the ranking value, the better the performance. As can be observed, α = 2.5 and F = 0.5 present the best results in most cases, thus this combination of parameter values is adopted in the following experiments.

2.4.3. Algorithm Steps

In this section, a comprehensive algorithm step of DHHO/M based multilevel color image thresholding technique is given in Algorithm 2.
Algorithm 2 Pseudocode of dynamic Harris hawks optimization with a mutation mechanism (DHHO/M) based multilevel color image thresholding
Input: The given color image.
Output: Segmentation thresholds.

Get   information   about   the   image
24. Read the given color image.
25. Extract the histogram of each color component (R, G, and B).

Harris   Hawks   Optimization
26. Initialize the position of the hawks X i and the rabbit X rabbit .
27. Initialize the fitness values of the hawks f i and the rabbit f rabbit .
28. Set population size N and maximum number of iterations t max .
29. Set the dimensions of the optimization problem d i m , namely the number of thresholds.
30. While (termination condition is not met ( t < t max ) ) do
31.  Check the boundary and evaluate the fitness value of each hawk f i using Equations (4)–(6).
32.  Update the location X rabbit and fitness value f rabbit of rabbit if there is a better one.
33.  For (each hawk ( i = 1 : N )) do
34.   Update the energy of the rabbit E using Equation (21). % Dynamic control parameter strategy
35.   If ( | E | 1 ) then
36.    Update the position using Equation (23).     % Mutation mechanism
37.   Else If ( | E | < 1 ) then
38.    If ( r 0.5   and   | E | 0.5 ) then
39.     Update the position using Equation (10).
40.    Else If ( r 0.5   and   | E | < 0.5 ) then
41.     Update the position using Equation (12).
42.    Else If ( r < 0.5   and   | E | 0.5 ) then
43.     Update the position using Equation (16).
44.    Else If ( r < 0.5   and   | E | < 0.5 ) then
45.     Update the position using Equation (17).
46.    End If
47.   End If
48.  End For
49. End While
Fitness function (Kapur’s entropy)
Input: Histogram of a color component, and segmentation thresholds X i .
Output: Fitness function value f i .
  • The histogram is divided into n + 1 parts by n thresholds.
  • Calculate the proportion of pixels in each gray level ( p j , j [ 0 , 255 ] ) to the total based on the histogram.
  • Compute the sum of the probabilities of the pixels ( ω k , k [ 0 , n ] ) contained in each part.
  • Calculate the Kapur’s entropy of each part ( H k , k [ 0 , n ] ) .
  • The sum of the entropies of all parts represents the fitness function value.
  • f i = H 0 + H 1 + + H n

2.4.4. Computational Complexity

The computational complexity of DHHO/M and HHO algorithms mainly depends on three processes: initialization, function evaluation, and position update. The improvement strategy of this paper only replaces the relevant formulations, and does not add additional judgments, loops or other commands. Therefore, the computational complexity of DHHO/M is approximately the same as that of HHO, which can be computed as follows:
O ( D H H O / M ) = O ( i n i t i a l i z a t i o n ) + O ( f u n c t i o n   e v a l u a t i o n ) + O ( p o s i t i o n   u p d a t e ) = O ( N ) + O ( t max × N × C o f ) + O ( t max × N × d i m ) = O ( N × ( 1 + t max × ( C o f + d i m ) ) ) ,
where N is the size of population, C o f is the cost of function evaluation, d i m is the dimension of the given optimization problem, and t max is the maximum number of the iterations [46].
In fact, the proposed DHHO/M runs slightly faster than the original HHO, which will be demonstrated and analyzed in subsequent experiments.

3. Discussions

3.1. Experimental Setup and Database

In order to assess the performance of the proposed approach, a series of experiments were carried out on various images. More specifically, the experiments can be divided into four parts. For the first part of the experiment, the performance influence of each strategy (dynamic control parameter and mutation mechanism) is investigated, and the relevant algorithms are denoted by DHHO and HHO/M, respectively. Two DHHO variants with different parameter values were also included to further evaluate the effect of the parameter ( α ) . The algorithm with parameter α equal to three is represented as DHHO+, and DHHO indicates α = 2 . For the second part of the experiment, the proposed method is compared with the traditional HHO and four advanced multilevel thresholding techniques, namely teaching-learning-based optimization (TLBO) [47], whale optimization algorithm with TH heuristic (WOA-TH) [48], improved differential search algorithm (IDSA) [49], and beta differential algorithm (BDE) [50]. The experiment is conducted on eight satellite images (see Figure 5) with low and high threshold levels [51]. For the third part of the experiment, DHHO/M is compared with other thresholding techniques based on different criteria to illustrate its feasibility and universality, such as Tsallis entropy based modified grasshopper optimization algorithm (MGOA) [20] and modified artificial bee colony (MABC) [13], as well as Otsu method based modified flower pollination algorithm (MFPA) [48] and grey wolf optimizer (GWO) [7]. For the fourth part of the experiment, four oil pollution images are used to further evaluate its practicality and feasibility on real engineering problems.
It is worth noting that the parameters of each thresholding method used are the same as the original literature, except for the population size N set to 30 and the number of iterations t max set to 500 for fair comparison. The parameter setting can be found in Table 2. Besides, the initial population of all algorithms is randomly generated in the solution space, and the termination condition in this work is to reach the maximum number of iterations t max . A number of 50 independent runs are conducted, and the best result is highlighted in boldface. All experiments are executed on a computer with “Microsoft Windows 10” system and “8GB” memory space under Matlab2017.

3.2. Performance Metrics.

Several performance metrics are briefed in this section, which can be observed in Table 3. Among the available measures, average fitness function value and standard deviation (Std) are used to evaluate the optimization capability of meta-heuristic algorithm, while peak signal to noise ratio (PSNR), mean squared error (MSE), structural similarity index (SSIM), and feature similarity index (FSIM) are utilized to assess the quality of the segmented image. Furthermore, Wilcoxon’s rank sum test and Friedman test were adopted for statistical analysis.

3.3. Experimental Series 1: Performance on Mathematical Functions

In this section, the proposed algorithm is evaluated on three mathematical test functions with different characteristics, such as unimodal, multi-modal, and composite. The test functions are taken from the CEC2005 special session, which are known as F1, F9, and F14. It can be found from Table 4 (bold is the best) that the proposed method outperforms others in most cases, showing its remarkable performance on mathematical test functions. This promising result motivates us to apply it to multilevel image thresholding domain.

3.4. Experimental Series 2: Influence of Dynamic Control Parameter Strategy and Mutation Mechanism

In this section, the influence of each improvement strategy (dynamic control parameter and mutation mechanism) is investigated, and the results can be found in Table 5 (bold is the best). As can be observed, DHHO/M presents the highest fitness function value in most cases, while HHO/M and DHHO are ranked second and third respectively. This promising result indicates that introduction of each strategy can improve the performance of the original algorithm, and the combination of them works better. Specifically speaking, the dynamic control parameter strategy prevents falling into local optimum, and the mutation mechanism effectively enhances the global search capability. Thus, DHHO/M exhibits superior performance compared to other algorithms.
Another comparison is made between DHHO, HHO, DHHO+, and DHHO. From the table, these algorithms can be arranged as DHHO > HHO ~ DHHO+ > DHHO, which strongly supports the discussion in the previous section. Smaller value of parameter ( α = 2 ) makes the amplitude of the disturbance term too large, although the diversity of the population is increased but the algorithm convergence ability is seriously affected. On the contrary, larger value of parameter ( α = 3 ) produces a small disturbance. Although the local optimum cannot be avoided, it does not affect the accuracy of the algorithm, so the result is similar to the original algorithm. However, the value of parameter α of DHHO is between the above two, which can enhance the ability to jump out of local optimum while ensuring convergence property, and obtain the best result in this comparison.
Furthermore, the Std values obtained by each algorithm are also given in Table 5. It can be found that DHHO/M presents the lowest value in general, which shows better stability than other algorithms. More detailed and comprehensive stability testing and analysis can be found in the next subsection.

3.5. Experimental Series 3: Comparison with Other Advanced Methods on Satellite Images

In this section, the proposed method is compared with other advanced multilevel image segmentation approaches on satellite images. The performance of each algorithm is comprehensively and objectively assessed in terms of precision, stability, statistical significance, scalability (high dimension), and convergence property.

3.5.1. Segmentation Accuracy

Accuracy is one of the main indicators for measuring image segmentation technique. In this experiment, the precision of each method is evaluated from five aspects, namely average fitness function value, MSE, PSNR, SSIM, and FSIM values. Firstly, the threshold values and segmented images obtained are given in Table 6 and Figure 6, respectively. The threshold levels selected (K = 2, 3, 4, and 5) are the same as [47,49]. As can be observed from Table 6, all methods give similar results when the threshold level is small; while the threshold values determined become different when K = 4 and 5. The reason for this phenomenon is that the difficulty of the current optimization problem increases with the number of thresholds, and some meta-heuristic algorithms with poor search ability are not competent, such as “Image3” at three threshold levels, “Image5” at three threshold levels, and “Image6” at five threshold levels. Considering the segmented images presented in Figure 6, it can be found that the images with high levels contain more information and details than that with low levels, because the of entropy of a given image reflects its average information content [18].
Secondly, the average fitness function values obtained by all algorithms are given in Table 7 (bold is the best), which reflects the optimization capability of each meta-heuristic algorithm. Higher Kapur’s entropy value indicates stronger search ability and higher precision. It can be observed that the values of all algorithms increase with the number of threshold levels, and DHHO/M outperforms other compared methods in general. This results strongly proves that high-quality segmented image with more information can be obtained at higher threshold level, such as K = 5. In DHHO/M, the introduction of dynamic control parameter strategy and mutation mechanism can improve the search efficiency of the algorithm and avoid premature convergence. Therefore, the superiority of the proposed algorithm is holistic rather than on a single image.
Thirdly, several performance metrics are introduced to assess the similarity between the segmented image and the original image, such as PSNR, MSE, SSIM, and FSIM. The experimental results are given in two tables. In Table 8, it can be found that the PSNR value increases with the number of threshold values, while MSE value is the opposite (bold is the best). Similarly, SSIM and FSIM values also increase with the number of threshold levels, as can be seen from Table 9 (bold is the best). These promising phenomena indicate that the quality of segmented image is improving gradually; while the results obtained by the proposed method are better than or at least comparable to other approaches, showing better feature and structure preserving ability.

3.5.2. Statistical Test

The experiments performed by each algorithm are the same, so the statistical test is necessary. Parametric statistical test is based on various assumptions, such as independence, normality, and synchronization [56]. However, these assumptions probably not valid when analyzing meta-heuristic algorithms with stochastic property. Therefore, two famous nonparametric statistical tests are adopted in this study, namely, Wilcoxon’s rank sum test [54] and Friedman test [55]. As a brief review, the null hypothesis ( H 0 ) assumes that there is no difference between the two compared algorithms (Wilcoxon’s rank sum test) or the whole (Friedman test); while the alternative hypothesis H 1 indicates the difference. More detailed discussion of these two statistical tests can be found in [55].
The experimental results are given in three tables. Table 10 presents the p-values and h-values of Wilcoxon’s rank sum test obtained by each comparison. It can be found that DHHO/M gives better results in 29 out of 32 cases (eight images and four threshold levels) for HHO, 30 cases for TLBO, 31 cases for WOA-TH, 30 cases for IDSA, and 30 cases for BDE, which shows significant difference between the proposed method and other approaches. As for the average rank of Friedman test given in Table 11, the proposed method ranks first in all cases, and the value of rank becomes smaller as the number of threshold level increases, indicating greater advantages over other compared methods (bold is the best). Table 12 gives the Chi-square (χ2) value and p-value of Friedman test at different threshold levels. According to the Chi-square distribution table, the critical value for 5 (6 algorithm - 1) degrees of freedom at 0.05 significant level is 11.07 [57]. It can be observed from the table that the Chi-square values obtained are much larger than the critical value and the p-values obtained are much smaller than the significant level. The results illustrate that the advantage of the proposed algorithm becomes more obvious as the dimension of the optimization problem increases, and there is a significant difference between the available methods.

3.5.3. Computational Time

Computational time was also one of the important indicators to measure the efficiency of the algorithm. The experimental results at different threshold levels can be found in Table 13 (bold is the best). Obviously, the IDSA algorithm is the fastest, although it cannot give the best results in most cases. DHHO/M is a little faster than the traditional HHO, because the greedy selection used in exploitation stage consumes a lot of time; while the dynamic control parameter strategy enables the exploration phase can be adopted in the later iteration, which reduces the complexity to some extent. It is worth noting that evaluating a meta-heuristic algorithm should be comprehensive and objective, not just from a single aspect. The proposed method presents better or at least competitive results in terms of average fitness value, PSNR, MSE, SSIM, and FSIM, despite not the fastest. Thus, DHHO/M can be considered as an efficient technique for multilevel color image segmentation.

3.5.4. Search Capability on High Dimensional Problems

In this section, all available thresholding approaches are compared at K=10, 15, and 20 to assess the performance on high dimensional tasks. Threshold levels selected are the same as [57]. The average fitness function values obtained are presented in Table 14 (bold is the best). It can be observed that DHHO/M gives the best results in 22 out of 24 cases (eight images and three thresholds); TLBO and WOA-TH only give the best results in one case respectively, while other algorithms give none. As discussed in above section, the proposed method has shown certain advantages when K = 5, but it is more obvious in higher dimensions. The reason for this phenomenon is that each image can be considered as a different optimization problem, and there is no algorithm can handle all tasks [27]. Even the meta-heuristic algorithm with strong search capability cannot show remarkable performance on all images. Therefore, DHHO/M is more competitive in the field of image segmentation than other algorithms, because it gives the best results in most cases.

3.5.5. Stability Analysis

The stability of each algorithm is analyzed and discussed in this section, including two performance measures, namely Std value and boxplot. The Std values obtained on satellite images can be found in Table 7, respectively. A lower Std value indicates better performance. As can be observed, the proposed approach presents the lowest Std values in most cases, showing remarkable stability, continuity and consistency. Furthermore, the boxplot with 20 threshold levels are given in Figure 7, because the difference between the algorithms is more obvious in high dimensions. It can be found that DHHO/M outperforms other approaches again. More specifically, the boxplot obtained is higher in position and more compact, indicating better data consistency. Compared with other algorithms, the proposed method has no bad points (in our study), which illustrates the superiority of DHHO/M based thresholding technique in satellite image segmentation.

3.5.6. Convergence Property

In this section, the convergence property of DHHO/M algorithm is evaluated. The convergence curves obtained by all algorithms at 20 threshold levels are presented in Figure 8. It is worth mentioning that these figures are drawn in a semi-logarithmic coordinate system for ease of observation. As can be seen, WOA-TH and IDSA are prone to premature convergence and local optimization, which is reflected in the shape of the curve. For example, under the circumstance of “Image3” segmentation, the fitness function values obtained by these two algorithms no longer change after 100 iterations, and obviously they do not get the best value at the end of iteration. Although the fast convergence speed is not a bad characteristic of meta-heuristic algorithm, local optimization seriously affects the accuracy as well as the quality of segmented image. On the contrary, the propose method can better balance the exploration and exploitation stages, neither premature convergence nor slow convergence. Therefore, the DHHO/M based technique can exhibit higher search efficiency in all selected satellite images.

3.6. Experimental Series 4: Performance Using Different Objective Functions

In this section, Tsallis entropy and Otsu between-class variance are served as objective function to assess the feasibility of DHHO/M. Eight different color images are selected from Figure 5. The experiment are conducted on high threshold levels and can be divided into two parts. For the former part, Tsallis entropy based thresholding techniques are selected for comparison, namely the MGOA and MABC. For the latter part, Otsu between-class variance based thresholding approaches are chosen for testing, namely the MFPA and GWO.
The experimental results of Tsallis entropy and Otsu’s method are presented in Table 15 and Table 16, respectively (bold is the best). As can be observed, DHHO/M gives the best objective function values in all cases, indicating the strong optimization capability. Considering other metrics, the proposed algorithm also produces better or at least competitive results, which shows that the segmented images with more information and details can be obtained. To sum up, the superior performance of the proposed algorithm is not limited to the objective function adopted in this paper, but has potential in other image fields.

3.7. Experimental Series 5: Oil Pollution Image Segmentation

In order to further verify the effectiveness of the proposed algorithm in solving practical engineering problems, four oil pollution images are selected for the experiment. These images were taken by the drone in the area of the eighth oil production plant, which can be found in Figure 9 [58]. As can be observed, the oil pollution area of (b) is relatively obvious; while the remaining three images all have strong interference, especially (d). This phenomenon undoubtedly increases the difficulty of segmentation operation, which can be considered as a challenging engineering problem.
The segmented results of the proposed method are presented in Figure 10. Because there is no absolute standard for a real engineering problem, the authors manually labeled the target region and separated it. Then took it as the ground truth for experimental comparison. It can be found from the figures that the oil pollution in (a)–(c) has been successfully separated from the complex background, which is similar to the ground truth. Considering figure (d), the ideal oil pollution area consists of two parts, while the DHHO/M based thresholding technique incorrectly classifies them as a whole. But in fact, it has achieved the desired goal, because the main area of oil pollution has been efficiently identified. To sum up, although the proposed method cannot achieve perfect results for complex images with mixed background and strong interference, it is competent for most cases and can still be considered as a competitive technique for the segmentation of oil pollution image.

4. Conclusions

An efficient satellite image segmentation technique based on DHHO/M is proposed in this paper. Dynamic control parameter strategy and mutation mechanism are used to avoid trapping into the local optimum and improve the search efficiency. In order to validate the superiority of the proposed method, a series of experiments are conducted on various images. For the first part of the experiment, the experimental results indicate that both the improvement strategies adopted can enhance the optimization capability. For the second part of the experiment (the most crucial part), it can be observed that the DHHO/M based technique gives better results for satellite images in terms of objective function value, Std, PSNR, SSIM, FSIM, and convergence property as well as the Wilcoxon’s rank sum test and Friedman test. For the third part of the experiment, the robustness of the proposed approach is assessed by the segmentation based on other criteria. For the last part of the experiment, the proposed method is applied to a real engineering problem, namely the segmentation of oil pollution image to further evaluate its practicality and feasibility. The experimental results strongly illustrate the remarkable performance of the DHHO/M based thresholding technique.
In the future, the authors will introduce more efficient technique to handle the images with mixed background and strong interference. Furthermore, due to the conflicts between different criteria, multi-objective optimization for image thresholding is also the research direction.

Author Contributions

H.J. and C.L. contributed to the idea of this paper; C.L. performed the experiments; C.L. wrote the paper; H.J., C.L., D.O., and X.P. contributed to the revision of this paper. W.S. provided fund support.

Funding

This research was supported by the Fundamental Research Funds for the Central Universities (No. 2572019BF04), the National Nature Science Foundation of China(No. 31470714), the Northeast Forestry University Horizontal Project (No. 43217002, No. 43217005, No. 43219002), and the Research Training Project for College Students of Northeast Forestry University.

Acknowledgments

The authors would like to thank the anonymous reviewers for their constructive comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Frongillo, M.; Gennarelli, G.; Riccio, G. Plane Wave Diffraction by Arbitrary-Angled Lossless Wedges: High Frequency and Time Domain Solutions. IEEE Trans. Antennas Propag. 2018, 66, 6646–6653. [Google Scholar] [CrossRef]
  2. Hinojosa, S.; Avalos, O.; Oliva, D.; Cuevas, E.; Pajares, G.; Zaldivar, D.; Gálvez, J. Unassisted thresholding based on multi-objective evolutionary algorithms. Knowl. Based Syst. 2018, 159, 221–232. [Google Scholar] [CrossRef]
  3. Guido, R.C.; Addison, P.; Walker, J. Introducing wavelets and time-frequency analysis. IEEE Eng. Biol. Med. Mag. 2009, 28, 13. [Google Scholar] [CrossRef]
  4. Guido, R.C. Practical and useful tips on discrete wavelet transforms. IEEE Signal Process. Mag. 2015, 32, 162–166. [Google Scholar] [CrossRef]
  5. Mallat, S.G. A Theory for Multiresolution Signal Decomposition: The Wavelet Representation. IEEE Trans. Pattern Anal. Mach. Intell. 1989, 2, 674–693. [Google Scholar] [CrossRef]
  6. Guariglia, E. Entropy and Fractal Antennas. Entropy 2016, 18, 84. [Google Scholar] [CrossRef]
  7. Khairuzzaman, A.K.M.; Chaudhury, S. Multilevel thresholding using grey wolf optimizer for image segmentation. Expert Syst. Appl. 2017, 86, 64–76. [Google Scholar] [CrossRef]
  8. Hinojosa, S.; Dhal, K.G.; Elaziz, M.A.; Oliva, D.; Cuevas, E. Entropy-based imagery segmentation for breast histology using the Stochastic Fractal Search. Neurocomputing 2018, 321, 201–215. [Google Scholar] [CrossRef]
  9. Oliva, D.; Cuevas, E.; Pajares, G.; Zaldivar, D.; Osuna, V. A Multilevel Thresholding algorithm using electromagnetism optimization. Neurocomputing 2014, 139, 357–381. [Google Scholar] [CrossRef] [Green Version]
  10. Aziz, M.A.E.; Ewees, A.A.; Hassanien, A.E. Whale Optimization Algorithm and Moth-Flame Optimization for multilevel thresholding image segmentation. Expert Syst. Appl. 2017, 83, 242–256. [Google Scholar] [CrossRef]
  11. He, L.; Huang, S. Modified firefly algorithm based multilevel thresholding for color image segmentation. Neurocomputing 2017, 240, 152–174. [Google Scholar] [CrossRef]
  12. Díaz-Cortés, M.; Ortega-Sánchez, N.; Hinojosa, S.; Oliva, D.; Cuevas, E.; Rojas, R.; Demin, A. A multi-level thresholding method for breast thermograms analysis using Dragonfly algorithm. Infrared Phys. Technol. 2018, 93, 346–361. [Google Scholar] [CrossRef]
  13. Bhandari, A.K.; Kumar, A.; Singh, G.K. Modified artificial bee colony based computationally efficient multilevel thresholding for satellite image segmentation using Kapur’s, Otsu and Tsallis functions. Expert Syst. Appl. 2015, 42, 1573–1601. [Google Scholar] [CrossRef]
  14. Kapura, J.N.; Sahoob, P.K.; Wongc, A.K.C. A new method for gray-level picture thresholding using the entropy of the histogram. Comput. Vis. Graph. Image Proc. 1985, 29, 273–285. [Google Scholar] [CrossRef]
  15. Tsallis, C. Possible generalization of Boltzmann–Gibbs statistics. J. Stat. Phys. 1988, 52, 479–487. [Google Scholar] [CrossRef]
  16. Otsu, N. A threshold selection method from gray-level histograms. IEEE Trans. Syst. Man Cybern. 1979, 9, 62–66. [Google Scholar] [CrossRef]
  17. Oliva, D.; Hinojosa, S.; Cuevas, E.; Pajares, G.; Avalos, O.; Gálvez, J. Cross entropy based thresholding for magnetic resonance brain images using Crow Search Algorithm. Expert Syst. Appl. 2017, 79, 164–180. [Google Scholar] [CrossRef]
  18. Lang, C.; Jia, H. Kapur’s Entropy for Color Image Segmentation Based on a Hybrid Whale Optimization Algorithm. Entropy 2019, 21, 318. [Google Scholar] [CrossRef]
  19. Jia, H.; Peng, X.; Song, W.; Lang, C.; Xing, Z.; Sun, K. Multiverse Optimization Algorithm Based on Lévy Flight Improvement for Multithreshold Color Image Segmentation. IEEE Access 2019. [Google Scholar] [CrossRef]
  20. Liang, H.; Jia, H.; Xing, Z.; Ma, J.; Peng, X. Modified Grasshopper Algorithm-Based Multilevel Thresholding for Color Image Segmentation. IEEE Access 2019, 7, 11258–11295. [Google Scholar] [CrossRef]
  21. Ouadfel, S.; Taleb-Ahmed, A. Social spiders optimization and flower pollination algorithm for multilevel image thresholding: A performance study. Expert Syst. Appl. 2016, 55, 566–584. [Google Scholar] [CrossRef]
  22. Beevi K., S.; Nair, M.S.; Bindu, G.R. Automatic segmentation of cell nuclei using Krill Herd optimization based multi-thresholding and Localized Active Contour Model. Biocybern. Biomed. Eng. 2016, 36, 584–596. [Google Scholar] [CrossRef]
  23. Pare, S.; Kumar, A.; Bajaj, V.; Singh, G.K. An efficient method for multilevel color image thresholding using cuckoo search algorithm based on minimum cross entropy. Appl. Soft Comput. 2017, 61, 570–592. [Google Scholar] [CrossRef]
  24. Pare, S.; Bhandari, A.K.; Kumar, A.; Singh, G.K. A new technique for multilevel color image thresholding based on modified fuzzy entropy and Lévy flight firefly algorithm. Comput. Electr. Eng. 2018, 70, 476–495. [Google Scholar] [CrossRef]
  25. Mlakar, U.; Potočnik, B.; Brest, J. A hybrid differential evolution for optimal multilevel image thresholding. Expert Syst. Appl. 2016, 65, 221–232. [Google Scholar] [CrossRef]
  26. Kotte, S.; Pullakura, R.K.; Injeti, S.K. Optimal multilevel thresholding selection for brain MRI image segmentation based on adaptive wind driven optimization. Measurement 2018, 130, 340–361. [Google Scholar] [CrossRef]
  27. Wolpert, D.H.; Macready, W.G. No free lunch theorems for optimization. Evolut. Comput. IEEE Trans. 1997, 1, 67–82. [Google Scholar] [CrossRef] [Green Version]
  28. Heidari, A.A.; Mirjalili, S.; Faris, H.; Aljarah, I.; Mafarja, M.; Chen, H. Harris hawks optimization: Algorithm and applications. Future Gener. Comput. Syst. 2019, 97, 849–872. [Google Scholar] [CrossRef]
  29. Das, S.; Mullick, S.S.; Suganthan, P.N. Recent advances in differential evolution—An updated survey. Swarm Evol. Comput. 2016, 27, 1–30. [Google Scholar] [CrossRef]
  30. Bhandari, A.K.; Kumar, A.; Singh, G.K. Tsallis entropy based multilevel thresholding for colored satellite image segmentation using evolutionary algorithms. Expert Syst. Appl. 2015, 42, 8707–8730. [Google Scholar] [CrossRef]
  31. Elaziz, M.A.; Oliva, D.; Ewees, A.A.; Xiong, S. Multi-level thresholding-based grey scale image segmentation using multi-objective multi-verse optimizer. Expert Syst. Appl. 2019, 125, 112–129. [Google Scholar] [CrossRef]
  32. Agrawal, S.; Panda, R.; Bhuyan, S.; Panigrahi, B.K. Tsallis entropy based optimal multilevel thresholding using cuckoo search algorithm. Swarm Evol. Comput. 2013, 11, 16–30. [Google Scholar] [CrossRef]
  33. Xiong, G.; Zhang, J.; Yuan, X.; Shi, D.; He, Y.; Yao, G. Parameter extraction of solar photovoltaic models by means of a hybrid differential evolution with whale optimization algorithm. Sol. Energy 2018, 176, 742–761. [Google Scholar] [CrossRef]
  34. Jadon, S.S.; Tiwari, R.; Sharma, H.; Bansal, J.C. Hybrid Artificial Bee Colony algorithm with Differential Evolution. Appl. Soft Comput. 2017, 58, 11–24. [Google Scholar] [CrossRef]
  35. Elaziz, M.A.; Xiong, S.; Jayasena, K.P.N.; Li, L. Task scheduling in cloud computing based on hybrid moth search algorithm and differential evolution. Knowl. Based Syst. 2019, 169, 39–52. [Google Scholar] [CrossRef]
  36. Xiong, G.; Zhang, J.; Shi, D.; He, Y. Parameter extraction of solar photovoltaic models using an improved whale optimization algorithm. Energy Convers. Manag. 2018, 174, 388–405. [Google Scholar] [CrossRef]
  37. Zhao, F.; Xue, F.; Zhang, Y.; Ma, W.; Zhang, C.; Song, H. A hybrid algorithm based on self-adaptive gravitational search algorithm and differential evolution. Expert Syst. Appl. 2018, 113, 515–530. [Google Scholar] [CrossRef]
  38. Xu, L.; Jia, H.; Lang, C.; Peng, X.; Sun, K. A Novel Method for Multilevel Color Image Segmentation Based on Dragonfly Algorithm and Differential Evolution. IEEE Access 2019, 7, 19502–19538. [Google Scholar] [CrossRef]
  39. Gao, W.; Huang, L.; Wang, J.; Liu, S.; Qin, C. Enhanced artificial bee colony algorithm through differential evolution. Appl. Soft Comput. 2016, 48, 137–150. [Google Scholar] [CrossRef]
  40. Ibrahim, R.A.; Elaziz, M.A.; Lu, S. Chaotic opposition-based grey-wolf optimization algorithm based on differential evolution and disruption operator for global optimization. Expert Syst. Appl. 2018, 108, 1–27. [Google Scholar] [CrossRef]
  41. Xu, Y.; Chen, H.; Luo, J.; Zhang, Q.; Jiao, S.; Zhang, X. Enhanced Moth-flame optimizer with mutation strategy for global optimization. Inf. Sci. 2019, 492, 181–203. [Google Scholar] [CrossRef]
  42. Hijime, K.; Isao, O.; Shingenobu, K. Theoretical Analysis of the Unimodal Normal Distribution Crossover for Real-coded Genetic Algorithms. Trans. Soc. Instrum. Control Eng. 2002, 2, 187–194. [Google Scholar]
  43. Akimoto, Y.; Nagata, Y.; Sakuma, J.; Ono, I.; Kobayashi, S. Proposal and evaluation of adaptive real-coded crossover AREX. Trans. Jpn. Soc. Artif. Intell. 2009, 24, 446–458. [Google Scholar] [CrossRef]
  44. Ariyarit, A.; Kanazaki, M. Multi-modal distribution crossover method based on two crossing segments bounded by selected parents applied to multi-objective design optimization. J. Mech. Sci. Technol. 2015, 29, 1443–1448. [Google Scholar] [CrossRef]
  45. Eshelman, L.J.; Schaffer, J.D. Real-Coded Genetic Algorithms and Interval-Schemata, Foundations of Genetic Algorithms 2; Morgan Kaufman Publisher: San Mateo, CA, USA, 1993; pp. 187–202. [Google Scholar] [CrossRef]
  46. Jia, H.; Lang, C.; Oliva, D.; Song, W.; Peng, X. Hybrid Grasshopper Optimization Algorithm and Differential Evolution for Multilevel Satellite Image Segmentation. Remote Sens. 2019, 11, 1134. [Google Scholar] [CrossRef]
  47. Gill, H.S.; Khehra, B.S.; Singh, A.; Kaur, L. Teaching-learning-based optimization algorithm to minimize cross entropy for Selecting multilevel threshold values. Egypt. Inf. J. 2019, 20, 11–25. [Google Scholar] [CrossRef]
  48. Bohat, V.K.; Arya, K.V. A new heuristic for multilevel thresholding of images. Expert Syst. Appl. 2019, 117, 176–203. [Google Scholar] [CrossRef]
  49. Kotte, S.; Kumar, P.R.; Injeti, S.K. An efficient approach for optimal multilevel thresholding selection for gray scale images based on improved differential search algorithm. Ain Shams Eng. J. 2018, 9, 1043–1067. [Google Scholar] [CrossRef] [Green Version]
  50. Bhandari, A.K. A novel beta differential evolution algorithm-based fast multilevel thresholding for color image segmentation. Neural Comput. Appl. 2018, 1–31. [Google Scholar] [CrossRef]
  51. Landsat Imagery Courtesy of NASA Goddard Space Flight Center and U.S. Geological Survey. Available online: https://landsat.visibleearth.nasa.gov/index.php?&p=1 (accessed on 7 October 2018).
  52. Ewees, A.A.; Elaziz, M.A.; Oliva, D. Image segmentation via multilevel thresholding using hybrid optimization algorithms. J. Electron. Imaging 2018, 27. [Google Scholar] [CrossRef]
  53. Shen, L.; Fan, C.; Huang, X. Multi-Level Image Thresholding Using Modified Flower Pollination Algorithm. IEEE Access 2018, 6, 30508–30519. [Google Scholar] [CrossRef]
  54. Frank, W. Individual Comparisons of Grouped Data by Ranking Methods. J. Econ. Entomol. 1946, 39, 269–270. [Google Scholar]
  55. Friedman, M. The use of ranks to avoid the assumption of normality implicit in the analysis of variance. J. Am. Stat. Assoc. 1937, 32, 676–701. [Google Scholar] [CrossRef]
  56. Derrac, J.; García, S.; Molina, D.; Herrera, F. A practical tutorial on the use of nonparametric statistical tests as a methodology for comparing evolutionary and swarm intelligence algorithms. Swarm Evol. Comput. 2011, 1, 3–18. [Google Scholar] [CrossRef]
  57. Mousavirad, S.J.; Ebrahimpour-Komleh, H. Multilevel image thresholding using entropy of histogram and recently developed population-based metaheuristic algorithms. Evol. Intell. 2017, 10, 45–75. [Google Scholar] [CrossRef]
  58. Jia, H.; Xing, Z.; Song, W. Three Dimensional Pulse Coupled Neural Network Based on Hybrid Optimization Algorithm for Oil Pollution Image Segmentation. Remote Sens. 2019, 11, 1046. [Google Scholar] [CrossRef]
Figure 1. Different stages of Harris hawks optimization (HHO).
Figure 1. Different stages of Harris hawks optimization (HHO).
Remotesensing 11 01421 g001
Figure 2. Variation of the disturbance term δ under different values of parameter α .
Figure 2. Variation of the disturbance term δ under different values of parameter α .
Remotesensing 11 01421 g002
Figure 3. Variation of the escaping energy of the rabbit E under different values of parameter α .
Figure 3. Variation of the escaping energy of the rabbit E under different values of parameter α .
Remotesensing 11 01421 g003
Figure 4. Variation of the escaping energy of the rabbit E using Levy distribution.
Figure 4. Variation of the escaping energy of the rabbit E using Levy distribution.
Remotesensing 11 01421 g004
Figure 5. Original satellite images and the corresponding histograms for each of color channels (red, green and blue).
Figure 5. Original satellite images and the corresponding histograms for each of color channels (red, green and blue).
Remotesensing 11 01421 g005aRemotesensing 11 01421 g005b
Figure 6. The segmented results of satellite images obtained by different algorithms.
Figure 6. The segmented results of satellite images obtained by different algorithms.
Remotesensing 11 01421 g006aRemotesensing 11 01421 g006bRemotesensing 11 01421 g006c
Figure 7. The boxplot of each method on satellite images (K = 20).
Figure 7. The boxplot of each method on satellite images (K = 20).
Remotesensing 11 01421 g007
Figure 8. The convergence curves for fitness function on satellite images (K = 20).
Figure 8. The convergence curves for fitness function on satellite images (K = 20).
Remotesensing 11 01421 g008aRemotesensing 11 01421 g008b
Figure 9. Original oil pollution images and the corresponding histograms for each of color channels (red, green and blue).
Figure 9. Original oil pollution images and the corresponding histograms for each of color channels (red, green and blue).
Remotesensing 11 01421 g009
Figure 10. The segmented results and ground truth of the oil pollution images. Panels (ad) indicate the segmented images, (eh) show the ground truth, (il) indicate the classification result.
Figure 10. The segmented results and ground truth of the oil pollution images. Panels (ad) indicate the segmented images, (eh) show the ground truth, (il) indicate the classification result.
Remotesensing 11 01421 g010
Table 1. Sum of ranks obtained by the proposed method (K = 5) under different of parameter values.
Table 1. Sum of ranks obtained by the proposed method (K = 5) under different of parameter values.
F = 0.1 F = 0.3 F = 0.5 F = 0.7 F = 0.9
α = 2 89412043109
α = 2.5 8019936112
α = 3 92656070115
Table 2. Parameters of the algorithms.
Table 2. Parameters of the algorithms.
No.AlgorithmParameter SettingYearReference
1.DHHO/M α = 2.5 , S F = 0.5
2.DHHO α = 2.5 ( control   parameter )
3.HHO/M S F = 0.5 ( scaling   factor )
4.HHO E [ 0 , 2 ] ( energy   of   a   rabbit ) 2019[28]
5.TLBO T F = 1 ( teaching   factor ) 2019[47]
6.WOA-TH G 0 = 40 , a 0 = 13 2019[48]
7.IDSA2018[49]
8.BDE a [ 0 , 1 ]   ( beta   distribution   parameter ) 2018[52]
9.MGOA β = 0.8   ( Levy   flight   parameter ) 2019[20]
10.MABC K = 300   ( chaotic   iteration ) 2015[13]
11.MFPA p [ 0.2 , 0.8 ]   ( switch   possibility ) 2018[53]
12.GWO a [ 0 , 2 ]   ( control   parameter ) 2017[7]
Table 3. Performance measures of the multilevel image segmentation methods.
Table 3. Performance measures of the multilevel image segmentation methods.
No.MeasuresFormulationReference
1.Average fitness function value Average = i = 1 N f i N [40]
2.Standard deviation (Std) Std = 1 N 1 i = 1 N ( f i Average ) 2 [46]
3.Peak signal to noise ratio (PSNR) PSNR = 10 log 10 ( 255 2 MSE ) [18]
4.Mean squared error (MSE) MSE = 1 M N i = 1 M j = 1 N [ I ( i , j ) K ( i , j ) ] 2 [18]
5.Structural similarity index (SSIM) SSIM = ( 2 μ x μ y + c 1 ) ( 2 σ x y + c 2 ) ( μ x 2 + μ y 2 + c 1 ) ( σ x 2 + σ y 2 + c 2 ) [53]
6.Feature similarity index (FSIM) FSIM = x Ω S L ( x ) · P C m ( x ) x Ω P C m ( x ) [23]
7.Average computation time T i m e = i = 1 N t i m e i N [17]
8.Wilcoxon’s rank sum test R + = d i > 0 r a n k ( d i ) + 1 2 d i = 0 r a n k ( d i )
R = d i < 0 r a n k ( d i ) + 1 2 d i = 0 r a n k ( d i )
[54]
9.Friedman test F f = 12 n k ( k + 1 ) [ j R j 2 k ( k + 1 ) 2 4 ] [55]
Table 4. The experimental results of different algorithms on mathematical test functions.
Table 4. The experimental results of different algorithms on mathematical test functions.
FDHHO/MHHOTLBOWOA-THIDSABDE
F1Best6.6218 × 101205.311 × 10−10654.91031.46393.4749 × 1032.7038 × 10−7
Mean2.7498 × 101051.7563 × 10−884.7014 × 1022.60781.0831 × 1049.2519 × 10−6
Worst5.2688 × 10975.2688 × 10−871.4603 × 1033.53832.1372 × 1044.5017 × 10−5
Std3.1671 × 10919.6195 × 10−883.1022 × 1020.456943.8111 × 1031.1476 × 10−5
F9Best002.2902 × 1021.5848 × 1021.7535 × 1023.1335 × 10−7
Mean002.7583 × 1022.1586 × 1022.0861 × 1022.1983 × 10−5
Worst003.1682 × 1022.8732 × 1022.4096 × 1021.597 × 10−4
Std0023.321336.192616.78573.7214 × 10−5
F14Best0.9980.9980.9980.9980.9981.0311
Mean1.0691.62541.36121.16242.23291.992
Worst2.80085.92885.92885.928819.23070.1815
Std0.330431.49811.59790.900243.42641.3037
Table 5. The experimental results of HHO variants on satellite images.
Table 5. The experimental results of HHO variants on satellite images.
ImageKDHHO/MDHHOHHO/MHHODHHO+DHHO
MeanStdMeanStdMeanStdMeanStdMeanStdMeanStd
1212.6254012.6253 1.00 × 10−412.6254012.6253 1.00 × 10−412.62538.46 × 10−512.6217 8.20 × 10−3
315.9977 1.40 × 10−315.99901.70 × 10−315.99901.40 × 10−315.9971 1.70 × 10−315.9984 1.80 × 10−315.9981 1.40 × 10−3
419.1537019.1537019.1537019.1537019.1537019.15370
522.02562.00 × 10−422.0254 1.00 × 10−422.0254 2.00 × 10−422.0253 5.60 × 10−522.0254 2.00 × 10−422.0254 2.00 × 10−4
3212.6999012.6999012.6999012.6999012.6999012.69990
315.84001.30 × 10−315.84001.30 × 10−315.84001.30 × 10−315.84001.60 × 10−315.84001.30 × 10−315.84001.30 × 10−3
418.9080 018.9080 8.45 × 10−518.90812.00 × 10−418.9080 1.00 × 10−418.9080 1.54 × 10−518.9080 2.00 × 10−4
521.69372.16 × 10−521.6589 4.75 × 10−221.6936 2.00 × 10−421.6936 3.07 × 10−521.6936 0.000121.6936 1.00 × 10−4
5212.2672012.2672012.2672012.2672012.2672012.26720
315.2240 1.99 × 10−215.23883.38 × 10−215.23881.00 × 10−315.2089 1.71 × 10−215.2314 1.63 × 10−215.2089 1.71 × 10−2
418.0756 9.90 × 10−518.0676 1.82 × 10−218.0675 1.81 × 10−218.0599 3.02 × 10−218.07571.00 × 10−418.06761.82 × 10−2
520.73871.88 × 10−220.7071 5.48 × 10−220.68175.05 × 10−220.7265 3.01 × 10−220.7277 3.17 × 10−220.6692 4.47 × 10−2
7212.4226012.4226012.4226012.4226012.4226012.42260
315.47921.99 × 10−915.47921.99 × 10−915.47901.99 × 10−915.47921.99 × 10−915.4790 1.99 × 10−915.4790 1.99 × 10−9
418.31821.95 × 10−518.3181 5.07 × 10−518.31811.95 × 10−518.3181 3.19 × 10−518.3181 2.38 × 10−518.3181 1.95 × 10−5
521.05013.37 × 10−521.04993.00 × 10−421.033163.74 × 10−221.0500 2.00 × 10−421.0500 1.00 × 10−420.9918 5.51 × 10−2
Rank1(12)3(8)2(9)4(6)4(6)6(5)
Table 6. Comparison of optimal thresholds for different algorithms using Kapur’s entropy at 2– 5 levels.
Table 6. Comparison of optimal thresholds for different algorithms using Kapur’s entropy at 2– 5 levels.
ImageKDHHO/MHHOTLBO
RGBRGBRGB
1258,171111,189117,18758,171111,189117,18758,171111,189117,187
356,123,18848,122,19242,128,19156,123,18848,122,19242,128,19156,123,18848,122,19285,143,203
453,110,163,20945,102,156,20741,95,149,20353,110,163,20945,102,156,20741,95,149,20353,110,163,20945,102,156,20741,95,149,203
545,86,130,173,21344,90,133,174,21637,85,129,170,21042,82,126,171,21344,90,133,174,21637,85,130,171,21042,82,126,171,21344,90,133,174,21637,85,129,170,210
2294,16498,16880,14694,16498,16880,14694,16498,16880,146
378,141,19966,132,19270,131,18578,141,19966,132,19270,131,18578,141,19966,132,19270,131,185
467,120,169,21757,107,152,19553,102,148,20067,120,169,21757,107,152,19556,107,154,20067,120,169,21757,107,152,19553,102,148,200
553,97,140,180,22055,101,145,186,22245,87,128,167,20253,97,140,180,22055,101,145,186,22244,85,126,167,20253,96,138,180,22055,101,145,186,22244,85,126,167,202
3286,16385,16280,16286,16385,16280,16286,16385,16280,162
384,152,20781,150,20765,117,17084,152,20781,150,20778,152,21084,152,20781,150,20765,117,170
470,117,165,21669,117,166,21665,114,165,21870,117,165,21667,115,165,21665,114,165,21870,117,165,21669,117,166,21665,114,165,218
536,82,125,168,21645,85,126,168,21657,93,132,171,21836,82,125,168,21645,85,126,168,21618,65,114,165,21836,82,125,168,21645,85,126,168,21658,95,134,172,218
4292,17996,17182,14692,17996,17182,14692,17996,17182,146
363,120,18272,118,17179,137,17663,120,18272,118,17179,137,17663,120,18272,118,17179,137,176
447,89,132,18270,114,159,19450,91,137,17647,89,132,18270,114,159,19450,91,137,17647,89,132,18270,114,159,19450,91,137,176
546,87,129,176,20952,88,124,164,20531,66,99,142,18346,87,129,176,20949,86,123,164,20549,90,134,162,19846,87,129,176,20948,85,122,164,20547,91,136,164,198
5297,15895,15681,14497,15895,15681,14497,15895,15681,144
397,148,18338,96,15683,138,18197,148,18338,96,15628,82,14497,148,18396,150,18783,138,180
474,109,150,18440,96,150,18728,83,138,18174,109,150,18438,96,150,18728,83,138,18178,112,150,18438,96,150,18728,83,138,181
564,93,121,153,18638,79,112,152,18824,52,88,138,18164,93,121,153,18638,82,114,152,18828,80,112,144,18320,74,109,148,18338,79,112,152,18828,80,112,144,183
6274,15077,15280,15474,15077,15280,15474,15077,15280,154
362,122,18368,129,18772,130,18762,122,18368,129,18772,130,18762,122,18368,129,18772,130,187
451,98,147,19461,105,151,19761,106,152,19751,98,147,19461,105,151,19761,106,152,19751,98,147,19461,105,151,19761,106,152,197
549,95,143,190,22861,103,148,191,22858,102,148,191,22749,95,143,190,22861,103,148,191,22858,102,148,191,22750,95,145,191,23061,103,148,191,22860,106,152,194,227
7296,16997,16064,12996,16997,16064,12996,16997,16064,129
380,135,19081,129,17755,105,15480,135,19081,129,17755,105,15480,135,19081,129,17755,105,154
467,113,159,20469,108,148,18854,101,148,19768,114,159,20469,108,148,18854,101,148,19767,113,159,20469,108,148,18854,101,148,197
557,94,132,171,20968,107,146,185,22648,87,125,163,19757,94,132,171,20924,70,110,149,18850,90,129,168,20615,68,114,159,20468,107,146,185,22648,87,125,163,197
8271,197111,203125,16871,197111,203126,17971,197111,203126,179
370,138,19796,150,20484,132,17970,138,19796,150,20426,86,14270,138,19796,150,20484,132,179
469,114,155,19756,106,152,20426,84,132,17969,114,155,19756,106,152,20427,84,132,17969,114,155,19754,105,152,20426,85,132,179
569,111,151,189,21954,97,138,171,20726,82,126,163,20369,111,151,189,21954,97,138,171,20726,82,126,163,19969,111,151,189,21956,105,148,189,21526,82,126,163,203
ImageKWOA-THIDSABDE
RGBRGBRGB
1258,171111,189117,18758,171111,189117,18758,171111,189117,187
356,123,18848,122,19284,142,20356,123,18848,122,19285,143,20356,123,18848,122,19285,143,203
453,110,163,20945,102,156,20741,95,149,20353,110,163,20945,102,156,20741,94,149,20353,111,162,20945,102,156,20741,95,149,203
542,82,126,171,21344,90,133,174,21637,85,129,170,21045,86,128,171,21144,90,133,174,21440,89,133,172,21042,82,126,171,21344,90,133,174,21637,85,129,170,210
2294,16498,16880,14694,16498,16880,14694,16498,16880,146
378,141,19966,132,19270,131,18578,141,19966,132,19270,131,18578,141,19966,132,19270,131,185
467,120,169,21757,107,152,19557,109,154,20060,109,155,19957,107,152,19556,107,154,20067,120,169,21757,107,152,19553,102,148,200
553,96,138,180,22055,101,145,186,22245,88,131,169,20253,95,137,180,22053,97,138,178,21545,87,128,167,20252,93,134,174,21755,101,145,186,22242,84,126,167,202
3286,16385,16280,16286,16385,16280,16286,16385,16280,162
384,152,20781,150,20765,117,17084,152,20781,150,20778,152,21084,152,20781,150,20765,117,170
470,117,165,21669,117,166,21665,114,165,21870,117,165,21669,117,166,21665,115,165,21870,117,165,21669,117,166,21665,114,165,218
536,82,125,168,21645,85,126,168,21658,95,134,172,21834,81,123,167,21645,85,126,168,21658,95,133,171,21836,82,125,168,21645,85,126,168,21658,95,134,172,218
4292,17996,17182,14692,17996,17182,14692,17996,17182,146
363,120,18272,118,17179,137,17662,120,18272,118,17179,137,17663,120,18272,118,17179,137,176
447,89,132,18271,115,164,20550,91,137,17662,118,176,20970,114,159,19450,91,137,17647,89,132,18270,113,158,19450,91,137,176
547,89,130,176,20949,86,123,164,20549,90,134,162,19842,84,124,176,20948,85,121,159,19449,90,134,162,19746,88,130,176,20952,87,122,160,19849,90,134,162,198
5297,15895,15681,14497,15895,15681,14497,15895,15781,144
397,148,18338,96,15683,138,18197,148,18338,96,15683,138,18197,148,18339,96,15683,138,181
474,109,150,18438,96,150,18728,83,138,18174,109,150,18438,96,150,18728,83,138,18074,109,150,18438,96,150,18728,83,138,182
521,78,112,150,18438,82,114,152,18828,80,112,144,18321,78,112,150,18438,79,113,151,18828,80,112,144,18121,74,109,150,18438,79,112,152,18830,83,113,156,197
6274,15077,15280,15474,15077,15280,15474,15077,15280,154
362,122,18368,129,18772,130,18762,122,18368,129,18772,131,18764,125,18468,129,18772,130,187
451,98,147,19461,105,151,19761,106,152,19752,99,147,19461,105,151,19662,107,152,19851,98,147,19461,105,151,19761,106,152,197
549,95,143,190,22861,103,148,191,22858,102,148,191,22755,99,146,190,22861,103,148,191,22851,86,123,161,20049,95,143,190,22861,103,148,191,22858,102,148,191,227
7296,16997,16064,12996,16997,16064,12996,16997,16064,129
380,135,19081,129,17755,105,15480,135,19081,129,17755,105,15480,135,19081,129,17755,105,154
468,114,159,20470,110,149,18854,101,148,19768,114,159,20469,108,148,18854,101,148,19767,113,159,20469,108,148,18854,101,148,197
557,94,132,171,20968,107,146,185,22650,90,129,168,20658,96,134,172,21168,107,146,185,22648,86,124,163,19757,94,132,171,20968,107,146,185,22648,87,125,163,197
8271,197111,203126,17971,197111,203126,17971,197111,203126,179
370,138,19796,150,20484,132,17970,138,19796,150,20484,132,17970,138,19796,150,20484,132,179
469,114,155,19756,106,152,20426,84,132,17969,115,157,19856,106,152,20479,125,163,20369,114,155,19756,106,152,20479,120,152,193
569,111,151,189,21956,104,146,189,21526,79,120,152,19368,111,151,189,21953,86,119,158,20426,79,120,152,19369,111,151,189,21954,103,148,190,21627,82,125,163,199
Table 7. The average fitness values and std values of different algorithms at 2, 3, 4, and 5 levels.
Table 7. The average fitness values and std values of different algorithms at 2, 3, 4, and 5 levels.
ImageKDHHO/MHHOTLBOWOA-THIDSABDE
MeanStdMeanStdMeanStdMeanStdMeanStdMeanStd
1212.62543.68 × 10−912.62544.73 × 10−312.62543.68 × 10−912.62543.68 × 10−912.62541.20 × 10−512.62543.68 × 10−9
315.99651.31 × 10−315.99961.61 × 10−315.99961.67 × 10−315.99961.31 × 10−315.99971.52 × 10−315.99971.61 × 10−3
419.15371.49 × 10−419.15373.86 × 10−519.15373.63 × 10−419.15373.68 × 10−919.1536.46 × 10−419.15235.59 × 10−3
522.02561.15 × 10−422.02512.11 × 10−422.02562.40 × 10−422.0253.54 × 10−422.02362.26 × 10−321.99031.42 × 10−2
2212.24463.68 × 10−912.24461.61 × 10−412.24463.68 × 10−912.24463.68 × 10−912.24463.68 × 10−912.24462.78 × 10−4
315.45531.51 × 10−415.45532.03 × 10−415.45531.92 × 10−515.45532.03 × 10−415.45539.03 × 10−415.45534.76 × 10−3
418.48391.84 × 10−418.48352.90 × 10−418.48394.84 × 10−418.48331.86 × 10−418.48132.30 × 10−318.48393.16 × 10−2
521.31811.64 × 10−421.31792.35 × 10−321.3182.65 × 10−321.31781.02 × 10−321.31251.27 × 10−221.31566.00 × 10−3
3212.69991.84 × 10−912.69991.84 × 10−912.69991.84 × 10−912.69991.84 × 10−912.69991.84 × 10−912.69995.73 × 10−5
315.84064.84 × 10−515.83771.00 × 10−315.84067.36 × 10−415.84061.18 × 10−315.83772.55 × 10−315.84061.00 × 10−3
418.90814.91 × 10−518.90796.67 × 10−418.9081018.90814.89 × 10−518.9082.96 × 10−318.90813.40 × 10−3
521.69372.00 × 10−521.6092.98 × 10−221.69343.81 × 10−521.69377.03 × 10−421.69351.19 × 10−321.69372.49 × 10−3
4212.4583.68 × 10−912.4583.68 × 10−912.4584.04 × 10−512.4583.68 × 10−912.4582.12 × 10−512.4581.21 × 10−3
315.58033.68 × 10−915.58033.68 × 10−915.58033.05 × 10−515.58033.68 × 10−915.58033.06 × 10−315.58033.07 × 10−4
418.48359.69 × 10−418.48341.39 × 10−318.48341.56 × 10−318.48348.01 × 10−518.48061.26 × 10−318.48147.04 × 10−4
521.22555.82 × 10−321.20898.52 × 10−321.22115.99 × 10−321.22546.81 × 10−321.21551.02 × 10−221.22116.13 × 10−3
5212.26721.84 × 10−912.26721.84 × 10−912.26721.84 × 10−912.26721.84 × 10−912.26721.84 × 10−912.26651.84 × 10−9
315.23931.12 × 10−315.20231.93 × 10−215.14981.67 × 10−215.23931.30 × 10−215.23931.30 × 10−215.23831.32 × 10−3
418.07581.63 × 10−218.07582.38 × 10−218.07561.24 × 10−418.07581.47 × 10−218.07573.41 × 10−218.07545.10 × 10−3
520.75971.78 × 10−220.69664.23 × 10−220.75062.79 × 10−220.66223.78 × 10−220.75883.35 × 10−220.67683.75 × 10−2
6212.64593.68 × 10−912.64593.68 × 10−912.64593.68 × 10−912.64593.68 × 10−912.64595.35 × 10−612.64593.16 × 10−5
315.83483.68 × 10−915.83483.68 × 10−915.83483.68 × 10−915.83483.68 × 10−915.83485.86 × 10−515.83467.96 × 10−4
418.80996.43 × 10−618.80996.43 × 10−618.80994.23 × 10−518.80993.03 × 10−518.80895.43 × 10−418.80991.48 × 10−3
521.58497.36 × 10−521.58491.47 × 10−421.58178.43 × 10−321.58497.23 × 10−521.55561.26 × 10−221.58491.02 × 10−2
7212.42263.68 × 10−912.42263.68 × 10−912.42263.68 × 10−912.42263.68 × 10−912.42263.68 × 10−912.42263.98 × 10−5
315.47925.52 × 10−915.47925.52 × 10−915.47925.52 × 10−915.47925.52 × 10−915.47925.37 × 10−515.47924.44 × 10−5
418.31822.12 × 10−518.31812.41 × 10−218.31826.18 × 10−518.31812.84 × 10−518.31814.27 × 10−318.31823.07 × 10−3
521.05027.07 × 10−520.96633.46 × 10−220.92776.17 × 10−421.05012.16 × 10−221.04922.60 × 10−221.05023.01 × 10−2
8212.08561.33 × 10−412.08532.78 × 10−312.08531.75 × 10−312.08531.33 × 10−512.08531.77 × 10−312.08538.05 × 10−4
315.33345.52 × 10−915.19373.60 × 10−215.33345.52 × 10−915.33345.52 × 10−915.33349.91 × 10−415.33344.08 × 10−4
418.25727.30 × 10−418.2565.24 × 10−318.25521.02 × 10−318.25624.88 × 10−418.25732.32 × 10−318.25472.94 × 10−3
521.01811.25 × 10−321.01553.78 × 10−221.01731.49 × 10−321.01423.70 × 10−221.01631.57 × 10−221.01362.38 × 10−3
Table 8. The PSNR values and MSE values of different algorithms at 2, 3, 4, and 5 levels.
Table 8. The PSNR values and MSE values of different algorithms at 2, 3, 4, and 5 levels.
ImageKDHHO/MHHOTLBOWOA-THIDSABDE
PSNRMSEPSNRMSEPSNRMSEPSNRMSEPSNRMSEPSNRMSE
1213.45683409.628813.50593409.628813.45683409.628813.45683409.628813.45683409.628813.45683409.6288
320.8043541.728118.7437541.728118.74371101.204818.74371079.856318.89121101.204818.86681101.2048
423.5379286.566823.5379288.250423.5379288.250423.5379288.250423.4893286.856323.4844288.2504
525.3432190.037825.2873190.426325.3432190.037825.2458190.037825.3091199.529924.8902192.3651
2213.7612793.635313.7612793.635313.7612793.635313.7612793.635313.7612793.635313.7612793.6353
317.16561274.882517.16561274.882517.16561274.882517.16561274.882517.16561274.882517.16561274.8825
419.8508690.969319.6349725.615519.8508690.969319.5626738.419120.126638.700519.8508690.9693
522.2059393.448521.9532415.659721.962414.815221.8581424.254822.073403.52521.8664423.4967
3214.60992250.56314.60992250.56314.60992250.56314.60992250.56314.60992250.56314.60992250.563
317.73431096.477217.07371316.394317.07371316.394317.07371316.394317.73431096.477217.07371316.3943
419.5277729.859619.4491743.221819.4491743.221819.4491743.221819.4482743.400519.4491743.2218
524.4587238.874222.7205395.943322.6857401.408422.6857401.408422.6956401.039322.6857401.4084
4214.15962518.004914.15962518.004914.15962518.004914.15962518.004914.15962518.004914.15962518.0049
317.30511316.175617.30511316.175617.30511316.175617.30511316.175617.34861308.801417.30511316.1756
420.3568652.828820.3524653.850120.3524653.850120.2761672.029919.4048764.848920.3524653.8501
524.0586266.30822.6141375.537922.8347354.150522.5577379.777122.9103356.313722.3481392.8026
5212.31493817.531612.31493817.531612.31493817.531612.31493817.531612.31493817.531612.32473809.1362
314.37532377.232813.93342642.142214.33882396.536414.37532377.232814.37532377.232814.37472377.5557
416.9041355.906916.85861368.813516.75361413.043216.85861368.813516.81581381.357416.85161370.5727
518.529941.672517.49321166.119417.61211131.20317.6611117.257517.56921142.552917.46111174.9531
6215.40961875.646515.40961875.646515.40961875.646515.40961875.646515.40961875.646515.40961875.6465
317.35861204.886917.35861204.886917.35861204.886917.35861204.886917.35761205.208917.32171213.4404
418.6145906.632518.6145906.632518.6145906.632518.6145906.632518.5634916.862318.6145906.6325
520.6141586.837420.6141586.837420.4926603.11820.6141586.837420.1165651.138420.6141586.8374
7212.96583316.678112.96583316.678112.96583316.678112.96583316.678112.96583316.678112.96583316.6781
316.25941542.330216.25941542.330216.25941542.330216.25941542.330216.25941542.330216.25941542.3302
419.9285710.12719.8763720.296519.9285710.12719.8208732.327319.8763720.296519.9285710.127
522.7765353.407921.725476.874822.5737406.434521.5981484.458421.6863481.355621.725476.8748
8214.24762953.528114.24522956.150414.24522956.150414.24522956.150414.24522956.150414.24522956.1504
318.4536854.952816.95441753.365716.95441753.365716.95441753.365716.95441753.365716.95441753.3657
424.6745222.572124.5916227.41724.5902227.479424.5902227.479420.78671039.720420.80961035.361
525.5344183.068925.2902193.252925.0615204.079725.1686199.83225.2937193.114725.1855198.1814
Table 9. The SSIM values and FSIM values of different algorithms at 2, 3, 4, and 5 levels.
Table 9. The SSIM values and FSIM values of different algorithms at 2, 3, 4, and 5 levels.
ImageKDHHO/MHHOTLBOWOA-THIDSABDE
SSIMFSIMSSIMFSIMSSIMFSIMSSIMFSIMSSIMFSIMSSIMFSIM
120.32880.54720.33190.5490.32880.54720.32880.54720.32880.54720.32880.5472
30.58930.73870.570.71480.570.71480.570.71480.5770.72030.57760.7205
40.74150.8270.74150.8270.74150.8270.74150.8270.74150.82690.74070.8261
50.79450.87630.79450.87610.79450.87610.79170.87380.7940.87550.78380.866
220.26960.55270.26960.55270.26960.55270.26960.55270.26960.55270.26960.5527
30.43760.64760.43760.64760.43760.64760.43760.64760.43760.64760.43760.6476
40.57370.74360.56340.73230.57370.73820.560.730.58290.73820.57370.7382
50.66690.79810.65770.79280.65840.79360.65360.790.6630.79650.65430.7902
320.36430.66270.36430.66270.36430.66270.36430.66270.36430.66270.36430.6627
30.42180.72390.42090.72170.42180.72170.42180.72170.42090.72390.42180.7217
40.48860.78080.48630.78010.48630.78010.48630.78010.48640.780.48630.7801
50.66930.84310.57950.8370.57850.8370.57850.8370.57880.8370.57850.837
420.40430.64630.40430.64630.40430.64630.40430.64630.40430.64630.40430.6463
30.57460.75840.57340.7580.57340.7580.57340.7580.57340.7580.57340.758
40.73280.83310.73230.83310.73230.83310.72910.8320.6980.8130.73230.8336
50.83360.88690.7920.86830.79660.86920.79090.86780.79780.87110.78830.8655
520.48060.63890.48060.63890.48060.63890.48060.63890.48060.63890.48020.6385
30.54320.670.6190.68650.46930.65370.54320.670.54320.670.54350.67
40.64590.71420.64540.71380.63980.71170.64540.71380.64590.71370.64560.714
50.6940.75150.6790.74170.71360.75060.70930.74740.71090.74850.71320.7504
620.40530.66280.40530.66280.40530.66280.40530.66280.40530.66280.40530.6628
30.47410.74090.47410.74090.47410.74090.47410.74090.4740.74080.47170.74
40.52830.78270.52830.78270.52830.78270.52830.78270.52590.78220.52830.7827
50.55950.80960.55890.80920.55480.80750.55890.80920.55890.80920.55890.8092
720.40610.61780.40610.61780.40610.61780.40610.61780.40610.61780.40610.6178
30.54560.71310.54560.71310.54560.71310.54560.71310.54560.71310.54560.7131
40.65690.78590.65570.78580.65690.78590.6540.78530.65570.78580.65690.7859
50.75570.83720.72340.8290.74340.83060.71940.82610.72190.82820.72340.829
820.22640.48410.22610.4840.22610.4840.22610.4840.22610.4840.22610.484
30.53190.6590.38640.61960.38640.61960.38640.61960.38640.61960.38640.6196
40.65470.74850.65390.74730.65320.74580.65390.74730.51220.70490.51520.7099
50.69610.79230.68290.7770.67240.76660.68040.77520.68290.77710.67440.7686
Table 10. Results of Wilcoxon’s rank sum test over all available satellite images at 2, 3, 4, and 5 levels.
Table 10. Results of Wilcoxon’s rank sum test over all available satellite images at 2, 3, 4, and 5 levels.
ImageKDHHO/M versus
HHOTLBOWOA-THIDSABDE
phphphphph
12<0.051<0.0510.41970<0.051<0.051
3<0.051<0.051<0.051<0.051<0.051
4<0.051<0.051<0.0510.09570<0.051
5<0.051<0.051<0.051<0.051<0.051
22<0.051<0.051<0.051<0.051<0.051
3<0.051<0.051<0.051<0.0510.07840
40.06920<0.051<0.051<0.051<0.051
5<0.051<0.051<0.051<0.051<0.051
32<0.051<0.051<0.051<0.051<0.051
3<0.051<0.051<0.051<0.051<0.051
4<0.0510.21170<0.051<0.051<0.051
5<0.051<0.051<0.051<0.051<0.051
42<0.051<0.051<0.051<0.051<0.051
3<0.051<0.051<0.051<0.051<0.051
4<0.051<0.051<0.051<0.051<0.051
5<0.051<0.051<0.051<0.051<0.051
52<0.051<0.051<0.051<0.051<0.051
3<0.051<0.051<0.051<0.051<0.051
4<0.051<0.051<0.0510.26240<0.051
5<0.051<0.051<0.051<0.051<0.051
62<0.051<0.051<0.051<0.051<0.051
3<0.0510.05440<0.051<0.051<0.051
4<0.051<0.051<0.051<0.0510.13420
50.27440<0.051<0.051<0.051<0.051
72<0.051<0.051<0.051<0.051<0.051
3<0.051<0.051<0.051<0.051<0.051
4<0.051<0.051<0.051<0.051<0.051
5<0.051<0.051<0.051<0.051<0.051
82<0.051<0.051<0.051<0.051<0.051
30.09630<0.051<0.051<0.051<0.051
4<0.051<0.051<0.051<0.051<0.051
5<0.051<0.051<0.051<0.051<0.051
Table 11. Friedman rank of different algorithms considering all experimental data at 2, 3, 4, and 5 levels.
Table 11. Friedman rank of different algorithms considering all experimental data at 2, 3, 4, and 5 levels.
KDHHO/MHHOTLBOWOA-THIDSABDE
23.104173.447923.531253.416673.635423.86458
32.395833.854173.833333.416673.604173.89583
41.937503.614583.364583.614584.677083.79167
51.291673.843753.375004.218753.802084.46875
Overall2.182293.690103.526043.666673.929694.00521
Table 12. Results of Friedman test over all available satellite images at 2, 3, 4, and 5 levels.
Table 12. Results of Friedman test over all available satellite images at 2, 3, 4, and 5 levels.
KChi-Square Valuep-Value
220.10958904109591.19191051533946 × 10−3
340.69339111592631.08206724960841 × 10−7
465.67723342939488.10948996472364 × 10−13
595.85659898477163.94253504591917 × 10−19
Overall194.2744590780814.76396506070144 × 10−40
Table 13. The average computation time (in seconds) of different algorithms at 2, 3, 4, and 5 levels.
Table 13. The average computation time (in seconds) of different algorithms at 2, 3, 4, and 5 levels.
KDHHO/MHHOTLBOWOA-THIDSABDE
22.1981332.2122331.7974381.8449331.5650332.057214
32.2233212.2295331.8001741.8525671.6547332.151333
42.2637672.3392671.8247331.8963331.6637672.314633
52.3424782.3900671.8286331.9271331.6815472.323733
Table 14. The fitness values obtained by different algorithms at higher threshold levels.
Table 14. The fitness values obtained by different algorithms at higher threshold levels.
ImageKDHHO/MHHOTLBOWOA-THIDSABDE
11034.210934.209834.210334.210734.178934.1621
1543.941743.932643.904943.933743.678743.8603
2052.0552.030551.438552.049951.616451.6291
21033.349733.333333.347933.349433.262833.333
1542.829142.809342.760842.826142.636342.7508
2050.719150.589950.38150.6750.238850.3818
31033.897133.894533.894933.895333.852633.8872
1543.734443.590743.565943.73543.540943.6546
2051.846351.749751.245751.804251.551.443
41033.024833.004932.998633.00632.871532.9839
1542.564342.555742.518242.549542.167942.4222
2050.40450.291149.888450.345349.996750.2819
51032.181632.139332.120432.108132.108332.0863
1541.595541.516241.52641.493641.236541.315
2049.363149.168848.577449.227948.330649.1217
61033.605533.597833.576233.564133.558933.5871
1543.456743.352343.295143.439343.299243.4059
2051.70651.625151.134751.680651.228351.4454
71032.823232.805532.82532.819232.691132.6068
1542.347542.311142.227342.296541.967742.3149
2050.358350.280249.674150.33149.563850.1
81032.885832.81332.883332.850532.814832.825
1542.378442.257442.338442.350941.962442.2916
2050.309550.224650.09350.237949.904650.0384
Table 15. Comparison of objective value, PSNR, SSIM, and FSIM values based on Tsallis entropy approaches.
Table 15. Comparison of objective value, PSNR, SSIM, and FSIM values based on Tsallis entropy approaches.
ImageKObjective ValuePSNR
DHHO/MMGOAMABCDHHO/MMGOAMABC
1103.66623.66493.666230.509327.770531.1043
155.33165.31095.330234.262628.630333.4524
206.9956.9486.986636.837133.858834.1846
3103.66613.6653.666125.738324.998926.0445
155.33135.31415.329933.364528.933430.1357
206.99476.95166.986335.318532.920533.9627
5103.66543.65713.665422.093520.886921.9164
155.32885.30845.326225.595727.706127.5186
206.98966.8376.974134.147533.088128.4798
7103.66593.66433.665828.028426.150627.273
155.33055.3185.328733.529129.808227.8041
206.99246.93116.983636.865332.212434.0468
ImageKSSIMFSIM
DHHO/MMGOAMABCDHHO/MMGOAMABC
1100.91570.88110.92250.95790.92360.9622
150.95790.89560.94950.98060.93510.9762
200.97560.95410.96380.98860.97720.9773
3100.69040.68830.69670.91560.87710.9157
150.88270.83220.80950.96360.92410.9442
200.92510.90180.89110.9760.95570.9633
5100.8160.78540.81130.85220.82170.8501
150.88440.85060.83680.88880.88440.9139
200.92240.91070.90770.95780.93440.9298
7100.90010.82940.88610.93760.88230.9324
150.9520.91270.89160.97250.9390.9371
200.97490.94680.95550.98420.95490.9733
Table 16. Comparison of objective value, PSNR, SSIM, and FSIM values based on Otsu’s methods.
Table 16. Comparison of objective value, PSNR, SSIM, and FSIM values based on Otsu’s methods.
ImageKObjective ValuePSNR
DHHO/MMFPAGWODHHO/MMFPAGWO
2101364.98841349.2641363.190227.446230.006627.6228
151376.51481365.3311375.160132.302933.634232.4534
201381.03971373.86391379.553836.596535.541936.5278
4101304.58921289.35951304.349729.095929.221629.4061
151315.03891302.80591314.467433.324731.41334.5001
201319.06921311.41331317.341436.107436.61836.1379
6105027.67375013.65055025.514926.592823.572625.7915
155043.22555034.35965041.603632.916827.959429.957
205049.67735040.98385047.562835.099833.31733.8307
810934.5128925.9764934.347331.604931.174430.8108
15943.6876934.3932942.269735.853732.972535.0404
20947.3744941.3551945.291637.374437.212237.3125
ImageKSSIMFSIM
DHHO/MMFPAGWODHHO/MMFPAGWO
2100.92180.89620.91060.93950.93520.9337
150.96510.94590.95770.97460.97020.9688
200.97960.96860.97790.9870.98240.9829
4100.95320.93430.94870.95610.95250.9533
150.98540.9540.97960.98530.9680.979
200.990.9850.98920.98990.9890.988
6100.81410.73690.78330.91510.88770.9125
150.94050.81230.86540.94960.93280.9619
200.96330.90280.92690.97520.96580.9723
8100.92710.90290.91550.95710.94330.9474
150.96550.93140.96110.9790.95790.9798
200.97640.96250.97450.98760.97990.9856

Share and Cite

MDPI and ACS Style

Jia, H.; Lang, C.; Oliva, D.; Song, W.; Peng, X. Dynamic Harris Hawks Optimization with Mutation Mechanism for Satellite Image Segmentation. Remote Sens. 2019, 11, 1421. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11121421

AMA Style

Jia H, Lang C, Oliva D, Song W, Peng X. Dynamic Harris Hawks Optimization with Mutation Mechanism for Satellite Image Segmentation. Remote Sensing. 2019; 11(12):1421. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11121421

Chicago/Turabian Style

Jia, Heming, Chunbo Lang, Diego Oliva, Wenlong Song, and Xiaoxu Peng. 2019. "Dynamic Harris Hawks Optimization with Mutation Mechanism for Satellite Image Segmentation" Remote Sensing 11, no. 12: 1421. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11121421

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop