Lidar Waveforms Research Articles

Integration of Structural Characteristics from GEDI Waveforms for Improved Forest Type Classification

Abstract: Forest types correspond to differences in structural characteristics and species composition that influence biomass and biodiversity values, which are essential measurements for ecological monitoring and management. However, differentiating forest types in tropical regions remains a challenge. This study aimed to improve forest type extent mapping by combining structural information from discrete full-waveform LiDAR returns with multitemporal images. This study was conducted in a tropical forest region over complex terrain in north-eastern Tanzania. First, structural classes were generated by applying time-series clustering algorithms. The results showed four different structural clusters corresponding to forest types, montane–humid forest, montane–dry forest, submontane forest, and non-forest, when using the Kshape algorithm. Kshape considers the shape of the full-sequence LiDAR waveform, requiring little preprocessing. Despite the overlap amongst the original clusters, the averages of structural characteristics were significantly different across all but five metrics. The labeled clusters were then further refined and used as training data to generate a wall-to-wall forest cover type map by classifying biannual images. The highest-performing model was a KNN model with 13 spectral and 3 terrain features achieving 81.7% accuracy. The patterns in the distributions of forest types provide better information from which to adapt forest management, particularly in forest–non-forest transitional zones.

Classification of compressed full-waveform airborne lidar data

ABSTRACT Airborne laser scanners produce 3D data that can be used for a range of applications, such as urban planning, facility monitoring, flood mapping, and forest management. Additional information on the surveyed area can be obtained from the backscattered waveforms recorded by modern light detection and ranging (lidar) sensors. However, the high-dimensional representation of full-waveform data has hindered progress in its use due to difficulties in processing and storage. This paper develops a quantized convolutional autoencoder network to compress lidar waveform data into a condensed feature representation, resulting in a compression rate of up to 20:1. This, together with height information, is fed into a U-net convolutional neural network that achieves an accuracy of 93.7% on six classes.

Forest 3D Radar Reflectivity Reconstruction at X-Band Using a Lidar Derived Polarimetric Coherence Tomography Basis

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.

Deep learning assisted exponential waveform decomposition for bathymetric LiDAR

Abstract. The processing of bathymetric LiDAR waveforms is an important task, as it provides range and radiometric information to determine the precise location of water surface and bottom, and other characteristics like amplitude. The exponential waveform decomposition proved to be an effective algorithm for bathymetric LiDAR waveforms processing, however, it heavily relies on the high-quality initial estimates of the model parameters. This paper proposes to make use of deep learning to obtain the initial values directly from the input received waveforms without any hand-crafted features and prior-knowledges. Additionally, to provide training samples, we presents a method to create the synthetic bathymetric LiDAR waveforms by simulating of the backscatter cross function returned from water bodies. Two networks with different sensitivities of weak signals were trained by these synthetic waveforms, and used to estimate the initial values of the model parameters, a least square optimization follows up to obtain the final waveform decomposition result. This deep learning assisted exponential waveform decomposition method is applied to the real waveforms acquired by RIEGL VQ-840-G. The results show that estimations with the help of deep learning is less influenced by the intermediate peaks backscattered from objects and particles in water, producing a cleaner point cloud with less isolated points below water surface than the original exponential waveform decomposition. Moreover, the proposed sensitive DL-XDC is even able to detect some very weak bottom returns with low SNR.

Characterizing Canopy Structure Variability in Amazonian Secondary Successions with Full-Waveform Airborne LiDAR

