An Unsupervised Method of Change Detection in Multi-Temporal PolSAR Data Using a Test Statistic and an Improved K & I Algorithm

In recent years, multi-temporal imagery from spaceborne sensors has provided a fast and practical means for surveying and assessing changes in terrain surfaces. Owing to the all-weather imaging capability, polarimetric synthetic aperture radar (PolSAR) has become a key tool for change detection. Change detection methods include both unsupervised and supervised methods. Supervised change detection, which needs some human intervention, is generally ineffective and impractical. Due to this limitation, unsupervised methods are widely used in change detection. The traditional unsupervised methods only use a part of the polarization information, and the required thresholding algorithms are independent of the multi-temporal data, which results in the change detection map being ineffective and inaccurate. To solve these problems, a novel method of change detection using a test statistic based on the likelihood ratio test and the improved Kittler and Illingworth (K&I) minimum-error thresholding algorithm is introduced in this paper. The test statistic is used to generate the comparison image (CI) of the multi-temporal PolSAR images, and improved K&I using a generalized Gaussian model simulates the distribution of the CI. As a result of these advantages, we can obtain the change detection map using an optimum threshold. The efficiency of the proposed method is demonstrated by the use of multi-temporal PolSAR images acquired by RADARSAT-2 over Wuhan, China. The experimental results show that the proposed method is effective and highly accurate.


Introduction
Change detection is the process of identifying the differences on the Earth's surface by multi-temporal images acquired in the same geographical area at different times [1,2].During recent years, because of the advantage of repetitive data acquisition, both satellite images and aerial photographs obtained by optical and synthetic aperture radar (SAR) sensors have been widely used in land-cover change detection.The change information is important because of its practical use in applications such as land-use and land-cover dynamic analysis [3,4] and environmental monitoring [5,6].
Although optical images have been widely applied in change detection [2,3,7], night-time and severe weather limit the use of optical images in practice.Because SAR is an active microwave sensor with all-weather, day-and-night operational imaging capability [8], the use of SAR sensors instead of optical sensors is attractive in change detection studies [9].Recently, many researchers have focused on change detection and time-series analysis based on SAR images [10][11][12].Several methods have been designed for single-channel SAR images [13][14][15][16], but the interpretation of the backscattering changes in land cover is limited [6].Compared with single-channel SAR images, polarimetric SAR (PolSAR) sensors acquire phase and amplitude information in different polarizations, and thus offer more scattering information [17,18].However, literature on change detection with PolSAR data is sparse [9].Because of the potential to provide improved detection capability and to extend the application of PolSAR data in change detection, we explore multi-temporal PolSAR data for change detection in this study.
Change detection based on multi-temporal PolSAR images can be classified into two categories: supervised and unsupervised [19].Although the supervised method, such as the widely used post-classification comparison (PCC), can provide more information on both the change extent and the type of change [18], it can be easily biased by the cumulative error caused by the classification of individual remote sensing image during change detection processing.Moreover, the supervised change detection technique depends on the classification accuracy of multi-temporal PolSAR imagery, which can be easily affected by manual intervention (such as the training sample selection).Compared with the supervised method, the unsupervised change detection approach is relatively straightforward, easy to implement and interpret without human intervention.Due to these reasons, this paper focuses on change detection based on an unsupervised method.
In general, unsupervised change detection in SAR images includes three steps: (1) preprocessing, (2) generating the comparison image (CI), and (3) making a decision based on the analysis of the CI in order to obtain the binary change map [9,18].The preprocessing of multi-temporal SAR images is a critical step in the change detection and involves radiometric calibration, image co-registration, and speckle filtering.Radiometric calibration is necessary for the comparison of SAR images acquired at different dates.After radiometric calibration, the pixel values of multi-temporal SAR images are related to the radar backscattering.Image co-registration [20,21] is aimed at reducing the errors caused by the misregistration of the images and aligns the images used in the change detection as precisely as possible.Speckle noise needs to be suppressed by speckle filtering [22] before change detection.In order to obtain a good result in unsupervised change detection, the second and third steps play an important role in processing.
The change detection result depends on both the quality of the CI (Step 2) and the quality of the threshold choice (Step 3) [11,14].In order to obtain the CI in the second step, there are many different approaches that can be used to compare the preprocessed images of the same geographical area at different dates.These methods can be divided into two main categories: (1) methods based on pixel radiometry and (2) methods based on local statistics [23].Each of these two categories has its own advantages and disadvantages.The approaches based on pixel radiometry are widely used and easy to design in multi-temporal SAR images, but they are easily influenced by speckle noise.The approaches based on local statistics combine the homogeneous pixels together, which effectively reduces the influence of speckle [24][25][26], but they do require many prior parameters.As a result, none of these methods can be regarded as optimal, and there is also no common standard at present [11].In this paper, we focus on change detection based on pixel radiometry.The methods based on pixel radiometry are widely used with multi-temporal SAR images.Such methods include the ratio or log-ratio operator of SAR amplitudes or intensities [13,15,16], change vector analysis (CVA) [27], principal component analysis (PCA) [28,29], the hidden Markov chain model [30,31], and the Kullback-Leibler divergence method [32].These methods are usually applied in multi-temporal single-channel SAR-based change detection.When a test statistic [5] was first used in change detection, it was assumed that the covariance (or coherency) matrices follow a complex Wishart distribution, but only the azimuthally symmetric case and the diagonal-element-only case were used to detect the changes.The method proposed in [4] uses all the elements of the covariance (or coherency) matrices to obtain the change detection map.In the third step of change detection analysis, the change detection map can be obtained by the decision threshold of the CI.Several algorithms have been widely used to automatically determine the threshold, such as the constant false alarm rate (CFAR) algorithm [33,34], Otsu's thresholding method [35], Kapur's entropy algorithm [36], and the Kittler and Illingworth (K&I) algorithm [37].
In this paper, we assume a complex Wishart distribution for the coherency matrix and apply a test statistic based on maximum likelihood estimation (MLE) and an improved K&I algorithm based on a generalized Gaussian model (GGM).The test statistic using more scattering information from the multi-temporal PolSAR images can obtain a high-precision CI.Moreover, the probability density function (PDF) of the CI can be simulated by the GGM.Because the traditional K&I algorithm using the Gaussian model to select the threshold automatically is not fit for CI from PolSAR data at different dates, we apply the GGM to the K&I algorithm to choose the threshold.After finishing the procedure, we can obtain the binary change detection map (changed and unchanged).
The remainder of this paper is organized as follows.Section 2 presents a brief introduction to the fundamental theory and details of the proposed change detection framework.Section 3 describes the experimental results.In Section 4, we draw our conclusion.

