Next Article in Journal
The Impact of Moss Species and Biomass on the Growth of Pinus sylvestris Tree Seedlings at Different Precipitation Frequencies
Next Article in Special Issue
Cross-Correlation of Diameter Measures for the Co-Registration of Forest Inventory Plots with Airborne Laser Scanning Data
Previous Article in Journal
Re-Greening Ethiopia: History, Challenges and Lessons
Previous Article in Special Issue
Estimation of the Timber Quality of Scots Pine with Terrestrial Laser Scanning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Correlating the Horizontal and Vertical Distribution of LiDAR Point Clouds with Components of Biomass in a Picea crassifolia Forest

1
The State Key Laboratory of Remote Sensing Science, Institute of Remote Sensing and Digital Earth, CAS Olympic S&T Park No. 20 Da Tun Road P.O. Box 9718, China
2
University of Chinese Academy of Sciences, Beijing 100039, China
3
Fujian Agriculture and Forestry University, Fuzhou 350002, China
*
Author to whom correspondence should be addressed.
Forests 2014, 5(8), 1910-1930; https://doi.org/10.3390/f5081910
Submission received: 18 December 2013 / Revised: 24 May 2014 / Accepted: 24 July 2014 / Published: 5 August 2014

Abstract

:
Light detection and ranging (LiDAR) has been widely used to estimate forest biomass. In this study, we aim to further explore this capability by correlating horizontal and vertical distribution of LiDAR data with components of biomass in a Picea crassifolia forest. Airborne small footprint full-waveform data were decomposed to acquire higher density point clouds. We calculated LiDAR metrics at the tree level and subplot level and correlated them to biomass components, including branch biomass (BB), leaf biomass (LB) and above-ground biomass (AGB), respectively. A new metric (Horizcv) describing the horizontal distribution of point clouds was proposed. This metric was found to be highly correlated with canopy biomass (BB and LB) at the tree level and subplot level. Correlation between AGB and Horizcv at the subplot level is much lower than that at tree level. AGB for subplot is highly correlated with the mean height metric (Hmean), canopy cover index (CCI) and the product of them. On the other hand, the relationship between the vertical distribution of LiDAR point and biomass was explored by developing two types of vertical profiles, including LiDAR distribution profiles and a biomass profile. Good relationships were found between these two types of vertical profiles and assessed by Pearson’s correlation coefficient (R) and the area of overlap index (AOI). These good correlations possess potential in predicting the vertical distribution of canopy biomass. Overall, it is concluded that not only the vertical, but also the horizontal distribution of LiDAR points should be taken into account in estimating components of biomass by LiDAR.

1. Introduction

Biomass is commonly used as an important parameter for the assessment of forest productivity and carbon sequestration rates [1]. It is generally estimated based on ground measurement [2], which is time-consuming and sometimes inefficient. Using remotely sensed data has become a preferable way to estimate the amount and spatial distribution of forest biomass [3]. Biomass shows significant correlation with the structure of forest, such as diameter at breast height (DBH) and tree height. Light detection and ranging (LiDAR), as an active remote sensing technology, exhibits immense potential in detecting the vertical and horizontal structure of forests. Therefore, LiDAR data has been widely used in the estimation of different types of forest biomass at various sites. The types of forest biomass cover tropical rainforest [4,5,6], deciduous and conifer forest in temperate area [3,7,8], and so on [9,10]. Height metrics, such as mean of height (Hmean), mode of height (Hmode) and percentile height, are the most commonly used LiDAR metrics to estimate biomass, because height metrics can be directly measured by LiDAR and can give insights into the vertical structure of individual trees and forest stands. Numerous allometric equations have been developed for different forest types in various areas when using LiDAR metrics to estimate biomass. However, most of the published literature has emphasized the effects of the vertical structure of forest stand on forest biomass. The horizontal structure of forest was described by the canopy cover index (CCI) or penetration index calculated using the simple ratio of canopy return or ground return [4,5,7,8,9]. Ni-Meister [11] indicated that the gap fraction combined with height metrics, which is closely related to the woody volume, can improve the aboveground biomass estimation at large scales by using large-footprint waveform LiDAR. The woody volume approximately equals the product of basal area and top 10% tree height [11]. Both the gap fraction and the woody volume take the vertical (tree height) and horizontal (basal area) distribution of canopy components into account.
For most studies, biomass was estimated at scales no smaller than plot sizes (approximately 0.01 hectare). Fewer studies estimated the biomass of individual trees [1,2,12,13]. However, calculating individual tree biomass is the first step for calculating biomass at larger scales in most studies. In these studies, height and crown diameter were identified from an interpolated raster surface (canopy height model, CHM) by some individual tree identification algorithms before biomass estimation. Not all of the trees could be identified, because some returns from the inner canopy under the raster surface were missing [14,15]. However, with the success of individual tree identification, it was possible to achieve a high resolution and to precisely map the biomass distribution, which can be used as a reference dataset instead of laborious ground measurements [13].
On the other hand, precisely estimated components of forest biomass improve the calculation of carbon storage. Branch biomass (BB) and leaf biomass (LB) are of vital importance to the estimation of canopy fuel load and monitoring of forest fire behavior. They are commonly used in the calculation of canopy fuel parameters, such as canopy base height, canopy bulk density and canopy fuel weight [16,17]. Good correlations were found between species-specific biomass components and LiDAR metrics [18]. Compared with aboveground biomass, canopy biomass, including BB and LB, mostly is distributed within the tree crown. Usually, the horizontal distribution range of the tree crown is comparable to or even wider than its crown depth. Therefore, LiDAR metrics that can reflect the horizontal distribution of the canopy possibly can give an indication on the canopy biomass. For instance, the canopy cover index (CCI) calculated using the simple ratio of canopy return has been found to be sensitive to the degree of canopy openness and shows good correlation with aboveground biomass at the plot level [4,5]. However, the relationship between CCI and canopy biomass has not been explored. To the best of our knowledge, the number of LiDAR metrics, which can give insight into the horizontal structure of the canopy, is quite limited among all of the traditional metrics. In addition, the horizontal distribution of the LiDAR point cloud has not been widely investigated. For instance, when tree heights reach a certain height, the canopy increases the speed of growth in the horizontal direction. This leads to the increment of canopy biomass (BB and LB). Different trees with similar vertical structures may have different horizontal structures. In this case, the variation of the horizontal distribution of LiDAR point clouds could be sensitive to the difference or variation of canopy biomass. We expected that when taking the horizontal characterization into account, the biomass estimating models could be improved by using LiDAR metrics. Therefore, investigation is needed to further explore the links between the distribution (both vertical and horizontal) of the LiDAR point cloud and biomass components.
In this paper, we explored the relationship between the distribution variation of LiDAR point cloud and biomass at the tree level and subplot level. Vertical and horizontal variations of the point distribution were investigated to assess their correlations with biomass of different canopy components. New LiDAR metrics were developed and used to estimate biomass components. The results were compared with those estimated using commonly selected metrics alone, which can be considered to deliver good supplementary information to previous studies in this study area [19,20].

