A TSVD-Based Method for Forest Height Inversion from Single-Baseline PolInSAR Data

The random volume over ground (RVoG) model associates vegetation vertical structure parameters with multiple complex interferometric coherence observables. In this paper, on the basis of the RVoG model, a truncated singular value decomposition (TSVD)-based method is proposed for forest height inversion from single-baseline polarimetric interferometric synthetic aperture radar (PolInSAR) data. In addition, in order to improve the applicability of TSVD for this issue, a new truncation method is proposed for TSVD. Differing from the traditional three-stage method, the TSVD-based inversion method estimates the pure volume coherence directly from the complex interferometric coherence, and estimates the forest height from the estimated pure volume coherence with a least-squares method. As a result, the TSVD-based method can adjust the contributions of the polarizations in the estimation of the model parameters and avoid the null ground-to-volume ratio assumption. The simulated experiments undertaken in this study confirmed that the TSVD-based method performs better than the three-stage method in forest height inversion. The TSVD-based method was also applied to E-SAR P-band data acquired over the Krycklan Catchment, Sweden, which is covered with mixed pine forest. The results showed that the TSVD-based method improves the root-mean-square error by 48.6% when compared to the three-stage method, which further validates the performance of the TSVD-based method.


Introduction
It is well known that vegetation height plays an important role in quantifying the terrestrial carbon cycle [1,2].Moreover, vegetation height is an essential factor for the estimation of the biomass stored in vegetation [3,4].Therefore, accurately extracting vegetation height at a large scale is an important task.Given the fact that polarimetric interferometric synthetic aperture radar (PolInSAR) can separate the scattering power of a single resolution cell into the contributions of surface, double-bounce, and volume scattering, it can be considered to be a viable remote sensing technique for estimating vegetation height in large-scale areas [4][5][6][7][8][9].
PolInSAR can be used to extract vegetation height through its sensitivity to the vegetation vertical structure [10][11][12][13][14].The complex interferometric coherence of the observed PolInSAR data has been related to the vertical distribution of the vegetation scattering [10][11][12].In a number of PolInSAR campaigns, the random volume over ground (RVoG) model [6,10] has been used to extract vegetation height from the complex interferometric coherence [15][16][17].The RVoG model is a physical PolInSAR model that integrates the complex coherence and the biophysical parameters.Based on the RVoG model, Papathanassiou [6] proposed six-dimensional nonlinear optimization method, which has been successfully evaluated with different types of PolInSAR data [6,15,16].However, for this method, the accuracy of the solution is greatly dependent on the selected initial value, and a poor choice of initial value can result in unstable parameter estimation.Furthermore, the iterative procedure used in this method consumes too much time and is unsuitable for the inversion of large-scale areas.To cope with this problem, Cloude [18] separated the six-dimensional nonlinear parameter optimization process into three stages.This solution is known as the "three-stage method", and has been widely used in vegetation height extraction for its simple, universal, and time-saving properties [15,19].However, the three-stage method assumes that there is at least one polarization channel without a ground scattering contribution (the null ground-to-volume ratio (GVR) assumption), and thus it is difficult to determine the polarization because the volume and ground scattering contributions are always mixed in all of the polarization channels due to the diverse penetration depths [20].As a result, the pure volume coherence estimated by the three-stage method is often inaccurate.Furthermore, with the three-stage method, it is not possible to adjust the contributions of the interferometric coherence observations in the estimation of the RVoG model parameters.
The aim of this work is to address the limitations described above and extract accurate vegetation heights from single-baseline PolInSAR data.Based on the RVoG model, linear observation equations are first developed from the complex interferometric coherence observations by Taylor expansion, in order to combine all the available interferometric coherence observations.Next, since the coefficient matrix derived from the linear equations is ill-conditioned, in order to overcome this ill-posed problem, the proposed truncated singular value decomposition (TSVD)-based method is used to estimate the pure volume coherence.Furthermore, due to the fact that the ordinary truncation method for TSVD is not suitable for this issue, a more adaptive truncation method is proposed so as to improve the accuracy of the estimated pure volume coherence.Finally, the forest height is extracted from the estimated pure volume coherence by the use of a least-squares method.
This paper is structured as follows.The principle of the RVoG model is introduced and discussed in Section 2. Section 3 presents the TSVD-based method for the estimation of vegetation height from complex interferometric coherence observations.The validation of the experiments are presented in Section 4. The discussions are presented in Section 5. Finally, the conclusions are drawn in Section 6.