The Model of PolSAR Data
Assuming that a p-dimensional random complex vector K = [k 1 , k 2 , . . ., k p ] T follows a complex multivariate normal distribution with mean 0 and dispersion matrix Σ K = E KK * T , the p × p matrix Σ is a Hermitian positive definite random matrix and follows a Wishart distribution.A PolSAR system measures the amplitude and phase of the backscattered signals in four combinations of linear receive and transmit polarizations: horizontal-horizontal (hh), horizontal-vertical (hv), vertical-horizontal (vh), and vertical-vertical (vv) [5,26].Assuming that the target reciprocity condition is satisfied [22], the polarimetric information can be expressed by a complex vector: where h and v denote the horizontal and vertical wave polarization states, T indicates the vector transposition, and S hv is the scattering element of the horizontal transmitting and vertical receiving polarizations.
For multi-look processed PolSAR data, the backscattered signal can be expressed as a covariance matrix: The covariance matrix C can be modeled by a complex Wishart distribution.The frequency function of covariance matrix C can be shown as follows: where tr(•) is the trace operator, n is the number of looks, Γ p (n) is a normalization factor, and

Test Statistic for the Equality of Two Covariance Matrices
We consider that the two covariance matrices (C 1 ,C 2 ) from the bi-temporal PolSAR images are independent and follow a Wishart distribution as follows: where p represents the dimensions of C 1 and C 2 , and m and n represent the number of looks of C 1 and C 2 , respectively.Σ C 1 and Σ C 2 represent the dispersion matrix of C 1 and C 2 , and the MLE of Σ C 1 and Σ C 2 can be shown as follows: Assuming that the null hypothesis H 0 : ∑ C 1 = ∑ C 2 means that the two matrices are equal and there is a strong possibility of non-change, the alternative hypothesis H 1 : ∑ C 1 = ∑ C 2 means that the two matrices are different and there is a strong possibility of change [5].
We suppose that the test statistics based on MLE have joint densities , where θ is the set of parameters of the probability function that has generated the data.Then H 0 states that θ ∈ H 0 , and the likelihood ratio of the test statistic is shown as follows [5]: where where is the likelihood function, and f (•) is the frequency function.
Putting Equation (3) into Equation ( 6) and assuming ) can be expressed as follows [5]: The MLE of Equation ( 6) can be simplified as follows [5]: Assuming m = n, the CI of the bi-temporal PolSAR images can be denoted by d: where ρ can be represented as follows:

Improved K&I
After obtaining the CI of the bi-temporal PolSAR data, an automatic method of threshold selection is introduced.The K&I algorithm for the automatic estimation of the optimal threshold is suitable for SAR data [4,14,15,26] and has been widely used to distinguish the changed and unchanged classes in the CI.The K&I thresholding method is an extension of Bayes' minimum-error probability theory and can be shown as follows: where h(d l ) and L represent the histogram and the number of possible gray levels of the CI, respectively, and c(d l , T) denotes the cost of classifying pixels by comparing the corresponding gray level d l and threshold T. P(ω i |d l , T) (i = u, c) is the posterior probability, which represents the unchanged (or changed) class under the condition of gray level d l and a specific value of the threshold T.
The optimal threshold corresponding to minimizing the classification error is the following cost function [37]: Using Bayes' theorem, the posterior probability can be expressed by the prior probability and class-conditional PDF, i.e., The traditional K&I method assumes that the class-conditional distribution follows a Gaussian distribution, but this assumption cannot accurately reflect the distribution of the CI.In order to improve the threshold selection, an improved model based on the GGM describing the statistical behaviors of the changed and unchanged classes in the CI can be used.The model of the class-conditional PDF is shown as follows [14]: where terms m i (T), σ i (T), and β i (i = u, c) are the means, the variances, and the shape parameters of the distribution based on the unchanged and changed classes, respectively.Different values of shape parameters can generate different styles of density function [38].In the proposed method, we choose the GGM to estimate the class-conditional PDF, and the criterion J(T) can be expressed as follows [14]: Associated with the unchanged and changed classes, p u (T) and p c (T) denote the prior probabilities, m u (T) and m c (T) denote the means, and σ u (T) and σ c (T) denote the variances.These parameters of the above formula are estimated by the gray level d l and the histogram h(d l ) in Equation ( 16): Appl.Sci.2017, 7, 1297 6 of 22

