A wealth of Full Waveform (FW) LiDAR data are available to the public from different sources, which is poised to boost the extensive application of FW LiDAR data. However, we lack a handy and open source tool that can be used by potential users for processing and analyzing FW LiDAR data. To this end, we introduce waveformlidar, an R package dedicated to FW LiDAR processing, analysis and visualization as a solution to the constraint. Specifically, this package provides several commonly used waveform processing methods such as Gaussian, adaptive Gaussian and Weibull decompositions, and deconvolution approaches (Gold and Richard-Lucy (RL)) with users’ customized settings. In addition, we also developed functions to derive commonly used waveform metrics for characterizing vegetation structure. Moreover, a new way to directly visualize FW LiDAR data is developed through converting waveforms into points to form the Hyper Point cloud (HPC), which can be easily adopted and subsequently analyzed with existing discrete-return LiDAR processing tools such as LAStools and FUSION. Basic explorations of the HPC such as 3D voxelization of the HPC and conversion from original waveforms to composite waveforms are also available in this package. All of these functions are developed based on small-footprint FW LiDAR data, but they can be easily transplanted to the large footprint FW LiDAR data such as Geoscience Laser Altimeter System (GLAS) and Global Ecosystem Dynamics Investigation (GEDI) data analysis. It is anticipated that these functions will facilitate the widespread use of FW LiDAR and be beneficial for better estimating biomass and characterizing vegetation structure at various scales. The package and code examples can be found at https://github.com/tankwin08/waveformlidar.

Full waveform (FW) LiDAR holds great potential for retrieving vegetation structure parameters at a high level of detail, but this prospect is constrained by practical factors such as lack of available handy processing tools and technical intricacy of waveform processing. This study introduces a new product, named the Hyper Point Cloud (HPC) derived from FW LiDAR data, and explore its potential applications such as tree crown delineation using the HPC-based intensity and percentile height (PH) surfaces, which show a promising solution to the constraints of using FW LiDAR data. Results of the HPC present a new direction to handle FW LiDAR data and offer prospects for studying the mid-story and understory of vegetation with high point density (~ 182 points/m2). The intensity-derived digital surface model (DSM) generated from the HPC shows that the ground region has larger maximum intensity (MAXI) and mean intensity (MI) than the vegetation region while having smaller total intensity (TI) and number of intensities (NI) at the given grid cell. Our analysis of intensity distribution contours at individual tree level exhibit similar patterns, indicating that the MAXI and MI are decreasing from the tree crown center to tree boundary while a rising trend is observed for TI and NI. These intensity variable contours provide a theoretical justification for using HPC-based intensity surfaces to segment tree crowns and exploit their potential for extracting tree attributes. The HPC-based intensity surfaces and the HPC-based PH Canopy Height Models (CHM) demonstrate promising tree segmentation results comparable to the LiDAR derived CHM for estimating tree attributes such as tree locations, crown widths and tree heights. We envision that products such as the HPC and the HPC-based intensity and height surfaces introduced in this study can open new perspectives to use FW LiDAR data and alleviate the technical barrier of exploring FW LiDAR data for detailed vegetation structure characterization.