Full-waveform LiDAR (FWF) offers a promising advantage over other technologies to represent the vertical canopy structure of secondary successions in the Amazon region, as the waveform encapsulates the properties of all elements intercepting the emitted beam. In this study, we investigated modifications in the vertical structure of the Amazonian secondary successions across the vegetation gradient from early to advanced stages of vegetation regrowth. The analysis was performed over two distinct climatic regions (Drier and Wetter), designated using the Maximum Cumulative Water Deficit (MCWD). The study area was covered by 309 sample plots distributed along 25 LiDAR transects. The plots were grouped into three successional stages (early—SS1; intermediate—SS2; advanced—SS3). Mature Forest (MF) was used as a reference of comparison. A total of 14 FWF LiDAR metrics from four categories of analysis (Height, Peaks, Understory and Gaussian Decomposition) were extracted using the Waveform LiDAR for Forestry eXtraction (WoLFeX) software (v1.1.1). In addition to examining the variation in these metrics across different successional stages, we calculated their Relative Recovery (RR) with vegetation regrowth, and evaluated their ability to discriminate successional stages using Random Forest (RF). The results showed significant differences in FWF metrics across the successional stages, and within and between sample plots and regions. The Drier region generally exhibited more pronounced differences between successional stages and lower FWF metric values compared to the Wetter region, mainly in the category of height, peaks, and Gaussian decomposition. Furthermore, the Drier region displayed a lower relative recovery of metrics in the early years of succession, compared to the areas of MF, eventually reaching rates akin to those of the Wetter region as succession progressed. Canopy height metrics such as Waveform distance (WD), and Gaussian Decomposition metrics such as Bottom of canopy (BC), Bottom of canopy distance (BCD) and Canopy distance (CD), related to the height of the lower forest stratum, were the most important attributes in discriminating successional stages in both analyzed regions. However, the Drier region exhibited superior discrimination between successional stages, achieving a weighted F1-score of 0.80 compared to 0.73 in the Wetter region. When comparing the metrics from SS in different stages to MF, our findings underscore that secondary forests achieve substantial relative recovery of FWF metrics within the initial 10 years after land abandonment. Regions with potentially slower relative recovery (e.g., Drier regions) may require longer-term planning to ensure success in providing full potential ecosystem services in the Amazon.

Estimating differential penetration of green (532 nm) laser light over sea ice with NASA's Airborne Topographic Mapper: observations and models

Abstract. Differential penetration of green laser light into snow and ice has long been considered a possible cause of range and thus elevation bias in laser altimeters. Over snow, ice, and water, green photons can penetrate the surface and experience multiple scattering events in the subsurface volume before being scattered back to the surface and subsequently the instrument's detector, therefore biasing the range of the measurement. Newly formed sea ice adjacent to open-water leads provides an opportunity to identify differential penetration without the need for an absolute reference surface or dual-color lidar data. We use co-located, coincident high-resolution natural-color imagery and airborne lidar data to identify surface and ice types and evaluate elevation differences between those surfaces. The lidar data reveals that apparent elevations of thin ice and finger-rafted thin ice can be several tens of centimeters below the water surface of surrounding leads, but not over dry snow. These lower elevations coincide with broadening of the laser pulse, suggesting that subsurface volume scattering is causing the pulse broadening and elevation shift. To complement our analysis of pulse shapes and help interpret the physical mechanism behind the observed elevation biases, we match the waveform shapes with a model of scattering of light in snow and ice that predicts the shape of lidar waveforms reflecting from snow and ice surfaces based on the shape of the transmitted pulse, the surface roughness, and the optical scattering properties of the medium. We parameterize the scattering in our model based on the scattering length Lscat, the mean distance a photon travels between isotropic scattering events. The largest scattering lengths are found for thin ice that exhibits the largest negative elevation biases, where scattering lengths of several centimeters allow photons to build up considerable range biases over multiple scattering events, indicating that biased elevations exist in lower-level Airborne Topographic Mapper (ATM) data products. Preliminary analysis of ICESat-2 ATL10 data shows that a similar relationship between subsurface elevations (restored negative freeboard) and “pulse width” is present in ICESat-2 data over sea ice, suggesting that biased elevations caused by differential penetration likely also exist in lower-level ICESat-2 data products. The spatial correlation of observed differential penetration in ATM data with surface and ice type suggests that elevation biases could also have a seasonal component, increasing the challenge of applying a simple bias correction.

Decomposition of full-waveform LiDAR data utilizing an adaptive B-spline-based model and particle swarm optimization

