The underlying topography is an important part of the three-dimensional structure of forests, and is used for a variety of applications, such as hydrology and water resource management, civil engineering projects, and forest resource surveying. Due to the three-dimensional imaging ability and strong penetration, the tomographic synthetic aperture radar (TomoSAR) with a long wavelength has been shown to be a useful tool to estimate the underlying topography. At present, most of the current methods use the local means method to estimate the sample covariance matrix, in which the vertical backscattering power is estimated. However, these methods cannot easily obtain high-precision underlying topography, and often lose some detailed information. In this paper, to solve this problem, a non-local means method is introduced to estimate the optimal covariance matrix by combining weighted neighborhood pixels. To validate the feasibility and effectiveness of this proposed method, a BioSAR 2008 campaign L-band dataset acquired from the northern forests of Sweden was used to inverse the underlying topography. The results show that the accuracy of the underlying topography retrieved by the proposed method is improved by more than 30% when compared with the traditional method.

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.