2. Experimental Section

2.1. Study Sites

The study area is a water conservation forest site located in Qilian Mountain Dayekou, Zhangye City, Gansu Province, China (Figure 1). It is a study area chosen by an airborne, satellite-borne and ground-based remote sensing experiment, Watershed Allied Telemetry Experimental Research (WATER) [21]. It covers three vegetative climate zones and significantly affects soil/water conservation and biodiversity protection. The area is the main body of the ecosystem in the Qilian Mountains, and the Hexi Corridor Oasis and is located at an elevation of 2500 m to 3800 m above sea level. The center coordinates of the area are 100°15’E and 38°32’N. The main species in the study area is P. crassifolia, occupying about 95% of the whole forestland over the test site [3].
Figure 1. Airborne ortho charged coupled device (CCD) image for the study area and the field plots.
Figure 1. Airborne ortho charged coupled device (CCD) image for the study area and the field plots.
Forests 05 01910 g001

2.2. LiDAR Data

Small footprint full-waveform LiDAR data were acquired for the Dayekou study area on June 23, 2008. The LiteMapper-5600 airborne laser scanning system was used and configured to emit laser pulses of 3.5 ns at 50 kHz in the near-infrared band (wavelength = 1550 nm) with a beam divergence fixed with 0.5 mrad. The average footprint diameter was 38 cm. Geographical coordinates (easting, northing and elevation) were recorded by a dual-frequency differential Global Positioning System with an inertial measurement unit. The horizontal and vertical accuracy of the instrument was 0.1 m and 0.3 m, respectively. Repetitive flights were carried out over the sample study area with a nominal height above ground of 700–800 m [22], leading to a pulse density of approximately 3.43 pts/m2. For each pulse, the transmitted pulse and the returned pulse were recorded in the full-waveform data. Point clouds were generated from the full-waveform data by the commercial software provided by the vendor. All of the points were geo-referenced to a projection system of universal transverse Mercator (UTM) Zone 47N/WGS-84. In order to increase the point density in this study, we produced another set of point clouds from the original waveform data by waveform decomposing. This is described in Section 2.5. The point cloud with the relatively lower point density described above was not used in this study.

2.3. Field Measurement

During the flight campaign, field measurements were conducted on two field plots distributing over the study area (Figure 1). The field plots were denoted as super-plot and line-plot. The super-plot was a 100 m × 100 m area with a slope of less than 20 degrees. It was divided into 16 subplots measuring 25 m × 25 m. The line-plot was a sample line consisting of 19 subplots measuring 20 m × 20 m distributed along the direction of flight. The distance between the subplots in the line-plot was 50 m. The total number of subplots was 35. The center location of each subplot was positioned using a differential global positioning system (DGPS) station. The DGPS mobile station was placed within the open space of the standing forest to ensure the positioning accuracy. The individual trees were located by the survey total stations with equipment model of TOPCON GTS-602 and TOPCON GPT-7002 (Figure 2). Tree height (H) and first branch height (FBH) were measured using a laser hypsometer; diameter at breast height (DBH) was measured with a tape measure, and crown size (R) in two orthotropic directions for each individual tree was measured using a ribbon tape. High correlations were found between the ground-measured features, as shown in Table 1 and Figure 3.
Figure 2. The spatial distribution of the field-measured individual trees in the super-plot with 16 subplots.
Figure 2. The spatial distribution of the field-measured individual trees in the super-plot with 16 subplots.
Forests 05 01910 g002
Table 1. Regression models of different ground-measured features.
Table 1. Regression models of different ground-measured features.
Dependent Variables (y)EquationR2Root Mean Squared Error (RMSE) (m)
Heighty = 0.5171 × DBH + 2.12680.831.24
Crown Sizey = 7.3496 × ln(DBH) − 8.5810.871.18
y = 0.2014 × Height + 1.37420.610.77
y = 0.1264 × DBH + 1.48940.740.71
Note: R2 = coefficient of determination.
Figure 3. Scatter plots for the field-measured attributes (a–d). Height is the tree height; DBH is the diameter at breast height; ln(DBH) is the log transformation of DBH; R2 is the coefficient of determination.
Figure 3. Scatter plots for the field-measured attributes (a–d). Height is the tree height; DBH is the diameter at breast height; ln(DBH) is the log transformation of DBH; R2 is the coefficient of determination.
Forests 05 01910 g003

2.4. Biomass Calculation

The biomass of different forest components was calculated based on their DBH and height according to the empirical relative growth equations [23] as follows:
stock biomass = 0.0478 × ((DBH)2 × H)0.8665
branch biomass = 0.0061 × ((DBH)2 × H)0.8905
leaf biomass = 0.2650 × ((DBH)2 × H)0.4701
fruit biomass = 0.0342 × ((DBH)2 × H)0.5779
The aboveground biomass (AGB) of each individual tree was calculated as the sum of biomass of all of the above components (stock, branch, leaf and fruit). Subplot biomass is calculated as the sum of ABG of each individual tree in each subplot.

2.5. Processing of LiDAR Data