High-precision waveform decomposition is crucial for LiDAR applications. Existing methods encounter challenges including poor target detection and low accuracy in extracting parameters of irregular components, especially in complex echoes. We introduce an adaptive B-spline-based decomposition (AdaptB-spline) method, which uses B-spline curves to adaptively adjust the shape and position of component through the particle swarm optimization (PSO); and proposes an initial parameter estimation method based on the B-spline and Richardson-Lucy (RL) deconvolution, which improves the noise immunity and component detection. Experiments were conducted on synthetic waveforms and satellite LiDAR waveforms by AdaptB-spline and other four methods (Gaussian (Gauss), B-spline-based (B-spline), skew-normal (SkewN), and multi-Gaussian (MultiGauss) decomposition). We concluded that AdaptB-spline exhibits superior performance in terms of component RMSE, CC, R2, component parameter error and range error metrics compared to the four methods. So AdaptB-spline can enhance component detection and accurately fit Gaussian or non-Gaussian waveforms, demonstrating outstanding target detection and ranging precision.

Using multiscale lidar to determine variation in canopy structure from African forest elephant trails

AbstractRecently classified as a unique species by the IUCN, African forest elephants (Loxodonta cyclotis) are critically endangered due to severe poaching. With limited knowledge about their ecological role due to the dense tropical forests they inhabit in central Africa, it is unclear how the Afrotropics are influenced by elephants. Although their role as seed dispersers is well known, they may also drive large‐scale processes that determine forest structure through the creation of elephant trails and browsing the understory, allowing larger, carbon‐dense trees to succeed. Multiple scales of lidar were collected by NASA in Lopé National Park, Gabon from 2015 to 2022. Utilizing two airborne lidar datasets in an African forest elephant stronghold, detailed canopy structural information was used in conjunction with elephant trail data to determine how forest structure varies on and off trails. Forest along elephant trails displayed different structural characteristics than forested areas off trails, with lower canopy height, canopy cover, and different vertical distribution of plant density. Less plant area density was found on trails at 1 m in height, while more vegetation was found at 12 m, compared to off trail locations. Trails in forest areas with previous logging history had lower plant area in the top of the canopy. Forest elephants can be considered as “logging light” ecosystem engineers, affecting canopy structure through browsing and movement. Both airborne lidar scales were able to capture elephant impact along trails, with the high‐resolution discrete return lidar performing higher than waveform lidar.

Retrieval of Tree Height Percentiles over Rugged Mountain Areas via Target Response Waveform of Satellite Lidar

The retrieval of tree height percentiles from satellite lidar waveforms observed over mountainous areas is greatly challenging due to the broadening and overlapping of the ground return and vegetation return. To accurately represent the shape distributions of the vegetation and ground returns, the target response waveform (TRW) is resolved using a Richardson–Lucy deconvolution algorithm with adaptive iteration. Meanwhile, the ground return is identified as the TRW component within a 4.6 m ground signal extent above the end point of the TRW. Based on the cumulative TRW distribution, the height metrics of the energy percentiles of 25%, 50%, 75%, and 95% are determined using their vertical distances relative to the ground elevation in this study. To validate the proposed algorithm, we select the received waveforms of the Global Ecosystem Dynamics Investigation (GEDI) lidar over the Pahvant Mountains of central Utah, USA. The results reveal that the resolved TRWs closely resemble the actual target response waveforms from the coincident airborne lidar data, with the mean values of the coefficient of correlation, total bias, and root-mean-square error (RMSE) taking values of 0.92, 0.0813, and 0.0016, respectively. In addition, the accuracies of the derived height percentiles from the proposed algorithm are greatly improved compared with the conventional Gaussian decomposition method and the slope-adaptive waveform metrics method. The mean bias and RMSE values decrease by the mean values of 1.68 m and 2.32 m and 1.96 m and 2.72 m, respectively. This demonstrates that the proposed algorithm can eliminate the broadening and overlapping of the ground return and vegetation return and presents good potential in the extraction of forest structure parameters over rugged mountainous areas.

Progressive Gaussian Decomposition of Airborne Bathymetric LiDAR Waveform for Improving Seafloor Point Extraction

