A deep learning solution for height estimation on a forested area based on Pol-TomoSAR data

Forest height and underlying terrain reconstruction is one of the main aims in dealing with forested areas. Theoretically, synthetic aperture radar tomography (TomoSAR) offers the possibility to solve the layover problem, making it possible to estimate the elevation of scatters located in the same resolution cell. This article describes a deep learning approach, named tomographic SAR neural network (TSNN), which aims at reconstructing forest and ground height using multipolarimetric multibaseline (MPMB) SAR data and light detection and ranging (LiDAR)-based data. The reconstruction of the forest and ground height is formulated as a classification problem, in which TSNN, a feedforward network, is trained using covariance matrix elements as input vectors and quantized LiDAR-based data as the reference. In our work, TSNN is trained and tested with P-band MPMB data acquired by ONERA over Paracou region of French Guiana in the frame of the European Space Agency's campaign TROPISAR and LiDAR-based data provided by the French Agricultural Research Center. The novelty of the proposed TSNN is related to its ability to estimate the height with a high agreement with LiDAR-based measurement and actual height with no requirement for phase calibration. Experimental results of different covariance window sizes are included to demonstrate that TSNN conducts height measurement with high spatial resolution and vertical accuracy outperforming the other two TomoSAR methods. Moreover, the conducted experiments on the effects of phase errors in different ranges show that TSNN has a good tolerance for small errors and is still able to precisely reconstruct forest heights

Estimation of Forest Biomass From Two-Level Model Inversion of Single-Pass InSAR Data

A model for aboveground biomass estimation from single-pass interferometric synthetic aperture radar (InSAR) data is presented. Forest height and canopy density estimates Delta h and eta(0), respectively, obtained from two-level model (TLM) inversion, are used as biomass predictors. Eighteen bistatic VV-polarized TanDEM-X (TDM) acquisitions are used, made over two Swedish test sites in the summers of 2011, 2012, and 2013 (nominal incidence angle: 41 degrees; height-of-ambiguity: 32-63 m). Remningstorp features a hemiboreal forest in southern Sweden, with flat topography and where 32 circular plots have been sampled between 2010 and 2011 (area: 0.5 ha; biomass: 42-242 t/ha; height: 14-32 m). Krycklan features a boreal forest in northern Sweden, 720-km north-northeast from Remningstorp, with significant topography and where 31 stands have been sampled in 2008 (area: 2.4-26.3 ha; biomass: 23-183 t/ha; height: 7-21 m). A high-resolution digital terrain model has been used as ground reference during InSAR processing. For the aforementioned plots and stands and if the same acquisition is used for model training and validation, the new model explains 65%-89% of the observed variance, with root-mean-square error (RMSE) of 12%-19% (median: 15%). By fixing two of the three model parameters, accurate biomass estimation can also be done when different acquisitions or different test sites are used for model training and validation, with RMSE of 12%-56% (median: 17%). Compared with a simple scaling model computing biomass from the phase center elevation above ground, the proposed model shows significantly better performance in Remningstorp, as it accounts for the large canopy density variations caused by active management. In Krycklan, the two models show similar performance

The NASA AfriSAR campaign: Airborne SAR and lidar measurements of tropical forest structure and biomass in support of current and future space missions

International audienceIn 2015 and 2016, the AfriSAR campaign was carried out as a collaborative effort among international space and National Park agencies (ESA, NASA, ONERA, DLR, ANPN and AGEOS) in support of the upcoming ESA BIOMASS, NASA-ISRO Synthetic Aperture Radar (NISAR) and NASA Global Ecosystem Dynamics Initiative (GEDI) missions. The NASA contribution to the campaign was conducted in 2016 with the NASA LVIS (Land Vegetation and Ice Sensor) Lidar, the NASA L-band UAVSAR (Uninhabited Aerial Vehicle Synthetic Aperture Radar). A central motivation for the AfriSAR deployment was the common AGBD estimation requirement for the three future spaceborne missions, the lack of sufficient airborne and ground calibration data covering the full range of ABGD in tropical forest systems, and the intercomparison and fusion of the technologies. During the campaign, over 7000 km2 of waveform Lidar data from LVIS and 30,000 km2 of UAVSAR data were collected over 10 key sites and transects. In addition, field measurements of forest structure and biomass were collected in sixteen 1-hectare sized plots. The campaign produced gridded Lidar canopy structure products, gridded aboveground biomass and associated uncertainties, Lidar based vegetation canopy cover profile products, Polarimetric Interferometric SAR and Tomographic SAR products and field measurements. Our results showcase the types of data products and scientific results expected from the spaceborne Lidar and SAR missions; we also expect that the AfriSAR campaign data will facilitate further analysis and use of waveform lidar and multiple baseline polarimetric SAR datasets for carbon cycle, biodiversity, water resources and more applications by the greater scientific community

