Forest canopy height estimation using ICESat/GLAS data and error factor analysis in Hokkaido,Japan

Spaceborne light detection and ranging (LiDAR) enables us to obtain information about vertical forest structure directly, and it has often been used to measure forest canopy height or above-ground biomass. However, little attention has been given to comparisons of the accuracy of the different estimation methods of canopy height or to the evaluation of the error factors in canopy height estimation. In this study, we tested three methods of estimating canopy height using the Geoscience Laser Altimeter System (GLAS) onboard NASA’s Ice, Cloud, and land Elevation Satellite (ICESat), and evaluated several factors that affected accuracy. Our study areas were Tomakomai and Kushiro, two forested areas on Hokkaido in Japan. The accuracy of the canopy height estimates was verified by ground-based measurements. We also conducted a multivariate analysis using quantification theory type I (multiple-regression analysis of qualitative data) and identified the observation conditions that had a large influence on estimation accuracy. The method using the digital elevation model was the most accurate, with a root-mean-square error (RMSE) of 3.2 m. However, GLAS data with a low signal-to-noise ratio (⩽10.0) and that taken from September to October 2009 had to be excluded from the analysis because the estimation accuracy of canopy height was remarkably low. After these data were excluded, the multivariate analysis showed that surface slope had the greatest effect on estimation accuracy, and the accuracy dropped the most in steeply sloped areas. We developed a second model with two equations to estimate canopy height depending on the surface slope, which improved estimation accuracy (RMSE = 2.8 m). These results should prove useful and provide practical suggestions for estimating forest canopy height using spaceborne LiDAR.     

Sub-footprint analysis to uncover tree height variation using ICESat/GLAS

Detailed forest height data are an indispensable prerequisite for many forestry and earth science applications. Existing research of using Geoscience Laser Altimeter System (GLAS) data mainly focuses on deriving average or maximum tree heights within a GLAS footprint, i.e. an ellipse with a diameter of 65 m. However, in most forests, it is likely that the tree heights within such ellipse are heterogeneous. Therefore, it is desired to uncover detailed tree height variation within a GLAS footprint. To the best of our knowledge, no such methods have been reported as of now. In this study, we aim to characterize tree heights’ variation within a GLAS footprint as different layers, each of which corresponds to trees with similar heights. As such, we developed a new method that embraces two steps: first, a refined Levenberg–Marquardt (LM) algorithm is proposed to decompose raw GLAS waveform into multiple Gaussian signals, within which it is hypothesized that each vegetation signal corresponds to a particular tree height layer. Second, for each layer, three parameters were first defined: Canopy Top Height (CTH), Crown Length (CL), and Cover Proportion (CP). Then we extracted the three parameters from each Gaussian signal through a defined model. In order to test our developed method, we set up a study site in Ejina, China where the dominant specie is Populus euphratica. Both simulated and field tree height data were adopted. With regard to the simulation data, results presented a very high agreement for the three predefined parameters between our results and simulation data. When our methods were applied to the field data, the respective R² become 0.78 (CTH), CL (R² = 0.76), CP (R² = 0.74). Overall, our studies revealed that large footprint GLAS waveform data have the potentials for obtaining detailed forest height variation.     

A new 500-m resolution map of canopy height for Amazon forest using spaceborne LiDAR and cloud-free MODIS imagery

In the present study, we aimed to map canopy heights in the Brazilian Amazon mainly on the basis of spaceborne LiDAR and cloud-free MODIS imagery with a new method (the Self-Organizing Relationships method) for spatial modeling of the LiDAR footprint. To evaluate the general versatility, we compared the created canopy height map with two different canopy height estimates on the basis of our original field study plots (799 plots located in eight study sites) and a previously developed canopy height map. The compared canopy height estimates were obtained by: (1) a stem diameter at breast height (D) – tree height (H) relationship specific to each site on the basis of our original field study, (2) a previously developed D–H model involving environmental and structural factors as explanatory variables (Feldpausch et al., 2011), and (3) a previously developed canopy height map derived from the spaceborne LiDAR data with different spatial modeling method and explanatory variables (Simard et al., 2011). As a result, our canopy height map successfully detected a spatial distribution pattern in canopy height estimates based on our original field study data (r = 0.845, p = 8.31 × 10⁻³) though our canopy height map showed a poor correlation (r = 0.563, p = 0.146) with the canopy height estimate based on a previously developed model by Feldpausch et al. (2011). We also confirmed that the created canopy height map showed a similar pattern with the previously developed canopy height map by Simard et al. (2011). It was concluded that the use of the spaceborne LiDAR data provides a sufficient accuracy in estimating the canopy height at regional scale.     