The Proposed Method
The entire procedure of the proposed method is as follows: Step (1): The bi-temporal PolSAR images need to be co-registered and filtered.Image registration is performed to align the images used in the change detection.Speckle filtering is commonly used to suppress speckle noise before the change detection and classification of PolSAR images.The preprocessing is important for change detection.In this study, the refined Lee filter [39] was used to remove the speckle noise.
Step (2): The CI is obtained by the test statistic using the covariance matrices of the bi-temporal images.In this step, bi-temporal fully PolSAR data are used to generate the CI.
Step (3): The improved K&I algorithm based on the GGM is used to simulate the distribution of the CI and obtain the optimum threshold for the CI.
Step (4): The unchanged and changed classes of position (i, j) are determined.If d ij < T, this means that the PolSAR data in the same position at different dates are similar; otherwise, the corresponding pixel positions are different.Finally, the change detection map is obtained.
The detailed processing flow of the proposed method is shown in Figure 1.

Evaluation Criterion
Quantitative evaluation is important to determine the result of the change detection.When the ground truth is available, a quantitative evaluation can be performed [18,40].We used the false alarm (FA) rate, the total errors (TEs), the overall accuracy (OA), and the Kappa coefficient [41] to verify the accuracy of the results based on the proposed method.These indicators are calculated as follows:

Evaluation Criterion
Quantitative evaluation is important to determine the result of the change detection.When the ground truth is available, a quantitative evaluation can be performed [18,40].We used the false alarm (FA) rate, the total errors (TEs), the overall accuracy (OA), and the Kappa coefficient [41] to verify the accuracy of the results based on the proposed method.These indicators are calculated as follows: where TP (true positives) means the number of changed points correctly detected; TN (true negatives) means the number of unchanged points correctly detected; FP (false positives) means the number of unchanged points incorrectly detected as changed (false alarms); FN (false negatives) means the number of changed points incorrectly detected as unchanged (missed detections).N u and N c are the number of unchanged points and changed points of the ground reference change map, respectively.

Study Area and Background
The city of Wuhan (Figure 2) is situated in the east of Hubei province, China and lies in the eastern Jianghan Plain, at the intersection of the middle reaches of the Yangtze River.Its climate is humid subtropical, with abundant rainfall and four distinct seasons.From 2011 to 2016, the construction of a new tunnel and its ancillary buildings took place on East Lake.After the new tunnel was completed, the ancillary buildings were removed in 2016.Moreover, in July 2016, a 50-year return period of rainfall occurred, and the rainfall was significantly higher than the average annual precipitation.Wuhan was affected by continuous heavy rain, and some areas were seriously flooded and dramatically changed.The use of optical sensors was limited by this severe weather.In this study, our aim was to monitor the changes associated with the tunnel construction on East Lake from three temporal PolSAR images (2011, 2015, and 2016) and to use bi-temporal PolSAR images (2015 and 2016) to detect the changes after the heavy rainfall.