The BIOMASS level 2 prototype processor : design and experimental results of above-ground biomass estimation

BIOMASS is ESA’s seventh Earth Explorer mission, scheduled for launch in 2022. The satellite will be the first P-band SAR sensor in space and will be operated in fully polarimetric interferometric and tomographic modes. The mission aim is to map forest above-ground biomass (AGB), forest height (FH) and severe forest disturbance (FD) globally with a particular focus on tropical forests. This paper presents the algorithms developed to estimate these biophysical parameters from the BIOMASS level 1 SAR measurements and their implementation in the BIOMASS level 2 prototype processor with a focus on the AGB product. The AGB product retrieval uses a physically-based inversion model, using ground-canceled level 1 data as input. The FH product retrieval applies a classical PolInSAR inversion, based on the Random Volume over Ground Model (RVOG). The FD product will provide an indication of where significant changes occurred within the forest, based on the statistical properties of SAR data. We test the AGB retrieval using modified airborne P-Band data from the AfriSAR and TropiSAR campaigns together with reference data from LiDAR-based AGB maps and plot-based ground measurements. For AGB estimation based on data from a single heading, comparison with reference data yields relative Root Mean Square Difference (RMSD) values mostly between 20% and 30%. Combining different headings in the estimation process significantly improves the AGB retrieval to slightly less than 20%. The experimental results indicate that the implemented retrieval scheme provides robust results that are within mission requirements

A Macroecological Analysis of SERA Derived Forest Heights and Implications for Forest Volume Remote Sensing

Individual trees have been shown to exhibit strong relationships between DBH, height and volume. Often such studies are cited as justification for forest volume or standing biomass estimation through remote sensing. With resolution of common satellite remote sensing systems generally too low to resolve individuals, and a need for larger coverage, these systems rely on descriptive heights, which account for tree collections in forests. For remote sensing and allometric applications, this height is not entirely understood in terms of its location. Here, a forest growth model (SERA) analyzes forest canopy height relationships with forest wood volume. Maximum height, mean, H100, and Lorey's height are examined for variability under plant number density, resource and species. Our findings, shown to be allometrically consistent with empirical measurements for forested communities world-wide, are analyzed for implications to forest remote sensing techniques such as LiDAR and RADAR. Traditional forestry measures of maximum height, and to a lesser extent H100 and Lorey's, exhibit little consistent correlation with forest volume across modeled conditions. The implication is that using forest height to infer volume or biomass from remote sensing requires species and community behavioral information to infer accurate estimates using height alone. SERA predicts mean height to provide the most consistent relationship with volume of the height classifications studied and overall across forest variations. This prediction agrees with empirical data collected from conifer and angiosperm forests with plant densities ranging between 102–106 plants/hectare and heights 6–49 m. Height classifications investigated are potentially linked to radar scattering centers with implications for allometry. These findings may be used to advance forest biomass estimation accuracy through remote sensing. Furthermore, Lorey's height with its specific relationship to remote sensing physics is recommended as a more universal indicator of volume when using remote sensing than achieved using either maximum height or H100

Features of Point Clouds Synthesized from Multi-View ALOS/PRISM Data and Comparisons with LiDAR Data in Forested Areas

LiDAR waveform data from airborne LiDAR scanners (ALS) e.g. the Land Vegetation and Ice Sensor (LVIS) havebeen successfully used for estimation of forest height and biomass at local scales and have become the preferredremote sensing dataset. However, regional and global applications are limited by the cost of the airborne LiDARdata acquisition and there are no available spaceborne LiDAR systems. Some researchers have demonstrated thepotential for mapping forest height using aerial or spaceborne stereo imagery with very high spatial resolutions.For stereo imageswith global coverage but coarse resolution newanalysis methods need to be used. Unlike mostresearch based on digital surface models, this study concentrated on analyzing the features of point cloud datagenerated from stereo imagery. The synthesizing of point cloud data from multi-view stereo imagery increasedthe point density of the data. The point cloud data over forested areas were analyzed and compared to small footprintLiDAR data and large-footprint LiDAR waveform data. The results showed that the synthesized point clouddata from ALOSPRISM triplets produce vertical distributions similar to LiDAR data and detected the verticalstructure of sparse and non-closed forests at 30mresolution. For dense forest canopies, the canopy could be capturedbut the ground surface could not be seen, so surface elevations from other sourceswould be needed to calculatethe height of the canopy. A canopy height map with 30 m pixels was produced by subtracting nationalelevation dataset (NED) fromthe averaged elevation of synthesized point clouds,which exhibited spatial featuresof roads, forest edges and patches. The linear regression showed that the canopy height map had a good correlationwith RH50 of LVIS data with a slope of 1.04 and R2 of 0.74 indicating that the canopy height derived fromPRISM triplets can be used to estimate forest biomass at 30 m resolution