RVoG Model
The RVoG model is a physical model that associates vegetation vertical structure parameters with multiple complex interferometric coherence observables [10].It is a basic and popular model for describing vegetation scenarios.The model depicts the vegetation layer as a volume with randomly oriented particles over an impenetrable ground surface.Without considering the temporal decorrelation, the complex interferometric coherence γ(ω) is expressed as [6]: where ω is the unitary polarization vector that defines the choice of scattering mechanism, ϕ 0 denotes the ground surface phase, µ(ω) represents the ground-to-volume ratio (GVR) accounting for the polarization diversity, which varies with ω, and γ v denotes the pure volume coherence, which is linked to the vegetation height and is expressed as [6]: where σ denotes the mean extinction coefficient, h v denotes the vegetation height, θ is the incidence angle, and k z is the vertical wave number, which depends on the imaging geometry and wavelength [21].k z is given as: where λ denotes the wavelength, and ∆θ represents the incidence angle difference between the master and slave images.
According to Equation (1), the theoretical loci of the complex coherence sets of different polarizations follow a straight line in the complex plane, as shown in Figure 1 [18].The figure also depicts the determination of the ground surface phase, which is applied in the three-stage method.However, the complex coherence observations no longer follows a common straight line in practice, due to the coherence fluctuations caused by all the possible decorrelations.In order to reconstruct the straight line accurately from noisy coherence sets, a line-fit approach based on least squares is proposed by the use of the following regression model [10,18]: where Re() and Im() denote the real and imaginary operations, respectively, c is the slope of the coherence line, and d is the intersection point to the imaginary axis.Once more than two complex coherence observations are provided, the coherence line can be determined by the line-fit approach.
The ground phase ϕ 0 can then be identified by the two intersections (ϕ 1 and ϕ 2 ) of the coherence line and the unit circle [18].However, the least-squares criterion may be unstable if the ellipticity of the coherence sets is high or if the complex coherence observations are too discrete, due to the noise of the polarizations.According to Equation (1), the theoretical loci of the complex coherence sets of different polarizations follow a straight line in the complex plane, as shown in Figure 1 [18].The figure also depicts the determination of the ground surface phase, which is applied in the three-stage method.However, the complex coherence observations no longer follows a common straight line in practice, due to the coherence fluctuations caused by all the possible decorrelations.In order to reconstruct the straight line accurately from noisy coherence sets, a line-fit approach based on least squares is proposed by the use of the following regression model [10,18]: where Re() and Im() denote the real and imaginary operations, respectively, is the slope of the coherence line, and is the intersection point to the imaginary axis.Once more than two complex coherence observations are provided, the coherence line can be determined by the line-fit approach.The ground phase can then be identified by the two intersections ( and ) of the coherence line and the unit circle [18].However, the least-squares criterion may be unstable if the ellipticity of the coherence sets is high or if the complex coherence observations are too discrete, due to the noise of the polarizations.
In addition, the three-stage method assumes that there is a polarization channel without ground scattering contribution, i.e., ( ) = 0, and then uses this assumption to further estimate the pure volume coherence.However, the volume and ground scattering contributions are always mixed in the polarization channels due to the diverse penetration depths.Therefore, it is difficult to find a polarization that can fit this hypothesis.As a result, when relying on this assumption, the estimated pure volume coherence will bias the extracted forest height.As shown in Equation ( 1), ( ) is polarization-dependent.The complex coherence observations of the different scattering mechanisms can be obtained by the use of different unitary polarization vectors .In other words, if the number of polarization vectors is m, then m equations like Equation ( 1) can be established.By using the least-squares criterion, the unknown parameters in ( ) can be estimated when the number of parameters is less than that of the equations.Therefore, the least-squares method is a viable way to estimate vegetation height from complex coherence observations.The least-squares method directly estimates the RVoG model parameters from the complex coherence observations without any assumptions, and it is able to adjust the contributions of the coherence observations in the estimation of the model parameters.

The TSVD-Based Method for the Estimation of Vegetation Height From Complex Interferometric Coherence Observations
In this section, a novel approach is introduced for the estimation of vegetation height from single-baseline PolInSAR data, on the basis of the least-squares criterion and the TSVD-based method.In addition, the three-stage method assumes that there is a polarization channel without ground scattering contribution, i.e., µ(ω) = 0, and then uses this assumption to further estimate the pure volume coherence.However, the volume and ground scattering contributions are always mixed in the polarization channels due to the diverse penetration depths.Therefore, it is difficult to find a polarization that can fit this hypothesis.As a result, when relying on this assumption, the estimated pure volume coherence will bias the extracted forest height.
As shown in Equation (1), γ(ω) is polarization-dependent.The complex coherence observations of the different scattering mechanisms can be obtained by the use of different unitary polarization vectors ω.In other words, if the number of polarization vectors is m, then m equations like Equation (1) can be established.By using the least-squares criterion, the unknown parameters in γ(ω) can be estimated when the number of parameters is less than that of the equations.Therefore, the least-squares method is a viable way to estimate vegetation height from complex coherence observations.The least-squares method directly estimates the RVoG model parameters from the complex coherence observations without any assumptions, and it is able to adjust the contributions of the coherence observations in the estimation of the model parameters.

The TSVD-Based Method for the Estimation of Vegetation Height From Complex Interferometric Coherence Observations
In this section, a novel approach is introduced for the estimation of vegetation height from single-baseline PolInSAR data, on the basis of the least-squares criterion and the TSVD-based method.