In order to acquire a higher point density, we used a waveform decomposition approach [24] to decompose the original waveform data into another point cloud dataset. This decomposition approach has been successfully applied to another LiDAR dataset collected in the urban area of Zhangye City around the same point in time during which the WATER experiment was conducted. It approximates the LiDAR waveform as a Gaussian mixture model. A stepwise strategy was used to implement waveform decomposition by defining the maximum iteration number and error expectations before decomposition [24]. The final point density was increased to 6.22 pts/m2 from 3.43 pts/m2. In the next step, ground points were filtered from the decomposed point clouds to generate a digital elevation model (DEM) at 0.5 m resolution. Heights of non-ground points were normalized by calculating the difference between the non-ground point heights and the corresponding interpolated DEM cell heights beneath the points. A digital surface model (DSM) at 0.5 m resolution was generated from the high density point clouds. The maximum height value of the points within each pixel was selected as the value of the interpolated DSM raster cell. Then, a canopy height model (CHM) was derived as the difference between the DEM and DSM.

2.6. Individual Tree Identification

The CHM was smoothed by a 3 m × 3 m Gaussian filter before being used to identify the individual trees. The search radius (3 m) was chosen based on the mean of all of the measured tree crown sizes. Then, the watershed algorithm [25] was applied to the smoothed CHM. The watershed algorithm firstly inverts the CHM to a raster converting peaks into depressions. Then, the canopy basins are delineated on the depression raster by watershed segmentation. Finally, the delineated basins raster is converted to segments. A height break (1.5 m) was set to extract the tree crown segments from the obtained segments. Then the crown segments were transferred into polygonally shaped files in Arcgis9.3. The highest pixel value obtained by zonal statistics within each canopy polygon was assigned as the tree height and its location as the tree location. This automatic delineation greatly helped to extract the individual trees, but manual intervention was still needed to remove some mistakes. Errors of omission and commission were unavoidable in individual tree delineation. It can clearly be seen that there are several field-measured trees clustering within one segment or that no field-measured trees are located in a particular segment. The center geographic location of each segment and its height were matched to the field measured trees accompanied by some visual interpretation. Before matching, visual interpretation and corrections were conducted on the segments containing more than one field-measured tree. Specifically, if the field-measured trees in the segment were too close to separate visually, the segment enclosing them was not selected in order to assure the quality of the crown segment. The difficulty in separating this kind of trees was caused either by tree clustering or the understory trees being sheltered by the tall trees. In the remaining cases, the field-measured trees were separated by splitting the segment manually according to the texture of CHM and the ortho CCD image. On the other hand, segments containing no field-measured trees were not selected. Finally, only the segments with heights and locations closest to the field-measured trees were selected as the final matched crown segments. The quality of tree delineation directly affects the accuracy of LiDAR metrics at tree level in this study. Differences in the locations and heights were limited to certain confidence intervals. Only the final matched crown segments were selected to clip the point clouds to calculate LiDAR metrics at tree level.

2.7. LiDAR Metrics Calculation

The decomposed point clouds within each extracted crown shape files were clipped from the original points. Each point contained attributes, including the geographic coordinates (X/Y), the normalized height (H), returned pulse amplitude, width, standard deviation and returns. The mean height of each individual tree point cloud was calculated using the conventional LiDAR height metrics. The coefficient of height variation (Hcv) is a commonly used LiDAR metric and is used to describe the height variation in the Z-axis [10]. The coefficient of variation (CV) is defined as the ratio of the standard deviation to the mean. The height Z is normalized by subtracting the DEM from DSM. Based on this, we expected to develop a similar, but slightly differing coefficient of variation in the X/Y-axis to describe the horizontal distribution of LiDAR points. Before calculating the new metric, the X/Y coordinates of the point clouds within each tree crown segment were normalized. The original X/Y coordinates were transformed to x/y as follows:
xi = (XiXmin) / (XmaxXmin)
yi = (YiYmin) / (YmaxYmin)
Then, the coefficient of variation for x/ y, xcv / ycv was calculated:
Forests 05 01910 i001
Forests 05 01910 i002
Where σx and σy are the standard deviations for x and y, respectively. Finally, the metric Horizcv was defined as:
Horizcv = xcv × ycv
The normalized x/y coordinates can be seen as the projections of the point clouds in the horizontal plane. They can describe the relative spatial locations of the point clouds within the crown segment in the horizontal plane. The coefficient of x/y variation can reflect how the point cloud distributions vary in the x/y direction. They actually are independent of the absolute locations. Theoretically, this works for tree point clouds (points within the tree crown segment) with coordinates characterized by different coordinate systems. For example, there are two LiDAR point clouds collected from two trees with the exact same structure and biomass located at different geographic positions. The interactions between the LiDAR pulse and the trees are assumed to be the same, and the only difference is the shift in x/y coordinates. There should be a transformation between the coordinates in the x/y plane. For example, we shift the x/y coordinates using the following formulas:
xshift_i = xi × a + b
yshift_i = yi × a + b
where xshift_i and yshift_i are the shifted coordinates, a and b are the shifting coefficients. The normalized coordinate x/y are the same as those before coordinate shifting, because the normalization process has removed the effect of the absolute locations.
Horizcv for each subplot was calculated by summing all of the Horizcv of the individual trees within the subplot. All of the first-returned echoes within each subplot were selected to obtain the canopy cover index (CCI). CCI was calculated by dividing the number of canopy first-returned echoes (H > 1.5 m) by the total number of first-returned echoes. The mean of the point cloud heights (Hmean) was also calculated, since biomass is closely related to tree height.

2.8. Vertical Biomass Profile Derivation