大光斑激光雷达数据的森林冠层高度反演

针对传统的LM波形分解算法在GLAS大光斑波形数据处理中容易陷于局部最优解,限制了GLAS大光斑激光雷达数据在森林结构参数反演方面应用的问题,该文结合GLAS大光斑数据特征,引进优化后的EM算法对大光斑全波形数据进行分解,获取波形参数最优值。结合波形前缘长度和波形后缘长度,建立树高反演模型,并与LM分解算法建立的模型进行对比分析。研究结果表明,通过改进的EM算法对GLAS大光斑激光雷达数据进行处理,波形特征参数的获取更为精确,达到了较高的树高反演精度。     

基于ICESatGLAS的云南省森林地上生物量反演

结合机载、星载激光雷达对GLAS(地球科学激光测高系统)光斑范围内的森林地上生物量进行估测,并利用MODIS植被产品以及MERIS土地覆盖产品进行了云南省森林地上生物量的连续制图。机载LiDAR扫描的260个训练样本用于构建星载GLAS的森林地上生物量估测模型,模型的决定系数(R2)为0.52,均方根误差(RMSE)为31Mg/ha。研究结果显示,云南省总森林地上生物量为12.72亿t,平均森林地上生物量为94Mg/ha。估测的森林地上生物量空间分布情况与实际情况相符,森林地上生物量总量与基于森林资源清查数据的估测结果相符,表明了利用机载LiDAR与星载ICESatGLAS结合进行大区域森林地上生物量估测的可靠性。     

多源遥感数据支持下的区域性森林冠层高度估测

结合机载LiDAR数据,提出了一种改进的GLAS光斑点冠层高度地形校正模型,以校正后的GLAS光斑点作为输入样本,结合MODIS遥感影像,利用支持向量回归（SVR）的方法对研究区森林冠层高度进行分生态区估测,并利用野外调查数据和机载LiDAR冠层高度结果对估测结果进行验证。结果显示:研究区的坡度等级直接影响GLAS光斑点森林冠层高度估测精度,改进的地形校正模型可以较好的减小坡度对GLAS光斑点森林冠层高度估测的影响,模型精度RMSE稳定在3.25~3.48 m;不同生态分区的SVR模型估测精度较为稳定,其RMSE=6.41~7.56 m;与算数平均高相比,样地的Lorey's高与制图结果拟合最好,不同生态分区平均估测精度为80.3%。机载LiDAR冠层高度结果的验证平均精度为79.5%,和Lorey's高验证结果呈现较好的一致性。     

基于ICESat-GLAS完整波形的坡地森林冠层高度反演研究——以吉林长白山林区为例

为了评价大光斑激光雷达在坡地反演森林冠层高度的能力,以长白山汪清林区为例,利用冰、云和陆地高程卫星(ICESat)上搭载的地学激光测高系统(GLAS)提供的完整波形数据,在不同地形坡度范围内,以激光波形长度和地形指数为变量建立森林冠层高度反演模型。研究结果显示,在0~5°坡度范围所建立的模型能够解释90%的冠层高度,随着坡度范围的扩大,相对应的模型反演精度显著下降。总体来讲,在地形坡度不大于15°的地区,GLAS波形能够比较合理地反演森林冠层高度,具有很大的应用潜力。     

