Temporal Coherence Estimators for GBSAR

<jats:p>Many Ground-Based Synthetic Aperture Radar (GBSAR) applications demand preliminary analysis to select areas with high-quality signal. That is, areas in which the phase can be processed to extract the desired information. The interferometric coherence and the amplitude dispersion index are important tools widely used in the literature to assess the quality of GBSAR images. So far, no direct relation has been found between the two. Indeed, they are parameters of different natures: amplitude dispersion index is calculated with only amplitude values, while coherence provides information also on the signal phase. The purpose of this article is to find a relation between the two parameters. Indeed, the amplitude dispersion index provides some practical advantages if compared to coherence estimators, especially to perform fast preliminary analysis. In this article, a theoretical relation between amplitude dispersion index and coherence is retrieved. GBSAR measurements acquired in different scenarios, at different working frequencies are presented and used to validate such a relation.</jats:p>

section_content: Ground-Based Synthetic Aperture Radar (GBSAR) [1] interferometry is a powerful tool that enables the measurement of target displacements with sub-wavelength precision. This popular technique is currently applied for monitoring natural and urban scenarios. Interferometry techniques exploit the signal phase information and give reliable results only when applied to signals with a sufficient phase quality. Therefore, preliminary studies to single out areas of the investigated scene characterized by high-quality signal are demanded. 

One of the main problems that affects radar interferometry is the temporal decorrelation of the radar echoes [2]. In fact, there are several factors that can degrade the phase information as changes in atmospheric conditions or modifications in the physical properties of the scene between acquisitions. This effect is particularly evident in vegetated scenarios, where (due to wind, vegetation growth, etc. ) targets are subject to non-coherent displacements between successive acquisitions. For instance, authors of [3] showed that the short-term variable reflectivity induced by wind, typical of vegetated scenarios, can be theoretically characterized and modeled. In [4] different theoretical models for the temporal decorrelation of GBSAR images in vegetated scenarios are examined. For the short-term signal variability, a model with dependence on the wind speed is presented and validated for experimental data acquired in different scenarios. 

To estimate the decorrelation degree of GBSAR images, the interferometric coherence [5] and the amplitude dispersion index [6] parameters are largely employed. 

The interferometric coherence of radar images is a benchmark of the achievable interferogram quality [5] and is widely used to select image pixels with a sufficiently stable phase.
section_name: Introduction

section_content: The amplitude dispersion index (D A ) [6, 7] was introduced for the determination of the so-called permanent scatterers, i. e., targets whose signal amplitude remains stable over time, and it is a popular tool currently employed in interferometric analysis. 

Whereas the coherence is calculated on complex image values, the D A is based only on the signal amplitudes and seems to provide no information on phase content. Nevertheless, the authors of [7] showed that the amplitude dispersion index can provide an estimate for the phase dispersion, under the hypothesis of high signal-to-noise ratio. 

As will be further discussed in the following section of this paper, the D A provides practical advantages with respect to the interferometric coherence in assessing the image quality. For instance, being based on amplitude values, it is less sensitive to changes in the atmospheric parameters. Thus, it does not require atmospheric phase compensation before calculation, which implies less computational effort. Furthermore, the amplitude dispersion index parameter gives reliable results if calculated on a limited measurement time series. 

Given the paramount importance and theoretical relevance of the coherence parameter on the one hand, and the practical advantages provided by amplitude dispersion index on the other, it would be interesting to find a way to relate them. Up to now, no direct relation has been established between the two parameters. The aim of this article is to find a theoretical relation between coherence and D A. 

In order to give the D A a theoretical formulation, the authors of this paper start from the theory of radar scattering statistics [8] [9] [10]. Indeed, as the amplitude dispersion index involves first and second momenta of the amplitude distribution, it is possible to retrieve an analytic expression of this parameter. 

The analytic expression retrieved, which relates D A and coherence, is tested using GBSAR experimental data resulting from three measurement campaigns performed in different scenarios.
section_name: of 16

