An Adaptive Cultural Algorithm with Improved Quantum-behaved Particle Swarm Optimization for Sonar Image Detection

This paper proposes an adaptive cultural algorithm with improved quantum-behaved particle swarm optimization (ACA-IQPSO) to detect the underwater sonar image. In the population space, to improve searching ability of particles, iterative times and the fitness value of particles are regarded as factors to adaptively adjust the contraction-expansion coefficient of the quantum-behaved particle swarm optimization algorithm (QPSO). The improved quantum-behaved particle swarm optimization algorithm (IQPSO) can make particles adjust their behaviours according to their quality. In the belief space, a new update strategy is adopted to update cultural individuals according to the idea of the update strategy in shuffled frog leaping algorithm (SFLA). Moreover, to enhance the utilization of information in the population space and belief space, accept function and influence function are redesigned in the new communication protocol. The experimental results show that ACA-IQPSO can obtain good clustering centres according to the grey distribution information of underwater sonar images, and accurately complete underwater objects detection. Compared with other algorithms, the proposed ACA-IQPSO has good effectiveness, excellent adaptability, a powerful searching ability and high convergence efficiency. Meanwhile, the experimental results of the benchmark functions can further demonstrate that the proposed ACA-IQPSO has better searching ability, convergence efficiency and stability.

Sonar imaging has the advantages of long-distance detection and strong penetration. Therefore, it is being extensively used for underwater inspections, hydrographic and bathymetric surveys, underwater positioning, dredging, cable-laying, pipe-line inspections, and numerous other marine applications. The underwater sonar image contains three kinds of regions, including object-highlight, shadow and background regions. The purpose of object detection is to segment the object-highlight and shadow regions from complex background region and preserve as much of the original edge information of underwater sonar image as possible 1 . A number of techniques have been proposed to detect underwater objects of the sonar image. A fuzzy clustering method on the sonar image has been proposed to solve the detection problem 2 . However, fuzzy clustering is very sensitive to speckle noise. A Markov segmentation algorithm was used for three-class sonar image detection 3 . Later, many improved image detection algorithms on Markov random field model were proposed by scholars 4 . Although their results are satisfactory, the processing procedures are quite complex and computationally costly. Maria Lianantonakis and Yvan R. Petillot developed active contours and level set methods, which were applied to the detection of regions like the seabed 5 . Subsequently, Implicit Active Contours were used in sonar image detection 6 . On this basis, Xiu-FenYe et al. proposed a new detection method of sonar images, Gauss-Markov random field model was integrated into level set energy function model to dynamically extract regions of interest 7 . Wang Xingmei et al. presented a narrowband Chan-Vese model by adaptive ladder initialization to precisely and quickly detect underwater objects of sonar images 8 . However, sometimes if the block mode k-means clustering algorithm cannot quickly and accurately complete initial segmentation, approximate position of object-highlight and shadow regions will not be obtained, and higher detection precision will not be possible.
In recent years, the cultural algorithm (CA) has gradually attracted more global attention. It can provide a powerful framework for the solution of complicated problems 9 . Youlin Lu et al. proposed a hybrid multi-objective cultural algorithm to solve short-term hydrothermal scheduling problems, which combined the differential evolutionary algorithm with CA framework 10 . It can obtain a more accurate solution. Later, Zhou Wei et al. presented the cultural particle swarm optimization algorithm (CPSO) to solve the partner selection problems of virtual enterprise 11 . The proposed algorithm has some feasibility and efficiency. To solve the optimization problems, Noor H. Awad et al. proposed a CA with an improved local searching algorithm 12 . Compared with the basic CA, the performance of this algorithm is greatly improved. Mostafa Z. Ali et al. combined the niche algorithm with the Tabu search algorithm in CA to solve engineering optimization problems, which was efficient and robust to some extent 13 . Tianyu Liu et al. introduced the QPSO into CA to solve multiobjective optimization problems 14 . The proposed algorithm has high efficiency. Through the comparative analysis, many intelligent optimization algorithms can be used as evolution strategy in the population space, which can increase population diversity, improve the searching ability and promote efficiency. Among intelligent optimization algorithms, QPSO is one of most commonly used algorithms in the population space of CA. This algorithm supposes that particles have quantum behaviour such that particles are attracted by a quantum potential well centered on its local attractor. It has fewer parameters and a relatively good searching ability 15 . The contraction-expansion coefficient is the only parameter in the QPSO algorithm, which plays an important role in balancing the global and local searching abilities. However, it is a random value, which easily leads to blindness in the searching process. To solve this problem, the revised QPSO algorithm regarded iterative times as an important factor to adjust the contraction-expansion coefficient 16 . The contraction-expansion coefficient can linearly decrease with the increase of iterative times in the interval [0. 5 1). Although this method is often used in practice, it only solves the linear problem and easily falls into the local optimal solution in the searching process of complex problems. Subsequently, in order to improve the performance of the QPSO algorithm, Jun Sun et al. further presented a diversity-maintained QPSO algorithm 17 . When the population diversity is lower than a set value, contraction-expansion coefficient is set as the boundary of convergence, or contraction-expansion coefficient linearly decreases. On this basis, Tian Jin constructed a new contraction-expansion coefficient using Sigmoid function to solve high-dimensional multimodal functions optimization problems in the QPSO algorithm 18 . This method can make the contraction-expansion coefficient decrease nonlinearly with the iterative times and increase the flexibility of the QPSO algorithm. Although these adjustments of the contraction-expansion coefficient in the QPSO algorithm can to some extent improve searching ability, these algorithms only regard iterative times as the factor to adjust the contraction-expansion coefficient. In fact, the quality of particles also influences their searching ability. Therefore, the fitness value is an important factor to adjust the contraction-expansion coefficient.
In addition, the communication protocol between population space and belief space can influence the performance of CA. Accept function is used to set the rate of accepted individuals in the population space, and the rate of accepted individuals usually decreases with iterative times. Ricardo Landa Becerra and Carlos A.Coello improved accept function by resetting the rate of accepted individuals when the best solution had not changed over the last several iterations 19 . Since this method can improve the performance of CA, it is widely used in CAs 20,21 . However, appropriate range of rate is different when solving different problems, the rate of accepted individuals still need to be reset. To enhance the flexibility, the fitness values of individuals were used in accept function 22 . Although this adjustment can somewhat enhance the flexibility, the utilization of evolution information was not adequately used. Influence function is used to guide the evolution of poor individuals in the population space by the knowledge in the belief space. Different types of knowledge are chosen according to different problems 23,24 . These types of knowledge can determine the searching step size and searching direction of the individuals in the population space.
In these regards, to obtain more accurate detection results, this paper presents an ACA-IQPSO to detect underwater sonar image. In the population space, the contraction-expansion coefficient of the IQPSO is adaptively adjusted according to iterative times and the fitness value of particles. In the belief space, a new update strategy is adopted to update the cultural individuals using the idea of SFLA. In addition, accept function and influence function are redesigned in the new communication protocol. The new communication protocol can make belief space with adequate evolutionary information that can more precisely guide the evolution of particles in the population space and further improve the searching ability of the algorithm. The new communication protocol can enhance convergence efficiency of the algorithm. The experimental results demonstrate that the ACA-IQPSO can locate good clustering centres according to the grey distribution information of underwater sonar images, and accurately complete underwater object detection. Through the analysis of benchmark functions, it can show that the proposed ACA-IQPSO is significantly better than other algorithms in searching ability, convergence efficiency and stability. Therefore, the proposed method has important theoretical and practical value.

CA and QPSO.
CA. CA is a model with double level evolutionary space. It defines both population space and belief space. The two spaces evolve respectively, and the communication protocol between these spaces is accomplished by accept function and influence function. The schematic diagram is shown in Fig. 1. In the population space, the intelligent optimization algorithm can be used as evolution strategy to achieve the evolution of individuals. These individuals in the population space can contribute their experience and evolution information to the belief space through accept function. In the belief space, the experience and evolutionary information are converted into the cultural individuals. And the knowledge is extracted from the cultural individuals to guide the evolution of individuals in the population space. There are five types of knowledge in the belief space, including situational knowledge, normative knowledge, historical knowledge, topographical knowledge and domain knowledge. These types of knowledge are used to influence the evolution of poor individuals in the population space by influence function, and guide individuals to generate elite offspring. QPSO. The QPSO supposes that particles have quantum behaviour, and they move in a quantum potential well centred on its local attractor. When one particle finds a local optimal state, some particles quickly converge to it, and others appear on any position of the whole space in a certain probability. Particles with quantum-behaved maintain high levels of population diversity, which can avoid falling into the local optimal solution to a certain extent. In addition, it only updates the position without speed. Therefore, the QPSO has relatively fewer parameters and a good searching ability.
In the D-dimensional space, N is the population size. The particle swarm is defined as X = [X 1 , X 2 , X 3 , ... X N ], and X i = [X i1 , X i2 , X i3 , ... X iD ] represents the current position of the i th particle. The position update of particles is shown as:  where μ is a random number in the interval [0,1], P m (t) is the local attractor, L(t) is the characteristic length of the wave function, and t is the current iterative times. P m (t) is defined as: where P i (t) is the personal best position of the particle, and G(t) is the global best position. L(t) is expressed as: i where β is the contraction-expansion coefficient, which can control the convergence speed of the algorithm. The contraction-expansion coefficient β has an important influence on the searching ability of particles. mbest (t) is the mean best position of all particles. mbest(t) is given by: According to Eq. (1) and Eq. (3), the position of the particles can be expressed as: The proposed ACA-IQPSO. Population space. The IQPSO is integrated into CA as the evolution strategy of the population space in the ACA-IQPSO. IQPSO: In the QPSO, the quality of particles influences their behaviour in the searching process. The contraction-expansion coefficient is an important factor to control the particles' behaviour. It only regards iterative times as the factor to generally adjust the contraction-expansion coefficient. However, this method easily falls into the local optimal solution in the searching process of complex problems. Furthermore, the quality of particles depends on information carried by themselves in each iteration, and all the information also has important influence on particles' behaviour in the searching process. Therefore, iterative times and the fitness value of particles are used to adaptively adjust the contraction-expansion coefficient in IQPSO. In each iteration, when the quality of particles is worse, the contraction-expansion coefficient is larger, and the global searching ability of the particles is relatively stronger. When the quality of particles is better, the contraction-expansion coefficient is relatively smaller, and the local searching ability of the particles is stronger.
The contraction-expansion coefficient β is defined as follows: where T is the maximum iterative times, f g is the fitness value of G(t), f is the average fitness value of particles in the population space, f i is the fitness value of the i th particle, and f w is the fitness value of the worst particle. k is a positive integer, that can adjust β to balance the relationship between iterative times and the fitness value of particles. The performance analysis of IQPSO: To verify the superiority of the IQPSO in searching ability, Sphere function and Griewank function are used in this paper to test position distribution of particles. Sphere function is unimodal and contains only one global optimal solution. Griewank function is multimodal and contains many local optimal solutions, but only one global optimal solution. Fig. 2 shows position distribution of particles in IQPSO and QPSO 16 . The relevant parameters are as follows. The dimension of the solution space is 2, the population size is 30, and the maximum number of iterations is 5.
As seen from Fig. 2, the solution scope is relatively larger in the IQPSO, and IQPSO more easily obtains the global optimal solution. Therefore, compared with the QPSO, the IQPSO can increase population diversity and improve searching ability.
To further verify the effectiveness of the IQPSO in searching ability, the fitness values are calculated by Sphere and Griewank functions in the IQPSO and the QPSO. The optimization results are shown in Fig. 3. The relevant parameters are as follows. The dimension of the solution space is 10, the population size is 30, the maximum number of iterations is 30, and the experiment runs for 30 times in each algorithm.  As depicted in Fig. 3, whether it is unimodal function or multimodal function, the fitness values of the IQPSO are always smaller than the QPSO in the iterative process. Therefore, the optimal solutions of the IQPSO are always better than the QPSO. The optimization results further show that the searching ability of the IQPSO is better than the QPSO.
Belief space. In the ACA-IQPSO, the cultural individuals in the belief space are updated by a new update strategy. Meanwhile, three types of knowledge are extracted from the cultural individuals. These types of knowledge are used to influence the evolution of particles in the population space by influence function, and guide particles to generate elite offspring.
Update strategy: Cultural individuals are updated by selection, crossover and mutation in many CAs 11 . However, this update strategy needs many parameters that make the algorithm not easy to control and lead to poor searching ability. In SFLA, the worst frog is updated by the step size that is the difference between the best frog and the worst frog. This update strategy has few parameters, powerful searching ability and is easy to implement 25 . Therefore, a new update strategy is adopted to update the cultural individuals in the belief space according to the idea of the update strategy in SFLA. The new update strategy adopts the difference between the situational knowledge and the cultural individuals as step size. New update strategy is expressed as follows: where Y j (t) is the position of the j th cultural individual, S(t) is situational knowledge of the t th iteration in belief space, and rand is a random number in the interval [0,1]. Knowledge structure: Situational knowledge, normative knowledge and domain knowledge are adopted in the belief space of ACA-IQPSO. In each iteration, these types of knowledge are updated by cultural individuals.
Situational knowledge contains the best cultural individual in the belief space. In each iteration, when the fitness value of the best cultural individual is larger than that of the current situational knowledge, the best cultural individual is used as situational knowledge in the update process. Update formula of situational knowledge is: Normative knowledge describes the feasible solution space of the problem in the belief space. It stores the boundary information that can guide particles in the population space to search for a better region. The update of normative knowledge reflects the change of the feasible search space. With the increase of iterative times, the searching scope can be concentrated in the dominant region. Therefore, when an excellent cultural individual is out of the current searching scope in each iteration, normative knowledge is updated. Update formula of normative knowledge is as follows: is lower limit of fitness value in the t th iteration, f u (t) is upper limit of fitness value in the t th iteration. Domain knowledge is used to statically or dynamically guide the particles in the population space to evolve along the predictive direction, and record a good evolutionary trend. In the belief space, the centre of gravity can reflect the overall distribution of the cultural individuals and guide the particles in the population space to search better solutions. The evolutionary direction of the particles in the population space can be predicted by updating the centre of gravity in each iteration, which can improve the searching efficiency. Therefore, domain knowledge stores the centre of gravity in the belief space to guide the evolution of particles in the population space in this paper. The centre of gravity GT(t) is defined as follows: Accept function: Accept function is used to set the rate of accepted particles in the population space. To further improve the flexibility of the rate, accept function is redesigned in this paper. The new accept function can adaptively change the rate of accepted particles according to the quality of particles in each iteration. Meanwhile, it can enhance utilization degree of information in the population space. When the quality of particles is better, the accept rate is larger. The new accept function is defined as follows: where f is the average fitness value of particles in the population space, f w is the fitness value of the worst particle in the population space, and f i is the fitness value of the i th particle.
Therefore, the number of accepted particles in population space is calculated according to Eq. (12). It is: where N is the population size.
For new accept function, Eq. (12) is analysed in detail. The fitness values of particles in population space are sorted in descending order. Therefore, if there is a position c where the fitness value f c is larger than f , and the fitness value f c + 1 is smaller than f . Af is: Therefore, Eq. (15) is as follows: Af As seen from Eq. (15), when c = N − 1, the maximum value of Af is 0.5. When c is relatively larger, there are more particles whose fitness values are larger than the average fitness value, and the quality of particles in the population space is better. Correspondingly, more particles are accepted. When c is relatively smaller, there are fewer particles whose fitness values are larger than the average fitness value, and the quality of particles in the population space is worse. Correspondingly, fewer particles are accepted. Therefore, new accept function is effective in this paper.
Influence function: Influence function is used to guide the evolution of poor particles in the population space by using the knowledge in the belief space. These types of knowledge determine searching step size and searching direction of the particles in the population space. To enhance utilization of information in the belief space and more precisely guide the evolution of poor particles in the population space, influence function is also redesigned in this paper. The new influence function adopts situational, normative and domain knowledge to guide the evolution of poor particles in the population space. Situational and domain knowledge are used to determine searching step size, and normative knowledge is used to determine searching direction. When the fitness values of the poor particles are less than the lower limit of fitness value in the belief space, domain knowledge is used to control searching step size. When the fitness values of poor particles are less than the upper limit of fitness value in the belief space, situational knowledge is used to control searching step size. Otherwise, the positions of poor particles are randomly generated in the solution space.
The influence function is defined as follows:  Fig. 4.The relevant parameters are as follows. The dimension of the solution space is 10, the population size is 30, the maximum number of iterations is 30, the size of cultural individuals is 8, and the experiment runs for 30 times in each algorithm.
As depicted in Fig. 4, whether it is unimodal function or multimodal function, the ACA-IQPSO can converge to global optimal solution after 10 iterations, convergence speed of ACA-IQPSO is faster in the iterative process. It can indicate that convergence efficiency of the ACA-IQPSO is superior to the CA-IQPSO. Meanwhile, the fitness values of the ACA-IQPSO are smaller than the CA-IQPSO. The new communication protocol can make belief space with adequate evolutionary information that can more precisely guide the evolution of particles in the population space. Therefore, the new communication protocol can further improve searching ability remarkably.

