Next Article in Journal
Integration of Remote Sensing and Mexican Water Quality Monitoring System Using an Extreme Learning Machine
Previous Article in Journal
Spatial Domain Terahertz Image Reconstruction Based on Dual Sparsity Constraints
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Underlying Topography Inversion Using Dual Polarimetric TomoSAR

1
School of Geography and Information Engineering, China University of Geosciences (Wuhan), Wuhan 430074, China
2
School of Geographical Sciences, Guangzhou University, Guangzhou 510006, China
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(12), 4117; https://doi.org/10.3390/s21124117
Submission received: 20 April 2021 / Revised: 1 June 2021 / Accepted: 10 June 2021 / Published: 15 June 2021
(This article belongs to the Section Remote Sensors)

Abstract

:
Underlying topography plays an important role in the national economic construction, military security, resource exploration and investigation. Since synthetic aperture radar tomography (TomoSAR) can achieve the three-dimensional imaging of forests, it has been widely used in underlying topography estimation. At present, there are two kinds of TomoSAR based on the applied datasets: single polarimetric TomoSAR (SP-TomoSAR) and fully polarimetric TomoSAR (FP-TomoSAR). However, SP-TomoSAR cannot obtain the underlying topography accurately due to the lack of enough observations. FP-TomoSAR can improve the estimation accuracy of underlying topography. However, it requires high-cost data acquisition for the large-scale application. Thus, this paper proposes the dual polarimetric TomoSAR (DP-TomoSAR) as another suitable candidate to estimate the underlying topography because of its wide swath and multiple polarimetric observations. Moreover, three frequently used spectral estimation algorithms, namely, Beamforming, Capon and MUSIC, are used in DP-TomoSAR. For validation, a series of simulated experiments was carried out, and the airborne P-band multiple polarimetric SAR data over the Lope, Gabon was also acquired to estimate the underlying topography. The results suggest that DP-TomoSAR in HH & HV combination is more suitable to estimate underlying topography over forest areas than other DP combinations. Moreover, the estimation accuracy of DP-TomoSAR is slightly lower than that of FP-TomoSAR but is higher than that of SP-TomoSAR.

1. Introduction

The digital terrain model (DTM) plays an important role in the national economic construction, military security, resource exploration and investigation and other fields [1,2,3]. However, over forest areas, the traditional aerial survey and optical remote sensing technique can only obtain forest surface elevation, rather than the true terrain elevation. The underlying topography is usually estimated by subtracting the average vegetation height from the forest surface elevation. This method has poor accuracy and low efficiency, particularly for large-scale (global, national and regional) underlying topography inversion [4,5].
Synthetic aperture radar tomography (TomoSAR) is a suitable candidate to solve this problem because of its three-dimensional imaging capability [6,7,8,9]. TomoSAR forms an additional aperture along the vertical direction, which can discriminate scatterers in different heights in the same resolution cell. In particular, the long wavelength (e.g., L and P-band) SAR system, compared with the short wavelength (e.g., X and C-band) SAR system, has a strong penetrating ability over forest areas [10,11,12,13]. It can penetrate the canopy to the ground, acquiring the scattering signal from the ground. Therefore, TomoSAR has been widely used in the inversion of underlying topography over forest areas.
Based on the available datasets, there are two kinds of TomoSAR to estimate the underlying topography: single polarimetric TomoSAR (SP-TomoSAR) and fully polarimetric TomoSAR (FP-TomoSAR). SP-TomoSAR primarily uses tomograms in HH polarization to estimate the ground height, as HH polarization is sensitive to ground scattering. Although the backscattering power of the HH polarization channel is primarily obtained from the ground, it still records the significant forest canopy scattering information by penetrating the forest canopy. Thus, with SP-TomoSAR, it is hard to accurately obtain the underlying topography [14,15,16,17]. Compared with the single polarimetric SAR, fully polarimetric SAR can distinguish different scattering mechanisms in the same resolution cell [13,18]. Consequently, FP TomoSAR can not only acquire scatterers’ heights, but also identify their scattering mechanisms [19,20]. This technique is helpful for the interpretation of tomograms and can improve the accuracy of underlying topography inversion [21]. However, the FP system also has some disadvantages compared with the dual polarimetric (DP) system, such as high pulse repetition frequency [22,23], narrow swath width [24,25,26], and complicated system design caused by additional channel gain [27]. In particular, as a special DP system, compact polarimetric SAR can simultaneously obtain a rich surface information and realize large width observation, which attracts more attention. This suggests that a DP system can better meet the demands of timely updating and complete coverage in large scale or global scale monitoring.
At present, to obtain the forest backscattering power along the vertical direction, a number of spectral estimation methods have been proposed, which can be divided into three kinds: nonparametric spectral estimation, parametric spectral estimation, and sparse spectral estimation. Among these methods, Beamforming [28,29], Capon [30] and MUSIC [31,32,33] are three frequently used estimators. These three methods have already been widely tested by single/fully polarimetric SAR data. However, it is still not yet fully clear how they perform in the dual polarimetric mode. The motivation of this research was, therefore, to fill this knowledge gap by evaluating the dual polarimetric TomSAR in underlying topography inversion. DP-TomoSAR combines the backscattering power of two different polarization channels to separate surface scattering contribution and forest canopy contribution. Beamforming, Capon and MUSIC were expanded in DP-TomoSAR. For validation, a series of simulated experiments were carried out, and the airborne P-band multiple polarimetric SAR data over the Lope, Gabon were also acquired to estimate the underlying topography. Moreover, the LiDAR data were applied to assess their performance under different polarimetric modes.

2. Materials and Methods

2.1. DP-TomoSAR Model

