[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Extending Nighttime Combustion Source Detection Limits with Short Wavelength VIIRS Data
Next Article in Special Issue
An Imaging Algorithm for Multireceiver Synthetic Aperture Sonar
Previous Article in Journal
Fusion of GNSS and Satellite Radar Interferometry: Determination of 3D Fine-Scale Map of Present-Day Surface Displacements in Italy as Expressions of Geodynamic Processes
Previous Article in Special Issue
Comparison of Computational Intelligence Methods Based on Fuzzy Sets and Game Theory in the Synthesis of Safe Ship Control Based on Information from a Radar ARPA System
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Adaptive Denoising and Detection Approach for Underwater Sonar Image

1
College of Computer Science and Technology, Harbin Engineering University, Harbin 150001, China
2
Acoustic science and Technology laboratory, Harbin Engineering University, Harbin 150001, China
3
College of Underwater Acoustic Engineering, Harbin Engineering University, Harbin 150001, China
4
Key Laboratory of Marine Information Acquisition and Security, Harbin Engineering University, Harbin 150001, China
5
Institute of Acoustics, Chinese Academy of Science, Beijing 100190, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(4), 396; https://doi.org/10.3390/rs11040396
Submission received: 14 January 2019 / Revised: 9 February 2019 / Accepted: 12 February 2019 / Published: 15 February 2019
(This article belongs to the Special Issue Radar and Sonar Imaging and Processing)
Graphical abstract
">
Figure 1
<p>The estimated values of the optimal threshold (image size: 277 × 325): (<b>a</b>) Original sonar image; (<b>b</b>) The estimated values of the optimal threshold.</p> ">
Figure 2
<p>The estimated values of the optimal threshold (image size: 173 × 167): (<b>a</b>) Original sonar image; (<b>b</b>) The estimated values of the optimal threshold.</p> ">
Figure 3
<p>The estimated values of the optimal threshold (image size: 197 × 211): (<b>a</b>) Original sonar image; (<b>b</b>) The estimated values of the optimal threshold.</p> ">
Figure 4
<p>The estimated values of the optimal threshold (image size: 205 × 201): (<b>a</b>) Original sonar image; (<b>b</b>) The estimated values of the optimal threshold.</p> ">
Figure 5
<p>The estimated values of the optimal threshold (image size: 130 × 106): (<b>a</b>) Original sonar image; (<b>b</b>) The estimated values of the optimal threshold.</p> ">
Figure 6
<p>The estimated values of the optimal threshold (image size: 146 × 135): (<b>a</b>) Original sonar image; (<b>b</b>) The estimated values of the optimal threshold.</p> ">
Figure 7
<p>The result of the linear fitting based on <a href="#remotesensing-11-00396-t001" class="html-table">Table 1</a>.</p> ">
Figure 8
<p>Schematic diagram of CA (cultural algorithm).</p> ">
Figure 9
<p>The denoising results based on <a href="#remotesensing-11-00396-f001" class="html-fig">Figure 1</a>a of the proposed denoising method in this paper and the previous denoising method: (<b>a</b>) The denoising result of the denoising method proposed in this paper; (<b>b</b>) The result of filtering the degree parameter is selected by two adaptive thresholds in the proposed denoising method; (<b>c</b>) The denoising result of the previous denoising method; (<b>d</b>) The result of the filtering degree parameter is selected by two previous thresholds in the previous denoising method.</p> ">
Figure 10
<p>The detection results of the fuzzy c-means (FCM) based on <a href="#remotesensing-11-00396-f009" class="html-fig">Figure 9</a>b,d: (<b>a</b>) FCM detection result when <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>0.01</mn> </mrow> </semantics></math>; (<b>b</b>) FCM detection result when <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>0.019</mn> </mrow> </semantics></math>; (<b>c</b>) FCM detection result when <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>0.032</mn> </mrow> </semantics></math>; (<b>d</b>) FCM detection result when <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>0.05</mn> </mrow> </semantics></math>.</p> ">
Figure 11
<p>Detection results of the image shown in <a href="#remotesensing-11-00396-f009" class="html-fig">Figure 9</a>a: (<b>a</b>) Detection result of NACA (new adaptive cultural algorithm); (<b>b</b>) Detection result of QSFLA-NSM (quantum-inspired shuffled frog leaping algorithm combining the new search mechanism); (<b>c</b>) Detection result of ACA-IQPSO (adaptive cultural algorithm with improved quantum-behaved particle swarm optimization); (<b>d</b>) Detection result of QSFLA (quantum-inspired shuffled frog leaping algorithm); (<b>e</b>) Detection result of CPSO (cultural particle swarm optimization algorithm); (<b>f</b>) Detection result of QPSO (quantum-behaved particle swarm optimization).</p> ">
Figure 12
<p>The chart of the best fitness function values in <a href="#remotesensing-11-00396-t002" class="html-table">Table 2</a>.</p> ">
Figure 13
<p>Detection results of the underwater sonar image (image size: 112×117): (<b>a</b>) Original sonar image; (<b>b</b>) The denoising result of the denoising method proposed in this paper; (<b>c</b>) Detection result of NACA; (<b>d</b>) Detection result of QSFLA-NSM; (<b>e</b>) Detection result of ACA-IQPSO; (<b>f</b>) Detection result of QSFLA; (<b>g</b>) Detection result of CPSO; (<b>h</b>) Detection result of QPSO.</p> ">
Figure 14
<p>The chart of the best fitness function values in <a href="#remotesensing-11-00396-t003" class="html-table">Table 3</a>.</p> ">
Figure 15
<p>Detection result after the first iteration in the detection process of <a href="#remotesensing-11-00396-f011" class="html-fig">Figure 11</a>: (<b>a</b>) Detection result of NACA; (<b>b</b>) Detection result of QSFLA-NSM; (<b>c</b>) Detection result of ACA-IQPSO; (<b>d</b>) Detection result of QSFLA; (<b>e</b>) Detection result of CPSO; (<b>f</b>) Detection result of QPSO.</p> ">
Figure 16
<p>The distribution of particle positions in the new update strategy and old update strategy: (<b>a</b>) Position distribution of particles on Sphere function; (<b>b</b>) Position distribution of particles on Griewank function.</p> ">
Figure 17
<p>The result of the fitness value calculation in the new update strategy and old update strategy: (<b>a</b>) The optimization results of the Sphere function; (<b>b</b>) The optimization results of the Griewank function.</p> ">
Figure 18
<p>Detection results of the underwater sonar image (image size: 259×368): (<b>a</b>) Original sonar image; (<b>b</b>) The denoising result of the denoising method proposed in this paper; (<b>c</b>) Detection result of NACA.</p> ">
Figure 19
<p>Detection results of the underwater sonar image (image size: 259×368): (<b>a</b>) Original sonar image; (<b>b</b>) The denoising result of the denoising method proposed in this paper; (<b>c</b>) Detection result of NACA.</p> ">
Figure 20
<p>Detection results of the underwater sonar image (image size: 203×101): (<b>a</b>) Original sonar image; (<b>b</b>) The denoising result of the denoising method proposed in this paper; (<b>c</b>) Detection result of NACA.</p> ">
Versions Notes

Abstract

:
An adaptive approach is proposed to denoise and detect the underwater sonar image in this paper. Firstly, to improve the denoising performance of non-local spatial information in the underwater sonar image, an adaptive non-local spatial information denoising method based on the golden ratio is proposed. Then, a new adaptive cultural algorithm (NACA) is proposed to accurately and quickly complete the underwater sonar image detection in this paper. Concretely, NACA has two improvements. In the first place, to obtain better initial clustering centres, an adaptive initialization algorithm based on data field (AIA-DF) is proposed in this paper. Secondly, in the belief space of NACA, a new update strategy is adopted to update cultural individuals in terms of the quantum-inspired shuffled frog leaping algorithm (QSFLA). The experimental results show that the proposed denoising method in this paper can effectively remove relatively large and small filtering degree parameters and improve the denoising performance to some extent. Compared with other comparison algorithms, the proposed NACA can converge to the global optimal solution within small epochs and accurately complete the object detection, having better effectiveness and adaptability.

Graphical Abstract">

Graphical Abstract

1. Introduction

With the development of sonar imaging technology, the underwater sonar image detection technology has been extensively used in marine exploration, research, and investigation [1]. The underwater sonar image must be detected in different regions including object-highlight, sea-bottom-reverberation, and shadow before the underwater object can be recognized. [2]. With the help of the underwater sonar image detection technology, we can accurately detect the above regions and effectively retain the underwater sonar image details [3]. However, because of the complexity of the underwater environment, the underwater sonar image is easily affected by the reverberation effect, strong speckle noise, fuzzy edge, and weak texture information [4]. Therefore, before the underwater sonar image detection, the denoising algorithm can be adopted to maximize noise removal.
The spatial-based method has been widely used in image processing. For cartoon components in the image fusion method, a proper spatial-based method is presented for morphological structure preservation [5]. To mitigate the boundary seams produced by spatial domain methods, a non-local means filtering is introduced on the final decision maps to generate the fusion weight maps for each of the source images [6]. In the image denoising method, the spatial information mainly includes local spatial information and non-local spatial information. Compared with the disadvantage that local spatial information only uses a single neighborhood information, non-local spatial information can utilize multiple similar neighborhood information in subregions. Therefore, the non-local spatial information makes full use of the image, which is more conducive to image denoising. The filtering degree parameter is a significant parameter. It is intensively related to denoising results [7]. The relatively large filtering degree parameter will cause the loss of some image detail information, especially the image edge information. On the contrary, the relatively small filtering degree parameter will mean that the image noise cannot be effectively eliminated. Therefore, how to select the appropriate filtering degree parameter is a challenging problem. To select the appropriate filtering degree parameter, a selection method based on the image noise level was proposed. In this method, filtering degree parameters were a single value [8]. In order to get over the disadvantage that the single filtering degree parameters are sensitive to the noise, another selection method on the basis of the statistical characteristic in the search window was proposed [9]. However, when the search window is small, the denoising ability of the non-local spatial information will be greatly reduced. Furthermore, the idea of threshold was introduced to remove the abnormal parameters and select the appropriate filtering degree parameters [7]. Nevertheless, the two thresholds are fixed values and only refer to one underwater sonar image. Different underwater sonar images have different filtering parameter distribution characteristics. Therefore, in some cases, the relatively large and small filtering degree parameters cannot be effectively removed.
Various techniques are widely used in the underwater sonar image detection. The Fuzzy clustering algorithm is used for underwater sonar image detection [10]. However, this detection algorithm is susceptible to noise. The Markov algorithm is proposed for underwater sonar image detection [11,12]. Although the detection results were satisfactory, the processing procedure was quite complicated and computationally costly. Later, a promising and multiphase detection framework related to the level set was developed [13], the key idea of which was to minimize the energy with the local mean. However, this framework was not suitable for the detection of underwater sonar images in the presence of speckle noise. To obtain better detection results, a large number of improved algorithms were proposed. Active contours were introduced in the level set method [14,15]. In this case, Markov random field and the level set are combined to extract the features from the underwater sonar image [2]. The morphological top-hat and bottom-hat transformation were used in the level set [16]. It is conducive to the extraction of underwater objects. Subsequently, an adaptive initialization narrowband Chan-Vese model is proposed and applied to underwater sonar image detection [1].
Over the years, as a framework to solve complex problems, cultural algorithm (CA) is attracting more and more attention [17]. Meanwhile, the intelligent optimization algorithm is also widely used in image processing [18,19] and has achieved good results. In order to further enhance the search ability, many intelligent optimization algorithms are used to improve the cultural algorithm. For accurate selection of business partners, a cultural particle swarm optimization algorithm (CPSO) was proposed can get better options [20]. Quantum-behaved particle swarm optimization (QPSO) was introduced into CA to figure out multi-objective optimization problems [21], which has high efficiency. Then, an adaptive cultural algorithm with improved quantum-behaved particle swarm optimization (ACA-IQPSO) was proposed to solve the problem of the underwater sonar image detection [22], which accurately completes underwater object detection by searching clustering centres. However, because the clustering centres are randomly initialized in the solution space, the situation when the clustering centres are very close can occur to some extent. Therefore, the convergence speed of the algorithm is too slow to complete the underwater sonar image detection accurately within small epochs.
A clustering algorithm is proposed to calculate the clustering centres by the potential entropy of the data field [23]. The comparative experiments have shown that this algorithm can get better clustering results. Nevertheless, it does not take the characteristics of underwater sonar images into account when calculating potential entropy, which is not conducive to the initialization of the clustering centres in the underwater sonar images.
In order to reduce the complexity of the algorithm and enhance the search ability in the process of detection, the idea of the shuffled frog leaping algorithm (SFLA) is used as the belief space update strategy in ACA-IQPSO. However, SFLA is easy to get into the local optimal solution. To make the population more diverse and improve the search ability, Quantum computing theory is used to improve SFLA [24]. A quantum-inspired shuffled frog leaping algorithm combining the new search mechanism (QSFLA-NSM) is proposed to the underwater sonar image detection [7]. The experimental results of QSFLA-NSM show that it can further enhance population diversity and search ability.
To overcome the drawback of two thresholds in work [7], in this work we proposed an adaptive non-local spatial information denoising method based on the golden ratio. Two thresholds are set adaptively according to the golden radio. Then, NACA is proposed in this paper to accurately detect underwater sonar image. The improvement of the detection algorithm focuses on two aspects: Firstly, inspired by the idea of data field and entropy [23], the AIA-DF is adopted to more accurately calculate the initialization of the clustering centres in this paper to enhance convergence efficiency. The AIA-DF can automatically extract the optimal value of threshold by using the potential entropy of data field from the underwater sonar images dataset. Subsequently, the threshold from the data field can adaptively initialize the clustering centres. Secondly, in the belief space, to improve the limitations of the update strategy, a new update strategy is used to update the cultural individuals in terms of QSFLA.

2. Non-Local Spatial Information Denoising Method

2.1. Non-Local Spatial Information

The underwater sonar image is represented as X = { x 1 , x 2 , x n } , n is the total number of pixels. x i is the i th pixel. x i ¯ is the non-local spatial information of x i and can be expressed as:
x i ¯ = p W i r w i p x p
where W i r represents the search window, its centre is x i and the radius is r . w i p and z i can be calculated by:
w i p = 1 Z i exp ( x ( N i ) x ( N p ) 2 , ρ 2 / h )
Z i = p W i r exp ( x ( N i ) x ( N p ) 2 , ρ 2 / h )
where N i is the neighborhood window, its centre is the x i and the radius is s . x ( N i ) is the vector and represents all the pixels in N i . h is the filtering degree parameter. x ( N i ) x ( N p ) 2 , ρ 2 is defined as:
x ( N i ) x ( N p ) 2 , ρ 2 = q = 1 ( 2 s + 1 ) 2 ρ ( q ) ( x ( q ) ( N i ) x ( q ) ( N p ) ) 2
where x ( q ) ( N i ) is the q th pixel in x ( N i ) .
ρ ( q ) can be presented as:
ρ ( q ) = t = max ( d , 1 ) s 1 ( 2 t + 1 ) 2 s
d = max ( | y q s 1 | , | z q s 1 | )
where y q = mod ( q , ( 2 s + 1 ) ) , z q = f l o o r ( q , ( 2 s + 1 ) ) + 1 . ( y q , z q ) represents the coordinates of the q th dimension in N i .

2.2. Our Proposed Denoising Method—An Adaptive Non-Local Spatial Information Denoising Algorithm Based on the Golden Ratio

h is a significant parameter. Its value can powerfully impact the weight of neighborhood configuration in Equation (2) and Equation (3). Therefore, h is closely related to the denoising result in the non-local spatial information denoising method. In the process of denoising, too large h will cause missing the image detail information. Meanwhile, too small h will reduce the denoising performance of non-local spatial information.
In our previous work [7], to select the appropriate h , two thresholds h m i n = 0.01 and h m a x = 0.05 were set to remove the inappropriate h . Experimental results showed that this method can effectively remove noise and had certain adaptability. However, it is found that different underwater sonar images have different h distribution characteristics. The two thresholds h m i n = 0.01 and h m a x = 0.05 are determined only based on one sonar image. In other words, we use the same thresholds for different underwater sonar images denoising. It obviously has limitations.
Although the distribution characteristics of h in the different underwater sonar images are different, further experiments show that the proportion of appropriate h in the total h is approximately the same, which satisfies the golden ratio within a certain error. Based on this discovery, in this paper, we propose an adaptive non-local spatial information denoising method based on the golden ratio which realizes two threshold adaptive settings and removes the underwater sonar image noise more effectively. The statistical law that the proportion satisfies is investigated on 20 different underwater sonar images by calculating the proportion of the h within two thresholds in the total h . Owing to space constraints, we only show partial results from 20 sonar images. Figure 1, Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6 show the estimated values of the optimal thresholds based on different underwater sonar images.
As seen from Figure 1b, Figure 2b, Figure 3b, Figure 4b, Figure 5b, and Figure 6b, different underwater sonar images have different distribution characteristics of h . This means that the estimated optimal thresholds are different. A linear fitting method is used to calculate the proportion satisfies based on different underwater sonar images.
For each underwater sonar image, the set of total h is defined as H = { h i | i N U } , h i is the filtering degree parameter corresponding to the i th pixel, N U is the number of the total pixels on the underwater sonar image. The set of h within two thresholds is defined as H = { h i | i N U h i h m i n h i h m a x } . The ratio of the cardinality of H to the cardinality of H is defined as:
η = c a r d   H c a r d   H
Table 1 shows the estimated optimal threshold and η of 20 different underwater sonar images.
Figure 7 shows the result of the linear fitting based on Table 1.
The fitting curve is y = 7.14 e 06 x + 0.613 in Figure 7. It demonstrates that η is approximately distributed in the vicinity of 0.613. The statistical law that the proportion satisfies is the golden radio within the range of error permitting. Two thresholds are set adaptively according to the golden ratio to remove too large and small h . On this basis, an adaptive non-local spatial information denoising method based on the golden ratio is proposed in this paper to effectively improve the denoising performance.
This is set that h less than h m i n is defined as H 1 = { h i | i N U h i h m i n } and that h larger than hmax is defined as H 1 = { h i | i N U h i h m a x } . On this basis, the ratio of the cardinality of H 1 to the cardinality of H is defined as η 1 = c a r d   H 1 c a r d   H . The ratio of the cardinality of H 2 to the cardinality of H is defined as η 2 = c a r d   H 2 c a r d   H . Comparing η 1 and η 2 with 1 η 2 , when the value exceeds the range of error permitting, the two thresholds are constantly updated by the dichotomy. The update process is not finished until the experiment result simultaneously satisfies | η 2 1 η 2 | ε and | η 1 1 η 2 | ε . ε is the range of error permitting. The result of two adaptive thresholds is respectively defined as A h m i n and A h m i n .
For the i th pixel in the sonar image, h i is defined as:
h i = { A h m i n h m i A h m i n A h m a x h m i A h m a x h m i o t h e r w i s e
h m i = mean p W i r { x ( N j ) x ( N p ) 2 , ρ 2 }
The adaptive non-local spatial information denoising method based on the golden ratio is described in Algorithm 1.
Algorithm 1: the proposed denoising method
Input: Original sonar image;
Output: Denoising image;
Initialization: r (the radius of search window), s (the radius of neighbourhood window);
Procedure:
1: For q = 1 to s 2 do:
2:  Calculate d using Equation (6);
3:     Calculate ρ ( q ) using Equation (5);
4: For every search window in the input image, do:
5:     For any two different neighbourhood windows in the search window, do:
6:      Calculate the weighted Euclidean distance using Equation (4);
7:     Calculate the filtering degree parameter h using Equation (9);
8: Calculate A h m i n and A h m a x when the result of Equation (8) satisfies the golden ratio;
9: For every search window in the input image, do:
10:   Calculate the filtering degree parameter h using Equation (8);
11:   Calculate Z using Equation (3);
12:   Calculate w using Equation (2);
13:   Calculate non-local spatial information x ¯ using Equation (1).

3. Culture Algorithm

3.1. Cultural Algorithm Framework

The bilevel evolutionary mechanism in the cultural algorithm mainly includes three elements: the population space, the belief space, and the communication protocol. The communication protocol consists of influence function and accept function. The population space and the belief space are two relatively independent evolutionary and update processes. However, they rely on the influence function and accept function to communicate with each other, manage relevant information, and guide the evolution and update of space. The important feature of the cultural algorithm is the introduction of belief space. The individual evolutionary experience in the population space can be transferred to the belief space by the accept function. The belief space transforms advanced experience into knowledge through certain rules. The knowledge in the belief space guides the evolution of the population space towards a more accurate direction through the influence function. Figure 8 shows the schematic diagram of CA.

3.2. Our Proposed Cultural Algorithm—A New Adaptive Culture Algorithm

In our previous work [22], we proposed an adaptive cultural algorithm with improved quantum-behaved particle swarm optimization (ACA-IQPSO) and applied it to the underwater sonar image detection. In the population space, to enhance the search ability of the ACA-IQPSO, we introduced the IQPSO as the evolution strategy. In the belief space, we adopted a new update strategy based on the idea of SFLA. Situational knowledge, normative knowledge, and domain knowledge were selected to form the knowledge structure. A new communication mechanism was designed to make better use of the information in the population space and belief space. Experimental results show that ACA-IQPSO has better searching ability and convergence efficiency.
However, we found that ACA-IQPSO could not obtain good initial clustering centres, and the update strategy of the belief space easily gets into the local optimal solution. Therefore, the convergence ability of ACA-IQPSO needs to be improved and the efficiency of underwater sonar image detection is reduced.
Therefore, in order to complete the underwater sonar image detection more accurately in shorter iterations, NACA is proposed for the underwater sonar image detection. NACA has two improvements over ACA-IQPSO. To begin with, inspired by the idea of data field and maximum entropy [23], we propose an adaptive initialization algorithm based on data field (AIA-DF) for underwater sonar image detection. AIA-DF is used in the population space to more accurately find the initial clustering centres. Secondly, in our previous work [7], we proposed QSFLA-NSM for the underwater sonar image detection. Experimental results show that compared with SFLA, QSFLA-NSM has stronger search ability. Therefore, we propose a new update strategy in the belief space based on the idea of QSFLA-NSM.

3.2.1. Adaptive Initialization Algorithm Based on Data Field

The dataset R is defined as R = { R 1 , R 2 , , R N } in the data space Ω. The ith data point is defined as Ri = (Ri1, Ri2, ⋯, RiD), N is the size of R. D is the dimension of a data point. Every data point in R interacts with the data points around it. In this way, a data field is formed according to the interaction of different data points. The interaction between different data points is measured by potential functions in the data field. Generally, the closer data points have stronger interaction and greater potential. In contrast, the further data points have weaker interaction and lower potential. The potential function is defined as:
φ R i 1 R i 2 = m R i e R i R * 2 σ 2
where φ R i 1 R i 2 is the potential of the data point Ri = (Ri1, Ri2) at the threshold R * , Ri1 is the grey value of the ith pixel in the underwater sonar image. Ri2 is the average of the grey values of pixels in the neighborhood window. The neighborhood window centre is the ith pixel. m R i is the frequency corresponding to the ith pixel in the two-dimensional histogram of the underwater sonar image. R i R * is the Euclidean distance and σ is the impact factor.
Furthermore, in the data field, the sum of the potential of all data points at the threshold R * is defined as:
φ = i = 1 N m R i e R i R * 2 σ 2
where N is the number of total pixels in the underwater sonar images.
The influence radius of every data point is 3 σ 2 in the data field [23]. Therefore, when calculating φ , the data points at work are in the region O, of which its centre is R * and the radius is 3 σ 2 . On this basis, φ is redefined as
φ = R i O m R i e R i R * 2 σ 2
The impact factor σ is a significant parameter. It is closely related to φ . The relatively large σ means that φ is relatively close, which makes the initialization algorithm easy to fall into the local optimal solution during the iteration. In contrast, the relatively small σ means that φ R * is relatively close to 0, which also affects the subsequent solution of optimal threshold. To accurately calculate the optimal threshold, σ is redesigned.
According to the criterion:
P ( μ 3 σ K μ + 3 σ ) = 0.997
where K is the possible grey levels in the underwater sonar image and μ = L 2 . L is grey levels in the underwater sonar image. σ is defined as:
σ = 6 L
According to the principle of the data field [23], when the uncertainty is the smallest, the potential can perfectly reflect the distribution of data points in Ω. The uncertainty is usually described as the entropy. The optimal threshold is calculated by entropy. The size of the data field is L × L, the scope of R * is 0 R * 1 L 1 , 0 R * 2 L 1 . Data points in the data field are divided into two classes: A = { ( 0 , 0 ) , , ( R * 1 , R * 2 ) } and B = { ( R * 1 + 1 , R * 2 + 1 ) , , ( L 1 , L 1 ) } . The entropy of A and B are respectively expressed as follows:
H ( A ) = i 1 = 0 R * 1 i 2 = 0 R * 2 φ R i 1 R i 2 Ψ A log φ R i 1 R i 2 Ψ A
H ( B ) = i 1 = R * 1 + 1 L 1 i 2 = R * 2 + 1 L 1 φ R i 1 R i 2 Ψ B log φ R i 1 R i 2 Ψ B
where Ψ A = i 1 = 0 R * 1 i 2 = 0 R * 2 φ R i 1 R i 2 , Ψ B = i 1 = R * 1 + 1 L 1 i 2 = R * 2 + 1 L 1 φ R i 1 R i 2 .
The sum of posterior entropy in the data field is defined as:
H ( R * ) = H ( A ) + H ( B )
When the maximum of H ( R * ) is obtained, R * is the optimal threshold.
AIA-DF is proposed in this paper according to Equations (12)–(17). The initialization algorithm adaptively calculates initial clustering centres in population space on the basis of the characteristics of the underwater sonar images. It improves the disadvantages that the initial cluster centers are randomly generated in the solution space and effectively avoids the prematurity of detection algorithm.
The influence region O of the optimal threshold is limited because of the influence radius. Besides, in the data field, the distribution of potentials has different characteristics in different underwater sonar images. Therefore, the optimal threshold is also different to some extent.
U = { ( R * 1 , R * 2 ) | 0 R * 1 L 1 , 0 R * 2 L 1 } is the initial scope of threshold R * . When the optimal threshold is obtained, the grey value is randomly selected in region O for initializing the clustering centres. The update formula for U is given as:
U ( n u m + 1 ) = U ( n u m ) U *
where n u m is the number of optimal thresholds currently obtained.
U* is defined as::
U * = { ( R ( n u m + 1 ) * 1 , ( R ( n u m + 1 ) * 2 | R n u m * 1 3 σ 2 R ( n u m + 1 ) * 1 R n u m * 1 + 3 σ 2 , R n u m * 2 3 σ 2 R ( n u m + 1 ) * 2 R n u m * 2 + 3 σ 2 }
When the value of R ( n u m + 1 ) * 1 and R ( n u m + 1 ) * 2 are less than 0 or larger than L − 1, the update process of U is over. AIA-DF is described in Algorithm 2.
Algorithm 2: the proposed AIA-DF
Input: Denoising image;
Output: Initialized image;
Initialization: s (the radius of the neighbourhood window), n u m (the number of the optimal threshold), U (the set of optimal threshold potential solutions);
Procedure:
1: For every pixel in the input image, do:
2:     Calculate the mean grey of its neighbourhood;
3:  Establish the two-dimensional grey histogram of the input image;
4: Establish the data field according to the grey histogram;
5: repeat:
6:      Calculate the entropy of H ( A ) using Equation (15);
7:      Calculate the entropy of H ( B ) using Equation (16);
8:      Calculate the optimal threshold when the maximum of Equation (17) is obtained;
9:    Calculate U * using Equation (19);
10:  Update U using Equation (18):
11:   Update n u m using n u m = n u m + 1 :
12: Until: the stopping criteria are met.

3.2.2. New Update Strategy

In the belief space, to complete the update process by using the new update strategy in terms of QSFLA-NSM. In the first place, cultural individuals are divided into the multiple subpopulations and sorted from large to small according to the fitness function value [25]. Each subpopulation executes the local search strategy. In the second place, with the subpopulation complete local search strategy, all subpopulations are mixed for global information exchange. Then, the belief space is divided into subpopulations again. Local search and global information exchange will alternate until the number of iterations is satisfied
Local search completes the update of the worst cultural individual in each subpopulation. The first step of the local search [5] is expressed as:
n e w Y w ( t + 1 ) = { P ( t ) + β | S b ( t ) Y w ( t ) | ln ( 1 μ 1 )     μ 1 < 0.5 P ( t ) β | S b ( t ) Y w ( t ) | ln ( 1 μ 1 )     otherwise
where t is the current local iterative times, μ1 is a random number in [0,1], P(t) is the local attractor, and β is the contraction–expansion coefficient. S b is the situational knowledge of the subpopulation on behalf of the local best cultural individual. Y w is the local worst cultural individual.
Comparing the fitness function values of n e w Y w and Y w , if the fitness function value of n e w Y w is larger, Y w is replaced by n e w Y w . Whereas, S g which represents the global best individual, will replace S b in Equation (20). n e w Y w is calculated by Equation (20) again. If the fitness function value of n e w Y w is still less than Y w , it will randomly generate an individual from the solution space and replace Y w . When the upper limit of the local iteration number is reached, the local search is over.

4. Experiments and Discussion

4.1. The Characteristics of the Underwater Sonar Images

The underwater sonar images include the object-highlight region, sea-bottom-reverberation region, and shadow region. However, because of the complexity of the underwater environment, the underwater sonar image is easily affected by the reverberation effect, strong speckle noise, fuzzy edge, and weak texture information. This means that the underwater sonar images contain a lot of noise.

4.2. The Effectiveness Verification of the Proposed Denoising Method

To prove the superiority of the proposed denoising method, Figure 9 shows the denoising results based on Figure 1a of the proposed denoising method and the previous denoising method [7]. The experiment environment is using Matlab R2012b with a 2.7 GHz Core processor and 8 GB of RAM.
As can be seen from Figure 9b,d, compared with the previous thresholds h m i n = 0.01 , the adaptive thresholds A h m i n = 0.019 can remove more fairly small filtering degree parameters and effectively remove image noise points. Meanwhile, the threshold A h m a x = 0.032 is smaller than h m a x = 0.05 which can remove more fairly large filtering degree parameters. It contributes to keeping the underwater sonar image details better when the image noise points are removed.
A relatively simple and effective the fuzzy c-means (FCM) [10] is used to further qualitatively verify the effectiveness of the denoising method proposed in this paper. Figure 10 shows the detection results of FCM based on the thresholds in Figure 9b,d.
As can be seen from Figure 10a,b, compared with the previous threshold h m i n , A h m i n removes smaller parameters in a larger range. Therefore, it can effectively improve the performance of sonar images. Meanwhile, as shown in Figure 10c,d, compared with the previous threshold h m a x , A h m a x can more accurately define the boundaries of relatively large h , which benefits preserving detail information in sonar images. The noise can be effectively removed on the basis of preserving image details by the proposed denoising method which contributes to the remainder of the image processing.

4.3. The Effectiveness Demonstration of the Proposed Detection Method

To testify the effectiveness of NACA and AIA-DF, the detection experiments based on different underwater sonar images are presented. The proposed non-local spatial information denoising method based on the golden ratio in this paper is used for removing noise in the underwater sonar images. On that basis, the proposed NACA is used for underwater sonar image detection. Meanwhile, the proposed NACA is compared with QSFLA-NSM [7], ACA-IQPSO [22], QSFLA [26,27], CPSO [20], and QPSO [28]. Except that the number of clustering centres is adaptively set according to AIA-DF in NACA, the number of clustering centres is 3 in other intelligent optimization algorithms. The population size is 20, the sub-populations size is 5 in proposed NACA, QSFLA-NSM, and QSFLA, the size of the belief space is 8 in NACA, ACA-IQPSO, and CPSO, the global maximum number of iterations is 10, the local maximum number of iterations is 2 in NACA, QSFLA-NSM, and QSFLA. The experiment environment is using Matlab R2012b with a 2.7 GHz Core processor and 8 GB of RAM.
To demonstrate the effectiveness of NACA, Figure 11 shows the detection results of the underwater sonar image shown in Figure 9a.
As is described in Figure 11, the regions of object-highlight and shadow are more accurately detected by NACA in Figure 11a. It is found that many noise points cannot be removed in Figure 11b,d. In Figure 11e,f, the detection results have poor integrity in the regions of object-highlight and shadow, respectively. It will cause serious difficulty to the subsequent underwater sonar image recognition. Therefore, compared with the detection results of QSFLA-NSM, ACA-IQPSO, QSFLA-NSM, QSFLA, CPSO, and QPSO, the proposed NACA can effectively remove the noise and accurately complete underwater sonar image detection.
Subsequently, the quantitative analysis of the detection results in Figure 12 is performed. The values of fitness function [25] after each iteration in the detection process of NACA and contrast algorithms is shown in Table 2.
Figure 12 shows the chart of the best fitness function values in Table 2.
As can be seen from Table 2 and Figure 12, the best fitness values of the proposed NACA is the largest after each iteration. Only NACA converges after four iterations. The results show that the proposed NACA has a better search ability and a faster convergence speed (the speed corresponds to the number of iterations). Meanwhile, the fitness function values after the first iteration are close to the fitness function values after 10 iterations in Figure 12. This means that AIA-DF can obtain better initial clustering centres in the proposed NACA. Through qualitative analysis and quantitative analysis of the underwater sonar detection results, the effectiveness of the proposed NACA is verified.
To further verify the effectiveness of the proposed NACA, Figure 13 shows the detection results of the underwater sonar image.
The values of fitness function [25] after each iteration in the detection process of NACA and contrast algorithms are shown Table 3. Figure 14 shows the chart of the best fitness function values in Table 3.
As can be seen in Figure 13, NACA more accurately completes the detection of different regions in Figure 13c. Figure 13d,f show over detection in the shadow region. The detection result fails to effectively remove the noise in Figure 13e,g,h. At the same time, the best fitness values of NACA are also the largest after each iteration in Table 3 and Figure 14. Comparing with other algorithms, the detection results demonstrate that NACA has the better search ability and a faster convergence speed (the speed corresponds to the number of iterations) based on the better initialization results that were obtained by AIA-DF in the paper. Therefore, through qualitative and quantitative analysis, the proposed NACA can effectively remove the noise and better complete the process of underwater sonar image detection.
To demonstrate the effectiveness of AIA-DF, Figure 15 shows the detection result after the first iteration in the detection process of Figure 11. The fitness function value after iteration is shown in the first line of Table 2.
As can be seen in Figure 15, only NACA can accurately detect the region of the object highlight region and shadow. It also effectively removes the noise in the underwater sonar image in Figure 15a. These results show that AIA-DF can find more accurate clustering centres. In addition, it can be seen from Figure 15a and Figure 11a that the detection result after the first iteration is almost identical with the final detection result. The experimental results demonstrate that AIA-DF contributes to improving convergence speed (the speed corresponds to the number of iterations), accuracy, and search ability in the proposed NACA. Through qualitative analysis and quantitative analysis of the underwater sonar image detection results after the first iteration, the effectiveness of AIA-DF is proved.

4.4. The Performance Analysis of the New Update Strategy in NACA

To demonstrate the superiority in the search ability of the new update strategy in NACA, in the benchmark functions, Sphere function and Griewank function are used to test the position distribution of particles in this paper. Sphere function is unimodal and only has one global optimal solution. Griewank function is multimodal and has many local optimal solutions, however only has one global optimal solution. Figure 16 shows the distribution of particle positions in the new update strategy and old update strategy [22]. The relevant parameters are as follows: the dimension of the solution space is 2, the size of the belief space is 30, the size of the subspaces is 5, the global maximum number of iterations is 10, and the local maximum number of iterations is 2.
It can be seen in Figure 16 that compared with the old update strategy, the distribution of particle positions in the new update strategy are more dispersed. The new update strategy is easier to obtain the global optimal solution and has a stronger search ability in this paper.
In order to further verify the effectiveness of the new update strategy in its search ability, the fitness function values are calculated by Sphere function and Griewank function in the new update strategy and old update strategy, and the result of the fitness function values are shown in Figure 17. The relevant parameters are as follows: the dimension of the solution space is 30, the size of the belief space is 10, the size of the subspaces is 2, the global maximum number of iterations is 100, and the local maximum number of iterations is 10.
As can be seen in Figure 17, the fitness function values of the old update strategy are always greater than that of the new update strategy in each iteration. Compared with the old update strategy, the new update strategy is more likely to jump out of the local optimization. The experiment results further demonstrate that the search ability of the new update strategy is better than the old update strategy.

4.5. The Adaptability Demonstration of the Proposed Denoising Method and Detection Method

To prove the adaptability of the proposed denoising method and NACA, Figure 18 shows the detection results of the structured seabed which is an object in sand ripples. Figure 19 shows the detection results of the starboard original sonar image with a ship. Figure 20 shows the detection results of a floating object. The experiment environment uses Matlab R2012b with a 2.7 GHz Core processor and 8 GB of RAM.
It can be seen from the detection results in Figure 18b, Figure 19b, and Figure 20b that the proposed denoising method can remove noise effectively in polytype underwater sonar images. The boundary information of the different regions is accurately detected in Figure 18c, Figure 19c, and Figure 20c. Moreover, the iteration times of the algorithm are only two times, which further proves the adaptability of the AIA-DF in this paper. Therefore, the adaptability of the proposed method has been effectively proved.

5. Conclusions

This paper proposed an adaptive approach to denoise and detect the underwater sonar image. The problem of the inappropriate filtering degree parameter means that it seriously affects the denoising performance in underwater sonar images, which is solved by the proposed adaptive non-local spatial information denoising method based on the golden ratio in the paper. NACA is proposed in this paper. Firstly, in the population space, AIA-DF is adopted to obtain good initial clustering centres by calculating the potential entropy of the data field from the underwater sonar images dataset. In the belief space, a new update strategy is adopted to update cultural individuals according to the idea of QSFLA. The new update strategy further improves the search ability of NACA.
We apply the proposed method to the different underwater sonar images. The experimental results demonstrate that the proposed denoising method can effectively remove noise and reduce the difficulty of the following underwater sonar image recognition. Compared to the comparison algorithms, the proposed NACA has more advantages in its search ability and convergence speed (the speed corresponds to the number of iterations). AIA-DF can better locate initial clustering centres and enhance convergence efficiency of NACA which contributes to accurately complete underwater sonar image detection and convergence to the global optimal solution within small epochs. This paper provides an effective and important method for underwater sonar image detection. It has much academic and practical importance.

Author Contributions

Conceptualization, X.W., J.Y., and Q.L.; Methodology, X.H. and W.H.; Software, W.H.; Investigation, X.H. and W.H.; Resources, X.W.; Writing—original draft preparation, Q.L. and W.H.; Writing—review and editing, X.W., Q.L., J.Y., and X.H.; Visualization, X.H.; Supervision, X.W. and J.Y.; Project administration, X.W.; Funding acquisition, X.W.

Funding

This research was supported in part by the National Natural Science Foundation of China grant number 41876110, the Fundamental Research Funds for the Central Universities of China grant number HEUCF180601, the Heilongjiang Province Outstanding Youth Science Fund of China grant number JC2017017, and in part by the Fok Ying-Tong Education Foundation of China grant number 151007.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, X.; Guo, L.; Yin, J.; Liu, Z.; Han, X. Narrowband Chan-Vese model of sonar image segmentation: An adaptive ladder initialization approach. Appl. Acoust. 2016, 113, 238–254. [Google Scholar] [CrossRef]
  2. Ye, X.; Zhang, Z.; Liu, P.; Guan, H. Sonar image segmentation based on GMRF and Level-set models. Ocean Eng. 2010, 37, 891–901. [Google Scholar] [CrossRef]
  3. Wang, X.; Liu, S.; Teng, X.; Sun, J.; Jiao, J. SFLA with PSO local search for detection sonar image. In Proceedings of the 2016 35th Chinese Control Conference, Chengdu, China, 27–29 July 2016; pp. 3852–3857. [Google Scholar]
  4. Wu, J.P.; Guo, H. A method for sonar image segmentation based on combination of MRF and region growing. In Proceedings of the 2015 5th International Conference on Communication Systems and Network Technologies, Gwalior, India, 4–6 April 2015; pp. 457–460. [Google Scholar]
  5. Zhu, Z.; Yin, H.; Chai, Y.; Li, Y.; Qi, G. A novel multi-modality image fusion method based on image decomposition and sparse representation. Inf. Sci. 2018, 432, 516–529. [Google Scholar] [CrossRef]
  6. Li, H.; Qiu, H.; Yu, Z.; Li, B. Multifocus image fusion via fixed window technique of multiscale images and non-local means filtering. Signal Process. 2017, 138, 71–85. [Google Scholar] [CrossRef]
  7. Wang, X.; Liu, S.; Liu, Z. Underwater sonar image detection: A combination of non-local spatial information and quantum-inspired shuffled frog leaping algorithm. PLoS ONE 2017, 12, e0177666. [Google Scholar] [CrossRef]
  8. Zhao, F.; Jiao, L.; Liu, H.; Gao, X. A novel fuzzy clustering algorithm with non-local adaptive spatial constraint for image segmentation. Signal Process. 2011, 91, 988–999. [Google Scholar] [CrossRef]
  9. Zhao, F. Fuzzy clustering algorithms with self-tuning non-local spatial information for image segmentation. Neurocomputing 2013, 106, 115–125. [Google Scholar] [CrossRef]
  10. Wang, L. Segmentation algorithm of fuzzy clustering on side scan sonar image. J. Huazhong Univ. Sci. Technol. 2012, 40, 25–29. [Google Scholar]
  11. Mignotte, M.; Collect, C.; Perez, P.; Bouthemy, P. Three-class markovian segmentation of high-resolution sonar image. Comput. Vis. Image Underst. 1999, 76, 191–204. [Google Scholar] [CrossRef]
  12. Mignotte, M.; Collect, C.; Perez, P.; Bouthemy, P. Sonar image segmentation using an unsupervised hierarchical MRF model. IEEE Trans. Signal Process. 2000, 9, 1216–1231. [Google Scholar] [CrossRef]
  13. Vese, L.; Chan, T. A multiphase level set framework for image segmentation using the mumford and shah model. Int. J. Comput. Vis. 2002, 50, 271–293. [Google Scholar] [CrossRef]
  14. Lianantonakis, M.; Petillot, Y. Sidescan sonar segmentation using active contours and level set methods. In Proceedings of the Oceans Europe 2005, Brest, France, France, 20–23 June 2005; pp. 719–724. [Google Scholar]
  15. Lianantonakis, M.; Petillot, Y. Sidescan sonar segmentation using texture descriptors and active contours. IEEE J. Ocean. Eng. 2007, 32, 744–752. [Google Scholar] [CrossRef]
  16. Liu, G.; Bian, H.; Shi, H. Sonar image segmentation based on an improved level set method. Phys. Procedia 2012, 33, 1168–1175. [Google Scholar] [CrossRef]
  17. Awad, N.; Ali, M.; Suganthan, P.; Reynolds, R. CADE: A Hybridization of cultural algorithm and differential evolution for numerical optimization. Inf. Sci. 2017, 378, 215–241. [Google Scholar] [CrossRef]
  18. Khatami, A.; Mirghasemi, S.; Khosravi, A.; Lim, C.P.; Nahavandi, S. A new PSO-based approach to fire flame detection using K-Medoids clustering. Expert Syst. Appl. 2017, 68, 69–80. [Google Scholar] [CrossRef]
  19. Morra, L.; Coccia, N.; Cerquitelli, T. Optimization of computer aided detection systems: An evolutionary approach. Expert Syst. Appl. 2018, 100, 45–156. [Google Scholar] [CrossRef]
  20. Wei, Z.; Bu, Y. Cultural particle swarm optimization algorithm and its application. In Proceedings of the 2012 24th Chinese Control and Decision Conference, Taiyuan, China, 23–25 May 2012; pp. 740–744. [Google Scholar]
  21. Liu, T.; Jiao, L.; Ma, W.; Ma, J.; Shang, R. A new quantum-behaved particle swarm optimization based on cultural evolution mechanism for multiobjective problems. Knowl. Based Syst. 2016, 101, 90–99. [Google Scholar] [CrossRef]
  22. Wang, X.; Hao, W.; Li, Q. An adaptive cultural algorithm with improved quantum-behaved particle swarm optimization for sonar image detection. Sci. Rep. 2017, 7, 17733. [Google Scholar] [CrossRef]
  23. Wang, S.; Wang, D.; Li, C.; Li, Y.; Ding, G. Clustering by fast search and find of density peaks with data field. Chin. J. Electron. 2016, 25, 397–402. [Google Scholar] [CrossRef]
  24. Zhang, Y.; Lu, K.; Gao, Y. Quantum algorithms and quantum-inspired algorithms. Chin. J. Comput. 2013, 36, 1835–1842. [Google Scholar] [CrossRef]
  25. Wang, X.; Liu, S.; Li, Q.; Liu, Z. Underwater sonar image detection: A novel quantum-inspired shuffled frog leaping algorithm. Chin. J. Electron. 2018, 27, 588–594. [Google Scholar] [CrossRef]
  26. Ding, W.; Wang, J.; Guan, Z.; Shi, Q. Enhanced minimum attribute reduction based on quantum-inspired shuffled frog leaping algorithm. J. Syst. Eng. Electron. 2013, 24, 426–434. [Google Scholar] [CrossRef]
  27. Ding, W.; Wang, J.; Guan, Z. A minimum attribute self-adaptive cooperative co-evolutionary reduction algorithm based on quantum elitist frogs. J. Comput. Res. Dev. 2014, 51, 743–753. [Google Scholar]
  28. Zhang, B.; Qi, H.; Sun, S.; Ruan, L.; Tan, H. Solving inverse problems of radiative heat transfer and phase change in semitransparent medium by using improved quantum particle swarm optimization. Int. J. Heat Mass Transf. 2015, 85, 300–310. [Google Scholar] [CrossRef]
Figure 1. The estimated values of the optimal threshold (image size: 277 × 325): (a) Original sonar image; (b) The estimated values of the optimal threshold.
Figure 1. The estimated values of the optimal threshold (image size: 277 × 325): (a) Original sonar image; (b) The estimated values of the optimal threshold.
Remotesensing 11 00396 g001
Figure 2. The estimated values of the optimal threshold (image size: 173 × 167): (a) Original sonar image; (b) The estimated values of the optimal threshold.
Figure 2. The estimated values of the optimal threshold (image size: 173 × 167): (a) Original sonar image; (b) The estimated values of the optimal threshold.
Remotesensing 11 00396 g002
Figure 3. The estimated values of the optimal threshold (image size: 197 × 211): (a) Original sonar image; (b) The estimated values of the optimal threshold.
Figure 3. The estimated values of the optimal threshold (image size: 197 × 211): (a) Original sonar image; (b) The estimated values of the optimal threshold.
Remotesensing 11 00396 g003
Figure 4. The estimated values of the optimal threshold (image size: 205 × 201): (a) Original sonar image; (b) The estimated values of the optimal threshold.
Figure 4. The estimated values of the optimal threshold (image size: 205 × 201): (a) Original sonar image; (b) The estimated values of the optimal threshold.
Remotesensing 11 00396 g004
Figure 5. The estimated values of the optimal threshold (image size: 130 × 106): (a) Original sonar image; (b) The estimated values of the optimal threshold.
Figure 5. The estimated values of the optimal threshold (image size: 130 × 106): (a) Original sonar image; (b) The estimated values of the optimal threshold.
Remotesensing 11 00396 g005
Figure 6. The estimated values of the optimal threshold (image size: 146 × 135): (a) Original sonar image; (b) The estimated values of the optimal threshold.
Figure 6. The estimated values of the optimal threshold (image size: 146 × 135): (a) Original sonar image; (b) The estimated values of the optimal threshold.
Remotesensing 11 00396 g006
Figure 7. The result of the linear fitting based on Table 1.
Figure 7. The result of the linear fitting based on Table 1.
Remotesensing 11 00396 g007
Figure 8. Schematic diagram of CA (cultural algorithm).
Figure 8. Schematic diagram of CA (cultural algorithm).
Remotesensing 11 00396 g008
Figure 9. The denoising results based on Figure 1a of the proposed denoising method in this paper and the previous denoising method: (a) The denoising result of the denoising method proposed in this paper; (b) The result of filtering the degree parameter is selected by two adaptive thresholds in the proposed denoising method; (c) The denoising result of the previous denoising method; (d) The result of the filtering degree parameter is selected by two previous thresholds in the previous denoising method.
Figure 9. The denoising results based on Figure 1a of the proposed denoising method in this paper and the previous denoising method: (a) The denoising result of the denoising method proposed in this paper; (b) The result of filtering the degree parameter is selected by two adaptive thresholds in the proposed denoising method; (c) The denoising result of the previous denoising method; (d) The result of the filtering degree parameter is selected by two previous thresholds in the previous denoising method.
Remotesensing 11 00396 g009
Figure 10. The detection results of the fuzzy c-means (FCM) based on Figure 9b,d: (a) FCM detection result when h = 0.01 ; (b) FCM detection result when h = 0.019 ; (c) FCM detection result when h = 0.032 ; (d) FCM detection result when h = 0.05 .
Figure 10. The detection results of the fuzzy c-means (FCM) based on Figure 9b,d: (a) FCM detection result when h = 0.01 ; (b) FCM detection result when h = 0.019 ; (c) FCM detection result when h = 0.032 ; (d) FCM detection result when h = 0.05 .
Remotesensing 11 00396 g010
Figure 11. Detection results of the image shown in Figure 9a: (a) Detection result of NACA (new adaptive cultural algorithm); (b) Detection result of QSFLA-NSM (quantum-inspired shuffled frog leaping algorithm combining the new search mechanism); (c) Detection result of ACA-IQPSO (adaptive cultural algorithm with improved quantum-behaved particle swarm optimization); (d) Detection result of QSFLA (quantum-inspired shuffled frog leaping algorithm); (e) Detection result of CPSO (cultural particle swarm optimization algorithm); (f) Detection result of QPSO (quantum-behaved particle swarm optimization).
Figure 11. Detection results of the image shown in Figure 9a: (a) Detection result of NACA (new adaptive cultural algorithm); (b) Detection result of QSFLA-NSM (quantum-inspired shuffled frog leaping algorithm combining the new search mechanism); (c) Detection result of ACA-IQPSO (adaptive cultural algorithm with improved quantum-behaved particle swarm optimization); (d) Detection result of QSFLA (quantum-inspired shuffled frog leaping algorithm); (e) Detection result of CPSO (cultural particle swarm optimization algorithm); (f) Detection result of QPSO (quantum-behaved particle swarm optimization).
Remotesensing 11 00396 g011
Figure 12. The chart of the best fitness function values in Table 2.
Figure 12. The chart of the best fitness function values in Table 2.
Remotesensing 11 00396 g012
Figure 13. Detection results of the underwater sonar image (image size: 112×117): (a) Original sonar image; (b) The denoising result of the denoising method proposed in this paper; (c) Detection result of NACA; (d) Detection result of QSFLA-NSM; (e) Detection result of ACA-IQPSO; (f) Detection result of QSFLA; (g) Detection result of CPSO; (h) Detection result of QPSO.
Figure 13. Detection results of the underwater sonar image (image size: 112×117): (a) Original sonar image; (b) The denoising result of the denoising method proposed in this paper; (c) Detection result of NACA; (d) Detection result of QSFLA-NSM; (e) Detection result of ACA-IQPSO; (f) Detection result of QSFLA; (g) Detection result of CPSO; (h) Detection result of QPSO.
Remotesensing 11 00396 g013
Figure 14. The chart of the best fitness function values in Table 3.
Figure 14. The chart of the best fitness function values in Table 3.
Remotesensing 11 00396 g014
Figure 15. Detection result after the first iteration in the detection process of Figure 11: (a) Detection result of NACA; (b) Detection result of QSFLA-NSM; (c) Detection result of ACA-IQPSO; (d) Detection result of QSFLA; (e) Detection result of CPSO; (f) Detection result of QPSO.
Figure 15. Detection result after the first iteration in the detection process of Figure 11: (a) Detection result of NACA; (b) Detection result of QSFLA-NSM; (c) Detection result of ACA-IQPSO; (d) Detection result of QSFLA; (e) Detection result of CPSO; (f) Detection result of QPSO.
Remotesensing 11 00396 g015
Figure 16. The distribution of particle positions in the new update strategy and old update strategy: (a) Position distribution of particles on Sphere function; (b) Position distribution of particles on Griewank function.
Figure 16. The distribution of particle positions in the new update strategy and old update strategy: (a) Position distribution of particles on Sphere function; (b) Position distribution of particles on Griewank function.
Remotesensing 11 00396 g016
Figure 17. The result of the fitness value calculation in the new update strategy and old update strategy: (a) The optimization results of the Sphere function; (b) The optimization results of the Griewank function.
Figure 17. The result of the fitness value calculation in the new update strategy and old update strategy: (a) The optimization results of the Sphere function; (b) The optimization results of the Griewank function.
Remotesensing 11 00396 g017
Figure 18. Detection results of the underwater sonar image (image size: 259×368): (a) Original sonar image; (b) The denoising result of the denoising method proposed in this paper; (c) Detection result of NACA.
Figure 18. Detection results of the underwater sonar image (image size: 259×368): (a) Original sonar image; (b) The denoising result of the denoising method proposed in this paper; (c) Detection result of NACA.
Remotesensing 11 00396 g018
Figure 19. Detection results of the underwater sonar image (image size: 259×368): (a) Original sonar image; (b) The denoising result of the denoising method proposed in this paper; (c) Detection result of NACA.
Figure 19. Detection results of the underwater sonar image (image size: 259×368): (a) Original sonar image; (b) The denoising result of the denoising method proposed in this paper; (c) Detection result of NACA.
Remotesensing 11 00396 g019
Figure 20. Detection results of the underwater sonar image (image size: 203×101): (a) Original sonar image; (b) The denoising result of the denoising method proposed in this paper; (c) Detection result of NACA.
Figure 20. Detection results of the underwater sonar image (image size: 203×101): (a) Original sonar image; (b) The denoising result of the denoising method proposed in this paper; (c) Detection result of NACA.
Remotesensing 11 00396 g020
Table 1. The estimated optimal threshold and η of 20 different underwater sonar images.
Table 1. The estimated optimal threshold and η of 20 different underwater sonar images.
- h m i n h m a x η
10.0300.0500.539
20.0400.2000.609
30.0150.0300.689
40.0200.0350.624
50.0100.0500.634
60.0100.0400.584
70.0050.0400.617
80.0100.0350.681
90.0050.0300.532
100.0100.0400.575
110.0020.0050.567
120.0300.1000.604
130.0050.0400.676
140.0100.0600.636
150.0050.0150.650
160.0050.0200.649
170.0100.0250.605
180.0200.035 0.571
190.0050.070 0.556
200.0300.050 0.539
Table 2. The values of fitness function after each iteration in the detection process.
Table 2. The values of fitness function after each iteration in the detection process.
Iterative TimesNACAQSFLA-NSMACA-IQPSOQSFLACPSOQPSO
13.0242.1882.0281.9691.7371.724
23.0242.4952.3202.5542.1002.024
33.0242.6952.5202.5542.2362.083
43.0242.6952.5202.5542.3152.083
53.0242.6952.5202.5542.4172.241
63.0872.6952.8062.5542.4172.241
73.0872.6952.8062.5542.4212.241
83.0872.8622.8062.7272.4752.410
93.0872.8622.8062.7272.4752.410
103.0872.8622.8062.7272.4752.410
Table 3. The values of fitness function after each iteration in the detection process.
Table 3. The values of fitness function after each iteration in the detection process.
Iterative TimesNACAQSFLA-NSMACA-IQPSOQSFLACPSOQPSO
12.3652.1842.0031.9561.7942.066
22.4772.1842.0032.0801.9512.066
32.4772.2122.0032.0802.1032.170
42.4772.2122.1332.3122.1032.170
52.5322.2122.1332.3372.1622.170
62.5322.4262.1332.3372.2262.170
72.5322.4262.2732.3372.2262.184
82.5322.4262.2732.3372.2532.209
92.5322.4262.2732.3372.2532.209
102.5322.4262.2732.3372.2532.209

Share and Cite

MDPI and ACS Style

Wang, X.; Li, Q.; Yin, J.; Han, X.; Hao, W. An Adaptive Denoising and Detection Approach for Underwater Sonar Image. Remote Sens. 2019, 11, 396. https://doi.org/10.3390/rs11040396

AMA Style

Wang X, Li Q, Yin J, Han X, Hao W. An Adaptive Denoising and Detection Approach for Underwater Sonar Image. Remote Sensing. 2019; 11(4):396. https://doi.org/10.3390/rs11040396

Chicago/Turabian Style

Wang, Xingmei, Qiming Li, Jingwei Yin, Xiao Han, and Wenqian Hao. 2019. "An Adaptive Denoising and Detection Approach for Underwater Sonar Image" Remote Sensing 11, no. 4: 396. https://doi.org/10.3390/rs11040396

APA Style

Wang, X., Li, Q., Yin, J., Han, X., & Hao, W. (2019). An Adaptive Denoising and Detection Approach for Underwater Sonar Image. Remote Sensing, 11(4), 396. https://doi.org/10.3390/rs11040396

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

Article Metrics

Back to TopTop