section_content: In this section, after having reviewed the theoretical definition of the signal coherence, some estimators commonly used to evaluate the phase quality of radar images are introduced and discussed. Then, to theoretically relate the amplitude dispersion index to the coherence, a statistical description of the radar backscattered signal is introduced, and an analytic formula that relates the two parameters is derived for the most common signal distributions. Finally, the theoretical relation is validated by simulations.
section_name: Materials and Methods

section_content: Given two zero-mean complex signals I n and I m, the complex coherence is defined as follows [5] :

where E[•] is the expectation value operator. The magnitude of this quantity is called the degree of coherence. The coherence is defined within [0, 1]. High values of γ nm correspond to high coherent signals. Generally, a signal is considered coherent when γ nm > 0. 75 ÷ 0. 9. 

In the context of GBSAR analysis, we are interested in evaluating the coherence between radar images of the same scenario, acquired at different times, say at time t 1, and t 2. Then, the above expression is called temporal or interferometric coherence and quantifies the signal quality of the image's pixel value I, between time t 1, and t 2. 

There are several factors that lead to a degradation of a radar image in terms of phase information, such as thermal effects, noise, atmospheric phase screen or physical changes in the scene [2, 3, 11, 12]. In GBSAR applications, we are particularly interested in the degradation related to physical changes in the scene. Indeed, these changes could fatally affect the interferometric image. 

Other sources of de-coherence, such as the atmospheric condition, could be compensated during the signal analysis. Indeed, these contributions are not related with the quality of the signal and can be removed by using specific techniques [13] [14] [15] [16] [17] [18] before coherence computing. 