An improved method for estimating forest canopy height using ICESat-GLAS full waveform data over sloping terrain: A case study in Changbai mountains,China

Light Detection And Ranging (LiDAR) has a unique capability for estimating forest canopy height, which has a direct relationship with, and can provide better understanding of the aboveground forest carbon storage. The full waveform data of the large-footprint LiDAR Geoscience Laser Altimeter System (GLAS) onboard the Ice, Cloud, and land Elevation Satellite (ICESat), combined with field measurements of forest canopy height, were employed to achieve improved estimates of forest canopy height over sloping terrain in the Changbai mountains region, China. With analyzing ground-truth experiments, the study proposed an improved model over Lefsky's model to predict maximum canopy height using the logarithmic transformation of waveform extent and elevation change as independent variables. While Lefsky's model explained 8–89% of maximum canopy height variation in the study area, the improved model explained 56–92% of variation within the 0–30° terrain slope category. The results reveal that the improved model can reduce the mixed effects caused by both sloping terrain and rough land surface, and make a significant improvement for accurately estimating maximum canopy height over sloping terrain.     

激光雷达回波模型辅助的坡地森林冠层高度反演

大光斑激光雷达数据已广泛应用于森林冠层高度提取,但通常仅限于地形坡度小于20°的平缓地区。在地形坡度大于20°的陡峭山区,地形引起的波形展宽使得地面回波和植被回波信息混合在一起,给森林冠层高度提取带来巨大挑战。本文利用激光雷达回波模型和地形信息,提出了一种模型辅助的坡地森林冠层高度反演算法。该方法以激光雷达回波信号截止点为参考,定义了波形高度指数H50和H75,使用激光雷达回波模型与已知地形信息模拟裸地的激光雷达回波,将裸地回波信号截止点与森林激光雷达回波信号截止点对齐,利用裸地回波计算常用的波形相对高度指数RH50和RH75,对森林冠层高度进行反演。并与高斯波形分解法和波形参数法的反演结果进行了比较。研究结果表明:(1)利用所提取的波形指数RH50和RH75对胸高断面积加权平均高(Lorey’s height)进行了估算,在坡度小于20°时,高斯波形分解法、波形参数法和模型辅助法的估算结果与实测值线性拟合的相关系数(R2)分别为0.70,0.78和0.98,对应的均方根误差(RMSE)分别为2.90 m,2.48 m和0.60 m,模型辅助法略优于其他两种方法;(2)在坡度大于20°时,高斯波形分解法、波形参数法和模型辅助法的R2分别为0.14,0.28和0.97,相应的RMSE分别为4.93 m,4.53 m和0.81 m,模型辅助法明显优于其他两种方法;(3)在0°—40°时,模型辅助法对Lorey’s height估算结果与实测值的R2为0.97,RMSE为0.80 m。本研究提出的模型辅助法具有更好的地形适应性,在0°—40°的坡度范围内具备对坡地森林冠层高度反演的潜力。     

Retrieving aboveground biomass of wetland Phragmites australis (common reed) using a combination of airborne discrete-return LiDAR and hyperspectral data

Wetland biomass is essential for monitoring the stability and productivity of wetland ecosystems. Conventional field methods to measure or estimate wetland biomass are accurate and reliable, but expensive, time consuming and labor intensive. This research explored the potential for estimating wetland reed biomass using a combination of airborne discrete-return Light Detection and Ranging (LiDAR) and hyperspectral data. To derive the optimal predictor variables of reed biomass, a range of LiDAR and hyperspectral metrics at different spatial scales were regressed against the field-observed biomasses. The results showed that the LiDAR-derived H_p99 (99th percentile of the LiDAR height) and hyperspectral-calculated modified soil-adjusted vegetation index (MSAVI) were the best metrics for estimating reed biomass using the single regression model. Although the LiDAR data yielded a higher estimation accuracy compared to the hyperspectral data, the combination of LiDAR and hyperspectral data produced a more accurate prediction model for reed biomass (R² = 0.648, RMSE = 167.546 g/m², RMSE_r = 20.71%) than LiDAR data alone. Thus, combining LiDAR data with hyperspectral data has a great potential for improving the accuracy of aboveground biomass estimation.     