SAR sensors have observed the same target area for N times, and N complex images can be obtained after imaging. The backscattering signal acquired by each pixel on the n-th image can be given by [7]:
g n ( x , r ) =   β ( x , r , z ) exp [ j 4 π b n z λ r sin θ ] d z =   β ( x , r , z ) exp [ j k z ( n ) z ] d z
where k z ( n ) = 4 π b n / λ r sin θ is the vertical wavenumber of the n-th image relative to the master image.   x , r and   z represent the azimuth direction, slant range direction and vertical direction, respectively. g n ( x , r ) indicates the pixel value at   ( x , r ) of the two-dimensional SAR complex image obtained from the n-th image, and β ( x , r , z ) indicates the continuous reflectivity equation. λ is the carrier wavelength; b n represents the length of the vertical baseline between the n-th image and the master image, and θ is the incidence angle.
Similar to FP SAR, DP SAR can not only acquire scatterers’ heights, but also identify their scattering mechanisms. There are three polarimetric combinations, including HH & HV, HH & VV and VV & VH. Here, taking HH & HV as an example, the continuous reflectivity function is discretised along the vertical direction Z, and the dual polarimetric SAR tomography model can be written as follows:
[ g HH 1 g HV 1   g HH N g HV N ] = [ k 1 exp ( j k z ( 1 ) z 1 ) k D exp ( j k z ( 1 ) z D ) k 1 exp ( j k z ( 2 ) z 1 ) k D exp ( j k z ( 2 ) z D ) k 1 exp ( j k z ( N ) z 1 ) k D exp ( j k z ( N ) z D ) ] [ β 1 β 2 β D ] + [ e HH 1 e HV 1   e HH N e HV N ]
where k d = [ k 1 , k 2 ] T , d = 1 , 2 , D   is the polarization measurement vector of two different polarization channels at the position d along the vertical direction and D is the number of discretised sampling.
Multi-look processing is performed in TomoSAR processing, and Formula (2) is written in matrix form:
G DP ( l ) = A ( z , K ) β ( l ) + e DP ( l ) , l = 1 , 2 , 3 , , L
In Equation (3), G DP = [ g HH 1 , g HV 1 , , g HH N , g HV N ] T is the observation vector with N images in DP mode; β = [ β 1 , , β D ] T represents the unknown parameter to be solved; e DP = [ e HH 1 , e HV 1 , , e HH N , e HV N ] T indicates the noise vector from N images in HH and HV channel, and L is the number of looks.
For the 2 N × D dimensional observation matrix A ( z , K ) , which is composed of D steering vector a ( z d , k d ) , the following equations are used:
a ( z d , k d ) = k d a ( z d ) A ( z , K ) = [ a ( z 1 , k 1 ) , , a ( z d , k d ) , a ( z D , k D ) ]
where a ( z d ) = [ exp ( j k z ( 1 ) z d ) , exp ( j k z ( 2 ) z d ) , exp ( j k z ( N ) z d ) ] , and K = [ k 1 , k 2 , , k D ] .

2.2. Classical Spectral Estimation Methods for DP-TomoSAR

At present, Beamforming, Capon and MUSIC are the three most common spectral estimation algorithms in the field of TomoSAR. Therefore, these three methods are used to explore the performance of DP-TomoSAR.

2.2.1. DP-Beamforming Estimator

The classical Beamforming method is regarded as a finite impulse response filter [34]. This ‘filter’ allows the signal of a certain spatial frequency to pass through without distortion and attenuates the signal of other frequencies.
Following the derivation of the classical Beamforming algorithm, the DP-Beamforming estimator is derived as follows. The optimisation of dual polarization Beamforming is given by:
min h h B H h B   s . t .   h B H a ( z d , k d ) = 1
where h B is the filter of the dual polarization Beamforming algorithm and it can be expressed as:
h B = a ( z d , k d ) a H ( z d , k d ) a ( z d , k d ) = a ( z d , k d ) N
The power of the filtered signal is as follows:
E { β D ( l ) β D H ( l ) } = a H ( z d , k d ) R a ( z d , k d ) N 2
It can be maximised according to the polarization state k d by:
max | | k d | | 2 = 1 a H ( z d , k d ) R a ( z d , k d ) N 2
where R is the covariance matrix of DP-TomoSAR model.
For the convenience of calculation, the steering vector a ( z d , k d ) is written as:
a ( z d , k d ) = ( I ( 2 × 2 ) a ( z d ) ) k d = B ( z d ) k d
Thus, the maximum problem in Equation (8) can be transformed into the maximum eigenvalue problem of the Hermitian matrix B H ( z d ) R B ( z d ) :
B H ( z d ) R B ( z d ) k d max = λ max k d max
Finally, the power by the DP-Beamforming is given by:
P ^ BF ( z d ) = λ max ( B H ( z d ) R ^ DP B ( z d ) ) N 2
where λ max ( ) denotes the maximum eigenvalue of the positive semidefinite matrix. R ^ DP = 1 L l = 1 L G DP ( l ) G DP H ( l ) is the sample covariance matrix.

2.2.2. DP-Capon Estimator

Some shortcomings of conventional Beamforming have been solved by the Capon algorithm, such as low height resolution, irregular sampling and serious sidelobe [28]. Thus, the principle of the Capon estimator is similar to that of Beamforming.
The optimisation problem using DP-TomoSAR can be written as follows:
min h h C H h C   s . t .   h C H a ( z d , k d ) = 1
By using the quadratic minimisation theorem, the filter of dual polarization Capon h C can be obtained:
h C = R 1 a H ( z d , k d ) a H ( z d , k d ) R 1 ( z d , k d )
Then, the filtered signal power can be expressed as follows:
E { β D ( l ) β D H ( l ) } = 1 a H ( z d , k d ) R - 1 a ( z d , k d )
Maximizing it as regards the polarizations k d yields:
max | | k d | | 2 = 1 1 a H ( z d , k d ) R 1 a ( z d , k d )
Combined with Equation (9), the problem of finding the maximum can be transformed into the problem of finding the minimum eigenvalue of Hermitian matrix B H ( z d ) R 1 B ( z d ) :
B H ( z d ) R 1 B ( z d ) k d min = λ min k d min
The power of DP-Capon can be obtained as follows:
P ^ CP ( z d ) = 1 λ min ( B H ( z d ) R ^ DP 1 B ( z d ) )
where λ min ( ) is the minimum eigenvalue operator.