Airborne bathymetric LiDAR (ABL) acquires waveform data with better accuracy and resolution and greater user control over data processing than discrete returns. The ABL waveform is a mixture of reflections from the water surface and bottom, water column backscattering, and noise, and it can be separated into individual components through waveform decomposition. Because the point density and positional accuracy of the point cloud are dependent on waveform decomposition, an effective decomposition technique is required to improve ABL measurement. In this study, a new progressive waveform decomposition technique based on Gaussian mixture models was proposed for universal applicability to various types of ABL waveforms and to maximize the observation of seafloor points. The proposed progressive Gaussian decomposition (PGD) estimates potential peaks that are not detected during the initial peak detection and progressively decomposes the waveform until the Gaussian mixture model sufficiently represents the individual waveforms. Its performance is improved by utilizing a termination criterion based on the time difference between the originally detected and estimated peaks of the approximated model. The PGD can be universally applied to various waveforms regardless of water depth or underwater environment. To evaluate the proposed approach, it was applied to the waveform data acquired from the Seahawk sensor developed in Korea. In validating the PGD through comparative evaluation with the conventional Gaussian decomposition method, the root mean square error was found to decrease by approximately 70%. In terms of point cloud extractability, the PGD extracted 14–18% more seafloor points than the Seahawk’s data processing software.

LAI estimation based on physical model combining airborne LiDAR waveform and Sentinel-2 imagery.

Leaf area index (LAI) is an important biophysical parameter of vegetation and serves as a significant indicator for assessing forest ecosystems. Multi-source remote sensing data enables large-scale and dynamic surface observations, providing effective data for quantifying various indices in forest and evaluating ecosystem changes. However, employing single-source remote sensing spectral or LiDAR waveform data poses limitations for LAI inversion, making the integration of multi-source remote sensing data a trend. Currently, the fusion of active and passive remote sensing data for LAI inversion primarily relies on empirical models, which are mainly constructed based on field measurements and do not provide a good explanation of the fusion mechanism. In this study, we aimed to estimate LAI based on physical model using both spectral imagery and LiDAR waveform, exploring whether data fusion improved the accuracy of LAI inversion. Specifically, based on the physical model geometric-optical and radiative transfer (GORT), a fusion strategy was designed for LAI inversion. To ensure inversion accuracy, we enhanced the data processing by introducing a constraint-based EM waveform decomposition method. Considering the spatial heterogeneity of canopy/ground reflectivity ratio in regional forests, calculation strategy was proposed to improve this parameter in inversion model. The results showed that the constraint-based EM waveform decomposition method improved the decomposition accuracy with an average 12% reduction in RMSE, yielding more accurate waveform energy parameters. The proposed calculation strategy for the canopy/ground reflectivity ratio, considering dynamic variation of parameter, effectively enhanced previous research that relied on a fixed value, thereby improving the inversion accuracy that increasing on the correlation by 5% to 10% and on R2 by 62.5% to 132.1%. Based on the inversion strategy we proposed, data fusion could effectively be used for LAI inversion. The inversion accuracy achieved using both spectral and LiDAR data (correlation=0.81, R2=0.65, RMSE=1.01) surpassed that of using spectral data or LiDAR alone. This study provides a new inversion strategy for large-scale and high-precision LAI inversion, supporting the field of LAI research.

Forest terrain and canopy height estimation using stereo images and spaceborne LiDAR data from GF-7 satellite

ABSTRACT Accurate estimation of forest terrain and canopy height is crucial for timely understanding of forest growth. Gao Fen-7 (GF-7) Satellite is China’s first sub-meter-level three-dimensional (3D) mapping satellite for civilian use, which was equipped with a two-line-array stereo mapping camera and a laser altimeter system that can provide stereo images and full waveform LiDAR data simultaneously. Most of the existing studies have concentrated on evaluating the accuracy of GF-7 for topographic survey in bare land, but few have in-depth studied its ability to measure forest terrain elevation and canopy height. The purpose of this study is to evaluate the potential of GF-7 LiDAR and stereo image for forest terrain and height measurement. The Airborne Laser Scanning (ALS) data were utilized to generate reference terrain and forest vertical information. The validation test was conducted in Pu’er City, Yunnan Province of China, and encouraging results have obtained. The GF-7 LiDAR data obtained the accuracy of forest terrain elevation with RMSE of 8.01 m when 21 available laser footprints were used for results verification; meanwhile, when it was used to calculate the forest height, R2 of 0.84 and RMSE of 3.2 m were obtained although only seven effective footprints were used for result verification. The canopy height values obtained from GF-7 stereo images have also been proven to have high accuracy with the resolution of 20 m × 20 m compared with ALS data (R2 = 0.88, RMSE = 2.98 m). When the results were verified at the forest sub-compartment scale that taking into account the forest types, further higher accuracy (R2 = 0.96, RMSE = 1.23 m) was obtained. These results show that GF-7 has considerable application potential in forest resources monitoring.