Study Area and Background
The city of Wuhan (Figure 2) is situated in the east of Hubei province, China and lies in the eastern Jianghan Plain, at the intersection of the middle reaches of the Yangtze River.Its climate is humid subtropical, with abundant rainfall and four distinct seasons.From 2011 to 2016, the construction of a new tunnel and its ancillary buildings took place on East Lake.After the new tunnel was completed, the ancillary buildings were removed in 2016.Moreover, in July 2016, a 50-year return period of rainfall occurred, and the rainfall was significantly higher than the average annual precipitation.Wuhan was affected by continuous heavy rain, and some areas were seriously flooded and dramatically changed.The use of optical sensors was limited by this severe weather.In this study, our aim was to monitor the changes associated with the tunnel construction on East Lake from three temporal PolSAR images (2011, 2015, and 2016) and to use bi-temporal PolSAR images (2015 and 2016) to detect the changes after the heavy rainfall.3a-c.Due to the large study areas, long intervals (2011-2015), and inclement weather conditions (heavy rainfall in July 2016), the regions labeled by the three red boxes are used to give a detailed assessment, to verify the efficiency of the proposed method.These labeled regions involve bridge, city, forest, and water bodies.Region 1 is East Lake, Region 2 is Yanxi Lake, and Region 3 is Yandong Lake.The ground reference maps over the three test sites were obtained based on the corresponding optical images obtained from Google Earth (version 7.18) and field surveys conducted by researchers.The images from the three different dates display the changes associated with the construction of the new tunnel on East Lake.In order to give a detailed assessment, we separated the dataset into two groups.The first group was made up of the data from 2011 and 2015, and the second group was made up of the data from 2015 and 2016.To verify the efficiency of the proposed method, comparative experiments were designed.We compared change detection methods based on log-ratio (HH) and improved K&I, log-ratio (HV) and improved K&I, log-ratio (VV) and improved K&I, CVA and improved K&I, a test statistic with a 5% significance level, a test statistic with a 1% significance level, a test statistic and K&I, a test statistic and CFAR, and a test statistic and improved K&I (the proposed method).The results of the change detection for the first and second groups of images are shown in Figures 4 and 5 The images from the three different dates display the changes associated with the construction of the new tunnel on East Lake.In order to give a detailed assessment, we separated the dataset into two groups.The first group was made up of the data from 2011 and 2015, and the second group was made up of the data from 2015 and 2016.To verify the efficiency of the proposed method, comparative experiments were designed.We compared change detection methods based on log-ratio (HH) and improved K&I, log-ratio (HV) and improved K&I, log-ratio (VV) and improved K&I, CVA and improved K&I, a test statistic with a 5% significance level, a test statistic with a 1% significance level, a test statistic and K&I, a test statistic and CFAR, and a test statistic and improved K&I (the proposed method).The results of the change detection for the first and second groups of images are shown in Figures 4 and 5, respectively.
The highlighted box labeled as Region 1 is used to give a detailed assessment.The land-cover types of this area are lake, bridge, city, and forest.The size of this area is 800 × 500 pixels, and the Pauli-RGB images are shown in Figure 6a-c.comparative experiments were designed.We compared change detection methods based on log-ratio (HH) and improved K&I, log-ratio (HV) and improved K&I, log-ratio (VV) and improved K&I, CVA and improved K&I, a test statistic with a 5% significance level, a test statistic with a 1% significance level, a test statistic and K&I, a test statistic and CFAR, and a test statistic and improved K&I (the proposed method).The results of the change detection for the first and second groups of images are shown in Figures 4 and 5, respectively.The highlighted box labeled as Region 1 is used to give a detailed assessment.The land-cover types of this area are lake, bridge, city, and forest.The size of this area is 800 × 500 pixels, and the Pauli-RGB images are shown in Figure 6a-c.The results of the CI maps and their 3-D plots are shown by linear mapping to [0, 1] in Figures 7 and 8, respectively.Although the CI maps of the log-ratio (HH) method (Figures 7a and 8a) show a strong response in the changed area, the contrast between the changed and unchanged areas is not obvious.Moreover, there is also a high response in the unchanged area because of the speckle noise, which impacts the accuracy of the change detection map.The CI maps of the CVA method (Figures 7b and 8b) show a weak response, which makes it difficult to distinguish the changed and unchanged areas.The proposed method (Figures 7c and 8c) using the fully polarimetric information also shows a strong response in the changed areas (the bridge and the edges of the water bodies).However, unlike the other methods, the proposed method shows a low response in the water bodies.Overall, the proposed method of generating the CI map is able to accurately distinguish the changed and unchanged areas.The results of the CI maps and their 3-D plots are shown by linear mapping to [0, 1] in Figures 7 and 8, respectively.Although the CI maps of the log-ratio (HH) method (Figures 7a and 8a) show a strong response in the changed area, the contrast between the changed and unchanged areas is not obvious.Moreover, there is also a high response in the unchanged area because of the speckle noise, which impacts the accuracy of the change detection map.The CI maps of the CVA method (Figures 7b and 8b) show a weak response, which makes it difficult to distinguish the changed and unchanged areas.The proposed method (Figures 7c and 8c) using the fully polarimetric information also shows a strong response in the changed areas (the bridge and the edges of the water bodies).However, unlike the other methods, the proposed method shows a low response in the water bodies.Overall, the proposed method of generating the CI map is able to accurately distinguish the changed and unchanged areas.
The results of the CI maps and their 3-D plots are shown by linear mapping to [0, 1] in Figures 7 and 8, respectively.Although the CI maps of the log-ratio (HH) method (Figures 7a and 8a) show a strong response in the changed area, the contrast between the changed and unchanged areas is not obvious.Moreover, there is also a high response in the unchanged area because of the speckle noise, which impacts the accuracy of the change detection map.The CI maps of the CVA method (Figures 7b and 8b) show a weak response, which makes it difficult to distinguish the changed and unchanged areas.The proposed method (Figures 7c and 8c) using the fully polarimetric information also shows a strong response in the changed areas (the bridge and the edges of the water bodies).However, unlike the other methods, the proposed method shows a low response in the water bodies.Overall, the proposed method of generating the CI map is able to accurately distinguish the changed and unchanged areas.After generating the CI maps, the change detection maps were obtained by the threshold methods of a 5% significance level, a 1% significance level, K&I, CFAR (Pfa = 0.005), and a test statistic and improved K&I (the proposed method).The results are shown in Figures 9 and 10, where it can be seen that the obvious changed areas can be easily detected by all methods.In order to further verify the efficiency of the proposed method in generating the CI, we used the same threshold strategy (improved K&I) in the methods of log-ratio of single-channel SAR, CVA, and test statistics, as shown in Figures 9a-d,i and 10a-d,i.As a result of the effect of speckle noise and the CI map with lower contrast, the results of log-ratio and improved K&I based on single-channel bi-temporal SAR images contain many wrong detections in some unchanged areas.The other change detection method of CVA and improved K&I results in many missed detections.Due to the high quality of the CI map After generating the CI maps, the change detection maps were obtained by the threshold methods of a 5% significance level, a 1% significance level, K&I, CFAR (Pfa = 0.005), and a test statistic and improved K&I (the proposed method).The results are shown in Figures 9 and 10, where it can be seen that the obvious changed areas can be easily detected by all methods.In order to further verify the efficiency of the proposed method in generating the CI, we used the same threshold strategy (improved K&I) in the methods of log-ratio of single-channel SAR, CVA, and test statistics, as shown in Figure 9a-d,i and Figure 10a-d,i.As a result of the effect of speckle noise and the CI map with lower contrast, the results of log-ratio and improved K&I based on single-channel bi-temporal SAR images contain many wrong detections in some unchanged areas.The other change detection method of CVA and improved K&I results in many missed detections.Due to the high quality of the CI map obtained using the fully polarimetric information based on the test statistic, the change detection map of the test statistic and improved K&I shows a better performance.These experiments show that using the multi-temporal full polarimetric data can achieve a better change detection result than using single-channel or multi-channel data.In order to compare the effect of the threshold choice, we used the CI map based on the test statistic and the threshold methods of 5% significance level, 1% significance level, K&I, CFAR (Pfa = 0.005), and improved K&I.Among these methods of choosing the threshold, 5% and 1% significance levels and K&I result in more false alarms in the region of city and forest than improved K&I.In addition, CFAR results in more false alarms in the region of water bodies than the proposed method.The quantitative comparisons of these detection schemes shown in Tables 1 and 2 also indicate that the proposed approach shows a better performance.The proposed approach of the test statistic and improved K&I achieves the best results in all four indicators (FA (%), TE (%), OA (%), and Kappa).This confirms that the proposed method is effective and shows a significant improvement over the other methods.(%), TE (%), OA (%), and Kappa).This confirms that the proposed method is effective and shows a significant improvement over the other methods.The highlighted box labeled as Region 2 (Yanxi Lake) is used to give a detailed assessment.This region contains grassland, city, and water bodies, and the size is 800 × 500 pixels.The main changes occurred in the water areas and the city.The RGB images in Pauli basis (|S hh − S vv | for red, |S hv | for green, and |S hh + S vv | for blue) and the ground reference of the changed areas are shown in Figure 11.
Test statistic with 5% significance level 6.17The highlighted box labeled as Region 2 (Yanxi Lake) is used to give a detailed assessment.This region contains grassland, city, and water bodies, and the size is 800 × 500 pixels.The main changes occurred in the water areas and the city.The results of the CI maps and their 3-D plots are shown by linear mapping to [0, 1] in Figure 12.Although the CI map of the log-ratio (HH) method (Figure 12a) shows a strong response in the changed area, the contrast between the changed and unchanged areas is not obvious.Moreover, there The results of the CI maps and their 3-D plots are shown by linear mapping to [0, 1] in Figure 12.Although the CI map of the log-ratio (HH) method (Figure 12a) shows a strong response in the changed area, the contrast between the changed and unchanged areas is not obvious.Moreover, there is also a high response in the unchanged area because of the speckle noise, which impacts the accuracy of the change detection map.The CI map of the CVA method (Figure 12b) shows a weak response, which makes it difficult to distinguish the changed and unchanged areas.The proposed method (Figure 12c) using the fully polarimetric information also shows a strong response in the changed areas.However, unlike the other methods, the CI values are almost homogeneous in the water bodies and forest areas.Overall, the proposed method of generating the CI map is effective at accurately distinguishing the changed and unchanged areas.
After generating the CI maps, the change detection maps were obtained by the threshold methods of 5% significance level, 1% significance level, K&I, CFAR (Pfa = 0.005), and the proposed method.The results are shown in Figure 13, where it can be seen that the obviously changed areas can be easily detected by all the methods.In order to further verify the efficiency of the proposed method in generating the CI map, we used the same threshold strategy (improved K&I) in the methods of log-ratio of single-channel SAR, CVA, and test statistics, as shown in Figure 13a-d,i.As a result of the speckle noise and the CI map with lower contrast, the result of log-ratio and improved K&I based on single-channel bi-temporal SAR images contains many wrong detections in some unchanged areas.The other change detection method of CVA and improved K&I results in many missed detections.Due to the high quality of the CI map obtained using the fully polarimetric information based on the test statistic, the change detection map of the test statistic and improved K&I shows a better performance.These experiments show that using the multi-temporal full polarimetric data can achieve a better change detection result than using the single-channel or multi-channel data.In order to compare the effect of the threshold choice, we used the CI map based on the test statistic and the threshold methods of 5% significance level, 1% significance level, K&I, CFAR (Pfa = 0.005), and improved K&I.Among these methods of choosing the threshold, the 5% and 1% significance levels and K&I result in more false alarms in the city and forest than the proposed method.In addition, CFAR results in more false alarms in the water bodies than the proposed method.The quantitative comparison of these detection schemes shown in Table 3 also indicates that the proposed approach shows a better performance.The proposed approach of the test statistic and improved K&I achieves the best results in all four indicators (FA (%), TE (%), OA (%), and Kappa).This confirms that the proposed method is effective and shows a significant improvement over the other methods.
The image of Yandong Lake (400 × 400 pixels) contains grassland, city, and water bodies.The main changes occurred in the water areas, because of the rain.The RGB images in Pauli basis (|S hh − S vv | for red, |S hv | for green, and |S hh + S vv | for blue), and the reference data of the changed areas are shown in Figure 14. is also a high response in the unchanged area because of the speckle noise, which impacts the accuracy of the change detection map.The CI map of the CVA method (Figure 12b) shows a weak response, which makes it difficult to distinguish the changed and unchanged areas.The proposed method (Figure 12c) using the fully polarimetric information also shows a strong response in the changed areas.However, unlike the other methods, the CI values are almost homogeneous in the water bodies and forest areas.Overall, the proposed method of generating the CI map is effective at accurately distinguishing the changed and unchanged areas.method.In addition, CFAR results in more false alarms in the water bodies than the proposed method.The quantitative comparison of these detection schemes shown in Table 3 also indicates that the proposed approach shows a better performance.The proposed approach of the test statistic and improved K&I achieves the best results in all four indicators (FA (%), TE (%), OA (%), and Kappa).This confirms that the proposed method is effective and shows a significant improvement over the other methods.The results of the CI maps and their 3-D plots are shown by linear mapping to [0, 1] in Figure 15.Although the CI map of the log-ratio (HH) method (Figure 15a) shows a strong response in the water bodies, the contrast between the changed and unchanged areas is not obvious.Moreover, there is also a high response in the unchanged area because of the speckle noise, which impacts the accuracy of the change detection map.The CI map of the CVA method (Figure 15b) shows a weak response, which makes it difficult to distinguish the changed and unchanged areas.The proposed method (Figure 15c) shows a strong response in the water bodies corresponding to the changed parts and low values in the water bodies corresponding to the unchanged parts.Overall, the proposed method of generating the CI map is effective at accurately distinguishing changed and unchanged areas.The results of the CI maps and their 3-D plots are shown by linear mapping to [0, 1] in Figure 15.Although the CI map of the log-ratio (HH) method (Figure 15a) shows a strong response in the water bodies, the contrast between the changed and unchanged areas is not obvious.Moreover, there is also a high response in the unchanged area because of the speckle noise, which impacts the accuracy of the change detection map.The CI map of the CVA method (Figure 15b) shows a weak response, which makes it difficult to distinguish the changed and unchanged areas.The proposed method (Figure 15c) shows a strong response in the water bodies corresponding to the changed parts and low values in the water bodies corresponding to the unchanged parts.Overall, the proposed method of generating the CI map is effective at accurately distinguishing changed and unchanged areas. non-change).
The results of the CI maps and their 3-D plots are shown by linear mapping to [0, 1] in Figure 15.Although the CI map of the log-ratio (HH) method (Figure 15a) shows a strong response in the water bodies, the contrast between the changed and unchanged areas is not obvious.Moreover, there is also a high response in the unchanged area because of the speckle noise, which impacts the accuracy of the change detection map.The CI map of the CVA method (Figure 15b) shows a weak response, which makes it difficult to distinguish the changed and unchanged areas.The proposed method (Figure 15c) shows a strong response in the water bodies corresponding to the changed parts and low values in the water bodies corresponding to the unchanged parts.Overall, the proposed method of generating the CI map is effective at accurately distinguishing changed and unchanged areas.After generating the CI maps, the change detection maps were obtained by the threshold methods of 5% significance level, 1% significance level, K&I, CFAR (Pfa = 0.005), and the proposed method.The results are shown in Figure 16, where it can be seen that the obviously changed areas can be easily detected by all the methods.In order to further verify the efficiency of the proposed method in generating the CI, we used the same threshold strategy (improved K&I) in the methods of log-ratio of single-channel SAR, CVA, and test statistics, as shown in Figure 16a-d,i.As a result of the speckle noise and CI map with lower contrast, the result of log-ratio and improved K&I based on single-channel bi-temporal SAR images contains many wrong detections in some unchanged areas.The other change detection method of CVA and improved K&I results in many missed detections.After generating the CI maps, the change detection maps were obtained by the threshold methods of 5% significance level, 1% significance level, K&I, CFAR (Pfa = 0.005), and the proposed method.The results are shown in Figure 16, where it can be seen that the obviously changed areas can be easily detected by all the methods.In order to further verify the efficiency of the proposed method in generating the CI, we used the same threshold strategy (improved K&I) in the methods of log-ratio of single-channel SAR, CVA, and test statistics, as shown in Figure 16a-d,i.As a result of the speckle noise and CI map with lower contrast, the result of log-ratio and improved K&I based on single-channel bi-temporal SAR images contains many wrong detections in some unchanged areas.The other change detection method of CVA and improved K&I results in many missed detections.Due to the high quality of the CI map obtained using the fully polarimetric information based on the test statistic, the change detection map of the test statistic and improved K&I shows a better performance.These experiments show that using the multi-temporal full polarimetric data can achieve a better change detection result than using single-channel or multi-channel data.In order to compare the effect of the threshold choice, we used the CI map based on the test statistic and the threshold methods of 5% significance level, 1% significance level, K&I, CFAR (Pfa = 0.005), and improved K&I.Among these methods of choosing the threshold, the 5% and 1% significance levels and K&I result in more false alarms in the city and forest than the proposed method.In addition, CFAR results in more false alarms in the water bodies than the proposed method.The quantitative comparison of these detection schemes shown in Table 4 also indicates that the proposed approach shows a better performance.The proposed approach of the test statistic and improved K&I achieves the best results in all four indicators (FA (%), TE (%), OA (%), Kappa).This confirms that the proposed method is effective and shows a significant improvement over the other methods.

