Next Article in Journal
Preparation and Characterization of Nano-Silver-Loaded Antibacterial Membrane via Coaxial Electrospinning
Next Article in Special Issue
Deep Convolutional Neural Network with Symbiotic Organism Search-Based Human Activity Recognition for Cognitive Health Assessment
Previous Article in Journal
RiPa-Net: Recognition of Rice Paddy Diseases with Duo-Layers of CNNs Fostered by Feature Transformation and Selection
Previous Article in Special Issue
SDE-YOLO: A Novel Method for Blood Cell Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Otsu Multi-Threshold Image Segmentation Based on Adaptive Double-Mutation Differential Evolution

1
Shandong Research Institute of Industrial Technology, Jinan 250100, China
2
School of Computer Science and Technology, Shandong University of Finance and Economics, Jinan 250014, China
*
Author to whom correspondence should be addressed.
Biomimetics 2023, 8(5), 418; https://doi.org/10.3390/biomimetics8050418
Submission received: 3 August 2023 / Revised: 5 September 2023 / Accepted: 6 September 2023 / Published: 8 September 2023
(This article belongs to the Special Issue Biomimetic and Bioinspired Computer Vision and Image Processing)

Abstract

:
A quick and effective way of segmenting images is the Otsu threshold method. However, the complexity of time grows exponentially as the number of thresolds rises. The aim of this study is to address the issues with the standard threshold image segmentation method’s low segmentation effect and high time complexity. The two mutations differential evolution based on adaptive control parameters is presented, and the twofold mutation approach and adaptive control parameter search mechanism are used. Superior double-mutation differential evolution views Otsu threshold picture segmentation as an optimization issue, uses the maximum interclass variance technique as the objective function, determines the ideal threshold, and then implements multi-threshold image segmentation. The experimental findings demonstrate the robustness of the enhanced double-mutation differential evolution with adaptive control parameters. Compared to other benchmark algorithms, our algorithm excels in both image segmentation accuracy and time complexity, offering superior performance.

1. Introduction

The technique of images segmentation involves breaking an image up into a number of distinct, non-overlapping parts and extracting the desired, human-interest-centered regions. From image processing to image analysis, it is a committed step that makes it easier for following computer vision, etc. There are other widely used image segmentation methods in addition to the Otsu threshold image segmentation method [1,2,3] in the field. The maximum entropy method, for instance [4]. A well-known edge detection strategy for extracting edge information from images is the Prewitt image segmentation technique [5]. This method identifies the edges in the image by detecting gradient changes and is based on gradient calculation of the grayscale values of image pixels. The watershed image segmentation method [6] is a region-based image segmentation method. In this method, the image is considered as a terrain map, where the height represents the grayscale value of the image. Image segmentation is achieved by finding watersheds, which divide the image into different regions or sets of different regions. The advantage of this method is that it can preserve areas with clear edges and rich textures in the image. However, there is also a drawback of over segmentation for images with large uniform areas. The Canny image segmentation method [7] is a classic edge detection algorithm that can accurately detect edges in an image. This algorithm determines image edges by detecting gradient changes, and improves the accuracy and continuity of edges through non maximum suppression, dual threshold processing, and edge connectivity. The Sobel image segmentation method [8] is a commonly used edge detection algorithm that is used to detect edges in images. The Sobel operator is used to calculate the gradient value of image pixels to determine the position and direction of the image. The Robert image segmentation method [9] is a classic edge detection algorithm that applies the Robert operator to convolution the image and calculates the gradient values between adjacent pixels to determine the position and direction of the edges. An image segmentation technique based on the Laplacian Gaussian filter is the log image segmentation approach [10]. By applying Gaussian and Laplacian filtering procedures to the image, this technique finds edges and other features in the image. The threshold image segmentation approach is a common one in image processing technology because of its benefits of simplicity, effectiveness, and strong robustness. Single-threshold image segmentation and multi-threshold image segmentation are two common threshold segmentation techniques. The threshold is determined using the Otsu threshold segmentation technique. The time complexity grows exponentially as the number of thresholds increases, as does the quantity of calculation.
One might think of the threshold image segmentation approach as an optimization issue. Researchers are increasingly merging threshold image segmentation with intelligent optimization algorithms as swarm intelligence optimization techniques grow over time. For instance, based on the enhanced Firefly algorithm [11], better particle swarm optimization algorithm [12], improved genetic algorithm [13], improved whale optimization algorithm [14], and improved cuckoo algorithm [15]. You can use these clever optimization methods to identify appropriate optimization objective functions. The Otsu method and the maximum entropy method are common objective functions used in image segmentation.
Differential evolution [16] (DE) was proposed by Storn and Price in 1997. Because of its small control parameters and strong robustness, differential evolution has been widely concerned by researchers. The basic operations of differential evolution include initial population, mutation, crossover and selection. As long as the population iteration conditions are met, the population will iterate continuously to find the optimal solution. Therefore, differential evolution is also a parallel search algorithm. In recent years, with the in-depth innovation of research, more and more researchers found that control parameters and mutation strategies are the main factors that affect the performance of differential evolution. The key control parameters in differential evolution mainly include population size NP, mutation operator F and crossover factor CR. Based on this, researchers continue to study and find variants of differential evolution. In reference [17], AC Sanderson et al. proposed an improved differential evolution (JADE) using a new mutation strategy of selecting mutation vectors from an external archived population, and the control parameters F and CR in the population obey Cauchy distribution and Normal distribution, respectively. Ryoji Tanabe et al. [18] developed an improved differential evolution algorithm (SHADE) based on JADE, which guides the selection of control parameters in subsequent iterations through the number of successful individuals with historical mutations. The experimental results showed that SHADE outperforms previous state-of-the-art DE algorithms on a large number of benchmark problems. Wu Deng et al. [19] proposed a differential evolution with a new hybrid mutation factor and adaptive control parameters. Experimental results show that compared with other DE variants, the algorithm has better performance on the above test functions.
The differential evolution algorithm (DE) is a population intelligent optimization algorithm that searches for the optimal solution by simulating the process of natural evolution. With the continuous deepening of research on differential evolution algorithms, their application fields are gradually expanding. More and more optimization problems are combined with differential evolution algorithms to solve the optimal solution, such as scheduling problems [20], traveling salesman problems [21], and image processing [22]. In scheduling problems, differential evolution algorithms optimize scheduling strategies and resource allocation to maximize work efficiency or minimize work time. Differential evolution combined with scheduling problems can handle complex constraint problems and quickly converge to find the optimal solution. The traveling salesman problem is a classic combinatorial optimization problem that requires finding the shortest path so that the traveling salesman can visit several cities and ultimately return to the starting city. The differential evolution algorithm searches for the optimal solution by adjusting the arrangement order of cities, with the goal of minimizing travel and finding the optimal solution for the shortest path. In this problem, differential evolution algorithm can efficiently search for the global optimal solution. In image processing problems, differential evolution algorithms are used for image registration, image segmentation, feature extraction, and image enhancement.
Bird swarm behavior is simulated via particle swarm optimization [23]. In the process of segmenting a picture, each particle represents a potential segmentation scheme, and the best segmentation outcome is sought by continual iteration. With the help of pheromone propagation and reinforcement techniques, the Ant Colony algorithm [24] eventually achieves the best segmentation outcomes by simulating the behavior of ants as they look for food. The genetic algorithm [25] mimics the evolution process in nature by repeatedly iterating the population through selection, crossover, and mutation operations, maintaining its diversity and exploratory capacity through mutation, and finally discovering the best answer. Due to its few control parameters, robustness, and difficulty finding optimal solutions, the differential evolution method [26] performs well and offers advantages in picture segmentation.
The differential evolution variant tested against pertinent benchmark functions performed really well. Researchers are increasingly using differential evolution to enhance and optimize image processing outcomes. To determine the ideal threshold and enhance the performance of image segmentation, Sushi L. Kumar et al. integrated differential evolution with the Otsu threshold segmentation method in reference [26]. The experimental results demonstrate that the application of differential evolution to image segmentation can significantly increase the segmentation quality when compared to the segmentation results obtained using the conventional Otsu threshold segmentation method. However, the overall performance of differential evolution depends on the ability to balance global search and local search. The reason for its greater impact is the selection of control parameters. However, the above differential evolution proposed by Sushi L Kumar et al. is applied to threshold image segmentation, with high time complexity and poor robustness. Unable to meet the processing objectives of multi-threshold image segmentation. In reference [27], Helen Vicente Humann Ayala et al. proposed an improved differential evolution (BDE) based on the generation of beta distribution adaptive control parameters F and CR values. In each iteration, the adaptive control parameters generate random F and CR values, enhancing the random search ability of the algorithm and improving the performance of the algorithm. Experimental results show that the improved differential evolution (BDE) is more effective than the FODPSO algorithm. However, due to the fact that the BDE algorithm randomly generates the values of control parameters F and CR, it cannot be ensured that the F and CR values are suitable for the entire iterative process of the algorithm. Based on the above algorithm. In this paper, a differential evolution based on adaptive double mutation is proposed for image segmentation. The population is divided according to the fitness value, and the double-mutation strategy is used to maximally balance the global search ability and local search ability of the algorithm. In addition, the values of control parameters F and CR are adaptively adjusted according to different stages of population evolution. The improved differential evolution in this paper is applied to Otsu multi-threshold image segmentation, and the maximum inter class variance method is regarded as the objective function to find the optimal threshold. The experimental results show that the algorithm is more accurate than other benchmark algorithms in image segmentation, with low time complexity and strong robustness.
Here, is the rest of the essay: The related work, such as the traditional differential evolution and Thresholding picture segmentation, are introduced in the second section. The third section introduces the TRDE algorithm proposed in this article. The fourth section introduces the experimental results, the fifth section is the analysis and research discussion of the experimental results, the sixth section is the references and future work, and the third part describes the improved double-mutation adaptive differential evolution algorithm (TRDE).