First Demonstration of Space-Borne Polarization Coherence Tomography for Characterizing Hyrcanian Forest Structural Diversity

Structural diversity is recognized as a complementary aspect of biological diversity and plays a fundamental role in forest management, conservation, and restoration. Hence, the assessment of structural diversity has become a major effort in the primary international processes, dealing with biodiversity and sustainable forest management. Because of prohibitive costs associated with the ground measurements of forest structure, despite their high accuracy, space-borne polarization coherence tomography (PCT) can introduce an alternative approach given its ability to provide a vertical reflectivity profile and spatiotemporal resolutions related to detecting forest structural changes. In this study, for the first time ever, the potential of space-borne PCT was evaluated in a broad-leaved Hyrcanian forest of Iran over 308 circular sample plots with an area of 0.1 ha. Two aspects of horizontal structure diversity, including standard deviation of diameter at breast height (σdbh) and the number of trees (N), were predicted as important characteristics in wood production and biomass estimation. In addition, the performance of prediction algorithms, including multiple linear regression (MLR), k-nearest neighbors (k-NN), random forest (RF), and support vector regression (SVR) were compared. We addressed the issue of temporal decorrelation in space-borne PCT utilizing the single-pass TanDEM-X interferometer. The data were acquired in standard DEM mode with single polarization of HH. Consequently, airborne laser scanning (ALS) was used to estimate initial values of height hv and ground phase φ0. The Fourier–Legendre series was used to approximate the relative reflectivity profile of each pixel. To link the relative reflectivity profile averaged within each plot with corresponding ground measurements of σdbh and N, thirteen geometrical and physical parameters were defined (P1−P13). Leave-one-out cross validation (LOOCV) showed a better performance of k-NN than the other algorithms in predicting σdbh and N. It resulted in a relative root mean square error (rRMSE) of 32.80%, mean absolute error (MAE) of 4.69 cm, and R2* of 0.25 for σdbh, whereas only 22% of the variation in N was explained using the PCT algorithm with an rRMSE of 41.56%. This study revealed promising results utilizing TanDEM-X data even though the accuracy is still limited. Hence, an entire assessment of the used framework in characterizing the reflectivity profile and the possible effect of the scale is necessary for future studies

First Demonstration of Space-Borne Polarization Coherence Tomography for Characterizing Hyrcanian Forest Structural Diversity

Structural diversity is recognized as a complementary aspect of biological diversity and plays a fundamental role in forest management, conservation, and restoration. Hence, the assessment of structural diversity has become a major effort in the primary international processes, dealing with biodiversity and sustainable forest management. Because of prohibitive costs associated with the ground measurements of forest structure, despite their high accuracy, space-borne polarization coherence tomography (PCT) can introduce an alternative approach given its ability to provide a vertical reflectivity profile and spatiotemporal resolutions related to detecting forest structural changes. In this study, for the first time ever, the potential of space-borne PCT was evaluated in a broad-leaved Hyrcanian forest of Iran over 308 circular sample plots with an area of 0.1 ha. Two aspects of horizontal structure diversity, including standard deviation of diameter at breast height (σdbh) and the number of trees (N), were predicted as important characteristics in wood production and biomass estimation. In addition, the performance of prediction algorithms, including multiple linear regression (MLR), k-nearest neighbors (k-NN), random forest (RF), and support vector regression (SVR) were compared. We addressed the issue of temporal decorrelation in space-borne PCT utilizing the single-pass TanDEM-X interferometer. The data were acquired in standard DEM mode with single polarization of HH. Consequently, airborne laser scanning (ALS) was used to estimate initial values of height hv and ground phase φ0. The Fourier–Legendre series was used to approximate the relative reflectivity profile of each pixel. To link the relative reflectivity profile averaged within each plot with corresponding ground measurements of σdbh and N, thirteen geometrical and physical parameters were defined (P1−P13). Leave-one-out cross validation (LOOCV) showed a better performance of k-NN than the other algorithms in predicting σdbh and N. It resulted in a relative root mean square error (rRMSE) of 32.80%, mean absolute error (MAE) of 4.69 cm, and R2* of 0.25 for σdbh, whereas only 22% of the variation in N was explained using the PCT algorithm with an rRMSE of 41.56%. This study revealed promising results utilizing TanDEM-X data even though the accuracy is still limited. Hence, an entire assessment of the used framework in characterizing the reflectivity profile and the possible effect of the scale is necessary for future studies