2.2.3. DP-MUSIC Estimator

The MUSIC estimator can identify multiple spatial signals. Specifically, the sample covariance matrix is decomposed into eigenvalues, and the eigenvalues are sorted in descending order [35]. When the signal-to-noise ratio (SNR) is high, the signal eigenvalue can be distinguished from the noise eigenvalue [36].
When the classical MUSIC algorithm is introduced to DP-TomoSAR, the true height values z d and their related reflection mechanisms k d are the solutions of:
a H ( z d , k d ) G G H a ( z d , k d ) = 0
The pseudospectrum of the dual polarization MUSIC algorithm is expressed as
P ^ MU ( z d ) = 1 λ min ( B H ( z d ) G ^ DP G ^ DP H B ( z d ) )
where λ min is the minimum eigenvalue; G ^ DP is the eigenvector matrix, which can be obtained from eigenvalue decomposition of the sample covariance matrix.

3. Numerical Simulation

In general, the distribution of forest backscattered signals along the vertical direction primarily consists of two parts. The first part is the surface scattering signal with small angular spread distribution, which has a scattering centre located on the ground. The other part is the forest canopy scattering contribution with a large angle spread distribution [37]. We carried out two groups of simulation experiments to verify the feasibility and effectiveness of these three DP-TomoSAR algorithms. The parameters such as wavelength, slant range and baseline used in the experiment are shown in Table 1 and Table 2, which are obtained from the F-SAR airborne system, as seen in Figure 1.

3.1. Forest Vertical Profile Reconstruction

The first group of simulation experiments is conducted from three aspects to investigate the performance of forest vertical profile reconstruction:
(1)
different SNR;
(2)
different number of looks (N_obs);
(3)
different height difference between the scattering centres of ground and canopy (Δh).
As shown in Figure 2, the accuracy of FP-TomoSAR, DP-TomoSAR and SP-TomoSAR is equal under low noise, a large number of looks and large height difference between the ground and canopy scattering centres. They can all separate surface scattering contribution and canopy scattering contribution, obtaining the locations of ground and canopy scattering centres accurately.
When decreasing SNR, as seen in Figure 3, the vertical profiles estimated by three methods based on FP data and DP data are also close to the true value, but the SP-TomoSAR estimation has some deviation from the true value. This finding indicates that in the case of strong noise, FP-TomoSAR and DP-TomoSAR have stronger abilities to obtain the ground scattering centres than SP-TomoSAR.
Then, we reduce the number of looks to a relatively low level. As shown in Figure 4, no evident effect is found on the Beamforming and MUSIC algorithm. However, the signal curve reconstructed by SP-Capon sharply deviates from the true value. Therefore, the results based on multi-polarization data have more advantages than those of SP-Capon algorithm.
Finally, we discuss the influence of height difference between the ground scattering centre and the canopy scattering centre (Figure 5). Given the weak ability of the SP model to distinguish ground scattering and canopy scattering mechanisms, the detected ground scattering centres of the SP model have a right-wing shift compared with the results of the multi-polarization model.

3.2. Statisitical Analysis of the Scatterer Separation

The second group of simulation experiments is conducted with a series of multiple repeated independent experiments to probe the resolution ability of three algorithms: Beamforming, Capon and MUSIC. As shown in Figure 6, Figure 7 and Figure 8, we calculate the RMSE of the estimated underlying topography versus SNR, N_obs and Δh.
Figure 6 shows that with the increase in SNR, the RMSE of underlying topography estimation indicates a convergent trend, and the single polarization data converge to a larger RMSE compared with multiple polarization data. As shown in Figure 7, the ability of DP-TomoSAR and FP-TomoSAR to detect the ground scattering centres is better than that of SP-TomoSAR. The convergent trend and convergent level of dual polarization are similar to those of the full polarization data. Figure 8 shows that the RMSE of underlying topography detected by DP-Beamforming and FP-Beamforming almost overlap, which are lower than the RMSE level of SP-Beamforming. The accuracy of dual polarization data based on Capon and MUSIC algorithms is between that of full polarization data and single polarization data.
Based on the two above-mentioned group experiments, under different conditions, the accuracy of underlying topography inversion obtained by DP-TomoSAR is better than that obtained by SP-TomoSAR, which is quite close to that of FP-TomoSAR. Therefore, using a dual polarimetric SAR tomography model to estimate the underlying topography is effective and feasible.

4. Experiments and Results

4.1. Study Area and Datasets

This paper uses 10 images at P-band obtained by the German F-SAR airborne system in the AfriSAR campaign to conduct experiments and verify the effectiveness of underlying topography inversion using DP-TomoSAR. Some basic parameters of the F-SAR system and baseline information are given in Table 1 and Table 2, respectively.
The data coverage is the Lope region, which is located in Gabon, Africa (Figure 9). Most of the experimental regions are hilly areas, primarily composed of savanna and dense forest areas [38], and the topography varies greatly, from 163 m to 583 m above sea level. In addition, the local slope of some areas can exceed 20°, and the maximum tree height is more than 50 m.

4.2. Results and Analysis

The tomogram of each pixel contains the backscattering power signals reconstructed by SAR tomography along the vertical direction. Hence, the analysis of the tomograms is a significant aspect for TomoSAR to obtain the 3D structure of forests. Tomograms are the key to underlying topography inversion.
We select two profiles along the north–south direction (red solid lines as shown in Figure 10) under geographic coordinates to perform SAR tomography imaging. In general, for airborne SAR, some differences in the vertical wavenumber and noise level from the near range to the far range are observed. Thus, we consider line aa′ and line bb′ as a representation of profiles at the near range and far range, respectively. The coordinate system of the SAR data is inconsistent with that of the LiDAR data. The SAR coordinate system is geocoded to the same geographic coordinate as LiDAR for the convenience of the follow-up experiments. Moreover, the window size used to estimate covariance matrix is 31 × 31, and the LiDAR DTM data are sampled to the same size.

4.2.1. Tomograms with Different Combinations