Impact of leaf phenology on estimates of aboveground biomass density in a deciduous broadleaf forest from simulated GEDI lidar

The Global Ecosystem Dynamics Investigation (GEDI) is a waveform lidar instrument on the International Space Station used to estimate aboveground biomass density (AGBD) in temperate and tropical forests. Algorithms to predict footprint AGBD from GEDI relative height (RH) metrics were developed from simulated waveforms with leaf-on (growing season) conditions. Leaf-off GEDI data with lower canopy cover are expected to have shorter RH metrics, and are therefore excluded from GEDI’s gridded AGBD products. However, the effects of leaf phenology on RH metric heights, and implications for GEDI footprint AGBD models that can include multiple nonlinear RH predictors, have not been quantified. Here, we test the sensitivity of GEDI data and AGBD predictions to leaf phenology. We simulated GEDI data using high-density drone lidar collected in a temperate mountain forest in the Czech Republic under leaf-off and leaf-on conditions, 51 d apart. We compared simulated GEDI RH metrics and footprint-level AGBD predictions from GEDI Level 4 A models from leaf-off and leaf-on datasets. Mean canopy cover increased by 31% from leaf-off to leaf-on conditions, from 57% to 88%. RH metrics < RH50 were more sensitive to changes in leaf phenology than RH metrics ⩾ RH50. Candidate AGBD models for the deciduous-broadleaf-trees prediction stratum in Europe that were trained using leaf-on measurements exhibited a systematic prediction difference of 0.6%–19% when applied to leaf-off data, as compared to leaf-on predictions. Models with the least systematic prediction difference contained only the highest RH metrics, or contained multiple predictor terms that contained both positive and negative coefficients, such that the difference from systematically shorter leaf-off RH metrics was partially offset among the multiple terms. These results suggest that, with consideration of model choice, leaf-off GEDI data can be suitable for AGBD prediction, which could increase data availability and reduce sampling error in some forests.

Influence of phenology on waveform features in deciduous and coniferous trees in airborne LiDAR

Information on forest structure is vital for sustainable forest management. Currently, airborne LiDAR remote sensing has been well established as an effective tool to characterize the structure of canopies and forest inventory variables. Radiometry and geometry are highly intertwined in LiDAR remote sensing of forest vegetation and phenology influences the geometric-optical properties of deciduous and evergreen trees causing seasonal variation in LiDAR observations. This variation may be considered as a nuisance or exploited in for example tree species identification. Airborne LiDAR data are also influenced by sensor functioning, acquisition settings, scan geometry and the atmosphere. Reliable estimation of subtle phenological effects calls for data in which the impact of the external factors is minimal. We experimented with such data and explored LIDAR waveforms (WFs) in boreal trees in winter, early summer and late summer. Our objectives were to i) assess the match of the multitemporal LiDAR data for observing true changes in vegetation; ii) quantify the influence of phenology in deciduous and evergreen trees; iii) study the effect of varying scan zenith angle (SZA) and canopy age on WF features in different phenostates; iv) assess the temporal feature correlation in individual living and dead standing trees. A WF-recording pulsed LiDAR sensor unit operating at the wavelength of 1550 nm was used in repeated acquisitions. WF attributes such as energy, peak amplitude and echo width were derived for each pulse and were localized vertically to crown, understory and ground components. Silver and downy birch, black alder, European aspen, Siberian larch, Scots pine, Norway spruce and dead standing spruce formed our strata. Results showed that phenology caused more variation in WF features of deciduous trees compared to evergreen conifers. Deciduous trees displayed substantial between-species variation that was linked with differences in branching pattern, leaf orientation and bark reflectance. Pine displayed a possible winter-early summer anomaly in canopy backscattering that may be linked with changes in foliage clumping or with the role of stamens in early summer trees. Trees displayed positive temporal correlation in WF features and correlations were the strongest in evergreen and deciduous conifers and decreased with time. SZA had minor influence on WF features whereas age exercised a strong effect on many features with parallel variation between species and phenostates. Structural changes following death, i.e. ‘aging’ changed the geometric WF features of dead standing trees. Our results provide new insights for enhancing tree species identification by using WF LiDAR and for LiDAR time-series analysis of vegetation.