Results and Discussion
Experimental results and discussion of sonar image detection. In this section, the original sonar image has serious noise. To remove some noise points and make the image smoother, Butterworth lower-pass  filter is used for the noise smoothing 26 . On this basis, numerical examples are shown to validate the effectiveness and adaptability of the proposed ACA-IQPSO for sonar image detection. Meanwhile, the ACA-IQPSO is compared with CA-IQPSO, CPSO 11 , IQPSO, QPSO 16 , and PSO 27 . In addition, the fitness function is mainly used to evaluate the quality of particles in the process of sonar image detection. Therefore, a fitness function combining intra-class difference with inter-class difference is adopted in these algorithms 28 .When the fitness value is larger, the detection result is better. The relevant parameters are as follows. The number of clustering centers is 4, the population size is 20, the maximum number of iterations is 20, the acceleration coefficient is c1 = c2 = 2, the inertia weight is w = 0.8 in CPSO and PSO, size of cultural individuals is M = 8 in the ACA-IQPSO, CA-IQPSO and CPSO, and the contraction-expansion coefficient β linearly decreases from 1.0 to 0.5 in the QPSO.   As depicted in Fig. 5, the proposed ACA-IQPSO can better detect object-highlight and shadow regions from complex background region. The edge information is preserved as much as possible in Fig. 5(c). While the detection result of the CA-IQPSO can only detect the shadow region in Fig. 5(d), it is not suitable for detecting the sonar image with floating objects. Although the CPSO can detect the object-highlight and shadow regions, it has serious information loss in the object-highlight region in Fig. 5(e). At the same time, IQPSO can also detect the object-highlight and shadow regions, but the detection result contains a lot of noise and the integrity of the underwater object is seriously weak in Fig. 5(f). The detection result of the QPSO has serious noise, which is not an ideal detection result in Fig. 5(g). PSO can only detect the object-highlight region in Fig. 5(h).Therefore, compared with the detection results of CA-IQPSO, CPSO, IQPSO, QPSO and PSO, the proposed ACA-IQPSO can locate good clustering centres according to the grey distribution information of the underwater sonar image with floating objects and accurately complete underwater object detection.
To further verify the effectiveness of the proposed ACA-IQPSO, Fig. 6 shows the detection results of the original sonar image with underwater stones on the bottom, which has a relatively weak contrast. Fig. 7 shows the detection results of structured seabed that is an object in sand ripples. Fig. 8 shows the detection results of larboard original sonar image including rocks, which are partly buried in the sand.  As seen from Figs 6, 7 and 8, the proposed ACA-IQPSO can obtain relatively better detection results in Figs 6(c), 7(c) and 8(c). While Fig. 6(d) has information loss in the object-highlight region. Although Fig. 6(e),(f),(g) and (h) can detect the object-highlight region and shadow region, the detection results have serious noise, especially in Fig. 6(g) and (h). Therefore, QPSO and PSO are failed to detect the sonar image with relatively weak contrast. Meanwhile, some information is lost in the object-highlight region and shadow region in Fig. 7(d),(e),(f),(g) and (h), which is not conducive to the subsequent feature extraction and underwater object recognition. In Fig. 8, the detection results of Fig. 8(d),(f),(g) and (h) show information loss in the shadow region and over detection in object-highlight region. At the same time, the detection results have different degrees of noise. Fig. 8(e) has serious information loss in object-highlight region. They are not ideal detection results.
Through the above comparative experiments, the proposed ACA-IQPSO can obtain relatively accurate results in sonar image detection. Moreover, the detection results of the CA-IQPSO are not better than ACA-IQPSO, which further verifies the effectiveness of the new communication protocol in this paper. Meanwhile, IQPSO can relatively obtain better detection results than QPSO, this indicates that IQPSO can improve searching ability of particles.
To demonstrate the advantages of proposed ACA-IQPSO more clearly, Table 1 shows the best fitness values after 20 iterations for the ACA-IQPSO and other intelligent optimization algorithms. Fig. 9 shows the variation of the best fitness values in each iteration.
It can be seen from Table 1 and Fig. 9 that the best fitness values of the proposed ACA-IQPSO are larger than other intelligent optimization algorithms after 20 iterations, which shows the powerful searching ability and high convergence efficiency of ACA-IQPSO. Meanwhile, the best fitness values of CA-IQPSO are larger than IQPSO  and the best fitness values of CPSO are larger than PSO, which can demonstrate the effectiveness of CA. Among IQPSO, QPSO and PSO, the best fitness values of IQPSO and QPSO are close, but the IQPSO is superior to the QPSO, and they are larger than PSO. These results indicate that IQPSO has the merit of searching ability.
To further verify the adaptability of the proposed ACA-IQPSO in this paper, Fig. 10 shows the detection results of another original sonar image with floating objects. Fig. 11 shows the detection results of starboard original sonar image with ship. Fig. 12 shows the detection results of an original sonar image with bottom tire.
From the results shown in Figs 10, 11 and 12, the proposed ACA-IQPSO can better detect original sonar image with floating objects, partly buried objects and objects on the bottom in this paper, and it has a certain effectiveness and adaptability. Moreover, it provides better preconditions for the subsequent feature extraction and underwater object recognition.
Experimental results and discussion of benchmark functions. To further verify the performance of the proposed ACA-IQPSO, 10 benchmark functions are used to compare the ACA-IQPSO with CA-IQPSO, CPSO, IQPSO, QPSO, and PSO. Among these benchmark functions, unimodal functions and multimodal functions are used to test local searching ability and global searching ability respectively. The details of benchmark functions are shown in Table 2.  Minimum value, maximum value, mean value and standard deviation are obtained after each algorithm running 20 times. These values are used to evaluate the performance of algorithms in this paper. The relevant parameters of each run are as follows. The population size is 50, the maximum number of iterations is 50, and the dimensions are 10, 30 and 50 respectively. The acceleration coefficient is c1 = c2 = 2, the inertia weight is w = 0.8 in CPSO and PSO, size of cultural individuals is M = 15 in ACA-IQPSO, CA-IQPSO and CPSO, and the contraction-expansion coefficient linearly decreases from 1.0 to 0.5 in QPSO.
The performance of the algorithms is related to the dimension. When the dimension increases, the performance differences will be more significant. To analyse the performance of the algorithms in different dimensions, Table 3 shows the minimum value, maximum value, mean value and standard deviation of the different algorithms when the dimension is 10. Table 4 shows the minimum value, maximum value, mean value and standard deviation of the different algorithms when the dimension is 30. Table 5 shows the minimum value, maximum value, mean value and standard deviation of the different algorithms when the dimension is 50.
From the comparative results in Table 3, Table 4 and Table 5, mean values of the IQPSO are smaller than the QPSO when the dimension is 10, 30 and 50, which indicates that the searching ability of the IQPSO is stronger than the QPSO. The standard deviation of the IQPSO and QPSO are close when dimension is 10, while standard    Similarly, Wilcoxon Signed-Rank Test and the Friedman test in nonparametric tests are adopted to further evaluate the performance of the algorithms by 10 benchmark functions in different dimensions 29,30 . Wilcoxon Signed-Rank Test is a paired comparison, which is used to compare the performance differences between two algorithms. If the p-values are less than or equal to the statistical significance value τ, the null hypothesis is rejected, which indicates that the two algorithms are different in performance. Friedman test is multiple comparison, which is used to compare the performance differences between more than two algorithms. When average rank is smaller, the performance of the algorithm is better. Table 6 shows the p-values of Wilcoxon Signed-Rank Test by 10 benchmark functions in different dimensions (the statistical significance value τ = 0.05). Table 7 shows the average rank of Friedman Test by 10 benchmark functions in different dimensions.
As seen in Table 6 and Table 7, the p-values between the ACA-IQPSO and other algorithms are smaller than the statistical significance value, which shows the significant performance differences between the ACA-IQPSO and other algorithms. Meanwhile, the average ranks of the ACA-IQPSO are smaller than the other algorithms. It can indicate that the performance of the ACA-IQPSO is superior to other algorithms. Therefore, through the analysis of benchmark functions, it can demonstrate that the proposed ACA-IQPSO is obviously better than other algorithms in searching ability, convergence efficiency and stability.