To assess how well the vertical distribution of point clouds correlated with the vertical distribution of biomass, three vertical profiles were calculated: canopy point density profile (CPD), canopy intensity profile (CIP) and canopy biomass density profile (CBP). The first two profiles are LiDAR-derived, and the third one is derived from ground-measured data. Trees in each subplot were sectioned into a series of small height bins (0.5 m). The number/intensity of LiDAR points in each height bin was summed to obtain CPD/CIP. The canopy biomass profiles (CBP) were derived during the derivation process of canopy bulk density (CBD). CBD is a measurement indicating the distribution of canopy biomass or fuel load. It is the dry weight of available canopy fuel per unit canopy volume [17], and it is a bulk property of the stand, not an individual tree. Canopy bulk density is derived following the steps [17]: (1) calculating the leaf and branch biomass for each tree in the subplot; (2) dividing that biomass equally into a 0.5 m bin from the base of the tree crown to the top of tree; (3) calculating the biomass contributed by each tree in the subplot for each 0.5 m bin (a vertical profile of biomass is created after this step, which we refer to as CBP); and (4) using a running mean to get the maximum value in the vertical profile (CBP), whereby the maximum value is selected as CBD. During the quantification of CBP, it is assumed that canopy biomass is distributed uniformly in the vertical direction within the tree crown [17]. Based on this, the vertical biomass distribution within the subplot or stand is approximately calculated by summing the canopy biomass contributed by each tree in each vertical bin. However, we have to admit that it is just an approximate way to describe the vertical distribution of biomass at the plot level or stand level rather than at the individual tree level. In fact, it is very hard or even unfeasible to directly measure the vertical distribution of tree biomass in the field.
Since height differences occur between LiDAR-derived heights and ground-measured heights, there should be height gaps between the CPD/CIP and CBP. Therefore, we limited the pairs for comparison to a height difference to 2 m in order to prevent potential errors. Different types of profiles possess different magnitudes and units. Before comparing and correlating different pairs of vertical profiles, all of the profiles were normalized by dividing the value for each vertical bin by the sum of those for all of the bins.
In order to correlate the LiDAR-derived profiles with the ground-measured profiles, we chose Pearson’s correlation coefficient (R) and the area of overlap index (AOI) to assess the correlation. Pearson’s correlation coefficient indicates how well two profiles match. The area of overlap index is defined as the fraction of area shared by two profiles, and it has been successfully used to compare the waveforms and profile data [5].

3. Results and Discussion

3.1. Individual Tree Identification

Ground-measured tree X/Y locations and their heights were compared to those delineated by the watershed segmentation method, as shown in Figure 4. Trees were matched if the difference between the coordinates and the difference between the tree heights were within certain confidence intervals, accompanied by visual interpretations. Errors of omission and commission for individual trees occurred as in the other studies [14,26,27]. These errors are highly related to the spatial distribution of the trees. Omission error is likely to occur when the heights of some trees are below the dominant canopy nearby [27]. Another important reason for errors of commission is that the only forest species, P. crassifolia, in this study is likely to cluster. In this study, we finally obtained 561 matched trees out of 1468 (38.22%) field-measured trees in the super-plot, as shown in Table 2. This ratio (38.22%) was higher than the 15.75% (117/743) obtained in [13] and lower than the 48% (84/175) obtained in [26]. The statistics of different extracted tree features are shown in Table 3.
Figure 4. Canopy height model overlaid with tree crown segments derived from watershed segmentation: (a) near the super-plot (within yellow boundary); (b) within the super-plot together with measured tree tops.
Figure 4. Canopy height model overlaid with tree crown segments derived from watershed segmentation: (a) near the super-plot (within yellow boundary); (b) within the super-plot together with measured tree tops.
Forests 05 01910 g004
Table 2. Statistics of watershed segmentation results for the 16 subplots in a super-plot.
Table 2. Statistics of watershed segmentation results for the 16 subplots in a super-plot.
SubplotTotal SegmentsMatched SegmentsMismatched SegmentsField MeasuredPercentage (%)
15338159838.78
263461711440.35
3292277728.57
45439158048.75
55038127054.29
6312479126.37
768551312045.83
892751715648.08
9272078722.99
102820810020.00
1155431211138.74
124939109441.49
133827116243.55
144533129435.11
15372897736.36
16181443737.84
Total7375611761,46838.22
Table 3. Statistics of light detection and ranging (LiDAR) extracted tree features matched to those from the ground measured trees.
Table 3. Statistics of light detection and ranging (LiDAR) extracted tree features matched to those from the ground measured trees.
StatisticsErrorX (m)Y (m)H (m)D (m)
Mean (abs)0.280.301.020.46
SD (abs)0.380.411.190.53
Number of Trees561 (38.22%)
Note: X = error of X-coordinate; Y = error of Y-coordinate; H = error of tree height; D = horizontal distance between LiDAR extracted tree and ground-measured tree; Mean = mean of errors; SD = standard deviation of errors; abs = absolute value.

3.2. Correlating Horizontal Variation of Point Clouds with Biomass