Extrapolation of canopy height and cover metrics of GEDI LiDAR in tropical montane forest ecosystem

ABSTRACT Global Ecosystem Dynamics Investigation (GEDI); a full waveform lidar sensor, acquires samples of the terrain and vegetation structures. To derive a spatially continuous estimate, the canopy height and cover metrics retrieved from GEDI L2A and L2B, respectively, were extrapolated using support vector regression models and explanatory variables derived from multiple remotely sensed datasets. Explanatory variables extrapolated GEDI canopy height with RMSE = 3.83 m (R2 = 0.84) as the performance was validated using a subset of GEDI data. However, the accuracy decreased to RMSE = 7.98 m (R2 = 0.65) as validated with field-measured canopy heights. The accuracy of the prediction of forest canopy covers retrieved from GEDI L2B generally remains poor; RMSE = 0.14 (R2 = 0.53).

Airborne dual-wavelength waveform LiDAR improves species classification accuracy of boreal broadleaved and coniferous trees

Tree species identification constitutes a bottleneck in remote sensing applications. Waveform LiDAR has been shown to offer potential over discrete-return observations, and we assessed if the combination of two-wavelength waveform data can lead to further improvements. A total of 2532 trees representing seven living and dead conifer and deciduous species classes found in Hyytiälä forests in southern Finland were included in the experiments. LiDAR data was acquired by two single-wavelength sensors. The 1064-nm and 1550-nm data were radiometrically corrected to enable range-normalization using the radar equation. Pulses were traced through the canopy, and by applying 3D crown models, the return waveforms were assigned to individual trees. Crown models and a terrain model enabled a further split of the waveforms to strata representing the crown, understory and ground segments. Different geometric and radiometric waveform attributes were extracted per return pulse and aggregated to tree-level mean and standard deviation features. We analyzed the effect of tree size on the features, the correlation between features and the between-species differences of the waveform features. Feature importance for species classification was derived using F-test and the Random Forest algorithm. Classification tests showed significant improvement in overall accuracy (74→83% with 7 classes, 88→91% with 4 classes) when the 1064-nm and 1550-nm features were merged. Most features were not invariant to tree size, and the dependencies differed between species and LiDAR wavelength. The differences were likely driven by factors such as bark reflectance, height growth induced structural changes near the treetop as well as foliage density in old trees.

GEDI waveform metrics in vegetation mapping—a case study from a heterogeneous tropical forest landscape

The distribution of different vegetation types is important information for landscape management, especially in the context of tackling global environmental change. Vegetation types can be mapped using satellite and airborne passive remote sensing. However, spectrally similar yet structurally different vegetation types, like different tree-dominated land covers, are often challenging to map using spectral information alone. We examined the potential of vertical vegetation structure acquired in the global ecosystem dynamics investigation (GEDI) mission that harnesses a space-borne waveform lidar sensor in vegetation mapping across a heterogeneous tropical landscape in Cambodia. We extracted 121 waveform metrics from Level-1B and Level-2A data products at 1062 locations across five key vegetation types. After reducing the relative height variables’ dimensionality through simple linear regressions, we developed a Random Forest classifier to predict vegetation classes based on 23 GEDI metrics. We then used this model to classify the vegetation types across more than 77 000 GEDI footprints in the study area. GEDI metrics alone were useful in identifying vegetation types with 81% accuracy. Cropland/grassland class had the highest prediction accuracy (user’s accuracy [UA] = 89%; producer’s accuracy [PA] = 91%), while dry deciduous forest had the lowest accuracy (UA = 73%; PA = 69%). By comparing the GEDI-only classification with an optical-radar map, we found that structural and topographic information from GEDI Level-1B and Level-2A can complement the spectral information in assessing natural habitats that neighbor other vegetation types in a heterogeneous landscape. The highest classification accuracy at the footprint scale was obtained from the combination of GEDI, Sentinel-1, and Sentinel-2 (88.3%). We also demonstrated how wall-to-wall vegetation mapping is possible by combining the three data sources. These findings expand the potential use of GEDI waveform lidar data in supporting the development of policy-relevant maps that depict the distribution of forests together with other vegetation types.