GLAS星载激光雷达和Landsat/ETM+数据的森林生物量估算

基于大脚印激光雷达数据和野外观测数据,该文提出一种获取脚印点内森林生物量的新思路,并结合陆地卫星数据应用于长白山地区森林地上生物量估算。首先,基于3种森林类型(针叶林、阔叶林和针阔混交林),采用多元逐步回归方法建立激光雷达波形指数与脚印点内实测平均树高的回归模型,估算全部脚印点内的平均树高;然后根据脚印点内样方的野外观测数据(平均树高和平均胸径)以及它们与样方生物量的拟合方程估算没有野外调查数据对应的脚印点的生物量;最后对3种森林类型的脚印点森林生物量在各森林覆盖度条件下进行分层分区统计得到生物量等级图。验证比较遥感估算的生物量与野外调查数据推算的生物量,总体误差在0~30(t·hm~(-2))之间,均方根误差为14.66(t·hm~(-2))。     

Scaling estimates of vegetation structure in Amazonian tropical forests using multi-angle MODIS observations

Detailed knowledge of vegetation structure is required for accurate modelling of terrestrial ecosystems, but direct measurements of the three dimensional distribution of canopy elements, for instance from LiDAR, are not widely available. We investigate the potential for modelling vegetation roughness, a key parameter for climatological models, from directional scattering of visible and near-infrared (NIR) reflectance acquired from NASA’s Moderate Resolution Imaging Spectroradiometer (MODIS). We compare our estimates across different tropical forest types to independent measures obtained from: (1) airborne laser scanning (ALS), (2) spaceborne Geoscience Laser Altimeter System (GLAS)/ICESat, and (3) the spaceborne SeaWinds/QSCAT. Our results showed linear correlation between MODIS-derived anisotropy to ALS-derived entropy (r² = 0.54, RMSE = 0.11), even in high biomass regions. Significant relationships were also obtained between MODIS-derived anisotropy and GLAS-derived entropy (0.52 ≤ r² ≤ 0.61; p < 0.05), with similar slopes and offsets found throughout the season, and RMSE between 0.26 and 0.30 (units of entropy). The relationships between the MODIS-derived anisotropy and backscattering measurements (σ⁰) from SeaWinds/QuikSCAT presented an r² of 0.59 and a RMSE of 0.11. We conclude that multi-angular MODIS observations are suitable to extrapolate measures of canopy entropy across different forest types, providing additional estimates of vegetation structure in the Amazon.     

吉林长白山森林冠顶高度激光雷达与MERSI联合反演

将激光雷达与光学遥感相结合进行区域林分冠顶高度联合反演,提出了大脚印激光雷达GLAS脚点波形数据处理和不同地形条件下的森林冠顶高度反演算法,并建立了区域尺度不同森林类型林分冠顶高度GLAS＋MERSI联合反演模型,制作了长白山地区森林冠顶高度图。     

An accurate and computationally efficient algorithm for ground peak identification in large footprint waveform LiDAR data

Large footprint waveform LiDAR sensors have been widely used for numerous airborne studies. Ground peak identification in a large footprint waveform is a significant bottleneck in exploring full usage of the waveform datasets. In the current study, an accurate and computationally efficient algorithm was developed for ground peak identification, called Filtering and Clustering Algorithm (FICA). The method was evaluated on Land, Vegetation, and Ice Sensor (LVIS) waveform datasets acquired over Central NY. FICA incorporates a set of multi-scale second derivative filters and a k-means clustering algorithm in order to avoid detecting false ground peaks. FICA was tested in five different land cover types (deciduous trees, coniferous trees, shrub, grass and developed area) and showed more accurate results when compared to existing algorithms. More specifically, compared with Gaussian decomposition, the RMSE ground peak identification by FICA was 2.82 m (5.29 m for GD) in deciduous plots, 3.25 m (4.57 m for GD) in coniferous plots, 2.63 m (2.83 m for GD) in shrub plots, 0.82 m (0.93 m for GD) in grass plots, and 0.70 m (0.51 m for GD) in plots of developed areas. FICA performance was also relatively consistent under various slope and canopy coverage (CC) conditions. In addition, FICA showed better computational efficiency compared to existing methods. FICA’s major computational and accuracy advantage is a result of the adopted multi-scale signal processing procedures that concentrate on local portions of the signal as opposed to the Gaussian decomposition that uses a curve-fitting strategy applied in the entire signal. The FICA algorithm is a good candidate for large-scale implementation on future space-borne waveform LiDAR sensors.     