Linear regression models were developed to indicate the correlation between the group of LiDAR metrics and the forest components of biomass. The regression models were evaluated using the coefficient of determination (R2), root mean squared error (RMSE) and relative root mean squared error (rRMSE).
Forests 05 01910 i003
Forests 05 01910 i004
where Pi is the observed value and Forests 05 01910 i005 is the predicted value. The value of RMSE is related to the magnitude of the observed variables. Therefore, we also selected rRMSE, a relative value, to compare the performance of different regression models [18]. A lower rRMSE value often indicates a better regression performance. We first explored the relationships between the conventional height metrics and branch biomass (BB), leaf biomass (LB) and aboveground biomass (AGB) at the individual tree level, respectively. Results (Figure 5a–c) show that biomass for each component of individual trees initially increases with the increasing mean point heights (Hmean) and is saturated when Hmean continues increasing. This means that when the tree height is comparably low, tree height contributes greatly to the amount of biomass. Whereas when tree heights reach a certain height, the canopy increases the speed of growth in the horizontal direction, which also increases the canopy biomass (branch and leaf). Therefore, the horizontal distribution of branches and leaves directly affects the amount of canopy biomass. High correlations were found between the new metric (Horizcv) and the biomass components with high R2 and low RMSE and rRMSE, as shown in Figure 5 and Table 4. Horizcv and AGB for the tree show a good relationship with R2 equal to 0.72, but a higher rRMSE (0.53) than that of BB (0.36) and LB (0.40). Specifically, a higher value of Horizcv means more dispersion of the LiDAR points in the horizontal XY plane. Since branches and leaves rarely grow in only one direction, a greater dispersion probably reflects a wider horizontal distance between the branch and leaf growth and a larger weight of canopy components. Combined metrics were developed by combining the horizontal metrics (CCI and Horizcv) with the vertical metric (Hmean). Both CCI and Horizcv can give an insight into the horizontal distribution of point clouds. Correlations were improved when using the combined metrics (Hmean*CCI and Horizcv*CCI) compared to when using the single metric alone, as shown in Table 4 and Figure 6a–f. The RMSEs and rRMSEs also decreased. Although improvements for the combined metric Horizcv*CCI is smaller than that of Hmean*CCI, the regression models for Horizcv*CCI obtained higher R2 and lower RMSEs and rRMSEs. This is due to the high correlation between Horizcv and the biomass components. However, no significant improvements were obtained when combining Hmean, Horizcv and CCI, as shown in Table 4.
Table 4. Regression models developed from biomass components and LiDAR metrics at the tree level.
Table 4. Regression models developed from biomass components and LiDAR metrics at the tree level.
TreeMetricsRegression EquationsR2RMSErRMSE
BBHmeany = 0.0054e0.6842x0.459.780.97
CCIy = 45.674x − 5.70330.538.020.81
Hmean*CCIy = 4.114x − 4.34110.597.470.76
Horizcvy = 69.523x − 0.36260.736.110.36
Hmean*Horizcvy = 6.1408x + 0.39720.736.070.59
Hmean*Horizcv*CCIy = 0.9624x + 3.71160.686.640.67
LBHmeany = 0.2487e0.3612x0.455.940.54
CCIy = 31.734x + 0.2420.525.260.48
Hmean*CCIy = 2.8721x + 1.14080.604.820.49
Horizcvy = 43.369x + 4.56290.674.390.40
Hmean*Horizcvy = 3.8244x + 5.04670.674.360.38
Hmean*Horizcv*CCIy = 0.6429x + 6.94870.635.140.52
AGBHmeany = 0.3638e0.5356x0.4598.430.82
CCIy = 472.47x − 42.5360.5481.810.69
Hmean*CCIy = 42.693x − 28.9190.6175.650.64
Horizcvy = 711.61x + 13.8550.7263.430.53
Hmean*Horizcvy = 62.835x + 21.6630.7262.860.50
Hmean*Horizcv*CCIy = 9.689x + 56.5580.6571.070.59
Note: x = independent variable (LiDAR metrics); y = dependent variable (biomass); BB = branch biomass; LB = leaf biomass; AGB = aboveground biomass; CCI = canopy cover index; R2 = coefficient of determination; RMSE = root mean squared error; rRMSE = relative root mean squared error.
Figure 5. Scatter plots for biomass components (branch biomass (BB) (a, d, g), leaf biomass (LB) (b, e, h) and aboveground biomass (AGB) (c, f, i)) and LiDAR metrics (mean height (Hmean) canopy cover index (CCI) and horizontal distribution of point clouds (Horizcv)) at the individual tree level.
Figure 5. Scatter plots for biomass components (branch biomass (BB) (a, d, g), leaf biomass (LB) (b, e, h) and aboveground biomass (AGB) (c, f, i)) and LiDAR metrics (mean height (Hmean) canopy cover index (CCI) and horizontal distribution of point clouds (Horizcv)) at the individual tree level.
Forests 05 01910 g005
In this study, the accuracy of LiDAR metrics calculated at the tree level highly depends on the quality of individual tree delineation. We only selected the crown segments with a high delineation quality to clip the point cloud to calculate the metrics at tree level in order to reduce the influence from error of tree delineation. Therefore, the LiDAR metrics at tree level were based on the degree of quality of the (selected) crown delineation. Based on this, an irregular delineated crown segment indicated an irregular canopy shape. Vegetation hits from different parts of the crown segment were dispersed, leading to a relatively high coefficient of variation. The more regular the canopy shape, the lower the coefficient of variation for the point cloud within the crown segment.
Figure 6. Scatter plots for biomass components (BB (a, d, g), LB (b, e, h) and AGB (c, f, i)) and the combined LiDAR metrics (Hmean*CCI, Hmean*Horizcv and Hmean*Horizcv*CCI) at the individual tree level.
Figure 6. Scatter plots for biomass components (BB (a, d, g), LB (b, e, h) and AGB (c, f, i)) and the combined LiDAR metrics (Hmean*CCI, Hmean*Horizcv and Hmean*Horizcv*CCI) at the individual tree level.
Forests 05 01910 g006
By summing the biomass of different forest components for each individual tree, three types of biomass (BB, LB and AGB) were calculated at the subplot level. Their relationship with LiDAR metrics was explored as shown in Table 5. Results indicate that the mean height (Hmean) shows good correlation with the three types of biomass, especially with AGB (R2 = 0.65). The canopy cover index (CCI) shows a good relationship with AGB (R2 = 0.60). When considering the product of Hmean and CCI as one metric, R2 for the regression model developed from this combined metric and AGB improved (R2 = 0.71), and RMSE and rRMSE decreased. This improvement was a consequence of taking both horizontal and vertical effects of the forest structure on AGB into account, which, in the end, agreed with Ni-Meister’s study [11].
Table 5. Regression models developed from biomass components and LiDAR metrics at the subplot level.
Table 5. Regression models developed from biomass components and LiDAR metrics at the subplot level.
SubplotMetricsRegression EquationsR2RMSErRMSE
BBHmeany = 101.71x − 426.780.49178.860.31
CCIy = 1,188.4x − 411.790.35202.000.35
Hmean*CCIy = 74.479x − 44.8840.48179.270.32
Horizcvy = 76.639x +7.05690.75126.170.22
Hmean*Horizcvy = 6.4054x + 95.3740.8499.760.18
Hmean*Horizcv*CCIy = 6.4395x + 160.20.83103.420.18
LBHmeany = 70.079x − 62.3660.42227.660.37
CCIy = 1,034.2x − 229.790.25223.40.36
Hmean*CCIy = 56.211x + 160.450.26221.390.36
Horizcvy = 84.016x + 7.7920.84101.920.16
Hmean*Horizcvy = 6.4469x + 147.110.8114.540.18
Hmean*Horizcv*CCIy = 6.4181x + 216.370.78122.020.2
AGBHmeany = 523.44x − 1,120.80.65665.190.17
CCIy = 6,984.2x − 1,760.10.6709.830.18
Hmean*CCIy = 404.22x + 672.210.71597.160.15
Horizcvy = 238.08x + 2,257.80.36894.310.22
Hmean*Horizcvy = 22.58x + 2,3340.52772.420.19
Hmean*Horizcv*CCIy = 24.377x + 2,456.10.59712.30.18
The accumulated Horizcv for the subplot still showed significant correlation with the branch biomass and leaf biomass, while not obviously with AGB. The combined metric Hmean*Horizcv even improved the performance of the BB prediction model with R2 increasing from 0.75 to 0.84, RMSE decreasing from 126.17 to 99.76 and rRMSE decreasing from 0.22 to 0.18. This indicates that the horizontal distribution of branches and leaves still greatly affects the canopy biomass (BB, LB) at the subplot level. The poor correlation between Horizcv and AGB at the subplot level might be due to the following: (1) individual trees were divided into several groups according to the area of extent of each subplot (Horizcv for each plot is calculated as the sum of those area extents of the trees within the subplot); (2) the total biomass of the whole tree (AGB) not only benefited from the contribution of the canopy components (branches and leaves), but also from that of other tree components, such as stock, root and fruit; and (3) a vital important structure attribute has been out of discussion when we have assessed the effects from the horizontal and vertical structure variations on biomass: the diameter at breast height (DBH). According to the relative growth equation, DBH greatly contributes to the estimation of AGB. However, up to now, airborne LiDAR still cannot directly measure the size of DBH. Finally, for the calculation of Horizcv, either at the tree level or subplot level, the influence of the LiDAR data should be considered, for instance the influence of the LiDAR flight properties, i.e., overlapping strips, different flying heights, scan pattern, scan angles, etc.

