Change Detection from SAR Images Based on Convolutional Neural Networks Guided by Saliency Enhancement

: Change detection is an important task in identifying land cover change in different periods. In synthetic aperture radar (SAR) images, the inherent speckle noise leads to false changed points, and this affects the performance of change detection. To improve the accuracy of change detection, a novel automatic SAR image change detection algorithm based on saliency detection and convolutional-wavelet neural networks is proposed. The log-ratio operator is adopted to generate the difference image, and the speckle reducing anisotropic diffusion is used to enhance the original multitemporal SAR images and the difference image. To reduce the inﬂuence of speckle noise, the salient area that probably belongs to the changed object is obtained from the difference image. The saliency analysis step can remove small noise regions by thresholding the saliency map, and interest regions can be preserved. Then an enhanced difference image is generated by combing the binarized saliency map and two input images. A hierarchical fuzzy c-means model is applied to the enhanced difference image to classify pixels into the changed, unchanged, and intermediate regions. The convolutional-wavelet neural networks are used to generate the ﬁnal change map. Experimental results on ﬁve SAR data sets indicated the proposed approach provided good performance in change detection compared to state-of-the-art relative techniques, and the values of the metrics computed by the proposed method caused signiﬁcant improvement.


Introduction
The synthetic aperture radar (SAR) imaging process is not affected by sunlight, clouds, or the atmosphere because of the microwave imaging principle. In the field of SAR image processing, change detection is a very important topic. SAR images are an important information resource for change detection when studying disaster relief, agricultural detection, and urban planning, especially when evaluating the damage caused by natural disasters [1][2][3][4]. Because of the interference of scattering echo, speckle noise will inevitably be generated; it has the nature of multiplicative noise, and it seriously affects the interpretation of SAR images [5]. Therefore, noise suppression is an important task in the process of change detection. In recent decades, many change detection methods utilizing SAR imagery have been introduced, and these approaches can be divided into two categories: coherent and incoherent change detection methods, depending on whether interferometric phase information is used. In this paper, we mainly discuss incoherent change detection methods [6].
In terms of incoherent SAR image change detection algorithms, supervised and unsupervised methods are the most used. The main problem associated with supervised methods is the lack of ground reference data, which often involves manual marking processes that are labor-intensive and time-consuming. Unsupervised SAR change detection generally includes three basic steps: speckle noise reduction, difference map generation, and classification [7].
There are problems that influence the effect of SAR image change detection. First, the inherent speckle noise in the SAR images may cause false positives. If we use the denoising method directly, we can also remove useful information in the denoising process. Second, the difference image (DI) influences the detection results, and the changed information may be lost when obtaining the DI. Finally, because the prior information is used in the supervised models to train the classifier, the supervised models may generate a performance superior to the unsupervised models. Nevertheless, prior information is usually achieved by manual annotation, which requires much work and affects the generalization of the model. Hence, in terms of the SAR image change detection, effective denoising, feature extraction, and prior information acquisition ought to be considered.
As the result of these findings, for the sake of suppressing the speckle noise and preserving interest information, one saliency detection method is used to extract the interesting regions that probably pertain to the changed objects. To extract the changed information, the convolutional-wavelet neural networks (CWNNs) model is utilized to learn features from the denoised images and the difference image. The speckle reducing anisotropic diffusion (SRAD) is used to enhance the input images and the difference image. The saliency detection is performed on the enhanced difference image, and the interest regions are extracted. The hierarchical fuzzy c-means clustering (HFCM) model is used for pre-classification. Finally, the CWNN model is used to generate the final change map. Experiments on qualitative and quantitative comparisons demonstrate the advantages and innovations of the proposed SAR image change detection algorithm. The major contributions of the proposed method are concluded as follows.
(1) A saliency detection model is used in the proposed method, which aims to generate the salient regions that probably belong to the changed objects. The saliency detection model can extract attractive and compact salient areas from the difference image with a simple operation. It can remove background pixels and suppress noise. (2) A hierarchical fuzzy c-means clustering (HFCM) model is introduced in the proposed method and is used to select pixels with high probability of becoming changed or unchanged. The samples from the changed and unchanged parts are selected as the training set for the convolutional neural network. (3) A convolutional neural network based on dual-tree complex wavelet transform is constructed that aims to enhance the accuracy of change detection.

SAR Image Preprocessing
To cope with the speckle noise in SAR images, many noise reduction algorithms have been proposed in the past few years. The classical methods are probabilistic-patchbased (PPB) filter [8], speckle reducing anisotropic diffusion (SRAD) [9], non-local means filter [10,11], total variation model [12], block-matching 3D filtering [13], deep multi-scale recurrent network [14], multilook and refined Lee filtering [15], etc. Zheng et al. [16] introduced a SAR change detection method based on PPB filter and k-means clustering model where the PPB filter is used to suppress the speckle noise, and the two difference images are generated by subtraction and log-ratio operators, respectively; then the k-means clustering model is performed on the combined difference image to obtain the final change map. Because of the high computational complexity of the PPB algorithm, the whole change detection algorithm is relatively time-consuming. Su et al. [17] proposed a SAR change detection technique via non-local means with ratio similarity measurement; the enhanced non-local means model plays a role in suppressing speckle noise. Some other denoising algorithms have also generated good results in remote sensing image change detection. Lou et al. [18] introduced the ROF model semi-implicit denoising method for speckle reduction into the SAR change detection process, and this denoising algorithm played an important role in improving the accuracy of change detection. Wang et al. [19] proposed a new change detection technique based on logarithmic transformation and a total variation denoising approach. Liu et al. [20] introduced a change detection method based on mathematical morphology and a k-means clustering model, and the accuracy of the change detection improved. The nonsubsampled contourlet transform (NSCT) and nonsubsampled shearlet transform (NSST) are widely used in image fusion and denoising [21][22][23][24][25][26][27]. Chen et al. [28] introduced the NSCT-hidden Markov tree (NSCT-HMT) model to the remote sensing image change detection; Li et al. [29] proposed a multitemporal remote sensing image change detection algorithm based on the NSCT denoising model. Zhou et al. [30] introduced the image denoising model via NSST for SAR change detection; Shen et al. [31] proposed an image change detection approach using a Kalman filter and nonlocal means filter in the NSST domain. Because the different decomposition layers of NSCT and NSST have a certain impact on the image denoising effect, finding the appropriate decomposition layer is also an important problem. As an important part of SAR image processing, the denoising model is very important in the pre-processing of change detection.

DI Generation
The commonly used methods for generating the difference image (DI) are the subtraction operator and the ratio operator. Because the model of speckle noise inherent in SAR images is multiplicative random noise [5], the subtraction operator cannot suppress speckle-noise effectively. Compared with the subtraction operator, the ratio operator is the division between two SAR images. It can overcome the disadvantage of sensitivity to multiplicative noise. Based on the ratio operator, some improved models such as the log-ratio operator, mean-ratio operator, and neighborhood-based ratio operator are utilized. The log-ratio operator is robust to calibration and radiometric errors [32], it can reduce the influence of speckle noise, and it is widely used to generate the difference image, but the noisy regions can still be seen in the difference image. To create a different image with rich information, a combination of different operators has been widely used in recently published articles. Hou et al. [32] introduced the Gauss-log ratio operator as an improvement of the log-ratio operator, and this change detection method combines the difference images obtained by Gauss-log ratio and log-ratio operators in the discrete wavelet transform domain to generate a fused difference image. Zhang et al. [33] introduced the SAR change detection approach via fusing the Gauss-log ratio and subtraction operation in the shearlet transform domain. In the references [34,35], the difference images generated by mean-ratio and log-ratio operators, respectively, are effectively fused in a wavelet domain.
In recent years, some scholars have been working to obtain a difference map with more regions of interest and noise suppression, and the saliency detection was introduced. The saliency detection model can extract attractive and compact salient regions from the difference image with a simple operation [36]. It can remove the background pixels and reduce the influence of speckle noise. Geng et al. [36] introduced saliency-guided deep neural networks for SAR change detection. Zheng et al. [37] proposed an unsupervised saliency-guided SAR image change detection method where the context-aware saliency detection model is used to compute the difference image of the log-ratio operator. The thresholding model then combines the original SAR images to generate the new difference image; finally, the principal component analysis (PCA) and k-means clustering models are used to create the change map. This method offers a better performance in remote sensing image change detection compared to other state-of-the-art algorithms without saliency detection. Wang et al. [38] used the pattern and intensity distinctive analysis model to generate the saliency map of the difference image. Majidi et al. [39] proposed a saliency-guided neighborhood ratio model for automatic change detection. Li et al. [40] introduced an unsupervised SAR image change detection algorithm combining PCANet with saliency detection, and the accuracy of change detection was significantly improved.

Application of Classification to Change Detection
For the classification step of SAR image change detection, the pixels in the difference image are usually classified into changed and unchanged areas. These classification methods can be divided into unsupervised and supervised models. In terms of unsupervised change detection, the thresholding-, active contours-, and clustering-based methods are widely utilized in the classical change detection algorithms, and the prior information of the training set is not needed. Moser et al. [41] utilized generalized minimum-error thresholding for unsupervised SAR image change detection. Celik et al. [42] used a region-based active contour model for change detection. Hao et al. [43] used a superpixel-based active contour model for SAR image change detection. Meanwhile, clustering-based algorithms have been introduced for change detection. Jia et al. [44] introduced the multiple kernel k-means clustering with local-neighborhood information for remote sensing image change detection. In reference [45], the modified fuzzy c-means (FCM) algorithm was performed on the difference image to generate the change map. The FCM-based models cannot generate accurate change detection results, mainly because of the use of overall optimization objectives. Li et al. [46] introduced the fuzzy local-information c-means clustering model (FLICM) into the field of SAR change detection.
The deep learning-based classification model has had a rapid development and application in recent years, and these classification methods have generated an excellent performance and higher classification accuracy rate than the abovementioned approaches. Zhang et al. [47] proposed two-phase object-based deep learning for unsupervised SAR change detection. Gao et al. [48] introduced the channel weighting-based deep cascade network for unsupervised change detection. The PCANet was first proposed by Chan et al. [49] in 2015, and the PCANet model was introduced into the change detection field by Gao et al. [50] in 2016. This PCANet-based unsupervised change detection model achieves good performance, but it is quite time-consuming. Liu et al. [51] presented the convolutional neural network (CNN) model for unsupervised SAR change detection. Duan et al. [52] introduced a SAR image segmentation algorithm using convolutional-wavelet neural networks (CWNNs), and this technique achieves good segmentation results and suppresses the speckle noise. Gao et al. [53] introduced the CWNN model for unsupervised SAR image change detection, and it performs well in change detection.
Inspired by the SAR image change detection algorithms introduced previously, an effective and practical SAR image change detection utilizing saliency-guided convolutional neural networks is proposed in this paper.

Proposed SAR Image Change Detection Method
In this section, we elaborate on the proposed unsupervised SAR image change detection method. The proposed method can be divided into the following steps: difference image, extraction of salient regions, pre-classification, and classification by CWNN model. The log-ratio operator is used to generate the initial difference image, then the saliency detection is utilized to obtain the salient regions. The HFCM model is used for preclassification. Finally, the CWNN model is used to generate the final change map. The flow chart of the proposed automatic change detection algorithm is given in Figure 1.