Conclusions
In this paper, we have presented an improved multi-temporal PolSAR image change detection approach.A test statistic based on MLE and improved K&I were used to deal with the problem of change detection of a complex urban area containing bridges, water bodies, and buildings.Multitemporal data of Wuhan, China, were used to assess the performance of the proposed method.The

Conclusions
In this paper, we have presented an improved multi-temporal PolSAR image change detection approach.A test statistic based on MLE and improved K&I were used to deal with the problem of change detection of a complex urban area containing bridges, water bodies, and buildings.Multi-temporal data of Wuhan, China, were used to assess the performance of the proposed method.The main advantages of the proposed method are (1) the use of fully polarimetric information, (2) the use of the GGM to choose the PolSAR threshold, and (3) the computational simplicity and stability.Because obtaining the CI map is the key step in unsupervised change detection, the log-ratio method using single-channel information and CVA using linear computation of multi-channel information cannot generate a high-quality CI map, which directly influences the precision of the change detection map.The proposed method using fully polarimetric information is better able to detect slight differences and can effectively distinguish between changed and unchanged areas.The use of a test statistic is appropriate for PolSAR images, and it can be used to obtain an accurate CI map.The use of K&I to automatically choose the threshold under the assumption that the PDF of the CI follows a Gaussian distribution does not accurately match the aforementioned conditional densities of classes.We therefore use the improved K&I algorithm based on the GGM to choose the threshold.Using the proposed method, we were able to accurately detect the changes associated with the construction of a tunnel on East Lake from 2011 to 2016, as well as the changes in water bodies caused by the heavy rainfall in July 2016.The experimental results obtained on the multi-temporal PolSAR images of Wuhan confirmed the effectiveness of the proposed approach.

