Digital Terrain, Surface, and Canopy Height Models From InSAR Backscatter-Height Histograms

Abstract

This article demonstrates how 3-D vegetation structure can be approximated by interferometric synthetic aperture radar (InSAR) backscatter-height histograms. Single-look backscatter measurements are plotted against the InSAR phase height and are aggregated spatially over a forest patch to form a 3-D histogram, referred to as InSAR backscatter-height histogram or simply InSAR histogram. InSAR histograms resemble LiDAR waveforms, suggesting that existing algorithms used to retrieve canopy height and ground topography from radar tomograms or LiDAR waveforms can be applied to InSAR histograms. Three algorithms are evaluated to generate maps of digital terrain, surface, and canopy height models: Gaussian decomposition, quantile, and backscatter threshold. Full-polarimetric L-band uninhabited aerial vehicle synthetic aperture radar (UAVSAR) data collected over the Gabonese Lopé National Park during the 2016 AfriSAR campaign are used to illustrate and compare the performance of the algorithms for the HH, HV, VV, HH+VV, and HH-VV polarimetric channels. Results show that radar-derived maps using the InSAR histograms differ by 4 m (top-canopy), 5 m (terrain), and 6 m (forest height) in terms of average root-mean-square errors (RMSEs) from standard maps derived from full-waveform laser, vegetation, and ice sensor (LVIS) LiDAR measurements.