Estimation of Pure Volume Coherence from the Complex Interferometric Coherence Observations
From Equations ( 1) and ( 2), it can be seen that γ(ω) and γ v are complex numbers.Therefore, the equations can be separated into two parts: a real part and an imaginary part.In order to simplify the nonlinear equations, we parameterize the pure volume coherence by γ v = a + bi, and then the complex interferometric coherence can be given by: and the observational function for the least-squares criterion can be formulated as: where ω j is the j-th unitary polarization vector, f is the RVoG model function, as described in Equation ( 5), µ ω j is the j-th GVR associated with vector ω j , and m represents the number of polarization projection vectors.It is clear that there are four unknown model parameters in each function, and m + 3 unknown parameters in m functions.
To estimate the unknown model parameters, least-squares estimation is adopted to adjust the contributions of the complex coherence observations and suppress noise [20,22,23].The corresponding least-squares criterion can be formulated as: where ∑ represents the summation operation, |•| represents the modulus operation, and γ ω j is the estimation of γ ω j .Since the complex coherence can be separated into a real part and an imaginary part, Equation ( 7) can be equally converted to: In this way, the residual functions for the least-squares criterion are formulated as: where V j Re denotes the residual of the real part corresponding to the j-th unitary polarization vector, and V j Im denotes the residual of the imaginary part.From Equation (9), it can be seen that, if the number of polarization projection vectors is m, then 2m residual equations can be established.The number of unknown parameters is m + 3. Therefore, the unknowns can be estimated by least-squares estimation when the number of polarization projection vectors is more than three.
From Equation (5), it can be seen that the complex coherence function is highly nonlinear, which greatly affects the efficacy of the least-squares estimation.In view of this, a linearized strategy based on a Taylor series is adopted to convert the nonlinear function into a linear function [23].In this case, the residual functions for the least-squares criterion can be rewritten as: (10) where ∂ represents the partial derivative operation, dϕ 0 , da, db, and dµ ω j denote the corrections of the approximations of the unknown parameters, and l j Re and l j Im are the real and imaginary differences between the observed complex coherence and the predicted initial value.The vector V denotes the residuals, X denotes the unknown corrections, A denotes the coefficient matrix, and L denotes the real and imaginary differences.The residual functions for the least-squares criterion can then be expressed as [23,24]: The unitary polarization vector can be constructed by linear-basis polarization, Pauli-basis polarization, magnitude diversity optimization polarization [6], and phase diversity polarization [25].Therefore, more than three polarization projection vectors can be obtained for the estimation of the unknown parameters.Using the least-squares criterion, an estimate of the unknown corrections can be given by: where X denotes the estimation of the unknown corrections.The estimation of the model parameters can then be obtained by combining the approximations and the corrections.However, during the computation, we find that matrix A is seriously ill-conditioned.It is well known that an ill-posed problem is a great hindrance to obtaining accurate parameter estimation from observation equations.
For an ill-posed equation, a small amount of noise in the observations can often bring large uncertainties to the estimation.The detrimental effect of the ill-posed problem is reflected in the large variance of the least-squares estimations.The variance-covariance matrix of the least-squares estimations can be given by [23]: where Cov X represents the variance-covariance matrix of the estimations, and the diagonal elements of the matrix are the variances of the estimations, σ 2 0 denotes the unit weight variance.We perform singular value decomposition (SVD) on matrix A [24]: where U is a 2m × 2m orthogonal matrix of the left singular vectors, G is an n × n orthogonal matrix of the right singular vectors, and n denotes the number of unknowns, which is equivalent to m + 3, S is the matrix of singular values, and λ are non-negative real numbers that are conventionally listed in decreasing order, i.e., Using the trace operation, the sum of the variances of the estimations can be given by: where Vars X represents the sum of the variances of the estimations, and Trace denotes the trace operation.If the equation is ill-posed, the singular values gradually decrease to zero, and λ 1 is much larger than λ n (λ n is close to zero).Equation (16) shows that small singular values greatly magnify the estimation variances.Thus, the least-squares estimation becomes highly unreliable and is unable to obtain accurate estimations of the parameters.
In order to overcome the ill-posedness of the problem, truncated singular value decomposition (TSVD) [26][27][28] is adopted.TSVD truncates the small singular values that greatly enlarge the variances to improve the least-squares estimation.Using the TSVD-based method, the estimation of the unknown parameters is given by: where Xt denotes the improved estimation by TSVD, S t denotes the inverse singular value matrix which truncates the small singular values, S P is the singular value matrix which represents the preserved large singular values, and S T is a zero matrix that represents the truncated small singular values.We can then obtain the variance of the TSVD estimations as follows: where Vars Xt represents the variance of the TSVD estimations.k denotes the number of preserved large singular values.Equation (19) shows that TSVD greatly reduces the variance of the least-squares estimations through truncating the small singular values.However, it is well known that TSVD results in a biased estimation.Truncating the small singular values not only reduces the variance but also introduces bias into the estimation.TSVD improves the estimation by reducing the mean-square error (MSE) of the least-squares estimation.The MSE is expressed as [26]: where Mse Xt represents the MSE of the TSVD estimation, and Bias Xt represents the bias introduced by TSVD.The sum of squares of Bias Xt is given by: where X represents the true values of the unknown parameters, and G i denotes the i-th right singular vector that corresponds to the i-th singular value.Equation (21) shows that the bias is introduced by truncating the small singular values.The more singular values that are truncated, the more bias is introduced.Considering Equation (20), it is clear that TSVD reduces the MSE by truncating the small singular values.However, the reduction of the MSE relies on the reduced variance being more than the introduced bias.Therefore, the truncation parameter k which determines the preserved large singular values and truncated small singular values is a key factor for TSVD to reduce the MSE.The method that is commonly used for the determination of the truncation parameter is related to the condition number, i.e., λ 1 /λ i , as defined in Equation (15).If the condition number is larger than the given upper limit, the singular value λ i should be truncated [28][29][30].However, large-scale PolInSAR data usually consist of millions of pixels, and the singular values in ill-conditioned matrices of pixels are different from each other.Therefore, it is difficult to determine a reasonable upper limit for the condition number.In order to determine a reasonable truncation parameter for this issue, we need to develop a new approach.From Equations ( 19) and ( 21), it can be seen that if singular value λ i is truncated, the reduced variances are σ 2 0 /λ 2 i and the introduced bias is the singular value should be truncated, and if σ 2 0 /λ 2 i < X T G i G T i X, the singular value should be preserved.In order to determine the small singular values, the values of σ 2 0 /λ 2 i and X T G i G T i X need to be calculated accurately.From the least-squares estimation, the estimation of is given by [23]: where U denotes the 2m × n dimensional left singular vectors, σ2 0 is the estimation of σ 2 0 [23,26], and n is the number of unknowns.It is clear that the small singular values have no adverse effect on the estimation of σ 2 0 .Therefore, σ 2 0 /λ 2 i can be calculated by the least-squares solution if 2m > n [23,26].In this paper, 10 polarizations are selected for the estimation of the pure volume coherence.Therefore, σ 2 0 can be estimated.
Using the SVD matrices of A, the estimation of G T i X by least squares can be expressed as: It can be seen from Equation ( 23) that the singular values greatly affect the estimation of G T i X.Small singular values seriously enlarge the value of G T i X, whereas large singular values do not.The variance of the estimation of G T i X can be given by: where Var g represents the variance of estimation G T i X.The equation shows that the variance of G T i X is negatively correlated with the singular value.Therefore, the variance of G T i X which corresponds to a large singular value is small.Furthermore, the estimation of G T i X is reliable due to the small variance.However, the estimation of G T i X, which corresponds to a small singular value, is unreliable due to the large variance.
Generally, if the standard deviation of G T i X is less than 3 σ0 , the estimation is considered to be reliable.Therefore, the reliable estimations with a standard deviation of less than 3 σ0 can be given by Since the values of the reduced variances σ2 0 /λ 2 i change as σ2 i is bigger than 90% of the values of J, it can be considered that σ 2 0 /λ 2 i > X T G i G T i X, and singular value λ i needs to be truncated.The other small singular values λ r that are smaller than λ i can also be denoted as σ 2 0 /λ 2 r > X T G r G T r X and need to be truncated.The small singular values which should be truncated in the TSVD-based method are thus determined.The ill-posed problem can be well solved by TSVD with the proposed truncation method.Finally, the pure volume coherence can be estimated by the proposed TSVD-based method.

Extraction of Vegetation Height from the Pure Volume Coherence
Since the pure volume coherence is parameterized by γ v = a + bi, the pure volume coherence is estimated as γv = â + bi by the proposed TSVD-based method.Using the estimated parameters â and b, the pure volume coherence which links to the vegetation height can be expressed as: The equation can then be separated into a real part and an imaginary part: From Equation ( 25), it can be seen that θ and k z are the known parameters, and the unknown parameters are σ and h v .Therefore, σ and h v can be estimated by the least-squares estimation from Equation (26).The Taylor series is used to convert the nonlinear function into a linear function.The residual functions for the least-squares criterion can then be expressed as [23]: where V p Re denotes the residual of the real part, and V p Im denotes the residual of the imaginary part, dσ and dh v denote the corrections of the approximations of the unknown parameters, and l p Re and l p Im are the real and imaginary differences between the estimated pure volume coherence and the predicted initial value.Using the vector V p to denote the residuals, X p denotes the unknown corrections, A p denotes the coefficient matrix, and L p denotes the real and imaginary differences.The residual functions for the least-squares criterion can then be expressed as [23,24,29]: Since the number of residual equations is the same as that of the unknowns, the least-squares criterion can be used to estimate the unknowns.The estimation is given by: where Xp denotes the estimation of the unknowns.The vegetation height can then be obtained from the estimation of dh v and the initial value of h v by: where ĥv denotes the estimated vegetation height, and h v0 denotes the initial value of h v .

The Determination of Initial Values of Model Parameters
From Equations ( 10) and ( 27), it can be seen that the proposed method is plagued by the initial value.If the initial value cannot be well determined, it is difficult to get reliable estimation.The classical six-dimensional nonlinear optimization method [6] is also confronted by this problem since it is difficult to obtain the priori information of the forest parameters (forest height and extinction).However, it is easy to obtain reliable initial values of the ground phase, the pure volume coherence and the ground-to-volume ratio according to the RVoG assumption and its linear geometrical expression in the complex plane [10].Based on this, Equation ( 10) is used to estimate the pure volume coherence, which is important to invert the forest parameters.Compared to the six-dimensional nonlinear optimization method, this method can avoid significant biases caused by the unreliable initial values of forest parameters.Moreover, in comparison to the three-stage method, the proposed method can provide more accurate pure volume coherence because it is free from the assumption that there is one polarization whose ground-to-volume power ratio should be less than −10 dB [10], which cannot be fulfilled for the low-frequency PolInSAR data or for the sparse forest [31].In addition, only one polarimetric observation is used to calculate the pure volume coherence in the three-stage method.However, multi-polarization observations are used to estimate the pure volume coherence under the TSVD based least-squares estimation framework, which enhances the method's ability to alleviate the effect of observational errors in the estimations.As a result, the obtained pure volume coherence can support more accurate forest parameter estimation.
Then, with the estimated pure volume coherence, Equation ( 27) is used to estimate forest height.The initial values for the forest height and the mean extinction coefficient are determined by the three-stage method.Although the three-stage method cannot give satisfactory results, especially for the low-frequency PolInSAR data or for sparse forest, these estimated results can be regarded as the initial values for Equation (27).The least-squares method estimates the corrections of initial values from the reliable pure volume coherence to improve the accuracy of the estimated unknown parameters.The final experiment has shown that with the accurate pure volume coherence and reasonable initial values, the proposed method can attain good forest height results.

Simulated Experiments
In order to evaluate the performance of the proposed TSVD-based inversion method for the estimation of vegetation height from single-baseline PolInSAR data, we simulated single-baseline PolInSAR data through the PolSARpro tool released by European Space Agency (ESA), Rome, Italy, using the following forest scenario, as shown in Table 1: simulated broad-leaved forest with an average height of 18 m; and ground phase of 0 degrees, corresponding to a ground elevation of 0 m.The complex interferometric coherence of the linear-basis polarizations (transmitted polarization and received polarization are horizontal polarizations, HH, transmitted polarization and received polarization are vertical polarizations, VV) and the Pauli-basis polarization (transmitted polarization is horizontal polarization and received polarization is vertical polarization, HV) could then be obtained.The simulated forest area and the Pauli RGB composite image are shown in Figure 2. Furthermore, using Pauli-basis polarization, magnitude diversity optimization polarization [6], and phase diversity polarization [25], the complex interferometric coherence observations of HH+VV, HH-VV, opt1, opt2, opt3, phase diversity (PD)-high, and phase diversity (PD)-low could be obtained.Therefore, 20 residual equations could be established from the above 10 complex interferometric coherence observations based on Equation (10) with 13 unknown parameters.The unknown parameters were then estimated by the TSVD-based method.From the estimated parameters and the pure volume coherence model which links to the forest height, the forest height could be extracted by Equation (29).The extraction results are shown in Figure 3b.For comparison, the extraction results of the three-stage method are shown in Figure 3a.Furthermore, using Pauli-basis polarization, magnitude diversity optimization polarization [6], and phase diversity polarization [25], the complex interferometric coherence observations of HH+VV, HH-VV, opt1, opt2, opt3, phase diversity (PD)-high, and phase diversity (PD)-low could be obtained.Therefore, 20 residual equations could be established from the above 10 complex interferometric coherence observations based on Equation (10) with 13 unknown parameters.The unknown parameters were then estimated by the TSVD-based method.From the estimated parameters and the pure volume coherence model which links to the forest height, the forest height could be extracted by Equation (29).The extraction results are shown in Figure 3b.For comparison, the extraction results of the three-stage method are shown in Figure 3a.
From Figure 3, it can be seen that the TSVD-based method performs much better than the three-stage method in forest height inversion in this test.Clearly, from Figure 3a, the three-stage method fails to invert the forest heights of the rectangular areas of the figure, but the TSVD-based method effectively improves the inversion of the rectangular areas, as can be seen in Figure 3b.The mean values of the extracted forest heights by the three-stage method and the TSVD-based method are 13.383 m and 15.4440 m.This demonstrates that the forest height estimated by TSVD is more accurate than that estimated by the three-stage method.For a further comparison, 16 forest stands were selected from the simulated forest area, and the root-mean-square error (RMSE) of each stand was adopted to compare the performance of the two methods.The RMSEs of each stand are shown in Figure 4. Furthermore, using Pauli-basis polarization, magnitude diversity optimization polarization [6], and phase diversity polarization [25], the complex interferometric coherence observations of HH+VV, HH-VV, opt1, opt2, opt3, phase diversity (PD)-high, and phase diversity (PD)-low could be obtained.Therefore, 20 residual equations could be established from the above 10 complex interferometric coherence observations based on Equation (10) with 13 unknown parameters.The unknown parameters were then estimated by the TSVD-based method.From the estimated parameters and the pure volume coherence model which links to the forest height, the forest height could be extracted by Equation (29).The extraction results are shown in Figure 3b.For comparison, the extraction results of the three-stage method are shown in Figure 3a.From Figure 3, it can be seen that the TSVD-based method performs much better than the threestage method in forest height inversion in this test.Clearly, from Figure 3a, the three-stage method fails to invert the forest heights of the rectangular areas of the figure, but the TSVD-based method effectively improves the inversion of the rectangular areas, as can be seen in Figure 3b.The mean values of the extracted forest heights by the three-stage method and the TSVD-based method are 13.383 m and 15.4440 m.This demonstrates that the forest height estimated by TSVD is more accurate than that estimated by the three-stage method.For a further comparison, 16 forest stands were selected from the simulated forest area, and the root-mean-square error (RMSE) of each stand was adopted to compare the performance of the two methods.The RMSEs of each stand are shown in Figure 4.   From Figure 3, it can be seen that the TSVD-based method performs much better than the threestage method in forest height inversion in this test.Clearly, from Figure 3a, the three-stage method fails to invert the forest heights of the rectangular areas of the figure, but the TSVD-based method effectively improves the inversion of the rectangular areas, as can be seen in Figure 3b.The mean values of the extracted forest heights by the three-stage method and the TSVD-based method are 13.383 m and 15.4440 m.This demonstrates that the forest height estimated by TSVD is more accurate than that estimated by the three-stage method.For a further comparison, 16 forest stands were selected from the simulated forest area, and the root-mean-square error (RMSE) of each stand was adopted to compare the performance of the two methods.The RMSEs of each stand are shown in Figure 4.For the 16 forest stands, the RMSE of the TSVD-based method is consistently less than that of the three-stage method.This indicates that the inverted forest height obtained by TSVD is closer to the true height than the height obtained by the three-stage method.This further validates the performance of the TSVD-based method.

Study Area and Data Sets
The proposed TSVD-based method was also applied to E-SAR P-band PolInSAR data, which were collected under the framework of the BioSAR 2008 campaign by the German Aerospace Center, Munich, Germany.The test site is a forest area within the Krycklan River catchment in Northern Sweden, and is mainly covered by mixed boreal forest with heights ranging from 0 to 35 m.The topography elevation is between 150 to 380 m above mean sea level (AMSL).The baseline PolInSAR data were acquired in the repeat-pass configuration.Moreover, as part of the BioSAR2008 campaign, a light detection and ranging (LiDAR) measurement was also obtained by the Swedish Defense Research Agency (FOI).The derived forest height is regarded as the reference in the following analysis.
In this experiment, the forest heights were extracted from single-baseline data.The temporal and spatial baselines were 70 min and 32 m, respectively.The vertical wavenumber ranged from 0.051 to 0.181.The Pauli-basis RGB composite intensity image for the test site is shown in Figure 5.
campaign, a light detection and ranging (LiDAR) measurement was also obtained by the Swedish Defense Research Agency (FOI).The derived forest height is regarded as the reference in the following analysis.
In this experiment, the forest heights were extracted from single-baseline data.The temporal and spatial baselines were 70 min and 32 m, respectively.The vertical wavenumber ranged from 0.051 to 0.181.The Pauli-basis RGB composite intensity image for the test site is shown in Figure 5.

Forest Height Inversion
Following the steps of the TSVD-based inversion method, the complex interferometric coherence observations of HH, VV, HV, HH+VV, HH-VV, opt1, opt2, opt3, PDhigh, and PDlow were used to establish the residual equations.The pure volume coherence was then estimated by the TSVD-based method from the residual equations and, finally, the forest heights were extracted from the pure volume coherence model using the estimated model parameters.The extracted forest heights are shown in Figure 6b.As in the simulated experiments, the three-stage method was also used to extract the forest heights, and the inversion results are shown in Figure 6a. Figure 6c shows the LiDAR forest heights used as a reference.

Forest Height Inversion
Following the steps of the TSVD-based inversion method, the complex interferometric coherence observations of HH, VV, HV, HH+VV, HH-VV, opt1, opt2, opt3, PDhigh, and PDlow were used to establish the residual equations.The pure volume coherence was then estimated by the TSVD-based method from the residual equations and, finally, the forest heights were extracted from the pure volume coherence model using the estimated model parameters.The extracted forest heights are shown in Figure 6b.As in the simulated experiments, the three-stage method was also used to extract the forest heights, and the inversion results are shown in Figure 6a. Figure 6c shows the LiDAR forest heights used as a reference.It can be seen from Figure 6 that the inversion results of the three-stage method and the TSVD-based method follow a similar spatial trend, but significant differences are also apparent.Compared to the LiDAR results, the inverted forest heights obtained by TSVD are clearly more accurate than those of the three-stage method.In order to analyze the differences, 272 forest stands characterized by nearly uniform tree heights were selected from the LiDAR results.We took the estimated forest height average for every stand and computed the difference between that and the LiDAR forest height.The RMSE and correlation coefficient (R 2 ) were calculated to validate the performance.The validated stand-level plots are displayed in Figure 7.It can be seen from Figure 6 that the inversion results of the three-stage method and the TSVD-based method follow a similar spatial trend, but significant differences are also apparent.Compared to the LiDAR results, the inverted forest heights obtained by TSVD are clearly more accurate than those of the three-stage method.In order to analyze the differences, 272 forest stands characterized by nearly uniform tree heights were selected from the LiDAR results.We took the estimated forest height average for every stand and computed the difference between that and the LiDAR forest height.The RMSE and correlation coefficient (R 2 ) were calculated to validate the performance.The validated stand-level plots are displayed in Figure 7.
The validated plots of the three-stage method and the TSVD-based method are characterized by R 2 values of 0.2166 and 0.5824, respectively.This indicates that the forest heights inverted by TSVD are closer to the LiDAR forest heights.The RMSEs of the three-stage method and TSVD are 6.6351 and 3.4096, respectively.Clearly, the inversion accuracy of TSVD is higher than that of the three-stage method, showing an improvement of 48.6%.Therefore, it is possible to state that the TSVD-based method can improve the inversion of the forest height in this test site.
Compared to the LiDAR results, the inverted forest heights obtained by TSVD are clearly more accurate than those of the three-stage method.In order to analyze the differences, 272 forest stands characterized by nearly uniform tree heights were selected from the LiDAR results.We took the estimated forest height average for every stand and computed the difference between that and the LiDAR forest height.The RMSE and correlation coefficient (R 2 ) were calculated to validate the performance.The validated stand-level plots are displayed in Figure 7.The validated plots of the three-stage method and the TSVD-based method are characterized by R 2 values of 0.2166 and 0.5824, respectively.This indicates that the forest heights inverted by TSVD are closer to the LiDAR forest heights.The RMSEs of the three-stage method and TSVD are 6.6351 and 3.4096, respectively.Clearly, the inversion accuracy of TSVD is higher than that of the three-stage method, showing an improvement of 48.6%.Therefore, it is possible to state that the TSVD-based method can improve the inversion of the forest height in this test site.

The Extracted Ground Surface Phases by the Three-Stage Method and TSVD
From Equation ( 10), it can be seen that the ground surface phase can also be estimated by the TSVD-based method.Since the ground surface phase plays an important part in the estimation of the

The Extracted Ground Surface Phases by the Three-Stage Method and TSVD
From Equation ( 10), it can be seen that the ground surface phase can also be estimated by the TSVD-based method.Since the ground surface phase plays an important part in the estimation of the underlying topography [30,32,33], the extracted ground surface phases obtained by the three-stage method and TSVD are shown in Figure 8a,b, respectively.underlying topography [30,32,33], the extracted ground surface phases obtained by the three-stage method and TSVD are shown in Figure 8a,b, respectively.It is difficult to see any difference between Figure 8a,b.This indicates that the three-stage method and TSVD perform similarly in extracting the ground surface phase.For a more in-depth analysis, the ground surface phase obtained by TSVD subtracted from that obtained by the three-stage method is shown in Figure 8c.It can be seen that most of the values in Figure 8c are close to zero.This further validates that TSVD is unable to improve the accuracy of the estimation of the ground surface phase and the underlying digital elevation model (DEM) [30].
The line-fit approach is used in the three-stage method to determine the ground surface phase.In order to reconstruct the straight line accurately from noisy coherence sets, 10 polarizations are used by least-squares based line-fit approach.Similarly, the same 10 polarizations are used in TSVD for the estimation of model parameters.Therefore, the observations used by line-fit approach and TSVD have the same noise and information, and the basic criterion in the line-fit approach and TSVD both are least-squares.This is a possible reason for this result.