22 Figure 1 .
Figure 1.The flowchart of the proposed method.

Figure 1 .
Figure 1.The flowchart of the proposed method.

Figure 2 .
Figure 2. Location of the study area.

Figure 2 .
Figure 2. Location of the study area.

Figure 3 .
Figure3a-c.Due to the large study areas, long intervals (2011-2015), and inclement weather conditions (heavy rainfall in July 2016), the regions labeled by the three red boxes are used to give a detailed assessment, to verify the efficiency of the proposed method.These labeled regions involve bridge, city, forest, and water bodies.Region 1 is East Lake, Region 2 is Yanxi Lake, and Region 3 is Yandong Lake.The ground reference maps over the three test sites were obtained based on the corresponding optical images obtained from Google Earth (version 7.18) and field surveys conducted by researchers.

Figure 7 .
Figure 7. CI maps and 3-D plots of Region 1 based on (a) log-ratio (HH), (b) CVA, and (c) test statistics for the different dates of 7 December 2011 and 25 June 2015.Figure 7. CI maps and 3-D plots of Region 1 based on (a) log-ratio (HH), (b) CVA, and (c) test statistics for the different dates of 7 December 2011 and 25 June 2015.

Figure 7 .
Figure 7. CI maps and 3-D plots of Region 1 based on (a) log-ratio (HH), (b) CVA, and (c) test statistics for the different dates of 7 December 2011 and 25 June 2015.Figure 7. CI maps and 3-D plots of Region 1 based on (a) log-ratio (HH), (b) CVA, and (c) test statistics for the different dates of 7 December 2011 and 25 June 2015.

