The upcoming space-borne LiDAR satellite Ice, Cloud and land Elevation Satellite-2 (ICESat-2) is scheduled to launch in 2018. Different from the waveform LiDAR system onboard the ICESat, ICESat-2 will use a micro-pulse photon-counting LiDAR system. Thus new data processing algorithms are required to retrieve vegetation canopy height from photon-counting LiDAR data. The objective of this paper is to develop and validate an automated approach for better estimating vegetation canopy height. The new proposed method consists of three key steps: 1) filtering out the noise photons by an effective noise removal algorithm based on localized statistical analysis; 2) separating ground returns from canopy returns using an iterative photon classification algorithm, and then determining ground surface; 3) generating canopy-top surface and calculating vegetation canopy height based on canopy-top and ground surfaces. This automatic vegetation height estimation approach was tested to the simulated ICESat-2 data produced from Sigma Space LiDAR data and Multiple Altimeter Beam Experimental LiDAR (MABEL) data, and the retrieved vegetation canopy heights were validated by canopy height models (CHMs) derived from airborne discrete-return LiDAR data. Results indicated that the estimated vegetation canopy heights have a relatively strong correlation with the reference vegetation heights derived from airborne discrete-return LiDAR data (R² and RMSE values ranging from 0.639 to 0.810 and 4.08 m to 4.56 m respectively). This means our new proposed approach is appropriate for retrieving vegetation canopy height from micro-pulse photon-counting LiDAR data.