3.3. Correlating Vertical Distribution of Point Clouds with Biomass

Table 6 shows the correlation between two groups of vertical profiles for each subplot in the study area. The Pearson’s correlation coefficient (R) values for the correlation between CBP and CPD ranged between 0.39 and 0.93 (mean = 0.71, SD = 0.15), and 28 pairs of 35 had values ≥0.6. For pairs of CBP and CID, R varied from 0.19 to 0.92 (mean = 0.62, SD = 0.19), and 21 pairs of 35 obtain values ≥0.6. On the other hand, AOI values for pairs of CBP and CPD varied from 0.50 to 0.84 (mean = 0.64, SD = 0.08), and 21 pairs of 35 had values ≥0.6. CIP played a weaker role than CPD: only 17 of 35 CBP/CIP pairs had values above 0.6 (mean = 0.61, SD = 0.09). The results of the correlations indicate that CBP representing the vertical distribution of biomass is highly correlated with the vertical distribution of LiDAR points. It is likely that both CPD and CIP contain the information of tree height. For the two LiDAR-derived profiles, CPD performs better than CIP. CPDs and CIPs for all of the subplots show very similar curve shapes. The profiles for Subplot S7 are shown in Figure 7. The similarity in the curve shapes for all of the subplots is mainly attributed to the similar forest structure in this study area. According to the previous studies [28], the Picea crassifolia forest in this study area consists primarily of half-mature and mature stands, and the complexity of the canopy vertical distribution and the differentiation grade of canopy height were low. In order to increase the usefulness of the outcome, the approaches should be tested in different forest conditions in the future. The distances between the peaks of CPD (CIP) and CBP for different subplots differ. The height at which these two pairs of profiles start to overlap relates to the stem density and the canopy structure of the Picea crassifolia in the subplot. Specifically, the higher the stem density, the more LiDAR pulses are reflected by the canopy and the more precisely the CPD represents the distribution of the canopy components. In addition, when branches and leaves within the canopy clump at a certain height, the peak of accumulated biomass for all of the height bins and the peak of the LiDAR point density are very likely to overlap around this height. Calculating the vertical biomass profiles is one of the steps to calculate the canopy fuel parameter, canopy bulk density (CBD). It is assumed that the canopy biomass is distributed uniformly in the vertical direction within the tree crown [17]. Based on this, the vertical biomass distribution within the subplot or stand is approximately calculated by summing the canopy biomass contributed by each tree in each vertical bin. However, we have to admit that it is just an approximate way to describe the vertical distribution of biomass at the plot level or stand level, rather than at the individual tree level. According to [17], in fact, canopy fuels have never been measured directly. However, we believe that the potential of predicting CBD, especially when using the vertical detecting capacity of LiDAR, deserves to be tested.
Table 6. Correlation between different pairs of vertical profiles.
Table 6. Correlation between different pairs of vertical profiles.
PlotPearson’s Correlation Coefficient (R)Area of Overlap Index (AOI)
CBP/CPDCBP/CIPCBP/CPDCBP/CIP
S10.670.650.590.59
S20.870.860.690.67
S30.860.840.660.64
S40.820.810.580.58
S50.790.580.660.61
S60.80.780.570.54
S70.910.920.770.77
S80.80.740.620.58
S90.720.70.680.68
S100.830.710.790.77
S110.670.660.550.54
S120.690.660.570.55
S130.780.730.60.58
S140.620.560.60.57
S150.390.190.590.54
L10.770.740.620.6
L20.40.310.560.53
L30.750.640.720.69
L40.760.70.590.56
L50.790.580.760.7
L60.460.330.560.47
L70.840.820.670.66
L80.930.910.750.74
L90.740.550.620.56
L100.880.810.750.71
L110.450.270.660.61
L120.730.680.560.54
L130.490.380.50.45
L140.40.320.540.51
L150.560.440.580.53
L160.710.380.840.77
L170.860.810.770.75
L180.660.60.650.62
L190.730.610.650.62
L200.710.60.710.66
Mean0.710.620.640.61
SD0.150.190.080.09
Note: S = subplot in super-plot; L=subplot in line-plot; CPD = canopy point density profile; CBP = canopy biomass profile; CIP = canopy intensity profile; Mean = mean of R2 or AOI for all of the plots; SD = standard deviation; numbers in bold represent the maximum or minimum values.
Figure 7. Vertical profiles for Subplot S7. CBP = canopy biomass profile; CPD = canopy point density profile; CIP = canopy intensity profile; normalized crown amplitude means the normalized values for CBP, CBD and CIP.
Figure 7. Vertical profiles for Subplot S7. CBP = canopy biomass profile; CPD = canopy point density profile; CIP = canopy intensity profile; normalized crown amplitude means the normalized values for CBP, CBD and CIP.
Forests 05 01910 g007