As mentioned in Section 2, the DP-TomoSAR model is composed of three polarisation combinations. At present, the suitable combination method for underlying topography inversion remains unknown. For the convenience of comparison, the LiDAR DTM is also projected onto the tomograms.
For line aa′, Figure 11 shows the tomograms estimated by the three algorithms using different combinations of DP-TomoSAR. Profile aa′ is located in a position that has less noise disturbance in the whole image. The tomograms are relatively complete, and few discontinuous pixels are observed in the whole section. Amongst the three kinds of dual polarization combination, the combination of HH and HV polarizations has the best performance. Several evident discontinuities are found in the other two combinations, particularly in the VV&VH channel, as shown in white circles in these figures. Although the underlying topography information acquired by all combinations of DP-TomoSAR is consistent with LiDAR DTM data, we can still determine that for line aa′, the HH & HV is the most suitable dual polarization combination for underlying topography according to its most continuous tomograms.
Line bb′ is located in a position that has a high noise level in the whole image, which leads to many discontinuities in its tomograms. As shown in Figure 12, the differences amongst ground scattering estimated by different combinations are more significant than those of line aa′. In addition, based on this combination, the results of the HH & HV method have an advantage because each algorithm has obtained a relatively complete layer of underlying topography, which presents a similar trend to LiDAR data. Moreover, the subfigure (i) in Figure 12 indicates that the tomogram estimated by the MUSIC method based on VV&VH is greatly affected by the canopy information, resulting in a serious deviation between the terrain layer and LiDAR. Therefore, this combination is almost entirely not suitable for underlying topography inversion using DP-TomoSAR.
By analysing the imaging mechanism of each polarization channel, the causes of the above-mentioned phenomenon are summarised as follows. Firstly, the corresponding phase of the HH channel is sensitive to the ground scattering; therefore, the VV&VH combination obtained poor inversion results. Secondly, the phase centre of the VV channel is located in the area between the forest canopy and the surface, and the phase centre of the HV channel corresponds to the vegetation canopy. Hence, the combination of HH & HV is consistent with the scattering characteristics over forest covered areas. Theoretically, HH & HV is the optimal combination for retrieving underlying topography when using the DP-TomoSAR model.
Therefore, DP-TomoSAR combining two polarization channels, HH and HV, can effectively identify almost all the surface scattering contribution in the forest area. Using this to estimate the underlying topography is effective and feasible.

4.2.2. Tomograms with Different Polarization Mode

In verifying the performance of DP-TomoSAR, we use SP-TomoSAR and FP-TomoSAR to perform SAR tomography imaging on selected profiles line aa′ and line bb′ and obtain the tomograms (Figure 13 and Figure 14). The window size of the covariance estimation, height imaging range and height sampling interval is the same as those of dual polarization data.
As shown in the figures below, for line aa′, no significant difference is found amongst the single polarization tomograms, dual polarization tomograms and full polarization tomograms estimated by Beamforming, Capon and MUSIC. All the three estimators can detect the ground scattering centres well, which have a good agreement with the LiDAR DTM data. However, based on the tomograms of line bb′, the SP-TomoSAR could not obtain continuous ground scattering centres (particularly the white circles in Figure 14), whereas DP-TomoSAR and FP-TomoSAR can obtain the height of ground corresponding to the profile, probably because the line bb′ profile is located at the far range in the SAR coordinate system. Furthermore, the quality of the data is greatly affected by noise. In this paper, the simulation experiments in Section 3 have proven that the accuracy of DP-TomoSAR and FP-TomoSAR is higher than that of SP-TomoSAR when the noise is strong, and based on the real data experiments, the inversion results of line bb′ also confirm this point.

5. Discussion

Beamforming, Capon and MUSIC algorithms are used in this paper to perform 3D focusing of SP, DP and FP data for the whole study area. The underlying topography under the same geographic coordinate system as the LiDAR DTM estimated by DP-Beamforming, DP-Capon and DP-MUSIC is shown in Figure 15.
As shown in the figure below, the underlying topography obtained by different DP-TomoSAR algorithms has a certain deviation from LiDAR DTM. However, they maintain a high consistency in texture features. This finding may be due to the large elevation span of the study area, in which the height above sea level of the terrain is from 163 m to 583 m. In addition, point cloud data at the edge of the image obtained by LiDAR are lacking, owing to the irregularity of the study area, which affects the accuracy of DTM generated by interpolating point cloud data to a certain extent. Thus, the results obtained by LiDAR data and TomoSAR data are different in the edges and corners of the image.
Moreover, RMSEs among different algorithms using SP-TomoSAR, DP-TomoSAR and FP-TomoSAR with respect to LiDAR DTM are calculated. The results are listed in Table 3, which makes a comprehensive analysis on the performance of three algorithms under different polarization mode.
For the Beamforming estimator, the RMSEs of underlying topography are 8.07, 8.25 and 9.24 m for full polarization, dual polarization and single polarization data, respectively. The result of DP-Beamforming is only 0.18 m higher than that of FP data. However, the difference between DP and SP data reached nearly 1 m. Thus, the multi-polarization mode can enhance the inversion accuracy of the Beamforming method. For the Capon algorithm, the RMSEs estimated by FP-TomoSAR, DP-TomoSAR and SP-TomoSAR are 7.92, 8.09 and 9.20 m, respectively. Based on the above-mentioned observation, the RMSE of SP-Capon exceeds those of DP-Capon and FP-Capon by a margin of more than 1 m. Similar to the estimation of the Beamforming method, the RMSE of DP-Capon is closer to that of FP-Capon and is smaller than that of SP-Capon. This result suggests that the underlying topography retrieved by DP-Capon is reliable. Furthermore, the RMSEs of MUSIC algorithm are 8.01, 8.17 and 9.21 m. Similar to the other two algorithms, the MUSIC method also shows certain advantages in dual-polarization mode.
In particular, no matter which polarization mode is used, Capon has the best performance in underlying topography inversion. Compared with DP-Beamforming, DP-Capon and DP-MUSIC have obtained a wider range of improvement in inversion accuracy, possibly because Capon has improved the resolution along the vertical direction of Beamforming algorithm. Moreover, the Capon method can suit the noise and interference signal in factual image data because of its adaptability; thus, this method can obtain a precise inversion result. As for MUSIC method, in this paper, we set the hyperparameter as 2 according to the canopy volume scattering and the surface scattering phase centres in a pixel over forest areas. However, the forest canopy includes several scatterers; thus, the number of scatterers in a resolution unit is still a relatively complex research object. Therefore, the performance of MUSIC algorithm is inferior to Capon in underlying topography inversion.
The accuracy of the results retrieved by all the three algorithms using DP-TomoSAR and FP-TomoSAR is close, and both of them are superior to those of SP-TomoSAR. This result provides a strong theoretical and data support for the potential application of the dual polarization model in underlying topography.
For spaceborne SAR systems, the dual polarization system is easy to maintain. Moreover, DP-TomoSAR has a higher resolution, wider swath and more data acquisition modes. Thus, it has great importance in optimising working modes of spaceborne SAR systems and estimating the large-scale underlying topography. In the future, DP-TomoSAR may be a primary method to obtain large-scale or global-scale underlying topography with spaceborne SAR data.