The accurate estimation of crop biomass during the growing season is very important for crop growth monitoring and yield estimation. The objective of this paper was to explore the potential of hyperspectral and light detection and ranging (LiDAR) data for better estimation of the biomass of maize. First, we investigated the relationship between field-observed biomass with each metric, including vegetation indices (VIs) derived from hyperspectral data and LiDAR-derived metrics. Second, the partial least squares (PLS) regression was used to estimate the biomass of maize using VIs (only) and LiDAR-derived metrics (only), respectively. Third, the fusion of hyperspectral and LiDAR data was evaluated in estimating the biomass of maize. Finally, the biomass estimates were validated by a leave-one-out cross-validation (LOOCV) method. Results indicated that all VIs showed weak correlation with field-observed biomass and the highest correlation occurred when using the red edge-modified simple ratio index (ReMSR). Among all LiDAR-derived metrics, the strongest relationship was observed between coefficient of variation (${\mathrm{H}}_{CV}$of digital terrain model (DTM) normalized point elevations with field-observed biomass. The combination of VIs through PLS regression could not improve the biomass estimation accuracy of maize due to the high correlation between VIs. In contrast, the ${\mathrm{H}}_{CV}$ combined with ${\mathrm{H}}_{mean}$ performed better than one LiDAR-derived metric alone in biomass estimation (R² = 0.835, RMSE = 374.655 g/m², RMSE_CV = 393.573 g/m²). Additionally, our findings indicated that the fusion of hyperspectral and LiDAR data can provide better biomass estimates of maize (R² = 0.883, RMSE = 321.092 g/m², RMSE_CV = 337.653 g/m²) compared with LiDAR or hyperspectral data alone.

The Ice, Cloud and Land Elevation Satellite-2 (ICESat-2) mission employs a micro-pulse photon-counting LiDAR system for mapping and monitoring the biomass and carbon of terrestrial ecosystems over large areas. In preparation for ICESat-2 data processing and applications, this paper aimed to develop and validate an effective algorithm for better estimating ground elevation and vegetation height from photon-counting LiDAR data. Our new proposed algorithm consists of three key steps. Firstly, the noise photons were filtered out using a noise removal algorithm based on localized statistical analysis. Secondly, we classified the signal photons into canopy photons and ground photons by conducting a series of operations, including elevation frequency histogram building, empirical mode decomposition (EMD), and progressive densification. At the same time, we also identified the top of canopy (TOC) photons from canopy photons by percentile statistics method. Thereafter, the ground and TOC surfaces were generated from ground photons and TOC photons by cubic spline interpolation, respectively. Finally, the ground elevation and vegetation height were estimated by retrieved ground and TOC surfaces. The results indicate that the noise removal algorithm is effective in identifying background noise and preserving signal photons. The retrieved ground elevation is more accurate than the retrieved vegetation height, and the results of nighttime data are better than those of the corresponding daytime data. Specifically, the root-mean-square error (RMSE) values of ground elevation estimates range from 2.25 to 6.45 m for daytime data and 2.03 to 6.03 m for nighttime data. The RMSE values of vegetation height estimates range from 4.63 to 8.92 m for daytime data and 4.55 to 8.65 m for nighttime data. Our algorithm performs better than the previous algorithms in estimating ground elevation and vegetation height due to lower RMSE values. Additionally, the results also illuminate that the photon classification algorithm effectively reduces the negative effects of slope and vegetation coverage. Overall, our paper provides an effective solution for estimating ground elevation and vegetation height from micro-pulse photon-counting LiDAR data.

Stem characteristics of plants are of great importance to both ecology study and forest management. Terrestrial laser scanning (TLS) may provide an effective way to characterize the fine-scale structures of vegetation. However, clumping plants, dense foliage and thin structure could intensify the shadowing effect and pose a series of problems in identifying stems, distinguishing neighboring stems, and merging disconnected stem parts in point clouds. This paper presents a new method to automatically detect stems in dense and homogeneous forest using single-scan TLS data. Stem points are first identified with a two-scale classification method. Then a clustering approach is used to group the candidate stem points. Finally, a direction-growing algorithm based on a simple stem curve model is applied to merge stem points. Field experiments were carried out in two different bamboo plots with a stem density of about 7500 stems/ha. Overall accuracy of the stem detection is 88% and the quality of detected stems is mainly affected by the shadowing effect. Results indicate that the proposed method is feasible and effective in detection of bamboo stems using TLS data, and can be applied to other species of single-stem plants in dense forests.

A comparative study of water indices and image classification algorithms for mapping inland water bodies using Landsat imagery was carried out through obtaining 24 high-resolution (≤5 m) and cloud-free images archived in Google Earth with the same (or ±1 day) acquisition dates as the Landsat-8 OLI images over 24 selected lakes across the globe, and developing a method to generate the alternate ground truth data from the Google Earth images for properly evaluating the Landsat image classification results. In addition to the commonly used green band-based water indices, Landsat-8 OLI’s ultra-blue, blue, and red band-based water indices were also tested in this research. Two unsupervised (the zero-water index threshold H0 method and Otsu’s automatic threshold selection method) and one supervised (the k-nearest neighbor (KNN) method) image classification algorithms were employed for conducting the image classification. Through comparing a total of 2880 Landsat image classification results with the alternate ground truth data, this study showed that (1) it is not necessary to use some supervised image classification methods for extracting water bodies from Landsat imagery given the high computational cost associated with the supervised image classification algorithms; (2) the unsupervised classification algorithms such as the H0 and Otsu methods could achieve comparable accuracy as the KNN method, although the H0 method produced more large error outliers than the Otsu method, thus the Otsu method is better than the H0 method; and (3) the ultra-blue band-based AWEI_nsuB is the best water index for the H0 method, and the ultra-blue band-based MNDWI_2uB is the best water index for both the Otsu and KNN methods.

Accurate land cover classification information is a critical variable for many applications. This study presents a method to classify land cover using the fusion data of airborne discrete return LiDAR (Light Detection and Ranging) and CASI (Compact Airborne Spectrographic Imager) hyperspectral data. Four LiDAR-derived images (DTM, DSM, nDSM, and intensity) and CASI data (48 bands) with 1 m spatial resolution were spatially resampled to 2, 4, 8, 10, 20 and 30 m resolutions using the nearest neighbor resampling method. These data were thereafter fused using the layer stacking and principal components analysis (PCA) methods. Land cover was classified by commonly used supervised classifications in remote sensing images, i.e., the support vector machine (SVM) and maximum likelihood (MLC) classifiers. Each classifier was applied to four types of datasets (at seven different spatial resolutions): (1) the layer stacking fusion data; (2) the PCA fusion data; (3) the LiDAR data alone; and (4) the CASI data alone. In this study, the land cover category was classified into seven classes, i.e., buildings, road, water bodies, forests, grassland, cropland and barren land. A total of 56 classification results were produced, and the classification accuracies were assessed and compared. The results show that the classification accuracies produced from two fused datasets were higher than that of the single LiDAR and CASI data at all seven spatial resolutions. Moreover, we find that the layer stacking method produced higher overall classification accuracies than the PCA fusion method using both the SVM and MLC classifiers. The highest classification accuracy obtained (OA = 97.8%, kappa = 0.964) using the SVM classifier on the layer stacking fusion data at 1 m spatial resolution. Compared with the best classification results of the CASI and LiDAR data alone, the overall classification accuracies improved by 9.1% and 19.6%, respectively. Our findings also demonstrated that the SVM classifier generally performed better than the MLC when classifying multisource data; however, none of the classifiers consistently produced higher accuracies at all spatial resolutions.