4. Conclusions

This paper tried to explore the correlation between the distribution of LiDAR point clouds and the biomass of different forest components. A new metric (Horizcv) describing the horizontal variation of LiDAR points was proposed and proven to be highly correlated with the canopy branch biomass and leaf biomass at both the individual tree and the subplot levels. Aboveground biomass (AGB) at the subplot level was more strongly correlated with the conventional LiDAR metrics (Hmean, CCI). The correlation was improved by combining these two metrics (Hmean*CCI) into one in order to predict the biomass components. On the other hand, CPDs and CIPs can be good curves to represent the vertical distribution of LiDAR points. Good relationships between the LiDAR-derived profiles (CPD, CIP) and the ground-measured-derived biomass profile (CBP) were found. These good correlations can be used to predict the vertical distribution of the canopy biomass. If this prediction can be put into practice, a new and promising way to obtain canopy bulk density (CBD) will be available, which is of critical significance for forest fire management. Above all, this study demonstrates that when we estimate forest biomass by using LiDAR metrics, not only the vertical, but also the horizontal distribution of LiDAR points should be taken into account.

Acknowledgments

This research was funded by China’s Special Funds for Major State Basic Research Project (2013CB733405, 2010CB950603), and the project was supported by the National Natural Science Foundation of China (41201345). We thank all of the people working on the projects, especially for their efforts in the WATER experiment. We are deeply grateful to the four anonymous reviewers for the very constructive and insightful comments that greatly improved this paper.

Author Contributions