Savannah woody structure modelling and mapping using multi-frequency (X-, C- and L-band) synthetic aperture radar data

Structural parameters of the woody component in African savannahs provide estimates of carbon stocks that are vital to the understanding of fuelwood reserves, which is the primary source of energy for 90% of households in South Africa (80% in Sub-Saharan Africa) and are at risk of over utilisation. The woody component can be characterised by various quantifiable woody structural parameters, such as tree cover, tree height, above ground biomass (AGB) or canopy volume, each been useful for different purposes. In contrast to the limited spatial coverage of ground-based approaches, remote sensing has the ability to sense the high spatio-temporal variability of e.g. woody canopy height, cover and biomass, as well as species diversity and phenological status – a defining but challenging set of characteristics typical of African savannahs. Active remote sensing systems (e.g. Light Detection and Ranging – LiDAR; Synthetic Aperture Radar – SAR), on the other hand, may be more effective in quantifying the savannah woody component because of their ability to sense within-canopy properties of the vegetation and its insensitivity to atmosphere and clouds and shadows. Additionally, the various components of a particular target’s structure can be sensed differently with SAR depending on the frequency or wavelength of the sensor being utilised. This study sought to test and compare the accuracy of modelling, in a Random Forest machine learning environment, woody above ground biomass (AGB), canopy cover (CC) and total canopy volume (TCV) in South African savannahs using a combination of X-band (TerraSAR-X), C-band (RADARSAT-2) and L-band (ALOS PALSAR) radar datasets. Training and validation data were derived from airborne LiDAR data to evaluate the SAR modelling accuracies. It was concluded that the L-band SAR frequency was more effective in the modelling of the CC (coefficient of determination or R2 of 0.77), TCV (R2 of 0.79) and AGB (R2 of 0.78) metrics in Southern African savannahs than the shorter wavelengths (X- and C-band) both as individual and combined (X + C-band) datasets. The addition of the shortest wavelengths also did not assist in the overall reduction of prediction error across different vegetation conditions (e.g. dense forested conditions, the dense shrubby layer and sparsely vegetated conditions). Although the integration of all three frequencies (X + C + L-band) yielded the best overall results for all three metrics (R2 = 0.83 for CC and AGB and R2 = 0.85 for TCV), the improvements were noticeable but marginal in comparison to the L-band alone. The results, thus, do not warrant the acquisition of all three SAR frequency datasets for tree structure monitoring in this environment.Council for Scientific and Industrial Research (CSIR) – South Africa, the Department of Science and Technology, South Africa (Grant Agreement DST/CON 0119/2010, Earth Observation Application Development in Support of SAEOS) and the European Union’s Seventh Framework Programme (FP7/2007-2013, Grant Agreement No. 282621, AGRICAB) for funding this study. The Xband StripMap TerraSAR-X scenes were acquired under a proposal submitted to the TerraSAR-X Science Service of the German Aerospace Center (DLR). The C-band Quad-Pol RADARSAT-2 scenes were provided by MacDonald Dettwiler and Associates Ltd. – Geospatial Services Inc. (MDA GSI), the Canadian Space Agency (CSA), and the Natural Resources Canada’s Centre for Remote Sensing (CCRS) through the Science and Operational Applications Research (SOAR) programme. The L-band ALOS PALSAR FBD scenes were acquired under a K&C Phase 3 agreement with the Japanese Aerospace Exploration Agency (JAXA). The Carnegie Airborne Observatory is supported by the Avatar Alliance Foundation, John D. and Catherine T. MacArthur Foundation, Gordon and Betty Moore Foundation, W.M. Keck Foundation, the Margaret A. Cargill Foundation, Mary Anne Nyburg Baker and G. Leonard Baker Jr., and William R. Hearst III. The application of the CAO data in South Africa is made possible by the Andrew Mellon Foundation, Grantham Foundation for the Protection of the Environment, and the endowment of the Carnegie Institution for Science.http://www.elsevier.com/locate/isprsjprs2016-07-31hb201