Effects on Estimation of Phase Height
As mentioned in Section 3, the TSVD-based method first estimates the pure volume coherence from the complex interferometric coherence and then extracts the forest height.The phase height, which is very important for the inversion of forest height [34,35], was also computed from the It is difficult to see any difference between Figure 8a,b.This indicates that the three-stage method and TSVD perform similarly in extracting the ground surface phase.For a more in-depth analysis, the ground surface phase obtained by TSVD subtracted from that obtained by the three-stage method is shown in Figure 8c.It can be seen that most of the values in Figure 8c are close to zero.This further validates that TSVD is unable to improve the accuracy of the estimation of the ground surface phase and the underlying digital elevation model (DEM) [30].
The line-fit approach is used in the three-stage method to determine the ground surface phase.In order to reconstruct the straight line accurately from noisy coherence sets, 10 polarizations are used by least-squares based line-fit approach.Similarly, the same 10 polarizations are used in TSVD for the estimation of model parameters.Therefore, the observations used by line-fit approach and TSVD have the same noise and information, and the basic criterion in the line-fit approach and TSVD both are least-squares.This is a possible reason for this result.

Effects on Estimation of Phase Height
As mentioned in Section 3, the TSVD-based method first estimates the pure volume coherence from the complex interferometric coherence and then extracts the forest height.The phase height, which is very important for the inversion of forest height [34,35], was also computed from the estimated pure volume coherence and compared with the result of the three-stage method.The phase heights calculated from the estimated pure volume coherence obtained by the three-stage method and TSVD are shown in Figure 9a,b.From Figure 9a,b, it can be seen that the estimated phase heights in Figure 9a are higher than those in Figure 9b.We subtracted the phase heights in Figure 9b from those in Figure 9a and display the results in Figure 9c.Clearly, most of the values in Figure 9c are positive numbers.This indicates that the phase height extracted by TSVD is a better fit for the theoretical phase height of the pure volume coherence than the phase height extracted by the three-stage method.For a further comparison, we computed the average phase height of each forest stand and show the results in Figure 10. Figure 10 clearly indicates that the phase heights estimated by the TSVD-based method are consistently higher than those estimated by the three-stage method in each forest stand.Therefore, the forest heights inverted by TSVD are more accurate than those inverted by the three-stage method in the forest stands.To confirm this conclusion, we also used the least-squares method to extract the forest heights from the pure volume coherence estimated by the three-stage method, as mentioned in Section 3.2.The results are shown in Figure 11.From Figure 9a,b, it can be seen that the estimated phase heights in Figure 9a are higher than those in Figure 9b.We subtracted the phase heights in Figure 9b from those in Figure 9a and display the results in Figure 9c.Clearly, most of the values in Figure 9c are positive numbers.This indicates that the phase height extracted by TSVD is a better fit for the theoretical phase height of the pure volume coherence than the phase height extracted by the three-stage method.For a further comparison, we computed the average phase height of each forest stand and show the results in Figure 10.From Figure 9a,b, it can be seen that the estimated phase heights in Figure 9a are higher than those in Figure 9b.We subtracted the phase heights in Figure 9b from those in Figure 9a and display the results in Figure 9c.Clearly, most of the values in Figure 9c are positive numbers.This indicates that the phase height extracted by TSVD is a better fit for the theoretical phase height of the pure volume coherence than the phase height extracted by the three-stage method.For a further comparison, we computed the average phase height of each forest stand and show the results in Figure 10. Figure 10 clearly indicates that the phase heights estimated by the TSVD-based method are consistently higher than those estimated by the three-stage method in each forest stand.Therefore, the forest heights inverted by TSVD are more accurate than those inverted by the three-stage method in the forest stands.To confirm this conclusion, we also used the least-squares method to extract the forest heights from the pure volume coherence estimated by the three-stage method, as mentioned in Section 3.2.The results are shown in Figure 11. Figure 10 clearly indicates that the phase heights estimated by the TSVD-based method are consistently higher than those estimated by the three-stage method in each forest stand.Therefore, the forest heights inverted by TSVD are more accurate than those inverted by the three-stage method in the forest stands.To confirm this conclusion, we also used the least-squares method to extract the forest heights from the pure volume coherence estimated by the three-stage method, as mentioned in Section 3.2.The results are shown in Figure 11.Comparing Figures 11 and 5b, which displays the inversion results of TSVD, clearly, the forest heights in Figure 5b are closer to the LiDAR forest heights than those in Figure 11.Therefore, we can conclude that the TSVD-based method has the capacity to improve the estimation of the pure volume coherence.Based on the improved pure volume coherence, the forest height can be extracted more accurately.

Limitations of the TSVD-Based Method
Two problems are worth discussing.Firstly, TSVD plays an important role in the proposed solution.Since the ordinary truncation method is not suitable for this issue, a more adaptive truncation method is proposed in this paper.From Equation (23), it can be seen that sufficient polarizations are needed to estimate , i.e., 2 > [23,26].In this paper, 10 polarizations are selected for the forest height inversion.Therefore, 2 − = 7 confirms the accuracy of the estimation of [23].Secondly, due to the scattering mechanism, the inverted forest heights always follow a spatial trend, which can be seen in the results of both the three-stage method and TSVD.As a consequence, it is apparent that the far-range areas result in the overestimation of the forest height in Figure 6a,b, and especially in Figure 6b.