Wang Li was the lead author of the article and made the analyses with Zheng Niu and Shuai Gao. Ni Huang contributed expertise in statistical analysis and spatial analysis. Hanyue Chen provided expertise in forest structure and management. The article was improved by the contributions of all of the co-authors at various stages of the analysis and writing process.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Allouis, T.; Durrieu, S.; Vega, C.; Couteron, P. Stem Volume and Above-Ground Biomass Estimation of Individual Pine Trees From LiDAR Data: Contribution of Full-Waveform Signals. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2013, 6, 924–934. [Google Scholar] [CrossRef]
  2. Bortolot, Z.J.; Wynne, R.H. Estimating forest biomass using small footprint LiDAR data: An individual tree-based approach that incorporates training data. ISPRS J. Photogramm. Remote Sens. 2005, 59, 342–360. [Google Scholar] [CrossRef]
  3. He, Q.-S.; Cao, C.-X.; Chen, E.-X.; Sun, G.-Q.; Ling, F.-L.; Pang, Y.; Zhang, H.; Ni, W.-J.; Xu, M.; Li, Z.-Y.; et al. Forest stand biomass estimation using ALOS PALSAR data based on LiDAR-derived prior knowledge in the Qilian Mountain, western China. Int. J. Remote Sens. 2012, 33, 710–729. [Google Scholar] [CrossRef]
  4. Drake, J.B.; Dubayah, R.O.; Clark, D.B.; Knox, R.G.; Blair, J.B.; Hofton, M.A.; Chazdon, R.L.; Weishampel, J.F.; Prince, S.D. Estimation of tropical forest structural characteristics using large-footprint lidar. Remote Sens. Environ. 2002, 79, 305–319. [Google Scholar] [CrossRef]
  5. Drake, J.B.; Dubayah, R.O.; Knox, R.G.; Clark, D.B.; Blair, J.B. Sensitivity of large-footprint lidar to canopy structure and biomass in a neotropical rainforest. Remote Sens. Environ. 2002, 81, 378–392. [Google Scholar] [CrossRef]
  6. Clark, M.L.; Roberts, D.A.; Ewel, J.J.; Clark, D.B. Estimation of tropical rain forest aboveground biomass with small-footprint lidar and hyperspectral sensors. Remote Sens. Environ. 2011, 115, 2931–2942. [Google Scholar] [CrossRef]
  7. Lefsky, M.A.; Cohen, W.B.; Harding, D.J.; Parker, G.G.; Acker, S.A.; Gower, S.T. Lidar remote sensing of above-ground biomass in three biomes. Glob. Ecol. Biogeogr. 2002, 11, 393–399. [Google Scholar] [CrossRef]
  8. Ahmed, R.; Siqueira, P.; Hensley, S. A study of forest biomass estimates from lidar in the northern temperate forests of New England. Remote Sens. Environ. 2013, 130, 121–135. [Google Scholar] [CrossRef]
  9. Huang, W.; Sun, G.; Dubayah, R.; Cook, B.; Montesano, P.; Ni, W.; Zhang, Z. Mapping biomass change after forest disturbance: Applying LiDAR footprint-derived models at key map scales. Remote Sens. Environ. 2013, 134, 319–332. [Google Scholar] [CrossRef]
  10. Næsset, E.; Bollandsås, O.M.; Gobakken, T.; Gregoire, T.G.; Ståhl, G. Model-assisted estimation of change in forest biomass over an 11year period in a sample survey supported by airborne LiDAR: A case study with post-stratification to provide “activity data”. Remote Sens. Environ. 2013, 128, 299–314. [Google Scholar] [CrossRef]
  11. Ni-Meister, W.; Lee, S.; Strahler, A.H.; Woodcock, C.E.; Schaaf, C.; Yao, T.; Ranson, K.J.; Sun, G.; Blair, J.B. Assessing general relationships between aboveground biomass and vegetation structure parameters for improved carbon estimate from lidar remote sensing. J. Geophys. Res. 2010, 115. [Google Scholar]
  12. Popescu, S.C. Estimating biomass of individual pine trees using airborne lidar. Biomass. Bioenerg. 2007, 31, 646–655. [Google Scholar] [CrossRef]
  13. Zhao, K.; Popescu, S.; Nelson, R. Lidar remote sensing of forest biomass: A scale-invariant estimation approach using airborne lasers. Remote Sens Environ. 2009, 113, 182–196. [Google Scholar] [CrossRef]
  14. Popescu, S.C.; Wynne, R.H.; Nelson, R.F. Estimating plot-level tree heights with lidar: Local filtering with a canopy-height based variable window size. Comput. Electron. Agrc. 2002, 37, 71–95. [Google Scholar] [CrossRef]
  15. Barraud, J. The use of watershed segmentation and GIS software for textural analysis of thin sections. J. Volcanol. Geoth. Res. 2006, 154, 17–33. [Google Scholar] [CrossRef]
  16. Andersen, H.-E.; McGaughey, R.J.; Reutebuch, S.E. Estimating forest canopy fuel parameters using LIDAR data. Remote Sens. Environ. 2005, 94, 441–449. [Google Scholar] [CrossRef]
  17. Scott, J.H.R.; Elizabeth, D. Assessing Crown Fire Potential by Linking Models of Surface and Crown Fire Behavior. Res. Pap. RMRS-RP-29; U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station: Fort Collins, CO, USA, 2001; p. 59. [Google Scholar]
  18. Pang, Y.; Li, Z.-Y. Inversion of biomass components of the temperate forest using airborne Lidar technology in Xiaoxing’an Mountains, Northeastern of China. Chin. J. Plant Ecol. 2013, 36, 1095–1105. [Google Scholar] [CrossRef]
  19. He, Q.; Chen, E.; An, R.; Li, Y. Above-Ground Biomass and Biomass Components Estimation Using LiDAR Data in a Coniferous Forest. Forests 2013, 4, 984–1002. [Google Scholar] [CrossRef]
  20. Tian, X.; Su, Z.; Chen, E.; Li, Z.; van der Tol, C.; Guo, J.; He, Q. Estimation of forest above-ground biomass using multi-parameter remote sensing data over a cold and arid area. Int. J. Appl. Earth Observ. Geoinf. 2012, 14, 160–168. [Google Scholar] [CrossRef]
  21. Li, X.; Li, X.; Li, Z.; Ma, M.; Wang, J.; Xiao, Q.; Liu, Q.; Che, T.; Chen, E.; Yan, G.; et al. Watershed Allied Telemetry Experimental Research. J. Geophys. Res. 2009, 114. [Google Scholar] [CrossRef]
  22. Cao, C.; Bao, Y.; Xu, M.; Chen, W.; Zhang, H.; He, Q.; Li, Z.; Guo, H.; Li, J.; Li, X.; et al. Retrieval of forest canopy attributes based on a geometric-optical model using airborne LiDAR and optical remote-sensing data. Int. J. Remote Sens. 2012, 33, 692–709. [Google Scholar] [CrossRef]
  23. Wang, J.Y.; Ju, K.J.; Fu, H.E.; Chang, X.X.; He, H.Y. Study on biomass of water conservation forest on North Slope of Qilian Mountains. J. Fujian Coll. For. 1998, 18, 319–323. [Google Scholar]
  24. Qin, Y.; Li, B.; Niu, Z.; Huang, W.; Wang, C. Stepwise decomposition and relative radiometric normalization for small footprint LiDAR waveform. China Ser. D 2010, 54, 625–630. [Google Scholar] [CrossRef]
  25. Vincent, L.; Soille, P. Watersheds on digital spaces-An efficient algorithm based on immersion simulations. IEEE Trans. Patt. Anal. 1991, 13, 583–598. [Google Scholar] [CrossRef]
  26. Liu, Q.; Li, Z.; Chen, E.; Pang, Y.; Wu, H. Extracting height and crown of individual tree using airborne LIDAR data. J. Beijing For. Univ. 2008, 30, 83–89. [Google Scholar]
  27. Edson, C.; Wing, M.G. Airborne Light Detection and Ranging (LiDAR) for Individual Tree Stem Location, Height, and Biomass Measurements. Remote Sens. 2011, 3, 2494–2528. [Google Scholar] [CrossRef] [Green Version]
  28. Li, W.; Niu, Z.; Gao, S.; Qin, Y.C. Analyzing and retrieving structure information of Picea crassifolia based on airborne light detection and ranging data. J. Remote Sens. 2013, 6, 1612–1626. [Google Scholar]

Share and Cite

MDPI and ACS Style

Li, W.; Niu, Z.; Gao, S.; Huang, N.; Chen, H. Correlating the Horizontal and Vertical Distribution of LiDAR Point Clouds with Components of Biomass in a Picea crassifolia Forest. Forests 2014, 5, 1910-1930. https://doi.org/10.3390/f5081910

AMA Style

Li W, Niu Z, Gao S, Huang N, Chen H. Correlating the Horizontal and Vertical Distribution of LiDAR Point Clouds with Components of Biomass in a Picea crassifolia Forest. Forests. 2014; 5(8):1910-1930. https://doi.org/10.3390/f5081910

Chicago/Turabian Style

Li, Wang, Zheng Niu, Shuai Gao, Ni Huang, and Hanyue Chen. 2014. "Correlating the Horizontal and Vertical Distribution of LiDAR Point Clouds with Components of Biomass in a Picea crassifolia Forest" Forests 5, no. 8: 1910-1930. https://doi.org/10.3390/f5081910

APA Style

Li, W., Niu, Z., Gao, S., Huang, N., & Chen, H. (2014). Correlating the Horizontal and Vertical Distribution of LiDAR Point Clouds with Components of Biomass in a Picea crassifolia Forest. Forests, 5(8), 1910-1930. https://doi.org/10.3390/f5081910

Article Metrics

Back to TopTop