Difference Image Generated by Log-Ratio Operator
Given two co-registered multi-temporal SAR images I 1 and I 2 , which are obtained from the same region at different times, t 1 and t 2 , respectively, the purpose of SAR image change detection is to produce a difference image that reflects the change information between t 1 and t 2 .
The initial SAR images will be affected by noise, and the speckle reducing anisotropic diffusion (SRAD) [9] is used to suppress the noise in multi-temporal SAR images I 1 and I 2 , generating the denoised images X 1 and X 2 corresponding to I 1 and I 2 , respectively. Figure 2 shows the example of the Ottawa images denoised by SRAD, and the peak signal-to-noise ratio (PSNR) index is used to measure the denoising effect; we can see that the noise is suppressed effectively from the area with the red box.
In this section, the initial difference image (DI1) is generated by the log-ratio (LR) operator, then SRAD is performed on the DI1, and the denoised difference image DI2 is obtained. The corresponding equation is:

Extraction of Salient Regions
Visual saliency regions contain information for visual image processing. In this section, the saliency detection theory is adopted to guide the change detection of SAR images. The initial difference image has a strong contrast region, and this is the salient region. We utilize the saliency detection model to locate the similar-change areas and optimize the proposed change detection task.
Suppose the X p denotes the intensity value of one pixel p in the image X. The saliency value V(p) of the pixel p can be calculated by the following [54]: where N is the total number of pixels in X. When two pixels have the same value of intensity, their saliency values are equal. Equation (2) can be modified as follows: where j depicts the pixel intensity, M j shows the number of pixels with an intensity equal to j, and L denotes the number of gray levels.
The saliency map of the denoised difference image DI2 is calculated by Equations (2) and (3), and we define the saliency map as D S .

Preclassification
In this section, the automatic threshold Otsu model is used to generate the binarized saliency map D s , and it is calculated by the following [40]: where p(x, y) denotes the gray value of the pixel in the salient area D s . The parameter τ shows the threshold computed by the Otsu model. In terms of D s , the value 1 indicates the salient pixels, and the value 0 indicates the non-salient pixels. The regions corresponding to the multi-temporal SAR images are extracted with the D s , and the following equation is given: where X i (i = 1, 2) shows the matrix generated by the denoised multi-temporal SAR images. denotes the dot product operator.
The new salient difference image D 2 is generated by utilizing the log-ratio model, and it is calculated by the following: When the D 2 is generated, the hierarchical FCM clustering (HFCM) model [55] is used to classify the D 2 into three components: the changed class Ω c , the unchanged class Ω u , and the intermediate class Ω i . The Ω c and Ω u are selected as the training samples, and Ω i is further classified by CWNN. More details of the hierarchical FCM clustering are contained in reference [55].