Conclusions
A TSVD-based method has been proposed in this paper for forest height inversion from single-baseline PolInSAR data.Differing from the traditional three-stage method, the new method estimates the pure volume coherence intuitively from the complex interferometric coherence, and has the capacity to adjust the contributions of the polarizations in the estimation of the model parameters.The TSVD-based method was first applied in forest height inversion from simulated PolInSAR data generated in PolSARpro.The results demonstrated that the TSVD-based method significantly improves the inversion results when compared to the three-stage method.This was also confirmed with airborne E-SAR P-band data obtained over a mixed boreal forest.The inverted forest heights obtained by TSVD showed an improvement in RMSE of 48.6% when compared to the results of the three-stage method.The phase heights of the estimated pure volume coherence were also well improved when compared to the results of the three-stage method.Comparing Figures 11 and 5b, which displays the inversion results of TSVD, clearly, the forest heights in Figure 5b are closer to the LiDAR forest heights than those in Figure 11.Therefore, we can conclude that the TSVD-based method has the capacity to improve the estimation of the pure volume coherence.Based on the improved pure volume coherence, the forest height can be extracted more accurately.

Limitations of the TSVD-Based Method
Two problems are worth discussing.Firstly, TSVD plays an important role in the proposed solution.Since the ordinary truncation method is not suitable for this issue, a more adaptive truncation method is proposed in this paper.From Equation (23), it can be seen that sufficient polarizations are needed to estimate σ 2 0 , i.e., 2m > n [23,26].In this paper, 10 polarizations are selected for the forest height inversion.Therefore, 2m − n = 7 confirms the accuracy of the estimation of σ 2 0 [23].Secondly, due to the scattering mechanism, the inverted forest heights always follow a spatial trend, which can be seen in the results of both the three-stage method and TSVD.As a consequence, it is apparent that the far-range areas result in the overestimation of the forest height in Figure 6a,b, and especially in Figure 6b.

Conclusions
A TSVD-based method has been proposed in this paper for forest height inversion from single-baseline PolInSAR data.Differing from the traditional three-stage method, the new method estimates the pure volume coherence intuitively from the complex interferometric coherence, and has the capacity to adjust the contributions of the polarizations in the estimation of the model parameters.The TSVD-based method was first applied in forest height inversion from simulated PolInSAR data generated in PolSARpro.The results demonstrated that the TSVD-based method significantly improves the inversion results when compared to the three-stage method.This was also confirmed with airborne E-SAR P-band data obtained over a mixed boreal forest.The inverted forest heights obtained by TSVD showed an improvement in RMSE of 48.6% when compared to the results of the three-stage method.The phase heights of the estimated pure volume coherence were also well improved when compared to the results of the three-stage method.

Figure 1 .
Figure 1.Geometrical interpretation of the coherence loci.

Figure 1 .
Figure 1.Geometrical interpretation of the coherence loci.

Figure 3 .
Figure 3. (a) forest height inversion result of the three-stage method; and (b) forest height inversion result of the truncated singular value decomposition (TSVD)-based inversion method.

Figure 3 .
Figure 3. (a) forest height inversion result of the three-stage method; and (b) forest height inversion result of the truncated singular value decomposition (TSVD)-based inversion method.

Figure 3 .
Figure 3. (a) forest height inversion result of the three-stage method; and (b) forest height inversion result of the truncated singular value decomposition (TSVD)-based inversion method.

16 Figure 6 .
Figure 6.(a) inversion results of the three-stage method; (b) inversion results of the TSVD-based method; (c) forest heights derived by light detection and ranging (LiDAR).

Figure 6 .
Figure 6.(a) inversion results of the three-stage method; (b) inversion results of the TSVD-based method; (c) forest heights derived by light detection and ranging (LiDAR).

Figure 7 .
Figure 7. (a) forest heights estimated by the three-stage method; (b) forest heights estimated by the TSVD-based method.

Figure 7 .
Figure 7. (a) forest heights estimated by the three-stage method; (b) forest heights estimated by the TSVD-based method.

Figure 8 .
Figure 8.(a) ground surface phase estimated by the three-stage method; (b) ground surface phase estimated by TSVD; (c) the difference between the ground surface phases obtained by the three-stage method and TSVD.

Figure 8 .
Figure 8.(a) ground surface phase estimated by the three-stage method; (b) ground surface phase estimated by TSVD; (c) the difference between the ground surface phases obtained by the three-stage method and TSVD.

Appl. 2017, 7 , 435 13 of 16 Figure 9 .
Figure 9. (a) the phase heights extracted from the pure volume coherence obtained by the three-stage method; (b) the phase heights extracted from the pure volume coherence obtained by TSVD; (c) the differences between the phase heights obtained by the three-stage method and TSVD.

Figure 10 .
Figure 10.The estimated phase heights for each forest stand.

Figure 9 .
Figure 9. (a) the phase heights extracted from the pure volume coherence obtained by the three-stage method; (b) the phase heights extracted from the pure volume coherence obtained by TSVD; (c) the differences between the phase heights obtained by the three-stage method and TSVD.

Appl. 2017, 7 , 435 13 of 16 Figure 9 .
Figure 9. (a) the phase heights extracted from the pure volume coherence obtained by the three-stage method; (b) the phase heights extracted from the pure volume coherence obtained by TSVD; (c) the differences between the phase heights obtained by the three-stage method and TSVD.

Figure 10 .
Figure 10.The estimated phase heights for each forest stand.

Figure 10 .
Figure 10.The estimated phase heights for each forest stand.

Appl. 2017, 7 , 435 14 of 16 Figure 11 .
Figure 11.Forest heights extracted by the least-squares method from the pure volume coherence estimated by the three-stage method.

Figure 11 .
Figure 11.Forest heights extracted by the least-squares method from the pure volume coherence estimated by the three-stage method.