Figure 7 .Figure 8 .
Figure 7. CI maps and 3-D plots of Region 1 based on (a) log-ratio (HH), (b) CVA, and (c) test statistics for the different dates of 7 December 2011 and 25 June 2015.

Figure 8 .
Figure 8. CI maps and 3-D plots of Region 1 based on (a) log-ratio (HH), (b) CVA, and (c) test statistics for the different dates of 25 June 2015 and 6 July 2016.

Figure 9 .
Figure 9. Change detection results in Region 1 between 7 December 2011 and 25 June 2015 based on the following: (a) log-ratio (HH) and improved K&I; (b) log-ratio (HV) and improved K&I; (c) logratio (VV) and improved K&I; (d) CVA and improved K&I; (e) test statistic with 5% significance level; (f) test statistic with 1% significance level; (g) test statistic and K&I; (h) test statistic and CFAR; (i) test statistic and improved K&I (the proposed method).(j) Ground reference (white denotes the change and black denotes the non-change).

Figure 9 .
Figure 9. Change detection results in Region 1 between 7 December 2011 and 25 June 2015 based on the following: (a) log-ratio (HH) and improved K&I; (b) log-ratio (HV) and improved K&I; (c) log-ratio (VV) and improved K&I; (d) CVA and improved K&I; (e) test statistic with 5% significance level; (f) test statistic with 1% significance level; (g) test statistic and K&I; (h) test statistic and CFAR; (i) test statistic and improved K&I (the proposed method); (j) Ground reference (white denotes the change and black denotes the non-change).

