An Elevation Ambiguity Resolution Method Based on Segmentation and Reorganization of TomoSAR Point Cloud in 3D Mountain Reconstruction

Tomographic Synthetic Aperture Radar (TomoSAR) is a breakthrough of the traditional SAR, which has the three-dimentional (3D) observation ability of layover scenes such as buildings and high mountains. As an advanced system, the airborne array TomoSAR can effectively avoid temporal de-correlation caused by long revisit time, which has great application in high-precision mountain surveying and mapping. The 3D reconstruction using TomoSAR has mainly focused on low targets, while there are few literatures on 3D mountain reconstruction. Due to the layover phenomenon, surveying in high mountain areas remains a difficult task. Consequently, it is meaningful to carry out the research on 3D mountain reconstruction using the airborne array TomoSAR. However, the original TomoSAR mountain point cloud faces the problem of elevation ambiguity. Furthermore, for mountains with complex terrain, the points located in different elevation periods may intersect. This phenomenon increases the difficulty of solving the problem. In this paper, a novel elevation ambiguity resolution method is proposed. First, Density-Based Spatial Clustering of Applications with Noise (DBSCAN) and Gaussian Mixture Model (GMM) are combined for point cloud segmentation. The former ensures coarse segmentation based on density, and the latter allows fine segmentation of the abnormal categories caused by intersection. Subsequently, the segmentation results are reorganized in the elevation direction to reconstruct all possible point clouds. Finally, the real point cloud can be extracted automatically under the constraints of the boundary and elevation continuity. The performance of the proposed method is demonstrated by simulations and experiments. Based on the airborne array TomoSAR experiment in Leshan City, Sichuan Province, China in 2019, the 3D model of the surveyed mountain is presented. Moreover, three kinds of external data are applied to fully verify the validity of this method.