2. Related Work

2.1. The Classical DE Algorithm

Differential evolution is a population intelligent optimization algorithm. DE algorithm includes the following four steps: population initialization, mutation operation, crossover operation and selection operation. The following is the basic process of the DE algorithm:
(1)
Population initialization
Randomly generate the initial population, initialize and generate NP individual vectors in d-dimensional space, denoted as:
x i , g = ( x 1 , g , x 2 , g , , x d , g ) i = 1 , 2 , , N P
Among them, NP is the population size, and g is the current evolutionary algebra. The NP individual vectors generated in d-dimensional space are randomly generated, and the specific formula is as follows:
x i , j = x i , j L + r a n d ( 0 , 1 ) × ( x i , j U x i , j L )
In the above formula,  x i , j  represents the ith randomly generated individual,  x i , j U , x i , j L  represents the Upper and lower bounds of the ith individual vector,  r a n d ( 0 , 1 )  represents a random number generated within the range of (0,1).
(2)
Mutation
The DE algorithm achieves individual vector mutation through differential strategy, which is also the core idea of the DE algorithm. Randomly select several individual vectors from the population for differential operation to generate differential vectors. Common mutation strategies are as follows:
  • DE/rand/1:
    v i , g = x r 1 , g + F × ( x r 2 , g x r 3 , g )
  • DE/rand/2:
    v i , g = x r 1 , g + F × ( x r 2 , g x r 3 , g ) + F × ( x r 4 , g x r 5 , g )
  • DE/best/1:
    v i , g = x b e s t , g + F × ( x r 1 , g x r 2 , g )
  • DE/best/2:
    v i , g = x b e s t , g + F × ( x r 1 , g x r 2 , g ) + F × ( x r 3 , g x r 4 , g )
    wherein  x r 1 , g , x r 2 , g , x r 3 , g , x r 4 , g , x r 5 , g ( r 1 r 2 r 3 r 4 r 5 )  represent five randomly selected individual vectors from the population that are different from each other,  v i , g  represents the variation vector of the ith individual in the generated generation, and F is the scaling factor.