Canopy foliar nitrogen retrieved from airborne hyperspectral imagery by correcting for canopy structure effects

A statistical relationship between canopy mass-based foliar nitrogen concentration (%N) and canopy bidirectional reflectance factor (BRF) has been repeatedly demonstrated. However, the interaction between leaf properties and canopy structure confounds the estimation of foliar nitrogen. The canopy scattering coefficient (the ratio of BRF and the directional area scattering factor, DASF) has recently been suggested for estimating %N as it suppresses the canopy structural effects on BRF. However, estimation of %N using the scattering coefficient has not yet been investigated for longer spectral wavelengths (>855 nm). We retrieved the canopy scattering coefficient for wavelengths between 400 and 2500 nm from airborne hyperspectral imagery, and then applied a continuous wavelet analysis (CWA) to the scattering coefficient in order to estimate %N. Predictions of %N were also made using partial least squares regression (PLSR). We found that %N can be accurately retrieved using CWA (R² = 0.65, RMSE = 0.33) when four wavelet features are combined, with CWA yielding a more accurate estimation than PLSR (R² = 0.47, RMSE = 0.41). We also found that the wavelet features most sensitive to %N variation in the visible region relate to chlorophyll absorption, while wavelet features in the shortwave infrared regions relate to protein and dry matter absorption. Our results confirm that %N can be retrieved using the scattering coefficient after correcting for canopy structural effect. With the aid of high-fidelity airborne or upcoming space-borne hyperspectral imagery, large-scale foliar nitrogen maps can be generated to improve the modeling of ecosystem processes as well as ecosystem-climate feedbacks.     

面向森林植被的星载大光斑激光雷达回波仿真

利用星载激光雷达的大光斑全波形数据估测植被结构参数、监测森林生态已受到广泛关注。为了更准确地理解森林植被的结构参数和光学特性对激光雷达回波波形的影响,利用实测森林植被数据提取植被空间分布的统计规律,考虑地形坡度变化和植被冠层反射特性的影响,生成参数化的森林植被空间轮廓反射模型,结合星载激光雷达的回波理论,建立了面向植被的星载激光雷达波形仿真器。由大兴安岭地区的实测植被数据提取的统计规律生成的森林目标仿真波形与地球科学激光测高仪系统（Geoscience Laser Altimeter System,GLAS）真实回波波形具有较好的一致性,平均相关系数R2达到0.91。通过波形仿真分析发现,光斑尺寸减小有利于大坡度地形的森林信息反演,研究成果对中国未来研制星载激光雷达载荷的系统参数设计具有参考意义。     

Fine-resolution forest tree height estimation across the Sierra Nevada through the integration of spaceborne LiDAR,airborne LiDAR,and optical imagery

ABSTRACT

Forests of the Sierra Nevada (SN) mountain range are valuable natural heritages for the region and the country, and tree height is an important forest structure parameter for understanding the SN forest ecosystem. There is still a need in the accurate estimation of wall-to-wall SN tree height distribution at fine spatial resolution. In this study, we presented a method to map wall-to-wall forest tree height (defined as Lorey’s height) across the SN at 70-m resolution by fusing multi-source datasets, including over 1600 in situ tree height measurements and over 1600?km² airborne light detection and ranging (LiDAR) data. Accurate tree height estimates within these airborne LiDAR boundaries were first computed based on in situ measurements, and then these airborne LiDAR-derived tree heights were used as reference data to estimate tree heights at Geoscience Laser Altimeter System (GLAS) footprints. Finally, the random forest algorithm was used to model the SN tree height from these GLAS tree heights, optical imagery, topographic data, and climate data. The results show that our fine-resolution SN tree height product has a good correspondence with field measurements. The coefficient of determination between them is 0.60, and the root-mean-squared error is 5.45?m.     