6. Conclusions

This paper primarily explores the effectiveness of DP-TomoSAR to estimate the underlying topography and selects three common spectral estimation algorithms (Beamforming, Capon and MUSIC) in TomoSAR to estimate the underlying topography based on dual polarization data.
Firstly, the ability of DP-TomoSAR to detect ground scattering centres is evaluated on the simulation data. We found that when the SNR is lower, or the difference between surface scattering centres and canopy scattering centres is smaller, the inversion results of Beamforming, Capon and MUSIC algorithms based on DP-TomoSAR remarkably improved compared with SP-TomoSAR. Furthermore, the accuracy of the Capon algorithm based on dual polarization data is evidently higher than that of single polarization data when the number of multi-looks is fewer. In addition, the inversion results of DP-TomoSAR are highly consistent with those of FP-TomoSAR.
Secondly, we use 10 images at the P-band obtained by AfriSAR campaign in 2016 to estimate underlying topography. We identify the combination mode of dual polarization that performs best in retrieving underlying topography. Consequently, the DP-TomoSAR model composed of the HH and HV channels can obtain the most complete tomograms. Moreover, tomograms of the selected profiles aa′ and bb′ based on SP and FP data are used as a control group to evaluate the performance of DP-TomoSAR. The results indicate that in terms of tomograms, the inversion effect and precision of DP-TomoSAR is reliable.
Finally, 3D focusing is performed under different polarization data modes in the whole research area. Compared with the DTM obtained by LiDAR, the RMSEs estimated by the Beamforming, Capon and MUSIC algorithms based on DP data are 8.25, 8.09 and 8.17 m, respectively. The RMSEs of underlying topography based on the FP data are 8.07, 7.92 and 8.01 m for the three estimators, and those of the SP data are 9.24, 9.20 and 9.18 m. All the above analyses show that the ability of FP-TomoSAR is slightly better than DP-TomoSAR in retrieving underlying topography. However, they are significantly better than that of single polarization data.
The dual TomoSAR model can be an important method to meet the demand of large-scale underlying topography inversion.

Author Contributions