(3)
Crossover
Crossover operation is to cross the generated mutated individual vector with the parent individual vector to generate an experimental individual vector. The commonly used binomial cross operation in DE algorithm is as follows:
u i , j , g = { v i , j , g i f r a n d C R o r j = j r a n d x i , j , g o t h e r w i s e
In the above formula, CR is the crossover operator,  j r a n d  is a random integer selected from {1, 2, …, d} to ensure that the mutated individual has at least one individual component inherited to the next generation, avoiding the same vector as the parent individual.
(4)
Selection
The selection operation in the DE algorithm adopts a greedy strategy, which selects the individual vector with the best fitness value as the evolutionary offspring vector through comparison. The specific formula is as follows:
x i , g + 1 = { u i , g , i f f ( u i , g ) f ( x i , g ) x i , g , o t h e r w i s e

2.2. The Otsu Threshold Segmentation Method

The maximum inter class variance method [28] was proposed by Japanese scholar Otsu in 1997. The basic principle of Otsu for image threshold segmentation is to divide the image into background and target categories based on the grayscale characteristics of the image, calculate the number of pixels in the image, and solve for the maximum variance between the image target and background to achieve image threshold segmentation.
For a pair digital image of M × N, where i is the grayscale value, L is the grayscale level, and the range of grayscale values is [0, L − 1].  n i  represents the number of pixels with grayscale level i, and the probability of grayscale level i appearing in the image is  p i , then:
P i = n i N , I = 0 , 1 , 2 , , L 1
i = 0 L 1 P i = 1
Based on the grayscale characteristics of digital images, the pixel points in the image are divided into two categories: target and background using threshold s, Represented by  H 0  and  H 1 . The grayscale values between [0, s] are classified as  H 0 , and the remaining grayscale values between [s + 1, L − 1] are classified as  H 1 . For the entire image, the average grayscale value is:
Q s = i = 0 L 1 i P i
The mean of  H 0  and  H 1  is:
q 0 = i = 1 s i p i t 0
q 1 = i = s + 1 L 1 i p i t 1
In the above formula:
t 0 = i = 0 s p i
t 1 = i = s + 1 L 1 p i = 1 t 0
Based on the above formula, we can obtain:
q s = t 0 q 0 + t 1 q 1
The definition of inter class variance is as follows:
σ B 2 = t 0 ( q 0 q s ) 2 + t 1 ( q 1 q s ) 2 = t 0 q 0 2 + t 1 q 1 2 q s 2 = t 0 q 0 2 + t 1 q 1 2 ( t 0 q 0 + t 1 q 1 ) 2 = t 0 t 1 ( q 0 q 1 ) 2
Through the above formula, we can see that the s value corresponding to the maximum inter class variance is the optimal segmentation threshold.

2.3. Existing Problems in Research and Motivation for Improvement

The segmentation criterion for traditional Otsu is the largest variance between the target and background of the picture, which divides the target and backdrop of the image into two categories. However, the histogram of the segmented Grayscale typically does not exhibit a bimodal form, making it challenging to obtain precise image segmentation findings. To obtain the ideal threshold, many researchers treat threshold picture segmentation as an optimization problem and pair it with clever optimization methods. The control settings and mutation approach in the traditional DE algorithm have the most effects on the algorithm’s overall performance. In the DE algorithm, the important control parameters are the scaling factor F and the crossover factor CR. The scaling factor F controls the magnitude of differential evolution vector, which affects the convergence ability and Rate of convergence of the algorithm. The crossover factor CR affects the probability of generating a vector of experimental individuals in crossover operations to inherit genes from their parents or mutated individuals. The second part is the improvement of mutation strategy, balancing the global search ability and local search ability of the algorithm. To sum up, this paper proposes an adaptive double-mutation differential evolution algorithm (TRDE), which divides the population according to the fitness value of the population, adopts corresponding mutation strategies for different sub-populations, and adaptively adjusts the values of control parameters F and CR to adapt to different stages of population evolution.

3. The TRDE Algorithm

Researchers have enhanced and created variants of the DE algorithm based on it in order to significantly balance the global and local search capabilities. There are primarily two components in the DE algorithm that greatly influence how well it performs. The first step is to adjust the values of the parameters F and CR. The scaling factor F determines the size of the differential evolution vector and has an impact on the algorithm’s capacity for convergence and rate of convergence. The possibility of creating a vector of experimental individuals in crossover operations who would inherit genes from their parents or modified individuals is influenced by the crossover factor CR. a variety of adaptive parameter control DE algorithm versions, such as in reference [29]. A differential evolution algorithm based on dual mutation strategy is proposed for image segmentation. The experimental results show that this algorithm is competitive with other improved DE algorithms in terms of optimization accuracy and convergence speed. Vladimir Stanovov et al. proposed a method in reference [30] to adjust the values of control parameters F and CR based on fitness values, and increase the biased Lemmer mean and LBR to improve its performance. The experimental results show that the algorithm performance has been greatly improved. The second part is the improvement of mutation strategy, balancing the global search ability and local search ability of the algorithm. Regarding the improvement of mutation strategy, there are also many variants of DE algorithm. For example, in reference [31], Wan-li Xiang et al. proposed an enhanced DE algorithm (EDE), which uses a reverse learning strategy to increase the diversity of the initial population, improves the mutation strategy on the basis of JADE algorithm, and improves the overall performance of the algorithm. As mentioned in reference [32], Ali Wagdy Mohamed proposed a mutation strategy based on triangular mutation rules, which greatly balances the global detection ability and local search ability of the population.

3.1. Dual Mutation Strategy

The main component of the DE algorithm is the mutation strategy, which can make the population’s individual evolution vectors and orientations display diversity. Based on the fitness scores of various population vectors, this article separates the population into two subpopulations. Mutation strategy 1 is used for populations with high fitness values, while mutation strategy 2 is used for populations with low fitness values. In general, it is simpler to locate the ideal solution close to people with high fitness values or those with low principle fitness values. Consequently, the dual mutation technique can balance population diversity and local search capability, speeding convergence.
After initializing the population, calculate the population fitness value, and divide the population into two subpopulations based on the fitness value. Subpopulation 1 stores the vector of individuals with high fitness values, while Subpopulation 2 stores the vector of individuals with poor fitness values. Improved mutation strategy 1:
v i , g = x i , g + F i × ( x d b e s t , g x r 1 , g ) + F i × ( x r 2 , g x r 3 , g )
x d b e s t , g  is a random selection of one of the top d% individual vectors from Subpopulation1 as the basis vector. In this article, the selection of the d parameter adopts a linearly decreasing function. The specific changes in parameter d are as follows:
d i = d min d max ( 1 g e n g e n e r t a t i o n s ) + ( d max d min d max + d min )
This article set  d max = 0.8 , d min = 0.4 , gen is the current population evolution algebra, and generations is the total evolution algebra. The variation of d parameter shows a decreasing trend with the increase of iteration times. As the number of iterations in the population continues to increase, the fitness values of individual vectors in the population continuously concentrate in a certain local region. Therefore, in order to improve the local convergence ability of the population and find the optimal solution, in the later stage of iteration, the selection of the optimal d% individual vector gradually decreases, accelerating population convergence and enhancing convergence ability. Generally speaking, individuals with high fitness values or poor principle fitness values are more likely to find the optimal solution, while individuals with poor fitness values can enhance population diversity. Therefore, the mutation strategy 2 proposed in this article is as follows:
v i , g = x i , g + F i × ( x r 1 , g x r 3 , g ) + F i × ( x r 2 , g x r 4 , g )
Based on the fitness value, randomly select four distinct integers  r 1 , r 2 , r 3 , r 4  in Subpopulation2, and generate a difference vector for mutation operation to improve population diversity.

3.2. Parameter Adaptation

The fixed values of the control parameters F and CR in the traditional DE algorithm make it impossible for them to adapt to the algorithm’s needs at various phases of evolution. As a result, adaptive parameter control techniques are being used by more and more DE variations. In general, a DE algorithm should have strong global search capabilities, retain population variety as much as feasible in the early stages of population evolution, and look for the world’s best answer. In the later stage of population evolution, as the number of iterations keeps rising, it becomes increasingly important to have strong local convergence capacity and accelerate convergence.
The scaling factor F in the DE algorithm is closely related to the Rate of convergence, which controls the amplitude of the difference vector. In the early stages of population evolution, a larger value of the scaling factor F is beneficial for maintaining population diversity and ensuring that the algorithm searches throughout the entire solution space. As the number of iterations continues to increase, in the later stage of evolution, the optimal solution of the individual vector gradually concentrates in a certain local region. It is necessary to gradually reduce the value of F and strengthen the local convergence ability of the algorithm. The crossover factor CR affects the probability of experimental individuals generated by crossover operations inheriting genes from mutated individuals or parent individuals. Generally speaking, in the early stages of algorithm evolution, to ensure the algorithm’s global search ability, the CR value is small. In the later stage of algorithm evolution, in order to improve the local search ability of the algorithm with a larger CR value, combined with changes in the fitness value of individual vectors in the population, this paper proposes adaptive parameter control as follows:
F i = F min ln ( F max F min ) × ( g e n g e n e r a t i o n s ) 2
where gen represents the evolutionary algebra of the current population, and generations represents the total number of iterations,  F max  represents the maximum value for setting the scaling factor,  F min  represents the minimum value for setting the scaling factor, This article set  F max = 0.7 , F min = 0.2 .
The changes in the scaling factor F proposed above show a decreasing trend, which is consistent with the convergence characteristics of differential evolution.
C R i = ( 1 C R min C R max ) × ( g e n g e n e r a t i o n s ) 2 + ln C R max C R min × ( C R min C R max C R min )
where gen represents the evolutionary algebra of the current population, and generations represents the total number of iterations,  C R max  represents the maximum value for setting the crossover factor, and  C R min  represents the minimum value for setting the crossover factor. In this article, set  C R max = 0.9 , C R min = 0.1 . The change in CR in the above formula shows a monotonic increasing trend, which greatly balances the contradiction between population diversity and Rate of convergence. From the test results, it can be seen that the improved adaptive control parameter method has advantages over other benchmark algorithms.

3.3. Pseudo Code of TRDE

1:Initialize population p individuals and calculate their fitness values, NP = 30, gen = 1, generations = 100;
2:for gen = 1 to generations do
3: for i = 1 to NP do
4:   F = F min ln ( F max F min ) × ( g e n g e n e r a t i o n s ) 2
5:   C R = ( 1 C R min C R max ) × ( g e n g e n e r a t i o n s ) 2 + ln C R max C R min × ( C R min C R max C R min )
6:  Implement mutation in Equations (18) and (20);
7:  Implement crossover in Equation (7);
8:   Implement selection in Equation (8);
9:   end for
10:   gen = gen + 1;
11: end for

4. Experiment

4.1. Benchmarking Datasets and Benchmarking Algorithms

Better adaptive double mutation; in this work, the segmentation performance of the differential evolution (TRDE) method is evaluated using five thresholds between [2,6] and eight test images from the Berkeley database, along with the other five benchmark algorithms. The JADE method with archive operations, the SHADE algorithm, which is an improvement of JADE, the BDE algorithm, the HSDE algorithm with adaptive parameter control [33], and the SAF-DE algorithm [34] are some of the benchmark algorithms utilized in this article.
The population size NP in this article is universally fixed to 30 and the maximum number of iterations is set to 100 in order to guarantee the validity of the experimental results. The TRDE algorithm and five other benchmark algorithms are each run 30 times for the chosen 8 test images. Finally, the experimental impact of picture segmentation is assessed using evaluation indicators. Figure 1 illustrates the picture segmentation outcomes of the TRDE method for 8 photos and 5 thresholds within the range of [2,6]. Additionally, three comparing methods (BDE, JADE, and SAF-DE) were chosen at random for this study in order to display the segmented picture findings in Figure 2, Figure A1 and Figure A2.
This article not only compares the improved TRDE algorithm with five other benchmark algorithms applied to Otsu threshold segmentation, but also compares the TRDE algorithm with other image segmentation methods (Prewitt, Sobel, Canny, Robert) for image segmentation, and uses evaluation indicators for quantitative statistics.

4.2. Common Image Segmentation Quality Evaluation Standards

In this study, we used benchmark techniques for image segmentation as well as the improved adaptive double-mutation differential evolution. According to references [35,36], we discovered that the structural similarity (SSIM) and the peak signal-to-noise ratio (PSNR), as well as the standard deviation of fitness value and the running time of the algorithm, are the most frequently used image segmentation quality evaluation criteria.
(1)
Structural Similarity (SSIM)
Digital images have certain similar features in space, and each pixel has correlation between them. Namely, structural similarity, by comparing structural information, further examining the distortion of the image. SSIM is a measure of image distortion based on human visual characteristics, defined as follows:
S S I M ( x , y ) = ( 2 μ x μ y + C 1 ) ( 2 σ x , y + C 2 ) ( μ x 2 + μ y 2 + C 1 ) ( σ x 2 + σ y 2 + C 2 )
In the above formula,  μ x  and  μ y  represent the mean of x and y, respectively, while  σ x  and  σ y  represent the variance of x and y.  C 1  and  C 2  are two positive constants. The value range of SSIM is 0 to 1. In practical operation, it is necessary to first unify the image and also perform grayscale processing on the image. Usually, the larger the SSIM value, the higher the structural similarity between two images. The SSIM value of the same two images is equal to 1.
(2)
Peak signal-to-noise ratio (PSNR)
The peak signal-to-noise ratio (PSNR) is an objective standard for evaluating images, typically used to measure the quality of compressed reconstructed images, measured in decibels (DB). Due to the fact that PSNR is based on error sensitive image quality evaluation and does not take into account the visual characteristics of the human eye, there may be inconsistencies between measurement results and subjective perception. The definition of PSNR is as follows:
P S N R = 10 log 10 255 2 M S E
MSE refers to the mean squared error of the original image and the segmented image. The definition of MSE is as follows
M S E = i = 1 M j = 1 N ( x ( i , j ) y ( i , j ) ) 2 M × N
In the above formula, x(i, j) and y(i, j) represent the original image and the segmented image, respectively, while MSE represents the average of the sum of squares of the pixel values of the original image and the segmented image. Generally speaking, the higher the PSNR value, the higher the similarity between the original image and the segmented image, and the higher the segmentation quality of the resulting image.
The population size and number of iterations for the improved adaptive double-mutation differential evolution (TRDE) and other benchmark algorithms (JADE, SHADE, BDE, HSDE) in this study are set to the identical values. Run to get the segmented image and ideal threshold for each test image. To assess its segmentation effect, we employ the widely used evaluation criteria PSNR and SSIM for picture segmentation quality. Standard deviation and the algorithm’s average execution time should be used to assess the algorithm’s performance.
Figure 2 shows the segmentation results obtained by the TRDE image segmentation algorithm and four other image segmentation methods. It can be seen that the TRDE algorithm preserves as much structural information of the original image as possible by selecting the optimal threshold, resulting in better image segmentation results. The images segmented by other edge detection algorithms may have noise or rough edges, and the segmentation effect is shown in Figure 2.

5. Discussion and Analysis of Experimental Results

Result Analysis

This study selected test images from the Berkeley segmentation dataset with each image size of 200,200, applied the TRDE algorithm and other benchmark algorithms to image segmentation, and evaluated the segmented images using evaluation indicators PSNR and SSIM. This study then conducted testing and analyzed experimental results. Table 1 displays the obtained PSNR data, while Table 2 displays the SSIM results.
These two evaluation indicators need to be considered together to evaluate the segmentation effect of an image. The larger the PSNR value and SSIM value, the better the segmentation effect of the image. The smaller the value, the worse the segmentation effect. Table 1 shows the peak signal-to-noise ratio (PSNR) calculated by the TRDE, JADE, SHADE, BDE, and HSDE algorithms for 8 test images at different thresholds. From Table 1, we can see that the segmentation performance of the TRDE algorithm is overall superior to other benchmark algorithms, with the PSNR values of the vast majority of segmented images ranking first. For Baboon images, when the thresholds are 3, 4, and 6, the PSNR values of image segmentation rank first. For camera images, only when the threshold is 5, the PSNR values of image segmentation rank second, and the rest of the PSNR values rank first. For pepper images, when the thresholds are 2, 3, and 4, the PSNR values of image segmentation rank first. For Saturn images, when the thresholds are 2, 3, 5, and 6, the PSNR values of image segmentation rank first. However, there are very few cases where the quality of image segmentation is poor. For example, when the pepper image threshold is 5, 6, the PSNR value of image segmentation ranks third and fourth. This is because the essence of differential evolution is a random algorithm, which may lead to poor image segmentation results.
Table 2 compares the segmented picture’s Structural Similarity Index (SSIM) to the effects of other benchmark algorithms JADE, SHADE, BDE, and HSDE on image segmentation. All SSIM scores for the Barbara image rank top for the 8 test images that were chosen. At thresholds of 3, 4, and 6, the TRDE algorithm performs the best image segmentation for the Lena image. According to the experimental findings, the TRDE algorithm successfully segments the great majority of images.
In comparison to other picture segmentation techniques, Table 3 displays the TRDE algorithm’s findings for the Structural Similarity Index (SSIM). The results of the TRDE algorithm applied to picture segmentation are shown in Table 4 along with comparisons to other image segmentation techniques. In contrast to standard SSIM, MS-SSIM evaluates the final assessment result by comparing the similarity of images at several scales and weighting the similarity at various scales. This paper presents an improved differential evolution technique for threshold image segmentation, which chooses an optimal threshold to preserve as much of the structural information of the original image as possible after segmentation. Due to the possible presence of noise or rough edges in the images segmented by edge detection algorithms, the structural similarity between the segmentation results and the original image may decrease. Therefore, the experimental results of the SSIM and MS-SSIM indicators for the segmented images of the TRDE algorithm in Table 3 and Table 4 are higher than those of the other four edge detection methods (Prewitt, Sobel, Canny, Robert).
The similarity index (UQI) of the TRDE algorithm used for picture segmentation is displayed in Table 5 along with comparisons to other image segmentation techniques. The UQI score ranges from 0 to 1, and a greater number denotes a higher degree of similarity and picture quality between the two photos. The UQI indicator thoroughly takes into account the mean, variance, and covariance of each image in order to assess how similar and high-quality each image is. Table 5 displays the experimental outcomes of the UQI picture quality evaluation index. The table shows that the TRDE algorithm consistently performs well when compared to other picture segmentation techniques. The Canny algorithm, for instance, places first and the TRDE method, second in photos such as Barbara and Lena. In cameraman images, the TRDE algorithm ranks first. Overall, TRDE has strong stability when applied to image segmentation algorithms, resulting in good segmentation image quality.
Taking into account the experimental results in Table 3, Table 4 and Table 5, the quantitative results show that the improved differential evolution algorithm (TRDE) applied to image segmentation has better image segmentation quality, strong algorithm robustness, and fast convergence speed compared to other image segmentation methods.
The enhanced adaptive double-mutation differential evolution (TRDE) and the other five benchmark algorithms (JADE, SHADE, BDE, SAF-DE, HSDE) are shown in Table 6 along with their average running times. The TRDE algorithm’s average running time is generally faster than other algorithms’ for photos of peppers and women. When the threshold is 2 or 3, the TRDE algorithm’s average running time for plane images is marginally greater than that of other techniques. Overall experimental findings demonstrate that this method’s rate of convergence is quicker than that of the other five benchmark algorithms (JADE, SHADE, BDE, SAF-DE, HSDE).
Table 1, Table 2, Table 6 and Table 7 illustrate the overall effectiveness of the enhanced adaptive double-mutation differential evolution (TRDE) in this study. The usual segmentation quality evaluation standards PSNR and SSIM are used to assess the segmented image. The improved adaptive double-mutation differential evolution has greater performance than previous benchmark algorithms, as can be observed from the experimental data listed in the table. Table 1 and Table 2 demonstrate that, when used for image segmentation, the method has good segmentation capabilities. Table 6 demonstrates that the enhanced differential evolution (TRDE) is more effective and has a faster Rate of convergence than previous algorithms. Table 7 demonstrates the TRDE algorithm’s strong resistance to the Otsu threshold image segmentation fitness function.

6. Conclusions

The adaptive double-mutation differential evolution method (TRDE) for picture segmentation is proposed in this research and employs Otsu as the goal function. In addition, compare this approach to five other benchmark algorithms (JADE, SHADE, BDE, SAF-DE, and HSDE), compute PSNR and SSIM values for the segmented images, determine the average algorithm running time, and thoroughly evaluate the technique’s performance. According to the experimental findings, the TRDE algorithm performs the best overall inside the [2,6] threshold range. Finally, the robustness of the method to the fitness function of Otsu threshold image segmentation was tested using the standard deviation of fitness values.
In future work, we will carry out studies from two aspects, improving different mutation strategies and optimizing the impact of adaptive control parameters on threshold selection for differential evolution. On the other hand, different threshold segmentation techniques are used to determine the optimal threshold.

Author Contributions

Y.G. and Y.W.; methodology, K.M.; software, Z.Z.; validation. writing—original draft preparation, Y.G.; writing—review and editing, Y.G. and Y.W.; visualization, K.M.; su-pervision, Z.Z.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data openly available in a public repository. The data that support the findings of this study are openly available in Berkeley Split Dataset at https://www2.eecs.berkeley.edu/Research/Projects/CS/vision/bsds/ (accessed on 20 June 2023).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. (a) Represents the original image, (b) segmented image with two thresholds, (c) segmented image with three thresholds, (d) segmented image with four thresholds, (e) segmented image with 5 thresholds, (f) segmented image with six thresholds (BDE).
Figure A1. (a) Represents the original image, (b) segmented image with two thresholds, (c) segmented image with three thresholds, (d) segmented image with four thresholds, (e) segmented image with 5 thresholds, (f) segmented image with six thresholds (BDE).
Biomimetics 08 00418 g0a1
Figure A2. (a) Represents the original image, (b) segmented image with two thresholds, (c) segmented image with three thresholds, (d) segmented image with four thresholds, (e) segmented image with 5 thresholds, (f) segmented image with six thresholds (JADE).
Figure A2. (a) Represents the original image, (b) segmented image with two thresholds, (c) segmented image with three thresholds, (d) segmented image with four thresholds, (e) segmented image with 5 thresholds, (f) segmented image with six thresholds (JADE).
Biomimetics 08 00418 g0a2
Figure A3. (a) Represents the original image, (b) segmented image with two thresholds, (c) segmented image with three thresholds, (d) segmented image with four thresholds, (e) segmented image with five thresholds, (f) segmented image with six thresholds (SAF-DE).
Figure A3. (a) Represents the original image, (b) segmented image with two thresholds, (c) segmented image with three thresholds, (d) segmented image with four thresholds, (e) segmented image with five thresholds, (f) segmented image with six thresholds (SAF-DE).
Biomimetics 08 00418 g0a3

References

  1. Abualigah, L.; Almotairi, K.H.; Elaziz, M.A. Multilevel thresholding image segmentation using meta-heuristic optimization algorithms: Comparative analysis, open challenges and new trends. Appl. Intell. 2023, 53, 11654–11704. [Google Scholar] [CrossRef]
  2. Su, H.; Zhao, D.; Elmannai, H.; Heidari, A.A.; Bourouis, S.; Wu, Z.; Cai, Z.; Gui, W.; Chen, M. Multilevel threshold image segmentation for COVID-19 chest radiography: A framework using horizontal and vertical multiverse optimization. Comput. Biol. Med. 2022, 146, 105618. [Google Scholar] [CrossRef] [PubMed]
  3. Han, Y.; Chen, W.; Heidari, A.A.; Chen, H. Multi-verse optimizer with rosenbrock and diffusion mechanisms for multilevel threshold image segmentation from COVID-19 chest X-ray images. J. Bionic Eng. 2023, 20, 1198–1262. [Google Scholar] [CrossRef] [PubMed]
  4. Mahajan, S.; Mittal, N.; Salgotra, R.; Masud, M.; Alhumyani, H.A.; Pandit, A.K. An efficient adaptive salp swarm algorithm using type II fuzzy entropy for multilevel thresholding image segmentation. Comput. Math. Methods Med. 2022, 2022, 2794326. [Google Scholar] [CrossRef]
  5. Prewitt, J.M. Edges as ridge lines in images. Proc. Natl. Acad. Sci. USA 1979, 76, 698–702. [Google Scholar]
  6. Roerdink, J.B.; Meijster, A. The watershed transform: Definitions, algorithms, and parallelization strategies. Fundam. Informa. 2000, 41, 187–228. [Google Scholar] [CrossRef]
  7. Mousavi, M.R.; Kakooee, O.R. Optimization of the Canny Edge Detection Algorithm Using Genetic Algorithm for Edge Detection of Biological Images. Int. J. Comput. Sci. 2015, 12, 186–194. [Google Scholar]
  8. Zhou, P.; Deng, K.; Liu, Q.; Chen, Y.; Zhao, Q. A novel adaptive efficient image segmentation algorithm based on the Sobel operator. PLoS ONE 2019, 14, e0212679. [Google Scholar]
  9. Napieralski, A.; Napieralska, M. Comparison of Roberts and Prewitt operators applied to the detection of stress patterns. Comput. Intell. 2014, 30, 199–211. [Google Scholar]
  10. Liao, S.H.; Chen, C.H. Efficient fish image segmentation on complex background using adaptive Log-Transformed CIELAB Color Space and Dual-Threshold Strategy. Neurocomputing 2015, 153, 84–94. [Google Scholar]
  11. Krishnan, V.G.; Rao, P.V.; Raja, S.; Bhargav, S.; Balaji, S.; Divya, V. An improved firefly algorithm for task scheduling in cloud computing environment. IEEE Access 2022, 10, 2363–2372. [Google Scholar]
  12. Xu, L.; Song, B.; Cao, M. An improved particle swarm optimization algorithm with adaptive weighted delay velocity. Syst. Sci. Control Eng. Open Access J. 2021, 9, 188–197. [Google Scholar] [CrossRef]
  13. Liang, H.; Zou, J.; Zuo, K.; Khan, M.J. An improved genetic algorithm optimization fuzzy controller applied to the wellhead back pressure control system. Mech. Syst. Signal Process. 2020, 142, 106708. [Google Scholar] [CrossRef]
  14. Kashani, A.R.; Camp, C.V.; Armanfar, M.; Slowik, A. Whale Optimization Algorithm; Springer: Berlin/Heidelberg, Germany, 2020. [Google Scholar]
  15. Yang, X.-S.; Deb, S. Cuckoo search via Lévy flights. In Nature-Inspired Algorithms for Optimisation; Springer: Berlin/Heidelberg, Germany, 2022; pp. 123–144. [Google Scholar]
  16. Storn, R.; Price, K. Differential evolution—A simple and efficient heuristic for global optimization over continuous spaces. J. Glob. Optim. 1997, 11, 341–359. [Google Scholar] [CrossRef]
  17. Zhang, J.; Sanderson, A.C. JADE: Adaptive differential evolution with optional external archive. IEEE Trans. Evol. Comput. 2009, 13, 945–958. [Google Scholar] [CrossRef]
  18. Tanabe, R.; Fukunaga, A.S. Improving the search performance of SHADE using linear population size reduction. In Proceedings of the IEEE Congress on Evolutionary Computation (CEC), Beijing, China, 6–11 July 2014; pp. 1658–1665. [Google Scholar]
  19. Deng, W.; Shang, S.; Cai, X.; Zhao, H.; Song, Y.; Xu, J. An improved differential evolution algorithm and its application in optimization problem. Soft Comput. 2021, 25, 5277–5298. [Google Scholar] [CrossRef]
  20. Tong, S.; Hu, L.; Wang, L. A hybrid discrete differential evolution algorithm for project scheduling problem. J. Intell. Manuf. 2023, 34, 763–778. [Google Scholar]
  21. Azarbadegan, B.; Jam, T.; Mukhopadhyay, A. A modified differential evolution algorithm for solving the traveling salesman problem. Eng. Optim. 2023, 55, 2089–2114. [Google Scholar]
  22. PShidiqi, T.; Aji, R.N.; Purnomo, A. Optimization of MRI image enhancement using improved differential evolution algorithm. J. Eng. Appl. Sci. 2022, 17, 1436–1442. [Google Scholar]
  23. Li, Y.; Zhang, L.; Wang, S.; Yang, J. Image segmentation using particle swarm optimization and expectation conditional maximization algorithm. Soft Comput. 2022, 1–15, 123–145. [Google Scholar]
  24. Zhang, Y.; Sun, L.; Wang, J. Image segmentation based on ant colony optimization for cost-sensitive artificial bee colony. J. Intell. Fuzzy Syst. 2022, 43, 3097–3105. [Google Scholar]
  25. Guo, G.; Chen, X.; Chen, W. Image segmentation using genetic algorithm optimized local region-based active contour model. Front. Comput. Sci. 2023, 17, 57–76. [Google Scholar]
  26. Kumar, S.; Pant, M.; Ray, A.K. Differential evolution embedded Otsu’s method for optimized image thresholding. In Proceedings of the Information and Communication Technologies (WICT), 2011 World Congress on, Mumbai, India, 11–14 December 2011; pp. 1–6. [Google Scholar]
  27. Ayala, H.V.H.; Santos, F.M.D.; Mariani, V.C. Image thresholding segmentation based on a novel beta differential evolution approach. Expert Syst. Appl. 2015, 42, 2136–2142. [Google Scholar] [CrossRef]
  28. Otsu, N. A thresholding selection method from gray-level histograms. IEEE Trans. Syst. Man Cybern. 1979, 9, 62–66. [Google Scholar] [CrossRef]
  29. Sun, Y.; Yang, Y. An Adaptive Bi-Mutation-Based Differential Evolution Algorithm for Multi-Threshold Image Segmentation. Appl. Sci. 2022, 12, 5759. [Google Scholar] [CrossRef]
  30. Stanovov, V.; Akhmedova, S.; Semenkin, E. Biased Parameter Adaptation in Differential Evolution. Inf. Sci. 2021, 566, 215–238. [Google Scholar] [CrossRef]
  31. Xiang, W.-L.; Meng, X.-L.; An, M.-Q. An Enhanced Differential Evolution Algorithm Based on Multiple Mutation Strategies. Comput. Intell. Neurosci. 2015, 2015, 285730. [Google Scholar] [CrossRef]
  32. Mohamed, A.W. An improved differential evolution algorithm with triangular mutation for global numerical optimization. Comput. Ind. Eng. 2015, 85, 359–375. [Google Scholar] [CrossRef]
  33. Yi, W.; Gao, L.; Li, X.; Zhou, Y. A new differential evolution algorithm with a hybrid mutation operator and self-adapting control parameters for global optimization problems. Appl. Intell. 2015, 42, 642–660. [Google Scholar] [CrossRef]
  34. An, H.C.; Hao, W. A differential evolution SAF-DE algorithm which jumps out of local optimal. In Proceedings of the 16th International Conference on Computational Intelligence and Security (CIS), Guangxi, China, 27–30 November 2020; pp. 333–337. [Google Scholar]
  35. Guo, Y. Image segmentation based on a hybrid artificial bee colony algorithm with SSIM and PSNR evaluation. Eng. Comput. 2022, 52, 1–16. [Google Scholar]
  36. Sun, Y. Differential evolution algorithm with SSIM and PSNR metrics for image segmentation. Multimed. Tools Appl. 2022, 81, 15577–15594. [Google Scholar]
Figure 1. (a) Represents the original image, (b) segmented image with 2 thresholds, (c) segmented image with 3 thresholds, (d) segmented image with 4 thresholds, (e) segmented image with 5 thresholds, (f) segmented image with 6 thresholds (TRDE).
Figure 1. (a) Represents the original image, (b) segmented image with 2 thresholds, (c) segmented image with 3 thresholds, (d) segmented image with 4 thresholds, (e) segmented image with 5 thresholds, (f) segmented image with 6 thresholds (TRDE).
Biomimetics 08 00418 g001aBiomimetics 08 00418 g001b
Figure 2. (a) Represents the original image and (b) the segmented image with TRDE, (c) the Canny image segmentation method, (d) the Sobel image segmentation method, (e) the Robert image segmentation method, and (f) the Prewitt image segmentation method.
Figure 2. (a) Represents the original image and (b) the segmented image with TRDE, (c) the Canny image segmentation method, (d) the Sobel image segmentation method, (e) the Robert image segmentation method, and (f) the Prewitt image segmentation method.
Biomimetics 08 00418 g002aBiomimetics 08 00418 g002b
Table 1. PSNR metrics for the JADE, SHADE, TRDE, SAF-DE, BDE and HSDE.
Table 1. PSNR metrics for the JADE, SHADE, TRDE, SAF-DE, BDE and HSDE.
ImagesK Mean PSNR Value
JADESHADETRDESAF-DEBDEHSDERank
baboon216.737815.903916.439115.673815.545615.72542
316.378916.097617.893616.098616.878916.78361
419.258518.325219.356219.035819.232919.02311
521.254820.189922.824821.907822.980621.46782
622.358420.314523.025421.302322.175921.39851
barbara218.752616.596218.369518.203618.268418.32062
320.256219.365421.148520.369821.036921.05891
422.256219.258522.589521.148521.658421.85481
523.541220.628524.562122.654124.625422.15232
623.095620.659824.356423.458924.213623.69841
cameraman215.340514.469816.422315.424215.134415.11441
318.646417.456618.735217.454717.455418.54681
419.687618.464520.468720.456420.457819.45271
521.252420.454622.454022.867622.456220.74622
620.476620.786222.678621.454522.456421.45321
lena 219.875718.466620.456619.464619.453419.67831
319.458718.453420.478219.782220.878619.45722
420.867719.456820.877820.786421.456420.78622
522.786221.786723.767622.213522.782421.45541
622.786221.867223.786722.786223.427822.72401
peppers220.254219.453420.553820.123020.201320.48781
323.456422.454523.477623.146323.456423.45331
424.454823.582125.486725.453423.453424.45321
526.786324.453426.456427.464627.456425.13204
624.436423.453426.135526.536826.453325.34513
plane217.135517.135517.151517.123517.135117.24331
318.345218.345618.351318.315317.353118.56782
423.315319.135524.135323.315824.315321.45642
523.345921.464825.315324.353124.353423.45561
623.345323.348323.415523.345323.135423.68331
saturn217.432417.456418.345517.343117.315417.45451
319.345319.345320.456319.211219.213319.78961
421.564521.345322.254322.354322.154622.34533
522.454622.312223.412823.123423.343423.34561
622.454522.545623.343222.925622.153122.45471
woman215.312316.546716.456516.453216.213516.45342
319.343418.456419.536919.451219.315519.48291
421.345320.456521.335521.123321.354121.14423
522.346821.344622.123821.345321.124321.79782
621.878621.467822.343120.345321.215421.78281
Table 2. SSIM metrics for the JADE, SHADE, TRDE, SAF-DE, BDE and HSDE.
Table 2. SSIM metrics for the JADE, SHADE, TRDE, SAF-DE, BDE and HSDE.
ImagesK Mean SSIM Value
JADESHADETRDESAF-DEBDEHSDERank
baboon20.645530.653470.683480.634530.638770.678651
30.745340.758660.758770.712380.787670.783243
40.834560.812340.833450.834550.814370.815372
50.864550.884340.914330.904530.873730.822451
60.944830.856650.945640.976840.945680.897862
barbara20.387640.446780.454640.394480.384650.346461
30.546450.444550.578760.578640.578460.567761
40.704530.588860.713440.698740.710440.698541
50.844530.694560.872130.824650.853450.724531
60.804530.687640.894640.813450.886440.776461
cameraman20.404560.456450.445450.434540.413430.34642
30.568640.546150.594560.584560.545650.645352
40.615340.653120.646310.689460.664640.701345
50.721250.704560.745680.745610.745640.745621
60.745650.745640.795630.797350.745640.776452
lena20.645660.554450.653480.645560.645640.654342
30.845610.846760.864620.746780.845640.848961
40.868780.843450.873540.823450.868760.834551
50.904320.894560.884730.898450.898540.884565
60.884560.845620.903450.884560.864560.884551
peppers20.494560.445640.554640.545640.545620.532451
30.534550.545450.545670.587650.53450.534522
40.534530.545340.545440.554330.534330.545242
50.531530.553430.543540.545330.512340.578663
60.564560.545340.564530.564670.564560.564655
plane20.145640.186760.195640.154540.154640.145641
30.286660.245680.254660.245640.256450.223433
40.774560.426450.794650.745640.786780.478641
50.784560.746540.884560.856640.755660.744241
60.745680.746780.845680.787450.764560.745661
saturn20.845660.845640.854750.854660.815640.854762
30.845640.854640.887660.845640.845640.845641
40.987640.845640.845640.864560.864560.845644
50.954340.845430.954430.945640.845640.845341
60.953430.954340.971630.934540.954640.945641
woman20.364660.368660.364560.356460.364540.364563
30.465460.454640.465560.456460.454660.464563
40.456460.464560.486550.445680.445640.445661
50.556410.554640.554660.464660.564560.564664
60.565670.532450.743130.664660.676450.545661
Table 3. SSIM metrics for the TRDE, Canny, Prewitt, Robert and Sobel methods.
Table 3. SSIM metrics for the TRDE, Canny, Prewitt, Robert and Sobel methods.
ImagesAverage SSIM Values of TRDE and Other Segmentation Methods
TRDECannyPrewittRobertSobel
baboon0.293290.004310.000750.000130.00363
barbara0.332790.012410.011240.012310.00224
cameraman0.618590.122190.113200.132650.10808
lena0.424940.010500.005640.011300.00046
peppers0.354510.018300.004620.002320.00124
plane0.639980.005720.005250.003480.00839
saturn0.732600.530860.506500.506480.49680
woman0.399270.010210.004320.014360.00868
Table 4. MS-SSIM metrics for the TRDE, Canny, Prewitt, Robert and Sobel methods.
Table 4. MS-SSIM metrics for the TRDE, Canny, Prewitt, Robert and Sobel methods.
ImagesAverage MS-SSIM Values of TRDE and Other Segmentation Methods
TRDECannyPrewittRobertSobel
baboon0.299040.002520.001760.002630.00634
barbara0.339190.010560.010540.005470.31222
cameraman0.606840.123140.116370.138350.11504
lena0.423840.008610.004790.008350.00339
peppers0.353380.015370.000540.002640.00203
plane0.635900.003220.003840.005480.01415
saturn0.720520.544050.515200.519480.51394
woman0.386470.005350.012220.023720.01807
Table 5. UQI metrics for the TRDE, Canny, Prewitt, Robert and Sobel methods.
Table 5. UQI metrics for the TRDE, Canny, Prewitt, Robert and Sobel methods.
ImagesAverage UQI Values of TRDE and Other Segmentation Methods
TRDECannyPrewittRobertSobel
baboon0.00012460.00010070.00014740.00011756.2078237901 × 10−5
barbara0.00013170.00017520.00012409.5289199533 × 10−57.0011508332 × 10−5
cameraman0.00017960.00016110.00010069.8649176392e × 10−58.0869347378 × 10−5
lena0.00013800.00016810.00011968.9136782626e × 10−56.3427345646 × 10−5
peppers0.00012420.00019710.00011180.00012347.4385365087 × 10−5
plane0.00018190.00014470.00015327.4707859548 × 10−55.1055016436 × 10−5
saturn0.00016730.00015240.00023560.00012550.0001545
woman0.00013070.00012390.00015019.6613189674 × 10−56.2780501489 × 10−5
Table 6. Average execution time for JADE, SHADE, TRDE, SAF-DE, BDE and HSDE.
Table 6. Average execution time for JADE, SHADE, TRDE, SAF-DE, BDE and HSDE.
ImagesKExecutionTimes (s)
JADESHADETRDESAF-DEBDEHSDERank
baboon21.445641.464561.413451.676151.979351.456851
31.467661.487661.445641.843431.645641.445651
41.687861.446461.448661.648861.986762.245682
51.586751.445451.445641.778462.786782.867622
61.534361.453461.453341.454442.134542.345351
barbara21.453411.434821.423451.743321.787341.454341
31.578631.483121.445631.855361.978831.445322
41.514561.416511.416161.856872.045642.545541
51.515361.454651.454431.954342.054462.345351
61.545331.445351.424631.945642.145342.345651
cameraman21.484341.445351.445331.854331.854351.531551
31.556431.456451.545641.945641.943341.553342
41.545341.454651.456461.956452.356452.056452
51.764661.456661.584662.184662.156482.505462
61.664651.545661576462.145642.246452.545642
lena21.546461.556461.456451.856431.945661.545611
31.564481.546461.446651.748781.456481.456511
41.564681.564651.664661.956542.445642.445693
51.612351.545641.645662.34642.145642.464663
61.686451.534551.564642.456462.243452.453452
peppers21.956461.964561.846542.354652.864601.956401
31.464561.964651.864642.645642.364651.956432
42.054641.954651.456462.645652.644552.456621
52.315342.543441.945332.553452.415462.856061
62.534342.345332.135452.753452.453442.456471
plane21.646461.454561.745381.875462.054451.554664
31.846451.545561.645601.914541.945641.545653
41.615561.546461.664561.956452.646652.534553
51.734531.545341.453432.043442.254342.354531
61.655461.545671.545642.064642.245672.445641
saturn21.686451.546541.564881.986481.935431.554332
31.745451.684631.635432.135352.153431.835431
41.745641.645561.645642.164562.356452.554662
52.145641.754651.754632.165462.345662.648641
61.954651.846451.756462.256462.456432.754561
woman21.956461.945431.856452.245642.464552.046451
32.545642.356452.356422.664652.256462.656461
42.984681.615561.545641.756482.464572.064561
51.684661.564641.584682.186782.178672.484422
61.784981.648481.548462.184682.584682.245641
Table 7. Standard deviation of the fitness value (k = 2, 3, 4, 5, 6).
Table 7. Standard deviation of the fitness value (k = 2, 3, 4, 5, 6).
ImagesK Standard Deviation of Fitness Value
JADESHADETRDESAF-DEBDEHSDERank
baboon200.0354300.01235400.045641
30.014530.034530.004340.005640.004780.046981
40.045520.046780.001350.043450.045440.010451
500.045380.004560.008790.064530.034552
60.013450.043440.003450.004560.034550.001342
barbara20.034640.1456300.045380.054640.086871
30.054640.345750.045640.054560.356440.153121
40.045680.345650.034560.045640.245640.078751
50.054680.375450.034360.045650.453440.076451
60.045670.565420.013450.087850.046340.037851
cameraman20.054640.045640.001540.034560.014568.61 × 10−62
30.007850.045640.003250.004530.045340.003541
40.565460.456870.575440.648430.684540.694343
50.015870.034550.005340.000450.015430.025432
60.008540.023450.004540.034460.045660.004561
lena20.054640.345640.005460.546440.145450.059441
30.054660.456640.023650.146650.378690.667801
40.156460.245640.026540.065440.535650.066561
50.074560.146460.065450.045450.157650.056874
60.043830.113460.024560.007670.045640.015462
peppers20.016460.343770.013430.246680.054450.145461
30.045660.154640.006450.145660.145660.064561
40.054300.164640.034340.064580.146450.045551
50.044550.173540.015440.354340.045560.054621
60.046450.087830.023450.073760.553880.034531
plane20.043450.345640.008760.278660.386860.067781
30.033240.234530.004540.008760.345350.154341
40.085540.154640.035450.045350.124550.153442
50.044450.254350.045610.034560.254640.03453
60.045650.453450.067710.064520.045620.05455
saturn20.032450.0454500.0435500.034541
30.045640.054640.045540.045650.544340.045432
40.097540.145640.045440.154650.045640.054451
50.144550.097870.056440.045640.245640.154442
60.024530.245340.004540.056450.145340.074561
woman20.156540.157670.001120.156420.004560.145671
30.054640.156440.052440.356440.345640.000252
40.054560.456440.045670.254650.054640.234571
50.057840.156440.015440.000420.006450.054342
60.045440.154640.042440.154640.015440.097982
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guo, Y.; Wang, Y.; Meng, K.; Zhu, Z. Otsu Multi-Threshold Image Segmentation Based on Adaptive Double-Mutation Differential Evolution. Biomimetics 2023, 8, 418. https://doi.org/10.3390/biomimetics8050418

AMA Style

Guo Y, Wang Y, Meng K, Zhu Z. Otsu Multi-Threshold Image Segmentation Based on Adaptive Double-Mutation Differential Evolution. Biomimetics. 2023; 8(5):418. https://doi.org/10.3390/biomimetics8050418

Chicago/Turabian Style

Guo, Yanmin, Yu Wang, Kai Meng, and Zongna Zhu. 2023. "Otsu Multi-Threshold Image Segmentation Based on Adaptive Double-Mutation Differential Evolution" Biomimetics 8, no. 5: 418. https://doi.org/10.3390/biomimetics8050418

APA Style

Guo, Y., Wang, Y., Meng, K., & Zhu, Z. (2023). Otsu Multi-Threshold Image Segmentation Based on Adaptive Double-Mutation Differential Evolution. Biomimetics, 8(5), 418. https://doi.org/10.3390/biomimetics8050418

Article Metrics

Back to TopTop