Classification by CWNN
CWNN was developed from convolutional neural networks (CNNs); it consists of convolutional layers, max-pooling layers, and fully connected layers [52,53]. In the CWNN model, the dual-tree complex wavelet transform (DTCWT) is introduced into the CNN model to reduce the effect of speckle noise in the SAR images. DTCWT has the advantages of good direction selectivity, limited redundancy, and a good reconstruction effect. DTCWT can decompose the layer preceding the pooling layer into eight components, including two low-frequency sub-bands LL 1 and LL 2 , and the high-frequency sub-bands in six orientations, ±15 • , ±45 • , and ±75 • (given by LH 1 , LH 2 , HL 1 , HL 2 , HH 1 , and HH 2 ) [56]. We chose the average of the two low-frequency sub-bands as the output of the pooling layer. Firstly, low-frequency components maintain the structures of the input layer according to the specified rules to better represent the patch of the input image. Secondly, some noises are suppressed by losing the high-frequency components.
In the CWNN method, the input of a wavelet pooling layer presents the output of the previous convolutional layer. In terms of each input feature map x, the DTCWT is utilized to generate the sub-bands: {LL 1 , LL 2 , LH 1 , LH 2 , HL 1 , HL 2 , HH 1 , where f (·) denotes the DTCWT function. The average of the low-frequency components is adopted as the output of the wavelet pooling layer, and the corresponding formula is defined as follows: where LL mean represents the output of the wavelet pooling layer. Figure 3 shows an example of the wavelet pooling layer. x i denotes one feature map after the convolutional layer. The eight sub-bands are generated by the DTCWT performed on the feature map, including two low-frequency sub-bands and six high-frequency subbands. The output feature map is achieved by averaging the two low-frequency sub-bands. The structure of CWNN is depicted in Figure 4. C 2 and C 4 denote the two convolutional layers, W 3 and W 5 present the wavelet pooling layer. Hence, the network can be depicted as {I 1 , C 2 , W 3 , C 4 , W 5 , F 6 , O 7 }. I 1 shows the input layer, and all the input image patches are resampled to 28 × 14. C 2 presents the convolutional layer with six convolutional kernels of 5 × 3. This layer generates six feature maps sized 24 × 12. W 3 shows the wavelet pooling layer. In this layer, all the input feature maps are decomposed with one-level DTCWT. This wavelet pooling layer generates six feature maps with the size of 12 × 6. C 4 depicts the convolutional layer with 12 convolutional kernels of 5 × 3. This layer produces 12 feature maps sized 8 × 4. W 5 shows the wavelet pooling layer, generating 12 feature maps sized 4 × 2. F 6 represents the fully connected layer with 96 units. O 7 presents the output layer with two units. In this section, the real samples and virtual samples are used as the training samples. When the CWNN training is finished, the image patches from Ω i are classified into changed and unchanged components, and the final change map is generated by the results of preclassification and CWNN classification.

Data Set Descriptions
In this section, four real SAR image data sets and one simulated SAR image data set were used to demonstrate the effectiveness of the proposed SAR image change detection method. Figures 5-9 show the SAR image data sets. Figure 5 is the Ottawa data set of two SAR images sized 290 × 350, captured by Radarsat-1, and they were acquired in May and August 1997. Figure 6 is the Coastline data set of two SAR images in Dongying, China, sized 450 × 280, captured by Radarsat-2; they were obtained in June 2008 and June 2009. Figure 7 shows the De Gaulle Airport data set of two SAR images acquired by ERS-1; they were taken in July 1997 and October 1998. The size of each SAR image is 240 × 370. Figure 8 is the Wenchuan data set of two SAR images sized 442 × 301, captured by ESA/ASAR on 3 March 2008 and 16 June 2008. These data mainly reflect the change caused by earthquake. The fifth data set in Figure 9 is the simulated data set of SAR images of the size 335 × 470 related to the village of Feltwell in the U.K. An image captured by the Daedalus 1268 Airborne Thematic Mapper (ATM) multispectral scanner was used as the reference image, and this image was assumed to be the Time 1 image of the data set. The Time 2 image was artificially generated from the reference one, and the land cover change was simulated by inserting some changes into the Time 1 image. The ground-truth images in Figures 5-9 were produced by artificial tagging with prior information. A description of the five data sets used in the experiment is given in Table 1.

Experimental Settings
In this section, 10 closely related algorithms are compared: PCAKM [57], change detection using log-ratio and Otsu (LROtsu) [34], change detection using mean-ratio and Otsu (MROtsu) [34], change detection using log-ratio and FCM (LRFCM) [58], change detection utilizing Gabor wavelet and two-level clustering (GaborTLC) [59], LMT [60], PCANet [50], NRELM [55], change detection using neighborhood-based ratio and collaborative representation (NRCR) [61], convolutional-wavelet neural networks for change detection (CWNN) [53]. We used the relevant parameter values of the original articles proposed by the authors. In the proposed method, the parameter's value setting in the CWNN model is consistent with that in reference [53]. A total of 10,000 pixels were randomly selected from Ω u and Ω c as the real samples, and 10,000 virtual samples were generated based upon these real samples. Table 2 shows the parameter settings of different methods. None MROtsu [34] None LRFCM [58] None GaborTLC [59] The size of median filter is 3 * 3 PCANet [50] σ = 1.20, stage number : 2, k = 5 NRELM [55] σ = 1.20, neighborhood size : 3, k = 5 NRCR [61] neighborhood size: 3 CWNN [53] patch size ω = 7 Proposed patch size ω = 7 To help judge the results of the change detection, five objective indicators were used as measures, namely, false positives (FPs) [50,62], false negatives (FNs) [50,63], overall errors (OEs) [50,64], percentage correct classification (PCC) [50,65], and kappa coefficient (KC) [50,66]. In the binary ground-truth image, we calculated the actual number of pixels belonging to the unchanged class (Nu) and the changed class (Nc). FP depicts the number of pixels belonging to the unchanged class but are falsely classified as the changed class. FN shows the number of pixels belonging to the changed class but are falsely classified as the unchanged class. The OE is the sum of FP and FN, and it is defined as follows: PCC is calculated by the following: KC is defined as follows: The quantitative evaluations of the change results computed by the proposed and comparison algorithms are summarized in Tables 3-8.

Results and Discussions
The corresponding experimental results of the Ottawa data set are reported in Table 3 and Figure 10. In Figure 10, we notice that the change map obtained by the proposed method outperformed other algorithms, and yielded better local consistency and fewer isolated pixels. The change maps generated by LROtsu, MROtsu, and LRFCM had many error pixels, and the three algorithms had the worst capacity to suppress the noise. Therefore, the FP values of LROtsu, MROtsu, and LRFCM were relatively high, the values were twice the FP value obtained by the proposed algorithm, as shown in Table 3. The change maps computed by PCAKM, GaborTLC, and LMT missed some changed regions, but the detection of the unchanged area was more accurate. Therefore, the FN values of the corresponding three methods were higher, while the FP values were lower, which was consistent with the data shown in Table 3. The PCANet, NRELM, NRCR, and CWNN techniques missed some unchanged regions. The OE value of the proposed algorithm for the Ottawa data set was reduced by 521, 229, 949, and 285 over PCANet, NRELM, NRCR, and CWNN, respectively. The PCC value of the proposed technique was improved by 0.51, 0.22, 0.93, and 0.28% over PCANet, NRELM, NRCR, and CWNN, respectively. KC was a comprehensive evaluation index; the KC value of the proposed algorithm was improved by 1.97, 1.03, 3.36, and 1.10% over PCANet, NRELM, NRCR, and CWNN, respectively.
The simulation results for the Coastline data set are shown in Figure 11 and Table 4. From Figure 11, we see that the PCAKM, LROtsu, MROtsu, LRFCM, GaborTLC, and NRCR had many noise spots, and the visual effect of change detection was poor with large FP. The change map generated by LMT was an improvement compared to the previous algorithms, but it missed some changed regions with large FN. Compared to the PCANet, NRELM, and CWNN methods, the OE value of the proposed method for the Coastline data set was reduced by 13,920, 5819, and 13,330, respectively. The PCC value of the proposed algorithm was improved by 11.04, 4.62, and 10.58% over PCANet, NRELM, and CWNN, respectively. The KC value of the proposed method was improved by 69.53, 55.20, and 69.40% over PCANet, NRELM, and CWNN, respectively. From the analysis of the results for the Coastline data set, the change map computed by the proposed method had a better change detection effect.    The simulation results for the De Gaulle Airport data set are given in Figure 12 and Table 5. From Figure 12, we see that the PCAKM, LROtsu, MROtsu, and LRFCM methods yielded much noise with large FP. The change results computed by GaborTLC, PCANet, NRELM, NRCR, and CWNN were poor. The LMT missed more changed regions with the largest FN value. Compared with the change maps generated by other algorithms, our method had a better performance in change detection on the De Gaulle Airport data set. In Table 5, we can observe that the values of OE, PCC, and KC generated by the proposed method were the best. The PCC value computed by the proposed method was improved by 12.93, 7.36, 15.52, and 9.07% over PCANet, NRELM, NRCR, and CWNN, respectively. The KC value computed by the proposed method was improved by 39.15, 31.52, 43.26, and 33.52% over PCANet, NRELM, NRCR, and CWNN, respectively. This means that the detection accuracy of our algorithm performed on the De Gaulle Airport data set was the highest compared to other approaches.
The simulation results for the Wenchuan data set are given in Figure 13 and Table 6. In Figure 13, we see that LROtsu, MROtsu, and LRFCM generated poor change detection results with some noise points; these three methods detected more unchanged areas as changed, so the corresponding FP values were large. The change detection results computed by PCAKM, GaborTLC, LMT, and NRELM were similar, but the results were still not ideal. In Table 6, we see that the proposed method had the least number of overall errors (OEs), and the values of PCC and KC were the best. The OE value of the proposed method for the Wenchuan data set was reduced by 1406, 3341, and 4900 over PCANet, NRCR, and CWNN, respectively. The PCC value of the proposed technique was improved by 1.05, 2.51, and 3.68% over PCANet, NRCR, and CWNN, respectively. The KC value of the proposed approach was improved by 4.82, 12.03, and 17.75% over PCANet, NRCR, and CWNN, respectively. The display of these data was obviously consistent with the analysis of the above image results. The change map generated by the proposed technique had higher detection accuracy and was closer to the reference true value image.
Remote Sens. 2021, 13, x FOR PEER REVIEW 15 of 23    The simulation results for the simulated data set are shown in Table 7 and Figure 14. Some isolated noise points were generated by the PCAKM, LROtsu, MROtsu, LRFCM, and LMT algorithms, and the FP values computed by the corresponding five methods were high. GaborTLC performed better than the previously mentioned change detection approaches. PCANet, NRELM, NRCR, and CWNN missed some changed regions with a large FN, although the four methods were effective for noise suppression. The proposed method effectively suppressed the noise and achieved a better detection effect that was closer to the true value image. The PCC value computed by the proposed method was improved by 0.83, 1.45, 0.79, and 1.24% over PCANet, NRELM, NRCR, and CWNN, respectively. The KC value generated by the proposed method was improved by 20.32, 41.38, 19.41, and 33.18% over PCANet, NRELM, NRCR, and CWNN, respectively.
To evaluate the superiority of the algorithm more accurately, we averaged the experimental data of five groups of SAR images, as shown in Table 8. Figure 15 shows the objective performance of different change detection algorithms on five SAR data sets. For each metric, the scores computed by an approach to different SAR data sets were connected to obtain a curve, and the average score is given in the legend, so we can directly see the fluctuation of data. In Table 8, we see that the values of OE, PCC, and KC computed by the proposed method were the best. Although the MROtsu method had the smallest miss detection rate, it had a higher false detection rate, and the values of FP and OE were the highest. LMT had the lowest false detection rate, but it had the highest miss detection rate, that is to say, the FN was the highest. In qualitative and quantitative analyses, our algorithm had absolute advantages in detection efficiency for SAR images.

Conclusions
In this paper, an effective SAR image change detection algorithm based on saliencyguided convolutional neural networks was proposed. The saliency map was adopted to guide the search for the interest regions in the initial difference image computed by a log-ratio operator, and the noise in the saliency map could be suppressed to some extent with the Otsu method. Then an enhanced difference image was generated by using the binarized saliency map and denoised input images. The hierarchical fuzzy c-means was used for pre-classification, and the final change map was obtained by the CWNN model. The experimental results demonstrated the effectiveness of the proposed change detection technique. Because the traditional saliency detection and the clustering model methods were used in the proposed algorithm, it was not an end-to-end deep learning model for SAR image change detection, so in future work, an end-to-end CWNN model for image change detection is what we need to construct.