X.P. (Xing Peng) conceived the idea, designed the experiments and improved the manuscript. S.L. performed the experiments, produced the results and drafted the manuscript. Y.W. improved the visualization of results. X.P. (Xiong Pan), Y.D. and Q.X. contributed to the discussion of the results and revision of the article. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Fundamental Research Funds for the Central Universi-ties, China University of Geosciences (Wuhan, Grant No. CUG200619) and the National Natural Science Foundation of China (No. 41804003). The airborne F-SAR datasets were provided by the ESA (No. 38329).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Huang, J.; Xing, Y.; Qin, L.; Xia, T. Ground Elevation Accuracy Verification by Using ICESat-2 Data under Forest. Infrared Laser Eng. 2020, 49, 122–131. [Google Scholar]
  2. Cazcarra-Bes, V.; Pardini, M.; Tello, M.; Papathanassiou, K.P. Comparison of Tomographic SAR Reflectivity Reconstruction Algorithms for Forest Applications at L-Band. IEEE Trans. Geosci. Remote Sens. 2019, 58, 1–18. [Google Scholar] [CrossRef]
  3. Li, L. Forest Vertical Information Extraction Based on P-Band SAR Tomography. Ph.D. Thesis, Chineese Academy of Forestry, Beijing, China, 2016. [Google Scholar]
  4. Minh, D.H.T.; Le Toan, T.; Rocca, F.; Tebaldini, S.; Villard, L.; Réjou-Méchain, M.; Phillips, O.L.; Feldpausch, T.R.; Dubois-Fernandez, P.; Scipal, K.; et al. SAR Tomography for the Retrieval of Forest Biomass and Height: Cross-Validation at Two Tropical Forest Sites in French Guiana. Remote Sens. Environ. 2016, 175, 138–147. [Google Scholar] [CrossRef] [Green Version]
  5. Peng, X.; Li, X.; Wang, C.; Zhu, J.; Liang, L.; Fu, H. SPICE-Based SAR Tomography over Forest Areas Using a Small Number of P-Band Airborne F-SAR Images Characterized by Non-Uniformly Distributed Baselines. Remote Sens. 2019, 11, 975. [Google Scholar] [CrossRef] [Green Version]
  6. Aghababaee, H.; Sahebi, M.R. Model-Based Target Scattering Decomposition of Polarimetric SAR Tomography. IEEE Trans. Geosci. Remote Sens. 2018, 56, 972–983. [Google Scholar] [CrossRef]
  7. Peng, X.; Wang, C.; Li, X.; Du, Y.; Fu, H.; Yang, Z.; Xie, Q. Three-Dimensional Structure Inversion of Buildings with Nonparametric Iterative Adaptive Approach Using SAR Tomography. Remote Sens. 2018, 10, 1004. [Google Scholar] [CrossRef] [Green Version]
  8. Minh, D.H.T.; Le Toan, T.; Rocca, F.; Tebaldini, S.; d’Alessandro, M.M.; Villard, L. Relating P-Band Synthetic Aperture Radar Tomography to Tropical Forest Biomass. IEEE Trans. Geosci. Remote Sens. 2013, 52, 967–979. [Google Scholar] [CrossRef]
  9. Aghababaei, H.; Ferraioli, G.; Ferro-Famil, L.; Huang, Y.; d’Alessandro, M.M.; Pascazio, V.; Schirinzi, G.; Tebaldini, S. Forest SAR Tomography: Principles and Applications. IEEE Geosci. Remote Sens. Mag. 2020, 8, 30–45. [Google Scholar] [CrossRef]
  10. Lombardini, F.; Reigber, A. Adaptive Spectral Estimation for Multibaseline SAR Tomography with Airborne L-Band Data. In Proceedings of the IGARSS 2003 2003 IEEE International Geoscience and Remote Sensing Symposium, Proceedings (IEEE Cat. No. 03CH37477), Toulouse, France, 21–25 July 2003; Volume 3, pp. 2014–2016. [Google Scholar]
  11. Reigber, A.; Moreira, A.; Papathanassiou, K.P. First Demonstration of Airborne SAR Tomography Using Multibaseline L-Band Data. In Proceedings of the IEEE International Geoscience & Remote Sensing Symposium, Toronto, ON, Canada, 24–28 June 2002. [Google Scholar]
  12. Frey, O.; Morsdorf, F.; Meier, E. Tomographic Imaging of a Forested Area by Airborne Multi-Baseline P-Band SAR. Sensors 2008, 8, 5884–5896. [Google Scholar] [CrossRef] [PubMed]
  13. Sauer, S.; Ferro-Famil, L.; Reigber, A.; Pottier, E. Three-Dimensional Imaging and Scattering Mechanism Estimation over Urban Scenes Using Dual-Baseline Polarimetric InSAR Observations at L-Band. IEEE Trans. Geosci. Remote Sens. 2011, 49, 4616–4629. [Google Scholar] [CrossRef] [Green Version]
  14. Fornaro, G.; Serafino, F.; Soldovieri, F. Three-Dimensional Focusing with Multipass SAR Data. IEEE Trans. Geoence Remote Sens. 2003, 41, 7–517. [Google Scholar] [CrossRef]
  15. Frey, O.; Meier, E. Analyzing Tomographic SAR Data of a Forest with Respect to Frequency, Polarization, and Focusing Technique. IEEE Trans. Geosci. Remote Sens. 2011, 49, 3648–3659. [Google Scholar] [CrossRef] [Green Version]
  16. Nannini, M.; Scheiber, R.; Moreira, A. Estimation of the Minimum Number of Tracks for SAR Tomography. IEEE Trans. Geosci. Remote Sens. 2009, 47, 531–543. [Google Scholar] [CrossRef]
  17. Ertin, E.; Moses, R.; Potter, L. Interferometric Methods for Three-Dimensional Target Reconstruction with Multipass Circular SAR. IET Radar Sonar Navig. 2010, 4, 464–473. [Google Scholar] [CrossRef]
  18. Huang, Y.; Ferro-Famil, L.; Reigber, A. Under-Foliage Object Imaging Using SAR Tomography and Polarimetric Spectral Estimators. IEEE Trans. Geosci. Remote Sens. 2011, 50, 2213–2225. [Google Scholar] [CrossRef] [Green Version]
  19. Liang, L.; Guo, H.; Li, X. Three-Dimensional Structural Parameter Inversion of Buildings by Distributed Compressive Sensing-Based Polarimetric SAR Tomography Using a Small Number of Baselines. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 4218–4230. [Google Scholar] [CrossRef]
  20. Huang, Y. Tomographic Processing of Polarimetric and Interferometric SAR Data for Urban and Forestry Remote Sensing. Ph.D. Thesis, European University of Brittany, Rennes, France, 2011. [Google Scholar]
  21. Ferro-Famil, L.; Huang, Y.; Pottier, E. Principles and Applications of Polarimetric SAR Tomography for the Characterization of Complex Environments. In VIII Hotine-Marussi Symposium on Mathematical Geodesy; Springer: Cham, Switzerland, 2015; pp. 243–255. [Google Scholar]
  22. Jianjun, Z.; Zhiwei, L.; Jun, H. Research Progress and Methods of InSAR for Deformation Monitoring. Acta Geod. Cartogr. Sin. 2017, 46, 1717. [Google Scholar]
  23. Velotto, D.; Migliaccio, M.; Nunziata, F.; Lehner, S. Dual-Polarized TerraSAR-X Data for Oil-Spill Observation. IEEE Trans. Geosci. Remote Sens. 2011, 49, 4751–4762. [Google Scholar] [CrossRef]
  24. Kugler, F.; Hajnsek, I.; Papathanassiou, K. Dual Pol-InSAR Forest Height Estimation by Means of TanDEM-X Data. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Munich, Germany, 22–27 July 2012; pp. 1–4. [Google Scholar]
  25. Kugler, F.; Hajnsek, I. Forest Characterisation by Means of TerraSAR-X and TanDEM-X (Polarimetric and) Interferometric Data. In Proceedings of the 2011 IEEE International Geoscience and Remote Sensing Symposium, Vancouver, BC, Canada, 24–29 July 2011; pp. 2578–2581. [Google Scholar]
  26. Li, H.; Wang, Y.; Luo, X. Tree Height Estimation at Plateau Mountains, Northwestern Sichuan, China Using Dual Pol-InSAR Data. In Proceedings of the 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Beijing, China, 10–15 July 2016; pp. 4698–4701. [Google Scholar]
  27. Wenxue, F.; Huadong, G.; Xinwu, L.; Bangsen, T.; Zhongchang, S. Extended Three-Stage Polarimetric SAR Interferometry Algorithm by Dual-Polarization Data. IEEE Trans. Geosci. Remote Sens. 2015, 54, 2792–2802. [Google Scholar] [CrossRef]
  28. Jr, E.; Parks, T.M. Direction Finding with an Array of Antennas Having Diverse Polarizations. IEEE Trans. Antennas Propag. 1983, 31, 231–236. [Google Scholar]
  29. Homer, J.; Longstaff, I.D.; Callaghan, G. High Resolution 3-D SAR via Multi-Baseline Interferometry. In Proceedings of the IGARSS ’96. 1996 International Geoscience and Remote Sensing Symposium, Lincoln, NE, USA, 27–31 May 1996. [Google Scholar]
  30. Capon, J. High-Resolution Frequency-Wavenumber Spectrum Analysis. Proc. IEEE 1969, 57, 1408–1418. [Google Scholar] [CrossRef] [Green Version]
  31. Stoica, P.; Moses, R.L. Spectral Analysis of Signals; Prentice Hall: Hoboken, NJ, USA, 2005. [Google Scholar]
  32. Schmidt, R.; Schmidt, R.O. Multiple Emitter Location and Signal Parameter Estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef] [Green Version]
  33. Gini, F.; Lombardini, F.; Montanari, M. Layover Solution in Multibaseline SAR Interferometry. IEEE Trans. Aerosp. Electron. Syst. AES 2002, 38, 1344–1356. [Google Scholar] [CrossRef]
  34. Hofmann-Wellenhof, B.; Moritz, H. Introduction to Spectral Analysis. In Mathematical and Numerical Techniques in Physical Geodesy; Springer: Berlin, Germany, 1986; pp. 157–259. [Google Scholar]
  35. Gini, F.; Lombardini, F. Multibaseline Cross-Track SAR Interferometry: A Signal Processing Perspective. IEEE Aerosp. Electron. Syst. Mag. 2005, 20, 71–93. [Google Scholar] [CrossRef]
  36. Bienvenu, G. Influence of the Spatial Coherence of the Background Noise on High Resolution Passive MethodsSpeech, & Signal Processing. In Proceedings of the IEEE International Conference on Acoustics, Washington, DC, USA, 2–4 April 1979. [Google Scholar]
  37. Cao, K.; Tan, W.; Li, X.; Xu, W.; Peng, X. Monitoring Broadleaf Forest Pest Based on L-Band SAR Tomography. IOP Conf. Ser. Earth Environ. Sci. 2019, 237, 052004. [Google Scholar] [CrossRef]
  38. Cazcarra-Bes, V.; Tello-Alonso, M.; Papathanassiou, K. 3D Forest Structure Estimation from SAR Tomography by Means of a Full Rank Polarimetric Inversion Based on Compressive Sensing. In Proceedings of the ESA POLinSAR Workshop, Frascati, Italy, 26–30 January 2015; pp. 26–30. [Google Scholar]