Forest Biomass Mapping Using Continuous InSAR and Discrete Waveform Lidar Measurements: A TanDEM-X/GEDI Test Study

Forest Biomass Mapping Using Continuous InSAR and Discrete Waveform Lidar Measurements: A TanDEM-X/GEDI Test Study

Gaussian convolution decomposition for non-Gaussian shaped pulsed LiDAR waveform

The full waveform decomposition technique is significant for LiDAR ranging. It is challenging to extract the parameters from non-Gaussian shaped waveforms accurately. Many parametric models (e.g. the Gaussian distribution, the lognormal distribution, the generalized normal distribution, the Burr distribution, and the skew-normal distribution) were proposed to fit sharply-peaked, heavy-tailed, and negative-tailed waveforms. However, these models can constrain the shape of the waveform components. In this article, the Gaussian convolution model is established. Firstly, a set of Gaussian functions is calculated to characterize the system waveform so that asymmetric and non-Gaussian system waveforms can be included. The convolution result of the system waveform and the target response is used as the model for fitting the overlapped echo. Then a combination method of the Richardson–Lucy deconvolution, layered iterative, and Gaussian convolution is introduced to estimate the initial parameters. The Levenberg–Marquardt algorithm is used for the optimization fitting. Through experiments on synthetic data and practical recorded coding LiDAR data, we compare the proposed method with two decomposition approaches (Gaussian decomposition and skew-normal decomposition). The experiment results revealed that the proposed method could precisely decompose the overlapped non-Gaussian heavy-tailed waveforms and provide the best ranging accuracy, component fitting accuracy, and anti-noise performance. However, the traditional Gaussian and skew-normal decomposition methods can not fit the components well, resulting in inaccurate range estimates.

Modelling forest volume with small area estimation of forest inventory using GEDI footprints as auxiliary information

• A suitable GEDI RH metric was selected to represent NFI dominant height at plot level. • A method that pairs up NFI plots with nearby GEDI footprints up to a critical distance suggested by the semi-variance examination. • Small area estimation of forest volume was improved with the help of GEDI auxiliary data. The French National Forest Inventory provides detailed forest information up to large national and regional scales. Forest inventory for small areas of interest within a large population is equally important for decision making, such as for local forest planning and management purposes. However, sampling these small areas with sufficient ground plots is often not cost efficient. In response, small area estimation has gained increasing popularity in forest inventory. It consists of a set of techniques that enables predictions of forest attributes of subpopulation with the help of auxiliary information that compensates for the small field samples. Common sources of auxiliary information usually come from remote sensing technology, such as airborne laser scanning and satellite imagery. The newly launched NASA’s Global Ecosystem Dynamics Investigation (GEDI), a full waveform Lidar instrument, provides an unprecedented opportunity of collecting large-scale and dense forest sample plots given its sampling frequency and spatial coverage. However, the geolocation uncertainty associated with GEDI footprints create important challenges for their use for small area estimations. In this study, we designed a process that provides NFI measurements at plot level with GEDI auxiliary information from nearby footprints. We demonstrated that GEDI RH 98 is equivalent to NFI dominant height at plot level. We stressed the importance of pairing NFI plots with nearby GEDI footprints, based on not only the distance in between but also their similarities, i.e., forest heights and forest types. Subsequently, these NFI-GEDI pairs were used for small area estimations following a two-phase sampling scheme. We showcased that, with an adequate sample size, small area estimation with GEDI auxiliary data can improve the accuracy of forest volume estimates.