The present study addresses the development, implementation and validation of a forest height mapping scheme based on the combination of TanDEM-X interferometric coherence and GEDI waveform measurements. The very general case where only a single polarisation TanDEM-X interferogram, a set of spatially discrete GEDI waveform measurements and no DTM are available is assumed. The use of GEDI waveforms to invert the TanDEM-X interferometric measurements is described together with a set of performance criteria implemented to ensure a certain performance quality. The emphasis is set on developing a methodology able to invert forest height at large scales. Combining 595 TanDEM-X scenes and about 15 million GEDI waveforms, a spatially continuous 25 m resolution forest height map covering the whole of Tasmania Island is achieved. The derived forest height map is validated against an airborne lidar-derived canopy height map available across the whole island.

Model-based forest height inversion from single-and multi-baseline Pol-InSAR data is today well-established. One of the critical performance points is the parameterisation of the vertical reflectivity profile. In this paper, the parameterisation of the vertical reflectivity profile using a tomographic reconstruction is proposed. To mitigate the limitations in the tomographic reconstruction induced by the limited vertical resolution especially for the ground scattering component, the ground and volume contributions are separated using the Sum of Kronecker Products (SKP) decomposition. For the forest height inversion, the Pol-InSAR line inclination angle is proposed, which, unlike the absolute coherence and the phase center, is invariant to the presence of a residual Dirac-like ground scattering contribution. The inversion performance is addressed and compared for both single- and multi-baseline cases, in the absence and presence of temporal decorrelation. In all cases the ground elevation is estimated from the tomographically reconstructed vertical reflectivity profiles and its accuracy directly influences the achieved inversion performance. The proposed methodology establishes a link between interferometric and tomographic measurements and is therefore relevant for missions that allow the implementation of both techniques. ESA’s BIOMASS is one such mission, and as such is being used as a test case to discuss different implementation scenarios. The different inversion strategies are applied on P-band campaign data acquired in the frame of the AFRISAR 2016 experiment and validated against reference measurements.

The realistic parameterization of the underlying (vertical) radar reflectivity profile is critical for a model-based inversion of forest height from polarimetric interferometric (Pol-InSAR) data. Indeed, an appropriate parameterization not only affects the final estimation performance, but also enables the inversion from a reduced observation space in terms of number of baselines and / or polarizations. Here, we investigate the possibility of using a full tomographic profile to parameterize the inversion. The proposed methodology is demonstrated and validated using Pol-InSAR and tomographic data acquired in the framework of relevant airborne campaigns.

Tomographic Synthetic Aperture Radar (SAR) allows the reconstruction of the 3D radar reflectivity of forests from a large(r) number of multi-angular acquisitions. However, in most practical implementations it suffers from limited vertical resolution and/or reconstruction artefacts as the result of non-ideal acquisition setups. Polarisation Coherence Tomography (PCT) offers an alternative to traditional tomographic techniques that allow the reconstruction of the low-frequency 3D radar reflectivity components from a small(er) number of multi-angular SAR acquisitions. PCT formulates the tomographic reconstruction problem as a series expansion on a given function basis. The expansion coefficients are estimated from interferometric coherence measurements between acquisitions. In its original form, PCT uses the Legendre polynomial basis for the reconstruction of the 3D radar reflectivity. This paper investigates the use of new basis functions for the reconstruction of X-band 3D radar reflectivity of forests derived from available lidar waveforms. This approach enables an improved 3D radar reflectivity reconstruction with enhanced vertical resolution, tailored to individual forest conditions. It also allows the translation from sparse lidar waveform vertical reflectivity information into continuous vertical reflectivity estimates when combined with interferometric SAR measurements. This is especially relevant for exploring the synergy of actual missions such as GEDI and TanDEM-X. The quality of the reconstructed 3D radar reflectivity is assessed by comparing simulated InSAR coherences derived from the reconstructed 3D radar reflectivity against measured coherences at different spatial baselines. The assessment is performed and discussed for interferometric TanDEM-X acquisitions performed over two tropical Gabonese rainforest sites: Mondah and Lopé. The results demonstrate that the lidar-derived basis provides more physically realistic vertical reflectivity profiles, which also produce a smaller bias in the simulated coherence validation, compared to the conventional Legendre polynomial basis.

The TerraSAR-X add on for Digital Elevation Measurement (TanDEM-X) mission provides Interferometric Synthetic Aperture Radar (InSAR) wall-to-wall data (not sparse) at high resolution and at global scale. In addition, the NASA Global Ecosystem Dynamics Investigation (GEDI) is a new spaceborne system that provides (from 51.6°N and 51.6°S) sparse measurements (not images) through LiDAR waveforms. Both systems are sensitivity to the canopy structure such as the forest height but with their own limitations. The TanDEM-X single polarization (HH) interferometric coherence magnitude at X-band provides a continuous mapping of the forest while GEDI provides accurate (but sparse) measurements of the forest. In this paper a methodology of how to combine both systems to estimated forest height is presented and applied to more than 900 TanDEM-X scenes over Gabon in Africa. The forest height results over an area of 1° by 1° are shown and compared respect to GEDI. Finally, a wall-to-wall forest map over the entire country of Gabon is presented as an example of large scale mapping towards a potential global (entire earth) forest height map.