Figure 1. TomoSAR imaging under F-SAR system.
Figure 1. TomoSAR imaging under F-SAR system.
Sensors 21 04117 g001
Figure 2. Forest backscattering profile reconstructed by different methods in SP, DP and FP mode when SNR = 20 dB, N_ob s = 100, and Δh = 30 m: (a) Beamforming, (b) Capon and (c) MUSIC.
Figure 2. Forest backscattering profile reconstructed by different methods in SP, DP and FP mode when SNR = 20 dB, N_ob s = 100, and Δh = 30 m: (a) Beamforming, (b) Capon and (c) MUSIC.
Sensors 21 04117 g002
Figure 3. Forest backscattering profile reconstructed by different methods in SP, DP and FP mode when SNR = 5 dB, N_obs = 100, and Δh = 30 m: (a) Beamforming, (b) Capon and (c) MUSIC.
Figure 3. Forest backscattering profile reconstructed by different methods in SP, DP and FP mode when SNR = 5 dB, N_obs = 100, and Δh = 30 m: (a) Beamforming, (b) Capon and (c) MUSIC.
Sensors 21 04117 g003
Figure 4. Forest backscattering profile reconstructed by different methods in SP, DP and FP mode when SNR = 20 dB, N_obs = 30, and Δh = 30 m: (a) Beamforming, (b) Capon and (c) MUSIC.
Figure 4. Forest backscattering profile reconstructed by different methods in SP, DP and FP mode when SNR = 20 dB, N_obs = 30, and Δh = 30 m: (a) Beamforming, (b) Capon and (c) MUSIC.
Sensors 21 04117 g004
Figure 5. Forest backscattering profile reconstructed by different methods in SP, DP and FP mode when SNR = 20 dB, N_obs = 100, and Δh = 8 m: (a) Beamforming, (b) Capon and (c) MUSIC.
Figure 5. Forest backscattering profile reconstructed by different methods in SP, DP and FP mode when SNR = 20 dB, N_obs = 100, and Δh = 8 m: (a) Beamforming, (b) Capon and (c) MUSIC.
Sensors 21 04117 g005
Figure 6. RMSE of underlying topography estimated by different methods in SP, DP and FP mode versus SNR: (a) Beamforming, (b) Capon and (c) MUSIC.
Figure 6. RMSE of underlying topography estimated by different methods in SP, DP and FP mode versus SNR: (a) Beamforming, (b) Capon and (c) MUSIC.
Sensors 21 04117 g006
Figure 7. RMSE of underlying topography estimated by different methods in SP, DP and FP mode versus N_obs: (a) Beamforming, (b) Capon and (c) MUSIC.
Figure 7. RMSE of underlying topography estimated by different methods in SP, DP and FP mode versus N_obs: (a) Beamforming, (b) Capon and (c) MUSIC.
Sensors 21 04117 g007
Figure 8. RMSE of underlying topography estimated by different methods in SP, DP and FP mode versus Δh: (a) Beamforming, (b) Capon and (c) MUSIC.
Figure 8. RMSE of underlying topography estimated by different methods in SP, DP and FP mode versus Δh: (a) Beamforming, (b) Capon and (c) MUSIC.
Sensors 21 04117 g008
Figure 9. Location of the study area.
Figure 9. Location of the study area.
Sensors 21 04117 g009
Figure 10. Two selected profiles (red solid lines aa′ and bb′) on LiDAR DTM image.
Figure 10. Two selected profiles (red solid lines aa′ and bb′) on LiDAR DTM image.
Sensors 21 04117 g010
Figure 11. Tomograms of line aa′ estimated by different methods using each combination: (ac) Beamforming, respectively, in HH & HV, HH & VV, and VV & VH; (df) Capon, respectively, in HH & HV, HH & VV, and VV & VH; (gi) MUSIC, respectively, in HH & HV, HH & VV and VV & VH. The black solid lines represent LiDAR DTM, and the white circles mark the significant differences among tomograms.
Figure 11. Tomograms of line aa′ estimated by different methods using each combination: (ac) Beamforming, respectively, in HH & HV, HH & VV, and VV & VH; (df) Capon, respectively, in HH & HV, HH & VV, and VV & VH; (gi) MUSIC, respectively, in HH & HV, HH & VV and VV & VH. The black solid lines represent LiDAR DTM, and the white circles mark the significant differences among tomograms.
Sensors 21 04117 g011
Figure 12. Tomograms of line bb′ estimated by different methods using each combination: (ac) Beamforming, respectively, in HH & HV, HH & VV, and VV & VH; (df) Capon, respectively, in HH & HV, HH & VV, and VV & VH; (gi) MUSIC, respectively, in HH & HV, HH & VV and VV & VH. The black solid lines represent LiDAR DTM.
Figure 12. Tomograms of line bb′ estimated by different methods using each combination: (ac) Beamforming, respectively, in HH & HV, HH & VV, and VV & VH; (df) Capon, respectively, in HH & HV, HH & VV, and VV & VH; (gi) MUSIC, respectively, in HH & HV, HH & VV and VV & VH. The black solid lines represent LiDAR DTM.
Sensors 21 04117 g012
Figure 13. Tomograms of line aa′ estimated by different methods in SP, DP and FP mode: (a) Beamforming, SP (b) Beamforming, DP (c) Beamforming, FP (d) Capon, SP (e) Capon, DP (f) Capon, FP (g) MUSIC, SP (h) MUSIC, DP (i) MUSIC, FP.
Figure 13. Tomograms of line aa′ estimated by different methods in SP, DP and FP mode: (a) Beamforming, SP (b) Beamforming, DP (c) Beamforming, FP (d) Capon, SP (e) Capon, DP (f) Capon, FP (g) MUSIC, SP (h) MUSIC, DP (i) MUSIC, FP.
Sensors 21 04117 g013
Figure 14. Tomograms of line aa′ estimated by different methods in SP, DP and FP mode: (a) Beamforming, SP (b) Beamforming, DP (c) Beamforming, FP (d) Capon, SP (e) Capon, DP (f) Capon, FP (g) MUSIC, SP (h) MUSIC, DP (i) MUSIC, FP. The white circles mark the significant differences among tomograms.
Figure 14. Tomograms of line aa′ estimated by different methods in SP, DP and FP mode: (a) Beamforming, SP (b) Beamforming, DP (c) Beamforming, FP (d) Capon, SP (e) Capon, DP (f) Capon, FP (g) MUSIC, SP (h) MUSIC, DP (i) MUSIC, FP. The white circles mark the significant differences among tomograms.
Sensors 21 04117 g014
Figure 15. Underlying topography obtained by different methods of the study area: (a) DP-Beamforming (b) DP-Capon (c) DP-MUSIC and (d) LiDAR.
Figure 15. Underlying topography obtained by different methods of the study area: (a) DP-Beamforming (b) DP-Capon (c) DP-MUSIC and (d) LiDAR.
Sensors 21 04117 g015
Table 1. Parameters of TomoSAR data at the P-band obtained from F-SAR.
Table 1. Parameters of TomoSAR data at the P-band obtained from F-SAR.
Wavelength (m)Polarization ModePlatform Height (m)Range of Incidence Angle (°)Resolution of Range Direction (m)Resolution of Azimuth Direction (m)
0.6890HH + HV + VV609620−503.842
Table 2. Baseline information of TomoSAR data in the Lope region.
Table 2. Baseline information of TomoSAR data in the Lope region.
IdentifierObtaining DateBaseline (m)
FL06-PS0210 February 20160
FL06-PS0310
FL06-PS0420
FL06-PS0540
FL06-PS0660
FL06-PS0780
FL06-PS08−20
FL06-PS09−40
FL06-PS10−60
FL06-PS11−80
Table 3. RMSE of underlying topography estimated by different algorithms using SP-TomoSAR, DP-TomoSAR and FP-TomoSAR.
Table 3. RMSE of underlying topography estimated by different algorithms using SP-TomoSAR, DP-TomoSAR and FP-TomoSAR.
MethodData TypeRMSE (m)
BeamformingSP9.24
DP8.25
FP8.07
CaponSP9.20
DP8.09
FP7.92
MUSICSP9.18
DP8.17
FP8.01
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Peng, X.; Long, S.; Wang, Y.; Xie, Q.; Du, Y.; Pan, X. Underlying Topography Inversion Using Dual Polarimetric TomoSAR. Sensors 2021, 21, 4117. https://doi.org/10.3390/s21124117

AMA Style

Peng X, Long S, Wang Y, Xie Q, Du Y, Pan X. Underlying Topography Inversion Using Dual Polarimetric TomoSAR. Sensors. 2021; 21(12):4117. https://doi.org/10.3390/s21124117

Chicago/Turabian Style

Peng, Xing, Shilin Long, Youjun Wang, Qinghua Xie, Yanan Du, and Xiong Pan. 2021. "Underlying Topography Inversion Using Dual Polarimetric TomoSAR" Sensors 21, no. 12: 4117. https://doi.org/10.3390/s21124117

APA Style

Peng, X., Long, S., Wang, Y., Xie, Q., Du, Y., & Pan, X. (2021). Underlying Topography Inversion Using Dual Polarimetric TomoSAR. Sensors, 21(12), 4117. https://doi.org/10.3390/s21124117

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

Article Metrics

Back to TopTop