Figure 10 .
Figure 10.Change detection results in Region 1 between 25 June 2015 and 6 July 2016 based on the following: (a) log-ratio (HH) and improved K&I; (b) log-ratio (HV) and improved K&I; (c) log-ratio (VV) and improved K&I; (d) CVA and improved K&I; (e) test statistic with 5% significance level; (f) test statistic with 1% significance level; (g) test statistic and K&I; (h) test statistic and CFAR; (i) test statistic and improved K&I (the proposed method); (j) Ground reference (white denotes the change and black denotes the non-change).

Figure 11 .
Figure 11.RADARSAT-2 PolSAR images of Region 2 indicated in Figure 3, acquired on (a) 25 June 2015 and (b) 6 July 2016; (c) The ground reference map (white denotes change and black denotes non-change).

Figure 12 .
Figure 12.CI maps and 3-D plots of Region 2 based on (a) log-ratio, (b) CVA, and (c) test statistics at the different dates of 25 June 2015 and 6 July 2016.Figure 12. CI maps and 3-D plots of Region 2 based on (a) log-ratio, (b) CVA, and (c) test statistics at the different dates of 25 June 2015 and 6 July 2016.

Figure 12 .
Figure 12.CI maps and 3-D plots of Region 2 based on (a) log-ratio, (b) CVA, and (c) test statistics at the different dates of 25 June 2015 and 6 July 2016.Figure 12. CI maps and 3-D plots of Region 2 based on (a) log-ratio, (b) CVA, and (c) test statistics at the different dates of 25 June 2015 and 6 July 2016.

Figure 13 .
Figure 13.Change detection results in Region 2 between 25 June 2015 and 6 July 2016 based on the following: (a) log-ratio (HH) and improved K&I; (b) log-ratio (HV) and improved K&I; (c) log-ratio (VV) and improved K&I; (d) CVA and improved K&I; (e) test statistic with 5% significance level; (f) test statistic with 1% significance level; (g) test statistic and K&I; (h) test statistic and CFAR; (i) test statistic and improved K&I (the proposed method).(j) Ground reference (white denotes the change and black denotes the non-change).

Figure 13 .
Figure 13.Change detection results in Region 2 between 25 June 2015 and 6 July 2016 based on the following: (a) log-ratio (HH) and improved K&I; (b) log-ratio (HV) and improved K&I; (c) log-ratio (VV) and improved K&I; (d) CVA and improved K&I; (e) test statistic with 5% significance level; (f) test statistic with 1% significance level; (g) test statistic and K&I; (h) test statistic and CFAR; (i) test statistic and improved K&I (the proposed method); (j) Ground reference (white denotes the change and black denotes the non-change).
Yandong Lake (400 × 400 pixels) contains grassland, city, and water bodies.The main changes occurred in the water areas, because of the rain.The RGB images in Pauli basis ( blue), and the reference data of the changed areas are shown in Figure14.

Figure 14 .
Figure 14.RADARSAT-2 PolSAR images of Region 2 indicated in Figure 3, acquired on (a) 25 June 2015 and (b) 6 July 2016.(c) The ground-truth map (white denotes the change and black denotes the non-change).

Figure 14 .
Figure 14.RADARSAT-2 PolSAR images of Region 2 indicated in Figure 3, acquired on (a) 25 June 2015 and (b) 6 July 2016.(c) The ground-truth map (white denotes the change and black denotes the non-change).

Figure 15 .
Figure 15.CI maps and 3-D plots of Region 3 based on (a) log-ratio, (b) CVA, and (c) test statistics at the different dates of 25 June 2015 and 6 July 2016.

Figure 15 .
Figure 15.CI maps and 3-D plots of Region 3 based on (a) log-ratio, (b) CVA, and (c) test statistics at the different dates of 25 June 2015 and 6 July 2016.

Figure 16 .
Figure 16.Change detection results in Region 3 between 25 June 2015 and 6 July 2016 based on the following: (a) log-ratio (HH) and improved K&I; (b) log-ratio (HV) and improved K&I; (c) log-ratio (VV) and improved K&I; (d) CVA and improved K&I; (e) test statistic with 5% significance level; (f) test statistic with 1% significance level; (g) test statistic and K&I; (h) test statistic and CFAR; (i) test statistic and improved K&I (the proposed method).(j) Ground reference (white denotes the change and black denotes the non-change).

Figure 16 .
Figure 16.Change detection results in Region 3 between 25 June 2015 and 6 July 2016 based on the following: (a) log-ratio (HH) and improved K&I; (b) log-ratio (HV) and improved K&I; (c) log-ratio (VV) and improved K&I; (d) CVA and improved K&I; (e) test statistic with 5% significance level; (f) test statistic with 1% significance level; (g) test statistic and K&I; (h) test statistic and CFAR; (i) test statistic and improved K&I (the proposed method); (j) Ground reference (white denotes the change and black denotes the non-change).

Table 1 .
Performance evaluation of the change detection on East Lake between 2011 and 2015.

Table 2 .
Performance evaluation of the change detection on East Lake between 2015 and 2016.

Table 3 .
Performance evaluation of the change detection on Yanxi Lake between 2015 and 2016.

Table 3 .
Performance evaluation of the change detection on Yanxi Lake between 2015 and 2016.

Table 4 .
Performance evaluation of the change detection on Yandong Lake between 2015 and 2016.

Table 4 .
Performance evaluation of the change detection on Yandong Lake between 2015 and 2016.