Conclusions
Considering the growing requirements of underwater sonar image detection, this paper proposed the ACA-IQPSO to detect underwater sonar images In the population space, iterative times and the fitness value of particles are used as factors to adaptively adjust the contraction-expansion coefficient of the QPSO. IQPSO can make particles adjust their behaviour to improve searching ability. In belief space, a new update strategy is adopted to update the cultural individuals according to the update strategy in SFLA. Moreover, to enhance utilization of information in population and belief spaces, accept function and influence function are redesigned in the new communication protocol. The new communication protocol can make belief space with adequate evolutionary information that can more precisely guide the evolution of particles in the population space and further improve the searching ability of the algorithm. Furthermore, the new communication protocol can enhance convergence efficiency of algorithm.
The proposed ACA-IQPSO is based on a clustering model. The object-highlight, shadow and background regions in the sonar image are detected by clustering centres. The experimental results demonstrate that ACA-IQPSO can locate good clustering centers and accurately complete underwater objects detection. Compared with other algorithms, the proposed ACA-IQPSO has good effectiveness and adaptability, and it has powerful searching ability and high convergence efficiency. Meanwhile, the performance of the proposed ACA-IQPSO is further demonstrated by the analysis of benchmark functions, it can show that the proposed ACA-IQPSO is significantly better than the other algorithms in searching ability, convergence efficiency and stability. Therefore, the proposed method can provide better preconditions for the subsequent feature extraction and underwater object recognition. It has important theoretical and practical value. 1