卫星激光测高植被目标回波的半解析仿真模型

目前激光星载激光测高仪已被广泛应用于植被目标特征的提取,证明了激光雷达在林业行业的巨大应用潜力。植被目标的回波波形复杂,本文提出了一种基于半解析法植被目标的回波仿真模型,可以较好地仿真特定输入参数产生的波形。使用GLAS测高系统经过大兴安岭区域的激光波形和实地林木样地参数为依据,仿真波形与GLAS回波的相关系数R²均达到0.91以上。利用回波仿真模型定量控制如冠层几何形状、区域坡度、地面粗糙度和林下植被的变化并快速获取大量波形的优势,独立地分析每个因素对回波波形的影响,为植被目标反演的数据源选取及对回波波形展宽的分析提供指导意见。     

Combined use of LiDAR data and multispectral earth observation imagery for wetland habitat mapping

Although wetlands play a key role in controlling flooding and nonpoint source pollution, sequestering carbon and providing an abundance of ecological services, the inventory and characterization of wetland habitats are most often limited to small areas. This explains why the understanding of their ecological functioning is still insufficient for a reliable functional assessment on areas larger than a few hectares. While LiDAR data and multispectral Earth Observation (EO) images are often used separately to map wetland habitats, their combined use is currently being assessed for different habitat types. The aim of this study is to evaluate the combination of multispectral and multiseasonal imagery and LiDAR data to precisely map the distribution of wetland habitats. The image classification was performed combining an object-based approach and decision-tree modeling. Four multispectral images with high (SPOT-5) and very high spatial resolution (Quickbird, KOMPSAT-2, aerial photographs) were classified separately. Another classification was then applied integrating summer and winter multispectral image data and three layers derived from LiDAR data: vegetation height, microtopography and intensity return. The comparison of classification results shows that some habitats are better identified on the winter image and others on the summer image (overall accuracies = 58.5 and 57.6%). They also point out that classification accuracy is highly improved (overall accuracy = 86.5%) when combining LiDAR data and multispectral images. Moreover, this study highlights the advantage of integrating vegetation height, microtopography and intensity parameters in the classification process. This article demonstrates that information provided by the synergetic use of multispectral images and LiDAR data can help in wetland functional assessment     

Estimates of forest structure parameters from GLAS data and multi-angle imaging spectrometer data

Quantitative estimates of forest vertical and spatial distribution using remote sensing technology play an important role in better understanding forest ecosystem function, forest carbon storage and the global carbon cycle. Although most remote sensing systems can provide horizontal distribution of canopies, information concerning the vertical distribution of canopies cannot be detected. Fortunately, laser radars have become available, such as GLAS (Geoscience laser altimeter system). Because laser radar can penetrate foliage, it is superior to other remote sensing technologies for detecting vertical forest structure and has higher accuracy. GLAS waveform data were used in this study to retrieve average tree height and biomass in a GLAS footprint area in Heilongjiang Province. However, GLAS data are not spatially continuous. To fill the gaps, MISR (multi- angle imaging spectrometer) spectral radiance was chosen to predict the regional continuous tree height by developing a multivariate linear regression model. We compared tree height estimated by the regression model and GLAS data. The results confirmed that estimates of tree height and biomass based on GLAS data are considerably more accurate than estimates based on traditional methods. The accuracy is approximately 90%. MISR can be used to estimate tree height in continuous areas with a robust regression model. The R², precision and root mean square error of the regression model were 0.8, 83% and 1 m, respectively. This study provides an important reference for mapping forest vertical parameters.