Throughout the last decades, different models [2, [19] [20] [21] and estimators for the temporal coherence have been defined, mostly for SAR applications [4, 5, 22, 23]. 

An efficient tool for retrieving information on the signal autocorrelation and coherence is provided by the signal spectral analysis. In the frequency domain, the coherence is usually computed as the amplitude of normalized cross spectral density [24] [25] [26]. 

The signal power spectral density function can provide important information also about the characterization of GBSAR signals. For instance, in [27] the power spectral density is proposed in the context of synthetic aperture radar (SAR) processing as a tool for the statistical characterization of the different SAR modes and interferograms. 

To understand how the signal power spectral density is related to the coherence, let us write the autocorrelation function C(τ) for a certain time interval τ, of a signal time series I(t) acquired at generic times t, and its power spectral density S(f), respectively, as,

where f is the conjugate variable of τ. 

The Wiener-Kinchin theorem states that, under the assumption of a stationary random process, the Fourier transform (FT) of the autocorrelation signal ( 2 ) is equal to the signal power spectral density,

Under the assumption of stationary processes, the spectral density itself can provide information on the overall coherence degree of the time series images. Indeed, the autocorrelation, if properly normalized is a common coherence estimator. For a stationary process it is E |I n | 2 = constant, for any t n, and normalization factors in Equation (1) do not depend on time. Therefore, by applying the Fourier transform to Equation (1), and combining Equations ( 2 ) and (4), we have

Equation ( 5 ) provides us with an expression of the signal coherence, in the frequency domain. If only one frequency dominates the spectral density S(f), then, the peak value of the spectral density, properly normalized, is an estimator for the coherence degree, relative to that frequency in the spectrum. Indeed, if the investigated area is stationary, or subject to a constant uniform movement, the spectrum will be characterized by a single peak at zero or at a constant frequency, respectively. This method provides us with an estimate of the temporal coherence averaged over the time interval. 

It is worth noting that almost regular sampling times are required to properly calculate the Fourier transform, hence, the signal spectral density. The GBSAR systems usually acquire long time series of images of the same scenario, with regular return times in the order of minutes. Thus, Equation (5) may represent an excellent tool for estimating the signal coherence in GBSAR applications. 

Finally, it is worth saying that the coherence estimate of Equation ( 5 ) presents great sensitivity to changes in atmospheric conditions, since the variations of the signal phase values over time influence the evaluation of the Fourier transform.
section_name: The Coherence Estimator

section_content: The amplitude dispersion index is an important tool for evaluating the signal stability over time, and was introduced in the context of permanent scatterers determination for interferometric analysis [6, 7, 28, 29]. It is defined as

where A = |I| is the image amplitude, and STD A, m A are the standard deviation and the mean value of the amplitude time series, respectively. Low values of D A correspond to high coherent signals. Generally, the threshold below which a pixel is considered a permanent scatterer is set to 0. 25. Since the D A evaluation is based on amplitude values, it is not very sensitive on changes in atmospheric conditions. Indeed, a time variation of the air refractive index profoundly changes the atmospheric phase contribution, hence, the signal phase information. Conversely, amplitude values are almost insensitive to these effects, and the D A evaluation does not demand atmospheric phase corrections. 

As already said, despite being calculated with amplitude values, it has been shown that the D A provides a good estimate of the phase dispersion in the hypothesis of high signal-to-noise ratio [7]. This result was obtained investigating the statistical properties of the signal amplitude. As it suggests a connection between D A and coherence of the interferometric phase, in the following section a statistical description of radar signal amplitudes is reviewed.
section_name: The Amplitude Dispersion Index

section_content: Statistical descriptions of radar signals have been studied at the very beginning of this technology [8-10, [30] [31] [32] [33]. Depending on the radar signal characteristics or the physically investigated scenario, different statistics have been proposed for both the interferometric phase and the signal amplitude. In what follows, we review basic concepts related to the statistical properties of the backscattered signal from distributed targets. Indeed, the signal amplitude statistics can provide important insights on the amplitude dispersion index and on its relation to the other parameters of interest. 

Let us start modelling the response of the so-called Gaussian scatterers, i. e., scatterers that can be decomposed into a sufficiently high number of random independent scatterers within a resolution cell [33]. If no scatterer dominates the others, for the central limit theorem, the radar image pixel value I(t) is a zero mean complex circular Gaussian random variable, arising from the contribution of independent scatterers located in the same resolution cell. Thus, the joint probability density function for the real and imaginary part of

where σ is the standard deviation. By performing a change of variables, one can derive a joint probability distribution for the amplitude A and phase ϕ of the signal. Specifically, by substituting Re(I) = Acos(ϕ)

into Equation (7), and after integration over the variable ϕ, we get the following expression for the signal amplitude statistics,

Equation ( 9 ) is a Rayleigh distribution [8], with scale parameter σ. Given the assumptions made (large number of independent contributions and no scatterer that dominates the others), this model has been regarded in the literature as that describing the speckle pattern [8]. Good agreement with the Rayleigh distribution has been observed for the amplitude of SAR experimental data in homogeneous regions with coarse spatial resolution [9]. Conversely, it fails in describing the radar statistics of heterogeneous areas imaged with fine spatial resolution. 

Let us now consider the case when a coherent background signal is summed to the independent scatterers contributions [8]. That is, a part of the scatterers within the resolution cell remains stable over time. Indeed, this may be the case in natural scenarios, when investigating land partially covered by vegetation. In this case the signal can be modeled as the sum of a random circular Gaussian variable and a constant complex phasor. In the following, we assume that the Gaussian scatterers are completely independent from the coherent background. Without loss of generality, we also assume the coherent background image value to be real, equal to A S. As a result, the real part of the image pixel I(t) has mean value equal to A S. Under these hypotheses, the joint probability distribution for the real and imaginary part of the image pixel value I(t), reads pdf(Re(I),

By performing the change of variables defined in Equation ( 8 ), and after integration over the phase variable ϕ, we retrieve the following relation:

which is called the Rice distribution [8]. The Rician model was found appropriate to describe low resolution SAR images [10], where strong scatterers are embedded into the surrounding weak clutter environment. Our intent is to understand whether this simple model can account for GBSAR amplitude images, acquired in vegetated scenarios.
section_name: Amplitude Statistics of Radar Data

section_content: An a-priori expression of the amplitude dispersion index can be retrieved by considering the first (m A ) and the second (STD A ) moment of the proper amplitude radar statistic. 

For example, considering the Rayleigh distribution (9), the value of dispersion index is constant, equal to

Indeed, this value corresponds to the dispersion index of radar clutter, which typically follows the Rayleigh distribution. 

If the amplitude signal is instead Rice distributed, it depends on the parameters A S, and σ, according to

where L -2 1 2 (•) is the Laguerre polynomial of degree 1 2. Equation ( 13 ) relates the amplitude dispersion index values to the ratio

Under the same assumptions on the statistical distribution, it is possible to retrieve an expression for the coherence of radar images, as a function of the parameters A S, and σ. Let us recall the expressions of the first and second distribution moment for the variables Re(I) and Im(I),

where, without loss of generality, we assumed the coherent background image value to be real, equal to A S (see Section 2. 3). 

By writing Equation ( 1 ) in terms of real and imaginary parts of the image I, and using the expressions in ( 14 ), we end with

Then, by combining Equations ( 13 ) and ( 15 ), we obtain an expression that relates the coherence with the amplitude dispersion index. 

This equation establishes a direct relation between the D A and the coherence, under the hypothesis of Rician amplitude statistics. As will be shown in the next section, for GBSAR systems, in rural environments where low vegetation is present, this condition is met, and Equation ( 16 ) successfully reproduces the parameters relation. 

It is worth noting that the method used to retrieve Equation ( 16 ) is rather general and could also be used for signals characterized by different statistics.
section_name: Relation between Amplitude Dispersion Index and Coherence

section_content: With the aim to provide solid basis for the above formulation, and to assess the validity of Equation ( 16 ), the two models considered in the above section have been simulated. 

The echo of a continuous wave stepped frequency signal operative in Ku-band were simulated. For the sake of simplicity, we limited the analysis in one dimension, by simulating the signal of a monostatic radar. The echo reflected from 100 scatterers, uniformly distributed inside the same one-dimensional resolution cell, was simulated. The amplitude reflectivity of each scatterer was uniformly distributed between 0 and 1. 

All the scatterers were subject to random displacements inside the resolution cell, between each acquisition, in order to simulate a model, in which the signal amplitude mimics the Rayleigh distribution. The Rice statistics were instead reproduced by fixing the position of a certain number of scatterers, while moving randomly the others, always inside the resolution cell. In the following we refer to these two models to as the "Rayleigh" and "Rice" model, respectively. We simulated acquisitions of N t consecutive measurements. For the Rice model, the simulation was replicated for different values of the ratio of stable to total scatterers, which was tuned from 0 to 1. 

To collect statistics and averaging over the initial positions of the scatterers, in every run of the simulation, the time series was repeated 100 times; before starting any new time series, the scatterers position inside the resolution cell was changed randomly. 

Figure 1 shows the obtained results concerning the amplitude values, both for Rayleigh (a), and Rice model (b). Histograms in gray represent simulated data, while the red and green curves are the Rayleigh and Rice functions (Equations ( 9 ) and (11), respectively), that result from a fit of the simulated data. 

It is worth noting that, in Figure 1a, when only random scatterers are present, the Rice distribution that best fits the data, reduces to the Rayleigh one. On the other hand, if both fixed and random scatterers are present (Figure 1b ), the Rayleigh distribution no longer fits simulated data. Simulated models successfully reproduced the expected theoretical amplitude distributions. For each time series, the coherence and the D were computed. The coherence was calculated by using the spectral density estimator defined in Equation ( 5 ). Indeed, this expression provides a coherence value for the whole time series, as the amplitude dispersion index, thus enabling the parameter comparison. Hence, for each run of the simulation, a single value of D and coherence is calculated. Specifically, we evaluated the spectral density of the signal and set the coherence equal to the peak value, properly normalized. 

Figure 2 For each time series, the coherence and the D A were computed. The coherence was calculated by using the spectral density estimator defined in Equation ( 5 ). Indeed, this expression provides a coherence value for the whole time series, as the amplitude dispersion index, thus enabling the parameter comparison. Hence, for each run of the simulation, a single value of D A and coherence is calculated. Specifically, we evaluated the spectral density of the signal and set the coherence equal to the peak value, properly normalized. 

Figure 2 For each time series, the coherence and the D were computed. The coherence was calculated by using the spectral density estimator defined in Equation ( 5 ). Indeed, this expression provides a coherence value for the whole time series, as the amplitude dispersion index, thus enabling the parameter comparison. Hence, for each run of the simulation, a single value of D and coherence is calculated. Specifically, we evaluated the spectral density of the signal and set the coherence equal to the peak value, properly normalized. 

Figure 2 12 ) and ( 16 ) for the Rayleigh and for the Rice model, respectively.
section_name: Simulations

section_content: Given the promising results obtained with simulations, authors of this paper aimed to validate Equation (16) for experimental data acquired in natural scenarios. That is, determining whether the simplified model introduced in Section 2, could effectively describe the scattering properties of different physical scenarios. To give solid basis to the relation between amplitude dispersion index and coherence, data from three GBSAR measurement campaigns were analyzed. This enabled us to test the relation in different scenarios and with different working frequencies. 

The image was focused on a range-azimuth grid using a backpropagation algorithm [34]. Using this algorithm, each range-azimuth resolution cell was covered by several pixels. The amplitude statistics were evaluated by averaging over an ensemble of the pixels whose overall size is approximately equal to that of a single resolution cell. Indeed, this ensemble of pixels properly represents the backscattered signal of the physical illuminated target. For each measurement series, we selected only resolution cells where the scatterers were present. Then, the amplitude dispersion index and the coherence were calculated to test the theoretical relation (16).
section_name: Results

section_content: For our scope, we analyzed GBSAR measurements acquired in Formigal, Spain, from 13 22 November 2006. The purpose of this campaign was to monitor a vegetated slope, subject to a landslide. The radar used was a C-band GBSAR able to acquire an image about every hour. Figure 3 shows an amplitude radar image (a) and a picture (b) the illuminated scene. The area was characterized by the presence of rocks, grass, and low vegetation. 

Remote Sens. 2022, 14, 3039 8 of 16 theoretical curves correspond to the analytic expressions obtained in Equations ( 12 ) and ( 16 ) for the Rayleigh and for the Rice model, respectively.
section_name: Formigal Measurement Campaign

section_content: For our scope, we analyzed GBSAR measurements acquired in Formigal, Spain, from 13 to 22 November 2006. The purpose of this campaign was to monitor a vegetated slope, subject to a landslide. The radar used was a C-band GBSAR able to acquire an image about every hour. Figure 3 shows an amplitude radar image (a) and a picture (b) the illuminated scene. The area was characterized by the presence of rocks, grass, and low vegetation. Figure 4 shows the histogram of signal amplitude values, obtained for the resolution cells in the region highlighted by the red box in Figure 3a. This region was a vegetated area characterized by meadows. The colored curves correspond to the theoretical distributions that best fit the experimental data. One can notice a good agreement between the experimental data and the Rice distribution. On such areas, amplitude dispersion index and coherence were calculated for each resolution cell. As done with simulated data, coherence was computed by using the spectral density function (5). As already stressed, changes in the atmospheric parameters affect the coherence values. Thus, a processing to remove the atmospheric phase screen was implemented. Results concerning D and coherence values are shown in Figure 5a. Experimental data appear to be arranged according to the theoretical curve, thus, validating the relation (16) in this specific case. To quantify the agreement between experimental and theoretical data, we calculated the squared deviation between the theoretical ( 16 ) and experimental value of D, by averaging on a group of almost one hundred targets with similar coherence. Figure 6b shows the obtained results. On such areas, amplitude dispersion index and coherence were calculated for each resolution cell. As done with simulated data, coherence was computed by using the spectral density function (5). As already stressed, changes in the atmospheric parameters affect the coherence values. Thus, a processing to remove the atmospheric phase screen was implemented. Results concerning D A and coherence values are shown in Figure 5a. Experimental data appear to be arranged according to the theoretical curve, thus, validating the relation (16) in this specific case. To quantify the agreement between experimental and theoretical data, we calculated the squared deviation between the theoretical (16) and experimental value of D A, by averaging on a group of almost one hundred targets with similar coherence. Figure 6b shows the obtained results. On such areas, amplitude dispersion index and coherence were calculated for each resolution cell. As done with simulated data, coherence was computed by using the spectral density function (5). As already stressed, changes in the atmospheric parameters affect the coherence values. Thus, a processing to remove the atmospheric phase screen was implemented. Results concerning D and coherence values are shown in Figure 5a. Experimental data appear to be arranged according to the theoretical curve, thus, validating the relation (16) in this specific case. To quantify the agreement between experimental and theoretical data, we calculated the squared deviation between the theoretical ( 16 ) and experimental value of D, by averaging on a group of almost one hundred targets with similar coherence. Figure 6b shows the obtained results.
section_name: Formigal Measurement Campaign

section_content: To investigate the relation between D and coherence for different frequency bands, measurements of a multi-frequency GBSAR campaign, acquired in Aquilonia, Italy, were analyzed. During this measurement campaign, three GBSAR systems working in C, S, and Ku-band illuminated a vegetated tailing dam [35]. The acquisition rate of the radars in C and S-band were about 580 s. The acquisition rate of the Ku-band radar was 240 s. 

For the present analysis, measurements acquired from 12 March to 3 April 2022 were studied, for a total of almost 4000 radar images. 

In Figure 6 amplitude images acquired in the three bands are shown, together with a picture of the illuminated scene. As can be seen in Figure 6d, the dam body is covered by grass and low vegetation. 

We focused the analysis on the dam body highlighted in red in Figure 6d. First, we tested the Rician statistics of the amplitude distribution. Areas for which we obtained excellent agreement are those characterized by meadows, without higher vegetations. 

In order to test whether Equation ( 16 ) could also apply for the other parts of the dam body, amplitude dispersion index and coherence have been calculated for each resolution cell in the area contoured in red. Figures 7a, 8a and 9a show the results concerning the relation between the two parameters obtained for S, C, and Ku-band data, respectively. Each blue dot corresponds to a specific resolution cell. It is worth noting good agreement between experimental data and the theoretical curve (16), represented by the red line, at
section_name: Aquilonia Test Site: Validation with Different Frequencies

section_content: To investigate the relation between D A and coherence for different frequency bands, measurements of a multi-frequency GBSAR campaign, acquired in Aquilonia, Italy, were analyzed. During this measurement campaign, three GBSAR systems working in C, S, and Ku-band illuminated a vegetated tailing dam [35]. The acquisition rate of the radars in C and S-band were about 580 s. The acquisition rate of the Ku-band radar was 240 s. 

For the present analysis, measurements acquired from 12 March to 3 April 2022 were studied, for a total of almost 4000 radar images. 

In Figure 6 amplitude images acquired in the three bands are shown, together with a picture of the illuminated scene. As can be seen in Figure 6d, the dam body is covered by grass and low vegetation. 

We focused the analysis on the dam body highlighted in red in Figure 6d. First, we tested the Rician statistics of the amplitude distribution. Areas for which we obtained excellent agreement are those characterized by meadows, without higher vegetations. 

In order to test whether Equation ( 16 ) could also apply for the other parts of the dam body, amplitude dispersion index and coherence have been calculated for each resolution cell in the area contoured in red. Figures 7a, 8a
section_name: Aquilonia Test Site: Validation with Different Frequencies

section_content: Finally, to further test the relation in an environment with different characteristics, a measurement series acquired in an urban scenario in the premises of the University of Florence is presented. In this case, the same GBSAR systems employed in the test site at Aquilonia, acquired images for 5 days. Figure 10 shows an amplitude image (a) and a picture (b) the illuminated scene. As can be observed, various kinds of targets are present in the area; namely, grass, some bushes, a tree, and buildings. In this case, the previous analysis was computed for the area corresponding to the whole focused image in Figure 10a. Results obtained are shown in Figures 11 12 13 for the S, C, and Ku-band, respectively. For all the three bands, good agreement with the theoretical curve is found.
section_name: Measurements in Urban Scenario: Validation with Different Frequencies

section_content: Finally, to further test the relation in an environment with different characteristics, a measurement series acquired in an urban scenario in the premises of the University of Florence is presented. In this case, the same GBSAR systems employed in the test site at Aquilonia, acquired images for 5 days. Figure 10 shows an amplitude image (a) and a picture (b) the illuminated scene. As can be observed, various kinds of targets are present in the area; namely, grass, some bushes, a tree, and buildings.
section_name: Measurements in Urban Scenario: Validation with Different Frequencies

section_content: Finally, to further test the relation in an environment with different characteristics, a measurement series acquired in an urban scenario in the premises of the University of Florence is presented. In this case, the same GBSAR systems employed in the test site at Aquilonia, acquired images for 5 days. Figure 10 shows an amplitude image (a) and a picture (b) the illuminated scene. As can be observed, various kinds of targets are present in the area; namely, grass, some bushes, a tree, and buildings. In this case, the previous analysis was computed for the area corresponding to the whole focused image in Figure 10a. Results obtained are shown in Figures 11 12 13 for the S, C, and Ku-band, respectively. For all the three bands, good agreement with the theoretical curve is found. In this case, the previous analysis was computed for the area corresponding to the whole focused image in Figure 10a. Results obtained are shown in Figures 11 12 13 for the S, C, and Ku-band, respectively. For all the three bands, good agreement with the theoretical curve is found.
section_name: Measurements in Urban Scenario: Validation with Different Frequencies

section_content: Amplitude dispersion index and coherence are widely used in GBSAR interferometric analysis as tools for selecting high quality pixels in the image. Since the coherence evaluation demands preprocessing analysis to give reliable results, the amplitude dispersion index could represent a precious efficient tool. Starting from simple general assumptions, we found a theoretical relation between the two parameters. 

The purpose of the experimental data analysis was to determine whether such theoretical relation would apply for GBSAR images acquired in scenarios with different characteristics. The three measurement campaigns considered for the analysis were performed in the following environments: on a slope characterized by the presence of rocks and partially covered by grass; on a dam, covered by low vegetation; on an urban scenario, characterized by the presence of both buildings and vegetation. 

For each data set, the agreement between experimental data and Rice statistics was tested for several areas in the wide range of scenarios considered. Excellent agreement was found for areas characterized by the presence of meadows. Then, both amplitude dispersion index and coherence were evaluated for all these areas. Results reported in Figures 5, 7, 8, 9, 11, 12 and 13 seem to confirm the relation between amplitude dispersion index and coherence retrieved analytically in Equation ( 16 ). Specifically, good agreement is found for low values of coherence (high values of D ), that is, for targets characterized by high signal-to-noise ratio. 

Indeed, as can be noted in Figures 7 8 9 and Figures 11 12 13, some points close to γ = 0, overcome the Rayleigh threshold (D ≃ 0. 52). These results were obtained by performing the analysis for the whole dam body (Figures 7 8 9 ) and for the entire focused image in the urban scenario (Figures 11 12 13 ). The points with γ ≃ 0, probably correspond to the targets of wooded areas, or in any case where there is more incoherent scatterers, and low signalto-noise ratio. 

Anyway, results obtained using data acquired at Formigal (Figure 5 ) in C-band, also show good agreement with the theoretical ones for low values of coherence. This is probably because the analysis has been carried out for a homogeneous area characterized by high signal-to-noise ratio. 

The results shown suggest the use of this relation as a quantitative reference to relate coherence and amplitude dispersion index, for targets with D < 0. 5, that is, with a sufficiently high signal-to-noise ratio. This method can represent a powerful tool when
section_name: Discussion

section_content: Amplitude dispersion index and coherence are widely used in GBSAR interferometric analysis as tools for selecting high quality pixels in the image. Since the coherence evaluation demands preprocessing analysis to give reliable results, the amplitude dispersion index could represent a precious efficient tool. Starting from simple general assumptions, we found a theoretical relation between the two parameters. 

The purpose of the experimental data analysis was to determine whether such theoretical relation would apply for GBSAR images acquired in scenarios with different characteristics. The three measurement campaigns considered for the analysis were performed in the following environments: on a slope characterized by the presence of rocks and partially covered by grass; on a dam, covered by low vegetation; on an urban scenario, characterized by the presence of both buildings and vegetation. 

For each data set, the agreement between experimental data and Rice statistics was tested for several areas in the wide range of scenarios considered. Excellent agreement was found for areas characterized by the presence of meadows. Then, both amplitude dispersion index and coherence were evaluated for all these areas. Results reported in Figures 5, 7-9 and 11-13 seem to confirm the relation between amplitude dispersion index and coherence retrieved analytically in Equation ( 16 ). Specifically, good agreement is found for low values of coherence (high values of D A ), that is, for targets characterized by high signal-to-noise ratio. 

Indeed, as can be noted in Figures 7 8 9 and 11-13, some points close to γ = 0, overcome the Rayleigh threshold (D A 0. 52). These results were obtained by performing the analysis for the whole dam body (Figures 7 8 9 ) and for the entire focused image in the urban scenario (Figures 11 12 13 ). The points with γ 0, probably correspond to the targets of wooded areas, or in any case where there is more incoherent scatterers, and low signal-to-noise ratio. 

Anyway, results obtained using data acquired at Formigal (Figure 5 ) in C-band, also show good agreement with the theoretical ones for low values of coherence. This is probably because the analysis has been carried out for a homogeneous area characterized by high signal-to-noise ratio. 

The results shown suggest the use of this relation as a quantitative reference to relate coherence and amplitude dispersion index, for targets with D A < 0. 5, that is, with a sufficiently high signal-to-noise ratio. This method can represent a powerful tool when
section_name: Discussion

section_name: Acknowledgments:

section_content: The authors would like to thank DIAN s. r. l., Matera, Italy, especially Giovanni Nico and Olimpia Masci, for their valuable support and cooperation.
section_name: Acknowledgments:

section_content: Data Availability Statement: Data available on request.

section_name: Conclusions

section_content: searching for high-signal quality pixels in GBSAR images, in a wide range of scenarios. The validation with different frequency bands also contributes to give solid basis to the relation.

section_content: In this paper, the estimation of the coherence of GBSAR images has been addressed. First, we reviewed some common estimators used to evaluate the phase stability of pixels in the radar image. Then, a statistical description of the backscattered signal amplitude is presented, and a relationship between the so-called amplitude dispersion index and the interferometric coherence was analytically found. 

The preliminary estimation of the targets signal quality usually performed in interferometric analysis, could benefit from this relation. Indeed, the amplitude dispersion index provide practical advantages, and represent an efficient tool for estimating the signal quality. A theoretical quantitative relation between coherence and amplitude dispersion index can give solid basis to the D A parameter: it relates a threshold on coherence to thresholds on D A. Thus, such a relation allows the obtaining of information on the coherence of the image, greatly simplifying the processing. 

Some simulations have been implemented to test the reliability of such a function under specific hypotheses. The theoretical relation confirmed by experimental data acquired in different environments and with different working frequencies. Obtained results suggest the general character of this quantitative relation and its applicability in a wide range of situations.
section_name: Conclusions

section_content: The authors declare no conflict of interest.
section_name: Conflicts of Interest: