Next Article in Journal
Feature Selection Solution with High Dimensionality and Low-Sample Size for Land Cover Classification in Object-Based Image Analysis
Previous Article in Journal
Integration of Information Theory, K-Means Cluster Analysis and the Logistic Regression Model for Landslide Susceptibility Mapping in the Three Gorges Area, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimating Forest Structural Parameters Using Canopy Metrics Derived from Airborne LiDAR Data in Subtropical Forests

Co-Innovation Center for Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing 210037, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2017, 9(9), 940; https://doi.org/10.3390/rs9090940
Submission received: 14 June 2017 / Revised: 6 September 2017 / Accepted: 6 September 2017 / Published: 11 September 2017
(This article belongs to the Section Forest Remote Sensing)

Abstract

:
Accurate and timely estimation of forest structural parameters plays a key role in the management of forest resources, as well as studies on the carbon cycle and biodiversity. Light Detection and Ranging (LiDAR) is a promising active remote sensing technology capable of providing highly accurate three dimensional and wall-to-wall forest structural characteristics. In this study, we evaluated the utility of standard metrics and canopy metrics derived from airborne LiDAR data for estimating plot-level forest structural parameters individually and in combination, over a subtropical forest in Yushan forest farm, southeastern China. Standard metrics, i.e., height-based and density-based metrics, and canopy metrics extracted from canopy vertical profiles, i.e., canopy volume profile (CVP), canopy height distribution (CHD), and foliage profile (FP), were extracted from LiDAR point clouds. Then the standard metrics and canopy metrics were used for estimating forest structural parameters individually and in combination by multiple regression models, including forest type-specific (coniferous forest, broad-leaved forest, mixed forest) models and general models. Additionally, the synergy of standard metrics and canopy metrics for estimating structural parameters was evaluated using field measured data. Finally, the sensitivity of vertical and horizontal resolution of voxel size for estimating forest structural parameters was assessed. The results showed that, in general, the accuracies of forest type-specific models (Adj-R2 = 0.44–0.88) were relatively higher than general models (Adj-R2 = 0.39–0.77). For forest structural parameters, the estimation accuracies of Lorey’s mean height (Adj-R2 = 0.61–0.88) and aboveground biomass (Adj-R2 = 0.54–0.81) models were the highest, followed by volume (Adj-R2 = 0.42–0.78), DBH (Adj-R2 = 0.48–0.74), basal area (Adj-R2 = 0.41–0.69), whereas stem density (Adj-R2 = 0.39–0.64) models were relatively lower. The combination models (Adj-R2 = 0.45–0.88) had higher performance compared with models developed using standard metrics (only) (Adj-R2 = 0.42–0.84) and canopy metrics (only) (Adj-R2 = 0.39–0.83). The results also demonstrated that the optimal voxel size was 5 × 5 × 0.5 m3 for estimating most of the parameters. This study demonstrated that canopy metrics based on canopy vertical profiles can be effectively used to enhance the estimation accuracies of forest structural parameters in subtropical forests.

Graphical Abstract

1. Introduction

Forested ecosystems are spatially dynamic and continuously changing and therefore comprise complex and heterogeneous forest structures [1,2]. Forest structure, defined as the spatiotemporal arrangement of structural components in specific vertical and horizontal spatial patterns within a forest stand [3,4,5], is recognized as both a product and driver of forest biophysical processes [6] and represents important forest information, which is useful for guiding multi-functional forest management [7]. Forest structural parameters (e.g., tree height, volume, biomass or stem density etc.) provide considerable information on the spatial and temporal distribution of forests as well as structural properties, and are considered critical components of forest inventory [8] and reliable diversity indicators across forest successional stages [3]. So obtaining spatially continuous estimates of forest structural parameters is valuable for supporting long-term sustainable forest management [9].
Subtropical forests are distributed in a transition zone between tropical and temperate zones, i.e., the region lying largely from 23.5° to 40° latitude in the northern or southern hemispheres [10]. Subtropical forests consist of both subtropical humid and subtropical dry forests, which have unique ecological characteristics when compared to tropical and temperate regions [11]. Subtropical forests, which account for approximately 9% of the world’s forest area [12], are considered a carbon sink contributing to global forest carbon sequestration, and have high species richness, complex structure of forest, high biodiversity and high net ecosystem productivity (NEP) [11]. Quantitative measurements of forest structural parameters of subtropical forests are required to understand forest ecological mechanisms, promote regional ecological developments, maintain biodiversity and enhance regional carbon balance [13].
Traditionally, forest structural parameters are assessed by conventional field inventories, which is time-consuming, costly and limited in spatial extent [2,14]. As a promising earth observation technique, remote sensing has shown great potential for providing multi-scale, multi-dimensional and multi-temporal earth surface information [15] for instantaneous, quantificational and accurate measurements of spatially continuous wall-to-wall properties of forest structure over large-scale areas in lieu of time-consuming and labor-intensive inventory [16]. Furthermore, integrating information from remotely sensed data with a high level of precision and temporal consistency has been recognized as having the ability to describe forest biophysical properties and effectively enhance the performance of forest structure estimations [17].
Estimates of forest stand structural parameters have been derived from optical remote sensing data for several decades [17,18]. However, passive remote signals are generally reflected or absorbed in the uppermost canopy layers and tend to “saturate”, especially in dense forest(i.e., high canopy closure), limiting the ability to characterize vertical structure [19,20]. Similarly, Radar (Radio Detection and Ranging) technology also reveals the aforementioned data saturation problems, due to noise introduced by terrain, surface moisture and other factors [20,21]. Conversely, as a promising active remote sensing technology, Light Detection and Ranging (LiDAR) can be used to directly estimate a spatially explicit three-dimensional (3D) canopy structure with submeter accuracy by transmitting short laser pulses and receiving returned signals [22,23]. Furthermore, LiDAR systems have the ability to overcome the data saturation problems in optical or Radar remote sensing, as a laser beam can strongly penetrate through even dense and multilayered forest canopies to the earth’s surface [24].
Means et al. (2000) [25] estimated forest structural parameters, i.e., tree height, basal area, and volume, using airborne LiDAR data over a Douglas-fir-dominated temperate forest in the Western Cascades of Oregon. They found that the estimation of tree height predicted by the metrics of height percentiles and resulted in R2 values of 0.93–0.98. The R2 values were 0.94–0.95 and 0.95–0.97 for basal area and volume, which were predicted using the metrics of height percentiles and canopy densities as independent variables. Silva et al. (2016) [26] predicted and mapped volume using LiDAR metrics in Eucalyptus plantations in tropical forests (located in São Paulo, Brazil), and found that volume (Adj-R2 = 0.84) was well predicted by the coefficient of variation of return height and the 99th height percentile from LiDAR. Tesfamichael and Beech (2016) [27] used height-related metrics (e.g., height percentiles, maximum height) and canopy density metrics to estimate plot-level structural attributes (i.e., mean height, maximum height, crown diameter and aboveground biomass) over a savanna ecosystem region located in the south western part of Zambia, and resulted in R2 values of 0.48–0.95. However, these studies often include height and density predictors with little physical justification for model formulation. Moreover, they usually neglected a mechanism to summarize complex canopy characteristics into simple parameters, which can potentially be used for estimates of forest structural parameters in different forest conditions [14,28], and the standard metrics (i.e., height-based and density-based metrics) tend to be strongly inter-correlated, and depend on forest conditions, plot sizes, point cloud density and geometrical distributions of point clouds etc. [29,30,31,32,33], and a large subset of these metrics are linked to only a few forest stand characteristics. Thus, these metrics generally have a relatively low transferability and are limited in describing the vertical heterogeneity of forest structure [24].
Canopy vertical profiles, defined as the distributed curves about the characteristic of forest structural components as a function of height above ground, are intimately linked to the vertical distributions of forest structural elements (e.g., foliage, branches, trunks etc.), and have strong potential in enhancing the theoretical explanations of vertical forest structure [4,34,35]. Canopy vertical profiles are important descriptors of forest structure. Lefsky et al. (1999) [36] developed an approach termed “canopy volume models (CVM) ” to characterize three-dimensional (3D) volumetric structure of forest canopies by quantifying the differences in the total volume and spatial organization of the tree foliage, and this approach is beneficial to distinguish specific volumetric canopy architecture. They found relatively high accuracies (R2 = 0.52–0.91) for estimating forest structural parameters using the metrics derived from canopy vertical profiles (i.e., canopy volume profiles (CVP) and canopy height profiles (CHP)). Lovell et al. (2003) [37] used airborne and terrestrial LiDAR data to derive foliage profiles (FP) and estimated effective leaf area index (LAI) in temperate forests located in southern Australia. They found that results compared with LAI derived from classified hemispherical photographs with agreement within 8%. Coops et al. (2007) [38] refined the CVM approach to adapt discrete return LiDAR data. In addition, a Weibull fitting approach was conducted to fit FP profiles and further obtain relevant LiDAR metrics, and finally a number of forest structural parameters (i.e., mean height, basal area) (R2 = 0.65–0.85) were estimated. Hilker et al. (2010) [39] assessed and compared canopy metrics derived from canopy vertical profiles using airborne and terrestrial LiDAR data. The results showed that airborne and terrestrial LiDAR were both able to accurately determine canopy height (absolute error of height was less than 2.5 m) and LAI (R2 = 0.86–0.90). However, most previous studies that estimate forest structural parameters using canopy metrics derived from canopy vertical profiles were conducted in temperate, tropical and boreal forests, and published studies of the subtropical forests are few.
In this study, the standard metrics and canopy metrics derived from airborne LiDAR data are used to estimate plot-level forest structural parameters (i.e., mean diameter at breast height, Lorey’s mean height, stem density, basal area, volume, and aboveground biomass) individually and in combination over a north subtropical secondary forest in southern Jiangsu Province, China. The objectives of this study are: (1) to derive two suites of canopy metrics, i.e., canopy volume (CV) metrics and Weibull-fitted (WF) metrics, using voxel-based CVM and Weibull fitting approaches separately; and (2) to assess the capability of standard metrics and canopy metrics based models and combination models for estimating forest structural parameters and to evaluate the accuracies of the models; and (3) to explore the optimal horizontal and vertical resolution of voxels for the predictive models.

2. Materials and Methods

An overview of the workflow for calculating plot-level forest structural parameters estimation is shown in Figure 1. Firstly, a 1-m Digital terrain model (DTM) was created by the last return points. The data was filtered to remove the above-ground returns, and then the DTM was created by calculating the average elevation from the remaining (ground) LiDAR returns within a cell. Cells that contain no points are interpolated by neighboring cells. The point cloud was then normalized against the ground surface height. Secondly, standard metrics (i.e., height-based (HB) metrics and density-based (DB) metrics) and canopy metrics (i.e., canopy volume (CV) metrics and Weibull-fitted (WF) metrics) were extracted. The suite of canopy volume (CV) metrics was derived by a voxel-based CVM approach, and another suite of WF-metrics was derived by calculating α (scale) and β (shape) parameters of Weibull function fitting to canopy height distribution (CHD) and FP. Thirdly, the estimation capability of standard metrics-based (SM) models, canopy metrics-based (CM) models and combination models were examined for estimating forest structural parameters separately. Finally, the accuracies of the models were assessed and validated by field measured data.

2.1. Study Area

This study was conducted in Yushan Forest, a state-operated forest and national park located near the town of Changshu in Jiangsu Province, southeastern China (120°42′9.4″E, 31°40′4.1″N). The total site area is approximately 1260 ha, which covers approximately 1140 ha of forests. Topographically, the site’s mountain terrain extends from northwest to southeast and the ridge line is more than 6500 m, with the elevation range between approximately 20 and 261 m above sea level. This site is situated in the north-subtropical monsoon climatic region with an annual mean temperature of 15.4 °C, and precipitation of 1047.7 mm, and annual mean relative humidity of approximately 80%. The highest monthly precipitation occurs from June to September. The soil type in Yushan is composed mainly of mountain yellow-brown earth. The forest in Yushan belongs to the north-subtropical mixed secondary forest with three main forest types: conifer-dominated, broad-leaved dominated and mixed forests. The dominant broad-leaved tree species include Oriental oak (Quercus variabilis Bl.), Chinese sweet gum (Liquidambar formosana Hance) and Sawtooth oak (Quercus acutissima Carruth.) of deciduous broad-leaved trees species, mixed with other evergreen broad-leaved tree species including Camphorwood (Cinnamomum camphora (L.) Presl.) and Chinese holly (Ilex chinensis Sims.). The primary coniferous forests are dominated by evergreen coniferous tree species, including Masson pine (Pinus massoniana Lamb.), Chinese fir (Cunninghamia lanceolata (Lamb.) Hook.), slash pine (Pinus elliottii Engelm.) and Japanese Blackbark Pine (Pinus thunbergii Parl.). Figure 2 shows an overview of the study site and distribution of sample plots and Figure 3 shows the field photos of three forest types.

2.2. Data Acquisition and Pre-Processing

2.2.1. LiDAR Data

Small footprint airborne LiDAR data were acquired on 17 August 2013 using a Riegl LMS-Q680i sensor flown at 900 m above ground level, with a flight speed of 55 m·s−1 and a flight line side-lap of ≥60%. The sensor recorded returned waveforms of laser pulse with a temporal sample spacing of 1 ns (approximately 15 cm). The LiDAR system was configured to emit laser pulses in the near-infrared band (1550 nm) at a 360 kHz pulse repetition frequency and a 112 Hz scanning frequency, with a scanning angle of ±30° from nadir and a swath of 1040 m. The dataset had an average beam footprint size of 0.45 m (nadir) in diameter. The average ground point distances of the dataset were 0.49 m (flying direction) and 0.48 m (scanning direction) in a single strip, with pulse density of approximately 5.06 pulse m−2. The final extracted point clouds and associated waveforms were stored in LAS 1.3 format (American Society for Photogrammetry and Remote Sensing, Bethesda, MD, USA).
In order to obtain the relative height of trees, raw point cloud data were first filtered by removing outliers. The data were filtered to remove non-ground points using an algorithm adapted from Kraus and Pfeifer (1998) [40], which was based on a method of linear least-squares interpolation, and then the data were smoothed by the median filter (moving square windows of size 5 × 5 m). After filtering the non-ground points, a 1-m Digital terrain model (DTM) was created by calculating the average elevation from the ground points within a cell (cells that contain no points were filled by interpolation using neighboring cells). Then, the point cloud was then normalized against the ground surface height and extracted for each plot. Point clouds for all plots (n = 51) were finally extracted using the coordinates of the lower left and upper right corners.

2.2.2. Field Data

The field data for the study site were collected from June to August in 2012 and in August of 2013. Throughout the Yushan study region, a total of 51 square sample plots (30 × 30 m) were established, covering the forest type, dominant species compositions, age classes, and site indices, according to an historical forest resource inventory data (2012). All plots were divided into broad-leaved forest (n = 14), coniferous forest (n = 14), and mixed forest (n = 23). The centers of each plot and plot corners were located using Trimble GeoXH6000 Handheld GPS (Trimble, Sunnyvale, CA, USA) units equipped with a dual frequency GNSS antenna, and corrected with high precision real-time differential signals received from the Jiangsu Continuously Operating Reference Stations (JSCORS), resulting in a submeter positional of accuracy of less than 0.5 m [41]. The plot directions and inclined angles were recorded by forest compass, and the border lengths were measured by PI tape. For all live trees with a diameter at breast height (DBH) over 5 cm, tree type, diameter, height, height to crown base, crown width in both cardinal directions, crown class, and crown transparency were measured. DBH was measured on all trees using a diameter tape. Heights of all trees were measured using a Vertex IV hypsometer (Haglöf, Långsele, Sweden). Crown widths were obtained by measuring the average of two values measured along two perpendicular directions from the location of the tree top. In addition, small trees (DBH < 5 cm) and dead wood were also tallied for total stem density, but not used in biomass calculations.
Several forest structural parameters were assessed in this study, including mean DBH, Lorey’s mean height (i.e., the basal area weighted height), stem density, basal area, volume and aboveground biomass. In addition, aboveground biomass of each tree was calculated by means of the species specific allometric equations from local or nearby province [42,43,44,45,46,47] (Appendix A (Table A1)), and the tree-based calculation results were summed within each plot to determine plot-level forest aboveground biomass. Plot-level volume was similarly calculated using provincial species specific volume equations of individual trees, which were based on DBH as predictor variables. A summary of plot-level forest structural parameters data is presented in Table 1.

2.3. Derived Metrics

2.3.1. Canopy Volume Model Approach

A voxel-based CVM approach was applied for point cloud data to derive metrics in this study. The canopy spaces were first organized as a matrix composed of voxels (5 × 5 × 0.5 m3), and these voxels were classified as either “filled” or “empty” volume depending on the presence or absence of LiDAR points within each voxel. “Filled” voxels were further classified as either “euphotic“ zone, if they were located in the uppermost 65% of all filled voxels, or as “oligophotic” zone if they were located below the point, whereas “empty” voxels were located either below (“closed gap”) or above the canopy (“open gap”) [38]. Open gap, euphotic, oligophotic and closed gap were determined as four canopy structure classes, with units defined as the volume of each class per unit area. All volume elements (Open gap, Oligophotic, Euphotic, Closed gap, Filled, Empty) were derived as canopy volume (CV) metrics using the CVM method and canopy volume profile (CVP) was visualized. Figure 4 shows the illustration of voxel-based CVM approach. Point clouds of a plot (30 × 30 m2) were voxelized, and divided into 36 vertical columns of voxels, and each column was further stratified with four canopy structure classes. All columns of a plot were expanded in a panel and the canopy volume distribution (CVD) was presented (Figure 4c). Finally, the volume percentages of canopy structure classes of each height interval (0.5 m) were calculated, resulting in CVP (Figure 4d).
Notably, an appropriate voxel volume size for CVM in this study was been considered because various voxel sizes likely change the distributions and proportions of canopy structure classes. Thus, this study also investigated the influence of various voxel sizes on the accuracies of the models. Given the average beam footprint size of 0.45 m, average ground point distances of 0.49 m (flying direction) and 0.48 m (scanning direction) and pulse density of approximately 5.06 pulse·m−2, horizontal resolutions of 1 m to 10 m were chosen (which were multiples of the footprint size and average ground point distances). Vertical resolutions of 0.5 m and 1 m were chosen to correspond to roughly three and six sampling intervals of the returned waveform. A sensitivity analysis was performed using CV-metrics (i.e., Open gap, Oligophotic, Euphotic, Closed gap, Filled, Empty).

2.3.2. Weibull Fitting Approach

Canopy height distributions (CHD), which describe vertical distributions of foliage elements and non-photosynthetic tissues within canopy spaces, were used to measure the distribution of laser returns within the 0.3-m bins (i.e., a 30 × 30 × 0.3 m3 rectangular section) from the ground to canopy top [48,49]. In this study, a two-parameter Weibull density function (PDF) was used to describe CHD on each plot. As a Weibull model is highly adaptive, ranging from an inversed J-shape to unimodal skewed and unimodal symmetrical curve, the Weibull model has flexibility in characterizing distributions of a range of forest attributes [50,51]. The two parameters, i.e., Weibull scale (α1) and Weibull shape (β1), were derived by the maximum likelihood estimation method. Weibull scale determines the basic shape of the distribution density curve and Weibull shape controls the breadth of the distribution [52]. Foliage profile (FP) can delineate the vertical distribution of canopy phytoelement (e.g., leaf, stem, twig, etc.) density above the ground within a forest stand [37]. FP is defined as the total one-sided leaf area that is involved in photosynthesis per unit canopy volume at canopy height z, and describes changes in the leaf area distribution with increasing height [53]. FP is highly related to leaf area index (LAI), which was demonstrated in previous studies [35,54], and the relationship between FP and LAI is:
L ( z ) = z 1 z 2 F P ( z ) d z ,
where L(z) is the cumulative leaf area index (LAIc) from the ground to a given height z; FP(z) represents the foliage area volume density at height z (is the vertical foliage profile in a thin layer or “slice” through a canopy as a function of height z); z1 and z2 are different canopy height. A height interval or each vertical “slice” was 0.3 m. Meanwhile, we assumed that foliage elements in a thin “slice” were very small so that occlusion can be neglected, and leaves presented Poisson random distribution. Because airborne LiDAR is incapable of resolving foliage angle distribution, clumping and non-foliage elements, the foliage profiles derived from airborne LiDAR are referred to here as “apparent” foliage profiles and effective LAI [37]. In this study, LAI can be indirectly determined from LiDAR by estimating the derived gap probability in the canopy [37,38], and the gap probability be estimated as the total number of laser hits up to a height z relative to the total number of LiDAR shots as follows:
L ( z ) = l n ( P g a p ( z ) ) = l n ( 1 ( # z j | z j > z ) N ) ,
where Pgap (z) is a gap probability measurement at height z, #z is the number of hits down to a height z above the ground, and N is the total number of shots emitted up to the sky. Previous studies have showed that Weibull distribution function can also delineate vertical foliage profiles distributions [37,55]. In this study, the Weibull fitted scale parameter (α2) and shape parameter (β2) were derived from the apparent FP by linking Weibull cumulative function to cumulative projected foliage area index [37,38]:
L ( z ) = 1 ( e ( 1 z / H max α 2 ) β 2 ) ,
where α2 and β2 are fitted parameters, z is the height, and H is the maximum height in a plot.
Moreover, another suite of standard metrics were calculated, including height-based (HD) metrics (h25, h50, h75, h95, hmean, hcv, hskewness and hkurtosis) and density-based (DB) metrics (d1, d3, d5, d7, d9, CC2m). A summary of these metrics with corresponding descriptions is shown in Table 2.

2.4. Metrics Selection and Statistical Analysis

All of the LiDAR metrics in Table 2 were used to analyze pair-wise relationships among different forest structural parameters (DBH, Lorey’s mean height, stem density, basal area, volume and AGB) by Pearson’s correlations (r). Then the metrics with low correlations (r < 0.2) were excluded and candidate metrics were used in the regression analysis. In the multiple regression analysis, all of the dependent variables and independent variables were transformed using the natural logarithm to improve linearity and corrected for bias using a bias correction factor (BCF) [56]. Some studies have applied log transformations to both dependent variables and independent variables for estimations of forest parameters [57,58]. Multiple regression models including forest type-specific (coniferous forest, broad-leaved forest, and mixed forest) models and general models of all plots were then established. Both stepwise variable selection and the maximum coefficient of determination (R2) improvement variable selection techniques were applied to select the metrics to be included in the models [59]. Independent variables were left in the model using an F-test with a p < 0.05 significance level. The standard least-squares method was used [60].
To ensure that the independent variables were not highly correlated, multicollinearity was evaluated using Principal Component Analysis (PCA) based on the correlation matrix. Models with condition number (k) lower than 30 were accepted to ensure that there was no serious multicollinearity in the selected models [57]. The best fitting models were then selected based on the lowest Akaike information criterion value [61]. The accuracies of predictive models were evaluated using adjusted coefficient of determination (Adj-R2), Root-Mean-Square Error (RMSE), which has been transformed back to original scale, and relative RMSE (rRMSE), which are defined as the percentage of the ratio of RMSE and the observed mean values. In this study, dummy variables (or class variables) were added to the selected models as the dependent variables to assess whether these models differ between forest types [62]. Once the best models were chosen, leave-one-out cross-validation was performed to evaluate the predictive accuracies of the models [63].
Adj - R 2 = 1 n 1 n p 1 ( 1 R 2 )
R M S E = 1 n i = 1 n ( x i x ^ i ) 2
r R M S E = R M S E x ¯ × 100 % ,
where x i is the observed value for plot i, x ¯ is the observed mean value for plot i, x ^ i is the estimated value for plot i, n is the number of plots i, and p is the number of variables.

3. Results

3.1. Profile Analysis

The plots of each forest type were stratified into three groups (low, medium, and high), according to the Lorey’s mean height from low to high. In each group, three plots were selected, and a total of nine typical plots were selected. For the typical plots, CVD, CVP and FP were extracted, as shown in Figure 5, Figure 6 and Figure 7. In addition, Figure 8 shows the mean LAIc for plots in different forest types and mean CVD.
Figure 5 shows the spatial arrangements of four canopy structure classes for coniferous, broad-leaved, and mixed forest plots. Generally, Oligophotic zones were larger than euphotic zone in filled volume; coniferous forests had the largest open gap zone and the smallest closed gap zone, whereas broad-leaved forests plots had a larger and wider spread of closed gap zone than mixed forest. Similarly, the percentage of closed gap volume was larger in broad-leaved forests than in mixed forests, and the lowest percentage of closed gap volume was in coniferous forests (Figure 6). The mean CVPs (Figure 6d,h,l) show that the percentages of open gap volume were the highest in coniferous forests, and the differences were not significant between the percentages of open gap volume in broad-leaved forests and mixed forests. The percentages of filled volume in coniferous and mixed forests were significantly higher than in broad-leaved forest, and the differences for the percentage of filled volume between coniferous and mixed forest were not significant.
Weibull models were fitted to canopy foliage distribution and matched the shape of foliage profile relatively well (Figure 7). In general, the FP profiles first exhibited a strong increasing trend, followed by a decreasing trend. Particularly, the peaks of FP in coniferous and mixed forests occurred in the lower or middle portions of the canopies whereas the peaks of broad-leaved forests were distributed more toward middle and upper portions of the canopies. Comparing with the mean foliage profiles and Weibull curves of three forest types, the curve showing the spatial distribution of FP values was smoother in broad-leaved forests than those of coniferous or mixed forests. The Weibull shapes of mixed forest canopy were slightly steeper than those of coniferous forest stands, indicating a wider spread of foliage within the canopy (Figure 7d,h,l). This same trend can be seen in the mean CHDs (Figure 8b–d).
The mean LAIc values below the threshold of 12 m (approximately middle canopy) were relatively high for mixed forests, followed by coniferous forests and broad-leaved forests (Figure 8a). Above the tree height of 12 m, the increasing slope of the mean LAIc of the broad-leaved forests with increasing tree height was higher than that of coniferous forests, and maintained a relative high increasing trend, whereas the increasing trend of coniferous forests and mixed forests gradually tended to saturate. As a result, the mean LAIc value of broad-leaved forests was eventually higher than that of mixed forests, and lowest for coniferous forests.

3.2. Accuracy Assessments

The selected metrics and accuracy assessment results of all the multi-regression models (i.e., SM models, CM models, and combination models) are shown in Table A1, Table A2 and Table A3 and Table 3 summarizes their accuracies. All of the forest structural parameters were generally well estimated (Adj-R2 = 0.39–0.88, rRMSE = 5.13–29.86%). Overall, Lorey’s mean height (Adj-R2 = 0.61–0.88, rRMSE = 5.13–12.79%) and AGB (Adj-R2 = 0.54–0.81, rRMSE = 12.19–28.42%) was predicted most accurately. For volume, DBH and basal area, the R2 values were slightly lower and ranged from 0.42 to 0.78, 0.48 to 0.74 and 0.41 to 0.69, respectively. The lowest accuracy was found for stem density (Adj-R2 = 0.39–0.64, rRMSE = 18.68–29.86%). In comparison, most of forest structural parameters in type-specific models (Adj-R2 = 0.44–0.88, rRMSE = 5.13–28.42%) had higher accuracies than in general models (Adj-R2 = 0.39–0.77, rRMSE = 8.54–29.86%), indicating that the accuracies of forest type-specific models were generally improved rather than general models. Furthermore, the fitted models of the forest structural parameters were relatively more accurate for coniferous forests (Adj-R2 = 0.54–0.81, rRMSE = 8.59–26.55%) than broad-leaved forests (Adj-R2 = 0.50–0.88, rRMSE = 6.39–28.42%) and mixed forests (R2 = 0.44–0.84, rRMSE = 5.13–29.52%). Compared with canopy metrics based models (Adj-R2 = 0.39–0.83, rRMSE = 6.94–29.26%), standard metrics based models had a relatively higher performance (Adj-R2 = 0.42–0.84, rRMSE = 5.60–29.86%) and the combination models performed best (Adj-R2 = 0.45–0.88, rRMSE = 5.13–28.96%), indicating the inclusion of canopy metrics potentially improved the estimation performances of structural parameters.
For all of the general SM models, the standard metrics that were regressed against for fitting models included most of the standard metrics, indicating those had a relatively strong correlation with forest structural parameters. Overall, h95 (selected by four out of six models), d7 (selected by four out of six models), d3, hcv and d9 (each of them was selected by three out of six models) were the most frequently selected, indicating these metrics are more sensitive and representative to the forest structural parameters. For general CM models, all of CV metrics and WF metrics were selected for estimating forest structural parameters. Within CV metrics, the statistic of Oligophotic (all selected by six models), Empty (selected by four out of six models) and Open (selected by four out of six models) were sensitive to forest structural parameters and these metrics were selected both in the general models and forest type-specific models, suggesting that the three metrics have a strong ability to explain variations. Within WF metrics, α1 was relatively sensitive to structural parameters (selected two out of six models). In six general combination models, most of standard metrics (nine out of 14) and canopy metrics (four out of total 10) were used in combination for parameter estimations. The metrics of Oligophotic, Empty, h95 remained sensitive to structural parameters (selected by 2–4 out of six general combo models). Moreover, h75, d1 and β1 (selected by 2–3 out of 6) became more sensitive for DBH, Lorey’s mean height, and stem density in combination models than SM models.
Figure 9 shows the LiDAR estimated versus the field measured forest structural parameters as well as the results for cross-validation in all plots models based on standard metrics and canopy metrics. As indicated, Lorey’s mean height and AGB models were fitted best and resulted in R2 values of 0.79 and 0.66, followed by DBH (R2 = 0.60), volume (R2 = 0.60) and basal area (R2 = 0.52), whereas the accuracy of stem density model was the lowest (R2 = 0.49). For Lorey’s mean height, AGB, DBH, and volume estimations, their relationships were close to the 1:1 line whereas basal area and stem density had a relationship that deviated from the 1:1 line, with a slightly larger deviation in broad-leaved forests.

3.3. The Selection of Voxel Sizes

In this study, a sensitivity analysis was performed using different voxel sizes to derive CV metrics based on CVM approach and to quantify their influence on the results. As shown in Figure 10, a quantitative comparison of estimation accuracy for four main forest structural parameters (i.e., DBH, Lorey’s mean height, stem density, and basal area) was performed. In general, the R2 values of the models showed a trend of first increasing and then decreasing when horizontal resolutions of voxels were varied from 1 m to 10 m (Figure 10a,b), and the voxels in horizontal resolution of 5 m had the best performance. Figure 10a was subtracted from Figure 10b to calculate the result of Figure 10c, which demonstrated the difference of rRMSE values of forest structural parameters for various vertical resolutions (0.5 m and 1 m). The values presented were mostly positive, except for some of the differences were negative (e.g., G at 3 m horizontal resolution) (Figure 10c). In particular, DBH and stem density models had all positive values across 1 m to 10 m of horizontal resolutions, indicating the two parameters were strongly influenced by the vertical resolution of the voxels. As a result, the suitable voxel size in this study was 5 × 5 × 0.5 m3.

4. Discussion

4.1. Canopy Vertical Profiles

Canopy is an important constituent of forest structure [64], and canopy structure is critical for estimation of forest structural parameters [65]. Canopy vertical profile is one of the means to quantify and analyze complex forest canopy structure and further characterize the potential heterogeneity of forest spatial structure [66]. A wide range of forest structural parameters can be directly quantified from canopy vertical profiles such as canopy height and canopy vertical distribution [67]. Also, a set of forest structural parameters (aboveground biomass, basal area, volume, LAI, canopy cover, etc.) can be predicted by establishing empirical models from LiDAR data [68]. In this study, a voxel-based CVM and Weibull fitting approach were conducted to extract two key suites of metrics for estimating forest structural parameters and derive correlative canopy vertical profiles including CVD, CVP, CHD, FP, and LAIc. As mentioned above, the CVM approach provides a broad classification approach to categorize the canopy into photosynthetically active and less active zones [39]. Therefore, it can better reflect the spatial heterogeneity of forest structure, which is caused by the difference of light environment in the canopy. Furthermore, the CVP explicitly presented variation in the spatial arrangement of elements (i.e., open gap, euphotic, oligophotic, closed gap) within the vertical forest canopy [38]. As shown in Figure 5 and Figure 6, the broad-leaved forests had the largest closed gap volume and the smallest open gap volume when compared to coniferous forests and mixed forests. The explanations of these phenomena need to take into account the canopy geometry and tree architecture [36]. At our research site, coniferous forests are dominated by Masson pine and slash pine; these species usually consist of a regular and conical crown, demonstrating a heavily thinned upper canopy and a dense sub-canopy (Figure 3). Furthermore, more open upper canopies in coniferous stands allow more light to pass through to the lower canopy strata [69,70], so a shrubby understory may incrementally emerge, resulting in the most open gap and the lowest closed gap zones in coniferous forests. Conversely and notably, broadleaves with elliptical or spherical crown are very tall and have positively skewed canopies with a lower canopy transparency in this study area, as indicated by the large decrease in open gap zones. Additionally, the closed canopy volume generally increased with decreasing stand density [55], hence the broad-leaved forests with a lower stem density (1126.00 ha−1) also had a more closed canopy gap. Although with a much more shrubby understory, mixed conifer–broadleaf forests generally encompass median height broadleaved trees [65] with a high stem density (1431.78 ha−1) and canopy transparency, resulting in a higher amount of closed gap volume than coniferous forests and a slightly higher amount of open gap volume. On the other hand, as Yushan forest is in secondary succession, the forest canopy surface became more uneven, and the competitions among shade-intolerant species (e.g., Masson pine, Chinese sweet gum) were accelerated and further inhibited the establishment and growth of these species [71,72]. As a result, in late-successional stage, the shade-tolerant species (e.g., Oriental oak, camphorwood and Chinese holly) eventually dominated the canopy [69,72,73] and coexisted with other species. This process could cause the transmittance of light through the canopy to decline [74], which may result in an increase the spatial heterogeneity of the light environment [75,76] and a further enhancement of more microsite light availability in lower canopies [70,76,77,78,79]. Thus, for each forest type, the oligophotic zone, which represented a larger proportion of the total filled volume compared to the euphotic zone that represented photosynthetically active tissues (Figure 6). As mentioned above, the canopy architectures of the three forest types can help explain why the distributions of FP and CHD in coniferous forests and mixed forests inclined to the under canopy, whereas the curves of broad-leaved forests were distributed more towards the middle or upper canopy (Figure 7 and Figure 8b–d).
In general, due to a thinner upper canopy and dense under canopy for each forest type, the mean LAI increased rapidly and shifted to an infinitesimal increment from the ground up to the top of the canopy (Figure 8a). Below the threshold of 12 m (approximately middle canopy), dense foliage accumulated in the lower canopy of mixed forests and coniferous forests but mixed forests had more understory shrubs and slightly denser canopies than coniferous forests whereas broad-leaved forest had less shrubbery; therefore, there was a dramatically increased LAIc in mixed forests, followed by coniferous forests and broad-leaved forests. Along with still moderate density of foliage near the upper canopy in broad-leaved forests as well as thinned density of foliage in mixed forests and coniferous forests, the mean LAIc increased trend remained relatively stable in broad-leaved forests compared to other forest types. Eventually, broad-leaved forests had the highest mean LAIc, followed by mixed forests and coniferous forests, which is consistent with the findings of previous studies [80,81].

4.2. Predictive Models

In comparison, the forest type-specific models had higher accuracies (Adj-R2 = 0.44–0.88, rRMSE = 5.13–28.42%) than the general models (Adj-R2 = 0.39–0.77, rRMSE = 8.54–29.86%). Bouvier et al. (2015) [14] developed a separate model for coniferous, deciduous and mixed stands to estimate forest structural parameters in the Lorraine forests. The results demonstrated that the separate models reduce estimation errors (2.0–5.3%) compared to general models in some complex forests conditions, which was confirmed by our research results. Fu et al. (2011) [82] reported R2 values for AGB of 0.37 of the general model and 0.43–0.68 of forest type-specific models in subtropical forests (located in southern Yunnan province, China). In our study site, the multi-layered forest conditions in subtropical forests contained greater species diversity, making the effects of tree-species composition (classified as forest types) significant. Overall, the models of the forest structural parameters were relatively more accurate for coniferous forests than broad-leaved forests and mixed forests. The relationships between stand structure and the forest structural parameters are species-dependent, and coniferous forests are usually characterized by relatively simple stand structures when compared with broad-leaved or mixed stands. So it is likely that the model prediction accuracy may decrease in multispecies stands [14]. Xu et al. (2015) [83] estimated forest structural parameters (i.e., Lorey’s mean height, stem density, basal area and volume) in the subtropical deciduous mixed forests (on Purple Mountain, located in eastern Nanjing), using canopy height metrics (i.e., height percentile, mean height, maximum height and minimal height) and canopy density metrics. Compared with our results (rRMSE = 5.13–22.28%), theirs showed a relatively lower rRMSE for Lorey’s mean height (6.47%), stem density (27.04%), basal area (16.38%), and volume (6.93%). Compared with canopy metrics-based models (Adj-R2 = 0.39–0.83, rRMSE = 6.94–29.26%), standard metrics-based models had relatively higher performance (Adj-R2 = 0.42–0.84, rRMSE = 5.60–29.86%), except for the volume and AGB (both in forest type-specific models). The combination models performed best (Adj-R2 = 0.45–0.88, rRMSE = 5.13–28.96%), explaining a large amount of the variability for all forest structural parameters and indicating the increased utility of canopy metrics in capturing spatially explicit information describing a heterogeneous forest structure. For DBH, stem density, basal area, and AGB, Lefsky et al. (1999) [36] reported adjusted R2 values of 0.61, 0.52, 0.87, and 0.91 in boreal forests, markedly higher than reported in this study. The cause of the lower performance in this research is likely the complex structure of the subtropical forests, which are typified by multi-layered forests which that encompass some stands with considerable variability in tree height and stem density, especially in old-growth stands, whereas boreal forests have a much higher homogeneous composition and more discernible canopy architecture.
When considering the selected frequency of metrics in the fitted models, for CV metrics, Oligophotic, Empty metrics were mostly selected by the combination models. This may be explained by the higher proportions of canopy elements (Figure 5 and Figure 6), which were oligophotic, empty volume zones, revealing the strong sensitivity and representativeness of these two metrics to local forest structures. Previous studies [38,84] found that the Weibull scale and shape parameter were related to canopy attributes (e.g., crown depth and crown length), hence two Weibull parameters were both selected by the CM models of structural parameters (e.g., mean diameter, Lorey’s mean height, basal area, and volume models) linked to canopy attributes. In combination models, the most selected WF metrics were α1 and β1, indicating that both of them are suitable for estimating structure parameters in local forests. The capacity of the Weibull parameters to represent the key attributes of mean crown dimension is important, as it provides a mechanism to summarize complex canopy characteristics into simple parameters that can be empirically analyzed in relation to various forest stand characteristics. The two-parameter Weibull model was applied for characterizing many types of FPs and CHDs in this study. In general, these profiles of single layer canopies corresponded well (Figure 7). However, the unimodal Weibull distribution function applied to the profile is inadequate to describe properly multimodal structure, which may occur in multi-layered, multi-age, complex forest stands [52]. Thus, the relatively poor fit for multi-layered forests could result in errors in estimates of structural parameters, which may explain why the Weibull parameters are not statistically more significant predictors than CV metrics. In this regard, future work could focus on how to apply a multi curve fitting approach in order to further capture the full distributions of canopy vertical profiles. On the other hand, different plot selection strategies could influence the performance of predictive models [85]. The plot selection in this study was only according to forest type, thus our future work could also examine different plot selection strategies of field training plots (e.g., using LiDAR data or geographical factors as a prior information, etc.) and utilize a suitable strategy to improve the estimation accuracy of forest structural parameters.

4.3. The Selection of Voxel Sizes

Voxels representing canopy elements such as trunks and branches were abstracted by a volume grid and placed in a 3D grid [86]. As a method of volume visualization of LiDAR points, voxels have already been applied to airborne LiDAR data for improving calculations of forest attributes [87,88]. Voxel size is a key parameter pertaining to the scale of forest structural parameter estimates to the physical dimension of canopy components [89]. Thus, a sensitivity analysis was conducted to investigate the influence of various voxel sizes on forest structural estimations. As shown in Figure 10, in very low (i.e., 1 × 1 × 0.5 m3) or high (i.e., 10 × 10 × 1 m3) resolution conditions, the R2 values showed a relatively lower performance. If voxels are too small, a voxel-based CVM approach may produce redundant unfilled voxels of empty volume containing few tree canopy elements, which may lead to the underestimation of forest structural parameters; however, too large voxels may lead to too few voxels and result in statistically insignificant descriptions of canopy features [90]. In these conditions, the voxel approach could become ineffective at characterizing the vertical distribution of various canopy structures and the capability to capture 3D heterogeneity of canopy structure for CV metrics could be constrained, hence resulting in relatively lower performances of the models. After taking into account factors of plot size (30 × 30 m2), point cloud densities (3.74 pts·m−2), etc., Hilker et al. (2010) [39] used a voxel size of 6 × 6 × 1 m3 for discrete airborne LiDAR data to estimate the tree height and LAI in Douglas-fir-dominated forest stands with relatively high tree heights (30–35 m). Concerning a much higher point cloud density (5.06 pts·m−2) of LiDAR data and relatively lower tree heights (4.47–18.52 m) in this study site, a 1 m vertical resolution produced more coarse data than the vertical resolution of 0.5 m (approximately treble the temporal sample spacing of 1 ns (15 cm)), thus, constraining the ability of canopy volume metrics to describe the vertical variability of the forest canopy structures. Moreover, potential tree movement due to wind between laser acquisitions is also considered a source of uncertainty, as laser returns from the same target can be located in different voxels for different laser acquisitions. By using a voxel size larger than the pulse diameter, this issue can be slightly reduced [91]. Overall, the optimal voxel size is a key parameter to determine in order to improve characterizations of forest structure [92,93]. Consequently, the optimal voxel spatial resolution should be determined based on plot size, the characteristics of the LiDAR instrument used (e.g., beam diameter, footprint size, average point density and temporal sample spacing, etc.), and forest structure attributes (e.g., tree height, crown diameter, crown depth, etc.)

5. Conclusions

In this study, a set of canopy metrics derived from canopy vertical profiles, which has the potential to aid in our understanding of the physical characteristics of forest structure, was extracted. The capability of the standard metrics (extracted from the point cloud data) and canopy metrics for estimating forest structural parameters (i.e., DBH, Lorey’s mean height, stem density, basal area, volume, and AGB) was assessed, individually and in combination, over a subtropical forest in southeastern China. Moreover, a sensitive analysis of different voxel sizes was performed to investigate the optimal voxel size for estimating forest structural parameters.
The results demonstrated that the forest type-specific models had relatively higher accuracies (Adj-R2 = 0.44–0.88, rRMSE = 5.13–28.42%) compared with the general models (Adj-R2 = 0.39–0.77, rRMSE = 8.54–29.86%). The estimation accuracies of Lorey’s mean height and AGB were the highest, followed by volume, DBH and basal area, whereas stem density was relatively lower. Overall, metrics of Oligophotic, Empty, Open, α1 were the most frequently selected, indicating their potential capability for predicting forest structural parameters in the forest stands within the study site. The results demonstrated the synergistic use of standard metrics and canopy metrics for better predicting forest structural parameters (∆Adj-R2 = 0.01–0.20, ∆rRMSE = −5.71–1.39%), compared with models developed using standard metrics (only) and canopy metrics (only). In addition, the optimal voxel size for estimating forest structural parameters in this study is 5 × 5 × 0.5 m3, and the voxel vertical and horizontal resolutions should be determined based on plot size, the characteristics of the acquired LiDAR data (i.e., beam diameter, footprint size, average point density, and temporal sample spacing) and forest structure attributes (i.e., tree height, crown diameter, and crown depth).

Acknowledgments

The project was funded by the Natural Science Foundation of Jiangsu Province (No. BK20151515) and the National Natural Science Foundation of China (No. 31400492). This research was also supported by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD). Special thanks to Xin Shen, Kun Liu, and Ting Xu for field works. The authors gratefully acknowledge the foresters in Yushan forest for their assistance with data collection and sharing their rich knowledge and working experience of the local forest ecosystems.

Author Contributions

Zhengnan Zhang analyzed the data and wrote the paper. Lin Cao helped in project and study design, paper writing, and analysis. Guanghui She helped with data analysis and paper writing.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Allometric equations for aboveground biomass components of dominant tree species and species group in the study site.
Table A1. Allometric equations for aboveground biomass components of dominant tree species and species group in the study site.
Tree SpeciesComponentabR2References
Masson pineStem wood (Ws)0.1411.0920.9970Jiang et al. (1992) [42]
Live branches (Wb)0.0650.9910.9871
Foliage (Wf)0.1320.7450.9827
Chinese firStem wood (Ws)0.1240.6800.9704Ye and Jiang (1983) [43]
Live branches (Wb)0.2030.3850.7223
Foliage (Wf)0.8500.1890.6567
Slash pineStem wood (Ws)0.2350.9000.9523Wang and Shi (1990) [44]
Live branches (Wb)0.0801.0640.8520
Foliage (Wf)0.4560.6100.8802
Sawtooth oakStem wood (Ws)0.0181.0340.9864Xu et al. (2011) [46]
Live branches (Wb)0.000081.4680.9745
Foliage (Wf)0.0040.7690.8662
Sweet gumStem wood (Ws)0.0930.8010.9310Qian (2000) [45]
Live branches (Wb)0.0830.6490.9890
Foliage (Wf)1.0840.2170.6940
Other broadleaves aStem wood (Ws)0.0230.9850.9903Sun et al. (1992) [47]
Live branches (Wb)0.000043.7850.9623
Foliage (Wf)0.000031.3780.9456
Notes: The equation of W = a(D2H)b was used to calculate each biomass component. H = Tree height (m), D = DBH (cm) and a, b are the parameters. a The general equation of “Other broadleaves” includes tree species of Quercus variabilis, Quercus fabri, Quercus aliena, Quercus glandurifera var. brevipetiolata, Castanea sequinii, Liquidambar formasana and Pistacia chinensis.
Table A2. Predictive models and accuracy assessment results (by standard metrics).
Table A2. Predictive models and accuracy assessment results (by standard metrics).
VariablesPredictive ModelsAdj-R2RMSErRMSE %
All plots
DBH/cm exp ( 1.064 + 0.641 ln h 95 0.580 ln d 1 + 0.066 ln d 7 ) × 1.008 0.60 ***1.7212.33
hLorey/m exp ( 0.33 0.079 ln h cv + 1.012 ln h 95 0.028 ln d 7 0.001 ln d 9 ) × 1.006 0.75 ***0.979.15
N/(ha−1) exp ( 6.814 0.049 ln h cv + 2.124 ln d 1 0.296 ln d 7 0.049 ln d 9 ) × 1.052 0.42 ***423.7529.86
G/(m2·ha−1) exp ( 0.899 + 0.851 ln h 95 0.819 ln d 3 0.177 ln d 7 ) × 1.019 0.44 ***3.9917.23
V/(m3·ha−1) exp ( 3.445 0.137 ln h cv + 0.232 ln h 25 + 0.515 ln h 50 + 0.295 ln d 3 ) × 1.023 0.46 ***22.3417.46
AGB/(Mg·ha−1) exp ( 0.169 0.058 ln h kurtosis + 1.817 ln h 95 + 0.627 ln d 3 0.048 ln d 9 ) × 1.037 0.64 ***19.1722.47
Coniferous forests
DBH/cm exp ( 0.996 + 0.664 ln h 95 0.485 ln d 1 + 0.88 ln d 7 ) × 1.006 0.67 **1.209.50
hLorey/m exp ( 1.480 0.488 ln h cv + 1.211 ln h 95 0.102 ln d 7 0.008 ln d 9 ) × 1.013 0.661.0911.47
N/(ha−1) exp ( 6.39 0.851 ln h cv + 1.810 ln d 1 0.477 ln d 7 + 0.140 ln d 9 ) × 1.036 0.60315.7818.68
G/(m2·ha−1) exp ( 0.885 + 1.531 ln h 95 + 1.380 ln d 3 0.376 ln d 7 ) × 1.029 0.62 **4.5319.63
V/(m3·ha−1) exp ( 3.459 + 0.633 ln h cv + 3.732 ln h 25 2.335 ln h 50 + 0.491 ln d 3 ) × 1.029 0.69 **22.4019.22
AGB/(Mg·ha−1) exp ( 2.216 + 0.565 ln h kurtosis + 2.218 ln h 95 + 0.455 ln d 3 0.099 ln d 9 ) × 1.059 0.72 **16.8624.17
Broad-leaved forests
DBH/cm exp ( 0.949 + 0.680 ln h 95 0.654 ln d 1 + 0.028 ln d 7 ) × 1.009 0.61 **1.7011.12
hLorey/m exp ( 0.296 0.149 ln h cv + 0.981 ln h 95 0.082 ln d 7 0.031 ln d 9 ) × 1.004 0.84 ***0.786.91
N/(ha−1) exp ( 7.663 + 0.602 ln h cv + 2.500 ln d 1 0.144 ln d 7 0.071 ln d 9 ) × 1.056 0.60298.9926.55
G/(m2·ha−1) exp ( 3.045 + 0.046 ln h 95 + 0.516 ln d 3 0.001 ln d 7 ) × 1.010 0.542.6211.96
V/(m3·ha−1) exp ( 3.459 + 0.633 ln h cv + 3.732 ln h 25 2.335 ln h 50 + 0.491 n d 3 ) × 1.029 0.5619.4914.68
AGB/(Mg·ha−1) exp ( 1.958 0.060 ln h kurtosis + 1.150 ln h 95 + 0.579 ln d 3 + 0.057 ln d 9 ) × 1.063 0.5726.8028.42
Mixed forests
DBH/cm exp ( 1.184 + 0.360 ln h 95 1.116 ln d 1 + 0.110 ln d 7 ) × 1.008 0.48 **1.6611.94
hLorey/m exp ( 0.830 + 0.222 ln h cv + 0.700 ln h 95 + 0.135 ln d 7 0.067 ln d 9 ) × 1.003 0.81 ***0.605.60
N/(ha−1) exp ( 7.774 + 0.093 ln h cv + 4.386 ln d 1 0.391 ln d 7 + 0.157 ln d 9 ) × 1.031 0.48 **336.7328.52
G/(m2·ha−1) exp ( 1.308 + 0.772 ln h 95 + 1.036 ln d 3 0.046 ln d 7 ) × 1.009 0.45 **3.0812.86
V/(m3·ha−1) exp ( 3.144 + 0.006 ln h cv 2.519 ln h 25 + 3.238 ln h 50 + 2.386 n d 3 ) × 1.009 0.60 ***16.7612.70
AGB/(Mg·ha−1) exp ( 0.177 0.261 ln h kurtosis + 1.917 ln h 95 + 1.380 ln d 3 0.081 ln d 9 ) × 1.016 0.64 ***13.2014.77
Notes: Level of significance: NS = not significant (>0.05); ** <0.01; *** <0.001.
Table A3. Predictive models and accuracy assessment results (by canopy metrics).
Table A3. Predictive models and accuracy assessment results (by canopy metrics).
VariablesPredictive ModelsAdj-R2RMSErRMSE %
All plots
DBH/cm exp ( 2.852 + 0.169 ln O G 0.076 ln O l i g o + 0.081 ln E m p t y 0.576 ln β 2 ) × 1.010 0.50 ***1.8613.31
hLorey/m exp ( 0.235 0.421 ln O l i g o + 1.010 ln F illed 0.272 ln E m p t y 0.338 ln β 2 ) × 1.009 0.61 ***1.1811.13
N/(ha−1) exp ( 5.545 0.508 ln O G + 0.232 ln E u + 0.342 ln O l i g o + 1.332 ln β 1 ) × 1.055 0.39 ***415.1729.26
G/(m2·ha−1) exp ( 2.119 + 0.579 ln O l i g o + 0.579 ln E m p t y + 0.094 ln α 2 ) × 1.020 0.41 ***3.6715.82
V/(m3·ha−1) exp ( 3.063 + 0.249 ln O G + 0.237 ln E u + 0.423 ln O l i g o 0.694 ln α 1 ) × 1.025 0.42 ***22.3617.48
AGB/(Mg·ha−1) exp ( 1.650 0.292 ln O G + 1.021 ln O l i g o 0.071 ln C G + 0.717 ln E m p t y ) × 1.048 0.54 ***19.8423.25
Coniferous forests
DBH/cm exp ( 1.455 0.535 ln O G + 0.388 ln O l i g o + 0.622 ln E m p t y 0.052 ln β 2 ) × 1.009 0.541.4011.09
hLorey/m exp ( 4.620 4.700 ln O l i g o + 5.948 ln F i l l e d + 0.281 ln E m p t y 2.952 ln β 2 ) × 1.014 0.641.2112.79
N/(ha−1) exp ( 5.379 0.310 ln O G 0.726 ln E u + 1.019 ln O l i g o + 1.156 ln β 1 ) × 1.042 0.58431.6525.53
G/(m2·ha−1) exp ( 1.378 + 1.535 ln O l i g o + 0.195 ln E m p t y 2.715 ln α 2 ) × 1.043 0.554.7320.48
V/(m3·ha−1) exp ( 1.462 + 0.559 ln O G + 2.869 ln E u 0.765 ln O l i g o 6.184 ln α 1 ) × 1.028 0.72 **18.3215.72
AGB/(Mg·ha−1) exp ( 3.600 3.719 ln O G + 2.238 ln O l i g o 1.604 ln C G + 5.287 ln E m p t y ) × 1.043 0.74 **18.5126.55
Broad-leaved forests
DBH/cm exp ( 3.562 + 0.454 ln O G 0.053 ln O l i g o 0.099 ln E m p t y 0.958 ln β 2 ) × 1.010 0.511.8111.79
hLorey/m exp ( 1.030 0.149 ln O l i g o + 0.323 ln F illed + 0.250 ln E m p t y 0.613 ln β 2 ) × 1.004 0.83 ***0.887.72
N/(ha−1) exp ( 6.025 0.774 ln O G + 1.103 ln E u 0.015 ln O l i g o + 0.836 ln β 1 ) × 1.078 0.52299.8226.63
G/(m2·ha−1) exp ( 2.514 + 0.321 ln O l i g o 0.040 ln E m p t y 0.032 ln α 2 ) × 1.013 0.502.7412.48
V/(m3·ha−1) exp ( 5.339 0.234 ln O G 0.981 ln E u + 0.588 ln O l i g o + 0.715 ln α 1 ) × 1.013 0.5818.9714.28
AGB/(Mg·ha−1) exp ( 4.009 0.188 ln O G + 0.440 ln O l i g o + 0.615 ln C G 0.517 ln E m p t y ) × 1.055 0.6026.4528.05
Mixed forests
DBH/cm exp ( 3.172 + 0.169 ln O G + 0.060 ln O l i g o + 0.002 ln E m p t y 0.724 ln β 2 ) × 1.009 0.481.7812.79
hLorey/m exp ( 0.481 + 0.318 ln O l i g o + 0.390 ln F illed + 0.210 ln E m p t y 0.109 ln β 2 ) × 1.003 0.75 ***0.756.94
N/(ha−1) exp ( 4.939 0.262 ln O G 0.997 ln E u + 0.965 ln O l i g o + 1.804 ln β 1 ) × 1.034 0.44324.7322.68
G/(m2·ha−1) exp ( 1.909 + 0.626 ln O l i g o + 0.080 ln E m p t y + 0.153 n α 2 ) × 1.009 0.45 ***3.1213.01
V/(m3·ha−1) exp ( 0.871 + 0.626 ln O G + 0.080 ln E u + 0.088 ln O l i g o + 0.744 ln α 1 ) × 1.008 0.65 ***16.3112.36
AGB/(Mg·ha−1) exp ( 1.848 + 0.087 ln O G + 1.195 ln O l i g o + 0.087 ln C G + 0.119 ln E m p t y ) × 1.013 0.71 ***13.0014.55
Notes: Level of significance: NS = not significant (>0.05); ** <0.01; *** <0.001.
Table A4. Predictive models and accuracy assessment results (using both standard metrics and canopy metrics).
Table A4. Predictive models and accuracy assessment results (using both standard metrics and canopy metrics).
VariablesPredictive ModelsAdj-R2RMSErRMSE %
All forests
DBH/cm exp ( 0.802 0.060 ln h cv 0.095 ln h 75 + 0.7401 ln d 1 0.548 ln O l i g o ) × 1.008 0.61 ***1.6711.97
hLorey/m exp ( 0.091 0.053 ln h 50 + 0.154 ln d 1 + 0.971 ln E m p t y 0.315 ln β 1 ) × 1.009 0.77 ***0.908.54
N/(ha−1) exp ( 7.012 + 0.167 ln d 1 + 0.761 ln O l i g o + 0.741 ln α 1 + 1.643 ln β 1 ) × 1.049 0.45 ***410.0228.90
G/(m2·ha−1) exp ( 1.576 0.255 h kurtosis 0.304 ln h 25 + 1.031 ln h 95 + 0.073 ln E m p t y ) × 1.017 0.50 ***3.4714.96
V/(m3·ha−1) exp ( 0.201 0.463 h 75 0.502 ln O l i g o + 0.845 ln E m p t y + 2.343 ln α 1 ) × 1.018 0.58 ***21.0716.47
AGB/(Mg·ha−1) exp ( 0.519 + 0.193 ln h 95 + 1.458 ln d 3 + 1.300 ln C C 2 m 1.348 ln O l i g o ) × 1.036 0.65 ***18.2521.39
Coniferous forests
DBH/cm exp ( 1.041 0.164 ln h cv + 0.062 ln h 75 + 0.8131 ln d 1 0.182 ln O l i g o ) × 1.006 0.74 **1.088.59
hLorey/m exp ( 1.207 0.196 ln h 50 + 0.582 ln d 1 + 1.441 ln E m p t y 0.661 ln β 1 ) × 1.010 0.77**0.9910.43
N/(ha−1) exp ( 1.881 + 0.573 ln d 1 3.703 ln O l i g o + 2.036 ln α 1 + 0.029 ln β 1 ) × 1.033 0.64339.2920.07
G/(m2·ha−1) exp ( 0.242 0.323 h kurtosis 0.060 ln h 25 + 1.668 ln h 95 + 0.046 ln E m p t y ) × 1.029 0.69 **4.2318.32
V/(m3·ha−1) exp ( 1.585 0.537 h 75 0.599 ln O l i g o + 1.392 ln E m p t y + 3.014 ln α 1 ) × 1.024 0.78 **18.2115.63
AGB/(Mg·ha−1) exp ( 0.808 + 0.794 ln h 95 + 1.594 ln d 3 + 3.551 ln C C 2 m 4.560 ln O l i g o ) × 1.040 0.81 **14.5320.83
Broad-leaved forests
DBH/cm exp ( 0.876 0.095 ln h cv 0.325 ln h 75 + 0.5881 ln d 1 0.797 ln O l i g o ) × 1.008 0.68 1.5410.06
hLorey/m exp ( 1.020 + 0.047 ln h 50 0.226 ln d 1 + 0.633 ln E m p t y 0.255 ln β 1 ) × 1.003 0.88 ***0.726.39
N/(ha−1) exp ( 7.281 + 0.074 ln d 1 + 0.8081 ln O l i g o + 0.583 ln α 1 + 2.132 ln β 1 ) × 1.055 0.62273.4924.29
G/(m2·ha−1) exp ( 3.404 0.100 h kurtosis 0.166 ln h 25 + 1.109 ln h 95 0.902 ln E m p t y ) × 1.009 0.632.4911.34
V/( m3·ha−1) exp ( 9.533 + 2.028 h 75 + 1.027 ln O l i g o 1.599 ln E m p t y 3.683 ln α 1 ) × 1.011 0.6716.6512.54
AGB/(Mg·ha−1) exp ( 2.376 0.395 ln h 95 + 1.210 ln d 3 + 3.440 ln C C 2 m 4.200 ln O l i g o ) × 1.051 0.6626.6728.29
Mixed forests
DBH/cm exp ( 0.615 + 0.012 ln h cv 0.331 ln h 75 + 0.5901 ln d 1 1.911 ln O l i g o ) × 1.008 0.55 **1.5811.34
hLorey/m exp ( 0.089 0.136 ln h 50 + 0.295 ln d 1 + 1.043 ln E m p t y 0.969 ln β 1 ) × 1.002 0.84 ***0.555.13
N/(ha−1) exp ( 8.121 + 0.135 ln d 1 + 1.246 ln O l i g o + 0.120 ln α 1 + 4.921 ln β 1 ) × 1.029 0.50 ***319.0522.28
G/(m2·ha−1) exp ( 0.179 0.069 h kurtosis + 0.361 ln h 25 1.140 ln h 95 + 1.927 ln E m p t y ) × 1.009 0.56 **2.7711.56
V/(m3·ha−1) exp ( 0.201 0.463 h 75 0.502 ln O l i g o + 0.845 ln E m p t y + 2.343 ln α 1 ) × 1.018 0.71 ***15.8712.02
AGB/(Mg·ha−1) exp ( 0.519 + 0.630 ln h 95 + 1.113 ln d 3 + 0.073 ln C C 2 m + 1.788 ln O l i g o ) × 1.009 0.79 ***10.8912.19
Notes: Level of significance: NS = not significant (>0.05); ** <0.01; *** <0.001.

References

  1. Pan, Y.; Birdsey, R.A.; Phillips, O.L.; Jackson, R.B. The Structure, Distribution, and Biomass of the World’s Forests. Ann. Rev. Ecol. Evol. Syst. 2013, 44, 593–622. [Google Scholar] [CrossRef]
  2. Hill, S.; Lati, H.; Heurich, M.; Müller, J. Individual-tree-and stand-based development following natural disturbance in a heterogeneously structured forest: A LiDAR-based approach. Ecol. Inform. 2017, 38, 12–25. [Google Scholar] [CrossRef]
  3. Franklin, J.F.; Spies, T.A.; Pelt, R.V.; Carey, A.B.; Thornburgh, D.A.; Rae, D.; Lindenmayer, D.B.; Harmon, M.E.; Keeton, W.S.; Shaw, D.C.; et al. Disturbances and structural development of nsatural forest ecosystems with silvicultural implications, using Douglas-fir forests as an example. For. Ecol. Manag. 2002, 155, 399–423. [Google Scholar] [CrossRef]
  4. Palace, M.W.; Sullivan, F.B.; Ducey, M.J.; Treuhaft, R.N.; Herrick, C.; Shimbo, J.Z.; Mota-E-Silva, J. Estimating forest structure in a tropical forest using field measurements, a synthetic model and discrete return lidar data. Remote Sens. Environ. 2015, 161, 1–11. [Google Scholar] [CrossRef]
  5. McElhinny, C.; Gibbons, P.; Brack, C.; Bauhus, J. Forest and woodland stand structural complexity: Its definition and measurement. For. Ecol. Manag. 2005, 218, 1–24. [Google Scholar] [CrossRef]
  6. Spies, T.A. Forest Structure: A Key to the Ecosystem. Northwest Sci. 1998, 72, 34–39. [Google Scholar]
  7. Zimble, D.A.; Evans, D.L.; Carlson, G.C.; Parker, R.C.; Grado, S.C.; Gerard, P.D. Characterizing vertical forest structure using small-footprint airborne LiDAR. Remote Sens. Environ. 2003, 87, 171–182. [Google Scholar] [CrossRef]
  8. Miura, N.; Jones, S.D. Characterizing forest ecological structure using pulse types and heights of airborne laser scanning. Remote Sens. Environ. J. 2010, 114, 1069–1076. [Google Scholar] [CrossRef]
  9. Pasher, J.; King, D.J. Multivariate forest structure modelling and mapping using high resolution airborne imagery and topographic information. Remote Sens. Environ. 2010, 114, 1718–1732. [Google Scholar] [CrossRef]
  10. Corlett, R.T. Where are the Subtropics? Biotropica 2013, 45, 273–275. [Google Scholar] [CrossRef]
  11. Zhang, Y.J.; Cristiano, P.M.; Zhang, Y.F.; Campanello, P.I.; Tan, Z.H.; Zhang, Y.P.; Goldstein, G.; Cao, K.F. Carbon Economy of Subtropical Forests. Trop. Tree Physiol. 2016, 6, 337–355. [Google Scholar]
  12. Global Forest Resources Assessment 2015. Available online: http://www.fao.org/3/a-i4808e.pdf (accessed on 14 June 2017).
  13. Ouyang, S.; Xiang, W.; Wang, X.; Zeng, Y.; Lei, P.; Deng, X.; Peng, C. Significant effects of biodiversity on forest biomass during the succession of subtropical forest in south China. For. Ecol. Manag. 2016, 372, 291–302. [Google Scholar] [CrossRef]
  14. Bouvier, M.; De, M.; Renaud, J. Generalizing predictive models of forest inventory attributes using an area-based approach with airborne LiDAR data area-based approach with airborne LiDAR data. Remote Sens. Environ. 2015, 156, 322–334. [Google Scholar] [CrossRef]
  15. White, J.C.; Coops, N.C.; Wulder, M.A.; Vastaranta, M.; Hilker, T.; Tompalski, P.; White, J.C.; Coops, N.C.; Wulder, M.A.; Vastaranta, M.; et al. Remote Sensing Technologies for Enhancing Forest Inventories: A Review Remote Sensing Technologies for Enhancing Forest Inventories: A Review. Can. J. Remote Sens. 2016, 42, 619–641. [Google Scholar] [CrossRef]
  16. Mcroberts, R.E.; Tomppo, E.O. Remote sensing support for national forest inventories. Remote Sens. Environ. 2007, 110, 412–419. [Google Scholar] [CrossRef]
  17. Wulder, M.A.; White, J.C.; Nelson, R.F.; Næsset, E.; Ørka, H.O.; Coops, N.C.; Hilker, T.; Bater, C.W.; Gobakken, T. Lidar sampling for large-area forest characterization: A review. Remote Sens. Environ. 2012, 121, 196–209. [Google Scholar] [CrossRef]
  18. Wulder, M. Optical remote-sensing techniques for the assessment of forest inventory and biophysical parameters. Prog. Phys. Geogr. 1998, 22, 449–476. [Google Scholar] [CrossRef]
  19. Duncanson, L.I.; Niemann, K.O.; Wulder, M.A. Integration of GLAS and Landsat TM data for aboveground biomass estimation. Carbon 2010, 36, 129–141. [Google Scholar] [CrossRef]
  20. Lu, D.; Chen, Q.; Wang, G.; Moran, E.; Batistella, M.; Zhang, M.; Laurin, G.V.; Saah, D. Aboveground Forest Biomass Estimation with Landsat and LiDAR Data and Uncertainty Analysis of the Estimates. Int. J. For. Res. 2012, 2012. [Google Scholar] [CrossRef]
  21. Rosenqvist, Å.; Milne, A.; Lucas, R.; Imhoff, M.; Dobson, C. A review of remote sensing technology in support of the Kyoto Protocol. Environ. Sci. Policy 2003, 6, 441–455. [Google Scholar] [CrossRef]
  22. Yu, Y.; Yang, X.; Fan, W. Estimates of forest structure parameters from GLAS data and multi-angle imaging spectrometer data. Int. J. Appl. Earth Obs. Geoinf. 2015, 38, 65–71. [Google Scholar] [CrossRef]
  23. Li, W.; Niu, Z.; Chen, H.; Li, D.; Wu, M.; Zhao, W. Remote estimation of canopy height and aboveground biomass of maize using high-resolution stereo images from a low-cost unmanned aerial vehicle system. Ecol. Indic. 2016, 67, 637–648. [Google Scholar] [CrossRef]
  24. Véga, C.; Renaud, J.; Durrieu, S.; Bouvier, M. On the interest of penetration depth, canopy area and volume metrics to improve Lidar-based models of forest parameters. Remote Sens. Environ. 2016, 175, 32–42. [Google Scholar] [CrossRef]
  25. Means, J.E.; Acker, S.A.; Fitt, B.J.; Renslow, M.; Emerson, L.; Abstract, C.J.H. Predicting Forest Stand Characteristics with Airborne Scanning Lidar. Photogramm. Eng. Remote Sens. 2000, 66, 1367–1371. [Google Scholar]
  26. Silva, C.A.; Klauberg, C.; Hudak, A.T.; Vierling, L.A.; Liesenberg, V.; Carvalho, S.P.C.E.; Rodriguez, L.C.E. A principal component approach for predicting the stem volume in Eucalyptus plantations in Brazil using airborne LiDAR data. Forestry 2016, 89, 422–433. [Google Scholar] [CrossRef]
  27. Tesfamichael, S.G.; Beech, C. Combining Akaike’s Information Criterion and discrete return LiDAR data to estimate structural attributes of savanna woody vegetation. J. Arid Environ. 2016, 129, 25–34. [Google Scholar] [CrossRef]
  28. Sabol, J.; Procházka, D.; Patočka, Z. Development of models for forest variable estimation from airborne laser scanning data using an area-based approach at a plot level. J. For. Sci. 2016, 62, 137–142. [Google Scholar] [CrossRef]
  29. Frazer, G.W.; Magnussen, S.; Wulder, M.A.; Niemann, K.O. Simulated impact of sample plot size and co-registration error on the accuracy and uncertainty of LiDAR-derived estimates of forest stand biomass. Remote Sens. Environ. J. 2011, 115, 636–649. [Google Scholar] [CrossRef]
  30. Montealegre, A.L.; Riva, J.D.; García-martín, A. Use of low point density ALS data to estimate stand-level structural variables in Mediterranean Aleppo pine forest. Forestry 2016, 1–10. [Google Scholar] [CrossRef]
  31. Zolkos, S.G.; Goetz, S.J.; Dubayah, R. A meta-analysis of terrestrial aboveground biomass estimation using lidar remote sensing. Remote Sens. Environ. 2013, 128, 289–298. [Google Scholar] [CrossRef]
  32. Magnussen, S.; Næsset, E.; Gobakken, T. Reliability of LiDAR derived predictors of forest inventory attributes: A case study with Norway spruce. Remote Sens. Environ. 2010, 114, 700–712. [Google Scholar] [CrossRef]
  33. Magnussen, S.; Næsset, E.; Gobakken, T. Prediction of tree-size distributions and inventory variables from cumulants of canopy height distributions. Forestry 2013, 86, 583–595. [Google Scholar] [CrossRef]
  34. MacArthur, R.; MacArthur, J.W. On bird species diversity. Ecology 1961, 42, 594–598. [Google Scholar]
  35. Aber, J.D. A method for estimating foliage-height profiles in broad-leaved forests. Ecology 1979, 67, 35–40. [Google Scholar] [CrossRef]
  36. Lefsky, M.A.; Cohen, W.B.; Acker, S.A.; Parker, G.G. Lidar Remote Sensing of the Canopy Structure and Biophysical Properties of Douglas-Fir Western Hemlock Forests. Remote Sens. Environ. 1999, 361, 339–361. [Google Scholar] [CrossRef]
  37. Lovell, J.L.; Jupp, D.L.B.; Culvenor, D.S.; Coops, N.C.; Les, R. Using airborne and ground-based ranging lidar to measure canopy structure in Australian forests. Can. J. Remote Sens. 2003, 29, 607–622. [Google Scholar] [CrossRef]
  38. Coops, N.C.; Hilker, T.; Wulder, M.A.; St-Onge, B.; Newnham, G.; Siggins, A.; Trofymow, J.A. Estimating canopy structure of Douglas-fir forest stands from discrete-return LiDAR. Trees 2007, 21, 295–310. [Google Scholar] [CrossRef]
  39. Hilker, T.; Leeuwen, M.V.; Coops, N.C. Comparing canopy metrics derived from terrestrial and airborne laser scanning in a Douglas-fir dominated forest stand. Trees 2010, 24, 819–832. [Google Scholar] [CrossRef]
  40. Kraus, K.; Pfeifer, N. Determination of terrain models in wooded areas with airborne laser scanner data. ISPRS J. Photogramm. Remote Sens. 1998, 53, 193–203. [Google Scholar] [CrossRef]
  41. Song, Y.B.; Ding, Y.P. The development and latest progress of JSCORS. Bull. Surv. Mapp. 2009, 2, 73–74. [Google Scholar]
  42. Jiang, B.; Yuan, W.; Zhu, G. A preliminary study on the plantation biomass and produce structure of Pinus massoniana, Pinus elliottii and Pinus taeda. J. Zhejiang For. Sci. Technol. 1992, 12, 1–9. [Google Scholar]
  43. Ye, J.; Yang, Z. Biomass structure of planted Chinese fir in Southern Jiangsu province, China. Acta Ecol. Sin. 1983, 3, 7–14. [Google Scholar]
  44. Wang, Q.; Shi, Y. A preliminary study on the biomass and production of slash pine plantation in Jiangsu province. Acta Phytoecol. Geobot. Sin. 1990, 14, 2–12. [Google Scholar]
  45. Xu, J.; Wang, M.; Huang, Q.; Gong, S. Study on aboveground biomass model of natural individual trees of Quercus acutissima. Anhui For. Sci. Technol. 2011, 37, 3–6. [Google Scholar]
  46. Qian, G. Studies on the dynamic change of the net production quantity of liquidambar formosana plantation. Acta Agric. Univ. Jiangxiensis 2000, 22, 399–404. [Google Scholar]
  47. Sun, D.; Ruan, H.; Ye, J. Biomass Structure of Oak-Dominated Secondary Forest in Kongqingshan. Available online: http://refhub.elsevier.com/S0034-4257(16)30106-7/rf0260 (accessed on 14 June 2017).
  48. 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]
  49. Weishampel, J.F.; Drake, J.B.; Cooper, A.; Blair, J.B.; Hofton, M. Forest canopy recovery from the 1938 hurricane and subsequent salvage damage measured with airborne LiDAR. Remote Sens. Environ. 2007, 109, 142–153. [Google Scholar] [CrossRef]
  50. Bailey, R.; Dell, T. Quantifying diameter distributions with the Weibull function. For. Sci. 1973, 19, 97–104. [Google Scholar]
  51. Penner, M.; Woods, M.; Pitt, D.G. A comparison of airborne laser scanning and image point cloud derived tree size class distribution models in Boreal Ontario. Forests 2015, 6, 4034–4054. [Google Scholar] [CrossRef]
  52. Tompalski, P.; Coops, N.C.; White, J.C.; Wulder, M.A. Enriching ALS-derived area-based estimates of volume through tree-level downscaling. Forests 2015, 6, 2608–2630. [Google Scholar] [CrossRef]
  53. Zhao, J.; LI, J.; Liu, Q. Review of forest vertical structure parameter inversion based on remote sensing technology. J. Remote Sens. 2013, 17, 697. [Google Scholar]
  54. Jupp, D.L.B.; Culvenor, D.S.; Lovell, J.L.; Newnham, G.J.; Strahler, A.H.; Woodcock, C.E. Estimating forest LAI profiles and structural parameters using a ground-based laser called ‘Echidna®. Tree Physiol. 2008, 29, 171–181. [Google Scholar] [CrossRef] [PubMed]
  55. Hilker, T.; Coops, N.C.; Newnham, G.J.; Leeuwen, M.V.; Wulder, M.A.; Stewart, J.; Culvenor, D.S. Comparison of Terrestrial and Airborne LiDAR in Describing Stand Structure of a Thinned Lodgepole Pine Forest. J. For. 2012, 10, 97–104. [Google Scholar]
  56. Sprugel, D. Correcting for Bias in Log-Transformed Allometric Equations. Ecology 1983, 64, 209–210. [Google Scholar] [CrossRef]
  57. Næsset, E. Predicting forest stand characteristics with airborne scanning laser using a practical two-stage procedure and field data. Remote Sens. Environ. 2002, 80, 88–99. [Google Scholar] [CrossRef]
  58. Næsset, E.; Bollandsås, O.M.; Gobakken, T. Comparing regression methods in estimation of biophysical properties of forest stands from two different inventories using laser scanner data. Remote Sens. Environ. 2005, 94, 541–553. [Google Scholar] [CrossRef]
  59. García-gutiérrez, J.; Martínez-Álvarez, F.; Troncoso, A.; Riquelme, J.C. A comparison of machine learning regression techniques for LiDAR-derived estimation of forest variables. Neurocomputing 2015, 167, 24–31. [Google Scholar] [CrossRef]
  60. Naesset, E.; Okland, T. Estimating tree height and tree crown properties using airborne scanning laser in a boreal. Remote Sens. Environ. 2002, 79, 105–115. [Google Scholar] [CrossRef]
  61. Akaike, H. A new look at the statistical model identification. IEEE Trans. Autom. Control 1974, 19, 716–723. [Google Scholar] [CrossRef]
  62. Kutner, M.; Nachtsheim, C.; Neter, J.; Li, W. Applied Linear Statistical Models, 5th ed.; McGraw-Hill/Irwin: New York, NY, USA, 2004. [Google Scholar]
  63. Bengio, Y. No Unbiased Estimator of the Variance of K-Fold Cross-Validation. J. Mach. Learn. Res. 2004, 5, 1089–1105. [Google Scholar]
  64. Parker, G.G. Structure and Microclimate of Forest Canopies. In Forest Canopies; Lowman, M.D., Nadkarni, N.M., Eds.; Academic Press: San Diego, CA, USA, 1995; pp. 73–106. [Google Scholar]
  65. Wilkes, P.; Suarez, L.; Andrew, H.; Andrew, M.; William, W.; Mariela, S.-B.; Skidmore, A.K. Using discrete-return ALS to quantify number of canopy strata across diverse forest types Using discrete-return airborne laser scanning to quantify number of canopy strata across diverse forest types. Methods Ecol. Evol. 2015, 1–13. [Google Scholar] [CrossRef]
  66. Lefsky, M.A.; Harding, D.; Cohen, W.B.; Parker, G.; Shugart, H.H. Surface Lidar Remote Sensing of Basal Area and Biomass in Deciduous Forests of Eastern Maryland, USA. Remote Sens. Environ. 1999, 67, 83–98. [Google Scholar] [CrossRef]
  67. Leiterer, R.; Torabzadeh, H.; Furrer, R.; Schaepman, M.E.; Morsdorf, F. Towards Automated Characterization of Canopy Layering in Mixed Temperate Forests Using Airborne Laser Scanning. Forests 2015, 6, 4146–4167. [Google Scholar] [CrossRef]
  68. Van Leeuwen, M.; Nieuwenhuis, M. Retrieval of forest structural parameters using LiDAR remote sensing. Eur. J. For. Res. 2010, 129, 749–770. [Google Scholar] [CrossRef]
  69. Ishii, H. The role of crown architecture, leaf phenology and photosynthetic activity in promoting complementary use of light among coexisting species in temperate forests. Ecol. Res. 2010, 25, 715–722. [Google Scholar] [CrossRef]
  70. He, J.; Zhu, S. The Preliminary Study on the Law of the Crown Structure of Massons Pine in Mid-Subtropical Region. J. Mt. Agric. Biol. 1990, 9, 61–72. [Google Scholar]
  71. Liang, X.; Ye, W. Advances in Study of Forest Gaps. J. Trop. Subtrop. Bot. 2001, 9, 355–364. [Google Scholar]
  72. Kneeshaw, D.D.; Bergeron, Y. Canopy Gap Characteristics and Tree Replacement in the Southeastern Boreal Forest. Ecology 1998, 79, 783–794. [Google Scholar] [CrossRef]
  73. An, S.; Zhao, R. Analysis of Characteristics of Secondary Forest Vegetation in the North Subtropical Zone of China. J. Nanjing Univ. Sci. 1991, 27, 324–331. [Google Scholar]
  74. Spies, T.A.; Frankli, J.F. Gap Characteristics and Vegetation Response in Coniferous Forests of the Pacific Northwest. Ecology 1989, 70, 543–545. [Google Scholar] [CrossRef]
  75. Yang, J.; Liu, X.; Yang, X. Forest canopy structure, light environment and their effects on the vegetation pattern and growth of understory in forests. J. Sci. Teach. Univ. 2015, 35, 57–62. [Google Scholar]
  76. Montgomery, R.A.; Chazdon, R.L. Forest structure, canopy architectura, and light transmittance in old-growth and second-growth tropical rain forests. Ecology 2001, 82, 2707–2718. [Google Scholar] [CrossRef]
  77. Zhu, S.; He, J. The Preliminary Study on the Law of the Crown Structure of Masson’s pine in South-subtropcial. J. Mt. Agric. Biol. 1993, 12, 36–44. [Google Scholar]
  78. Messier, C.; Sylvain, P.; Bergeron, Y. Effects of Overstory and Understory Vegetation on the Understory Light Environment in Mixed Boreal Forests. J. Veg. Sci. 1998, 9, 511–520. [Google Scholar] [CrossRef]
  79. Parent, S.; Messier, C. A simple and efficient method to estimate microsite light availability under a forest canopy. Can. J. For. Res. 1996, 26, 151–154. [Google Scholar] [CrossRef]
  80. Zhao, F.; Yang, X.; Schull, M.A.; Román-colón, M.O.; Yao, T.; Wang, Z.; Zhang, Q.; Jupp, D.L.B.; Lovell, J.L.; Culvenor, D.S.; et al. Measuring effective leaf area index, foliage profile, and stand height in New England forest stands using a full-waveform ground-based lidar. Remote Sens. Environ. 2011, 115, 2954–2964. [Google Scholar] [CrossRef]
  81. Tang, H.; Dubayah, R.; Brolly, M.; Ganguly, S.; Zhang, G. Large-scale retrieval of leaf area index and vertical foliage profile from the spaceborne waveform lidar (GLAS/ICESat). Remote Sens. Environ. 2014, 154, 8–18. [Google Scholar] [CrossRef]
  82. Fu, T.; Pang, Y.; Huang, Q.F.; Liu, Q.W.; Xu, G.C. Prediction of Subtropical Forest Parameters Using Airborne Laser Scanner. J. Remote Sens. 2011, 15, 1092–1104. [Google Scholar]
  83. Xu, Z.; Cao, L.; Ruan, H.; Li, W.; Jiang, S. Inversion of subtropical forest stand characteristics by integrating very high resolution imagery acquired from UAV and LiDAR point-cloud. Chin. J. Plant Ecol. 2015, 39, 849–856. [Google Scholar]
  84. Wilkes, P.; Jones, S.; Suarez, L.; Haywood, A.; Mellora, A.; Soto-Berelov, M.; Woodgate, W. MAUP and LiDAR derived canopy structure. Int. Geosci. Remote Sens. Symp. 2013, 173–175. [Google Scholar]
  85. Maltamo, M.; Bollandsås, O.M.; Næsset, E.; Gobakken, T.; Packalén, P. Different plot selection strategies for field training data in ALS-assisted forest inventory. Forestry 2011, 84, 23–31. [Google Scholar] [CrossRef]
  86. Van der Zande, D.; Stuckens, J.; Verstraeten, W.W.; Mereu, S.; Muys, B.; Coppin, P. 3D modeling of light interception in heterogeneous forest canopies using ground-based LiDAR data. Int. J. Appl. Earth Obs. Geoinf. 2011, 13, 792–800. [Google Scholar] [CrossRef]
  87. Stoker, J. Volumetric Visualization of multiple-return Lidar Data: Using Voxels. Photogramm. Eng. Remote Sens. 2009, 75, 109–112. [Google Scholar]
  88. Zheng, G.; Moskal, L.M. Computational-Geometry-Based Retrieval of Computational-Geometry-Based Retrieval of Effective Leaf Area Index Using Terrestrial Laser Scanning. IEEE Geosci. Remote Sens. Soc. 2015, 50, 3958–3969. [Google Scholar] [CrossRef]
  89. Grau, E.; Durrieu, S.; Fournier, R.; Gastellu-etchegorry, J.; Yin, T. Estimation of 3D vegetation density with Terrestrial Laser Scanning data using voxels. A sensitivity analysis of influencing parameters. Remote Sens. Environ. 2017, 191, 373–388. [Google Scholar] [CrossRef]
  90. Li, J.; Hu, B.; Noland, T.L. Classification of tree species based on structural features derived from high density LiDAR data. Agric. For. Meteorol. 2013, 171–172, 104–114. [Google Scholar] [CrossRef]
  91. Kükenbrink, D.; Schneider, F.D.; Leiterer, R.; Schaepman, M.E.; Morsdorf, F. Quantification of hidden canopy volume of airborne laser scanning data using a voxel traversal algorithm. Remote Sens. Environ. 2017, 194, 424–436. [Google Scholar] [CrossRef]
  92. Oshio, H.; Asawa, T.; Hoyano, A.; Miyasaka, S. Estimation of the leaf area density distribution of individual trees using high-resolution and multi-return airborne LiDAR data. Remote Sens. Environ. 2015, 166, 116–125. [Google Scholar] [CrossRef]
  93. Béland, M.; Widlowski, J.; Fournier, R.A. A model for deriving voxel-level tree leaf area density estimates from ground-based LiDAR. Environ. Model. Softw. 2014, 51, 184–189. [Google Scholar] [CrossRef]
Figure 1. An overview of the workflow for forest structural parameters estimation. DTM: Digital Terrain Model.
Figure 1. An overview of the workflow for forest structural parameters estimation. DTM: Digital Terrain Model.
Remotesensing 09 00940 g001
Figure 2. Study site and distribution of sample plots.
Figure 2. Study site and distribution of sample plots.
Remotesensing 09 00940 g002
Figure 3. Examples of the three main forest types in study site. (a) Coniferous forest; (b) broad-leaved forest; (c) mixed forest.
Figure 3. Examples of the three main forest types in study site. (a) Coniferous forest; (b) broad-leaved forest; (c) mixed forest.
Remotesensing 09 00940 g003
Figure 4. The illustration of voxel-based canopy volume model. (a) A plot (30 × 30 m2) was stratified with voxelization and height bin is 0.5 m; (b) a voxel column was stratified in four structure classes (open gap, euphotic, oligophotic, closed gap) with canopy volume model approach; (c) canopy volume distribution, which shows the distribution of canopy structure classes after all columns were expanded in a panel; (d) the canopy volume profile, which was transformed from the canopy volume distribution diagram, shows the volume percentage of each class of total volume in each height interval.
Figure 4. The illustration of voxel-based canopy volume model. (a) A plot (30 × 30 m2) was stratified with voxelization and height bin is 0.5 m; (b) a voxel column was stratified in four structure classes (open gap, euphotic, oligophotic, closed gap) with canopy volume model approach; (c) canopy volume distribution, which shows the distribution of canopy structure classes after all columns were expanded in a panel; (d) the canopy volume profile, which was transformed from the canopy volume distribution diagram, shows the volume percentage of each class of total volume in each height interval.
Remotesensing 09 00940 g004
Figure 5. Canopy volume distributions for the plots in different forest types. (ac) Three typical plots of coniferous forest; (df) three typical plots of mixed forest; (gi) three typical plots of broad-leaved forest.
Figure 5. Canopy volume distributions for the plots in different forest types. (ac) Three typical plots of coniferous forest; (df) three typical plots of mixed forest; (gi) three typical plots of broad-leaved forest.
Remotesensing 09 00940 g005
Figure 6. Canopy volume profiles for the plots in different forest types. (ac) Three typical plots of coniferous forest; (eg) three typical plots of mixed forest; (ik) three typical plots of broad-leaved forest; (d,h,l) the mean canopy volume profiles in each forest.
Figure 6. Canopy volume profiles for the plots in different forest types. (ac) Three typical plots of coniferous forest; (eg) three typical plots of mixed forest; (ik) three typical plots of broad-leaved forest; (d,h,l) the mean canopy volume profiles in each forest.
Remotesensing 09 00940 g006
Figure 7. Foliage profiles for the plots in different forest types. (ac) Three typical plots of coniferous forest; (eg) three typical plots of mixed forest; (ik) three typical plots of broad-leaved forest; (d,h,l) the mean foliage profiles in each forest type.
Figure 7. Foliage profiles for the plots in different forest types. (ac) Three typical plots of coniferous forest; (eg) three typical plots of mixed forest; (ik) three typical plots of broad-leaved forest; (d,h,l) the mean foliage profiles in each forest type.
Remotesensing 09 00940 g007
Figure 8. The mean cumulative leaf area index profiles for plots in different forest types (a) and mean canopy height distribution: (b) coniferous forest; (c) mixed forest; (d) broad-leaved forest.
Figure 8. The mean cumulative leaf area index profiles for plots in different forest types (a) and mean canopy height distribution: (b) coniferous forest; (c) mixed forest; (d) broad-leaved forest.
Remotesensing 09 00940 g008
Figure 9. Scatter plots of field-measured and LiDAR estimated mean diameter at breast height (a), lorey’s mean height (b), stem density (c), basal area (d), volume (e) and aboveground biomass (f) for cross-validation in combination models.
Figure 9. Scatter plots of field-measured and LiDAR estimated mean diameter at breast height (a), lorey’s mean height (b), stem density (c), basal area (d), volume (e) and aboveground biomass (f) for cross-validation in combination models.
Remotesensing 09 00940 g009
Figure 10. Comparison of the forest structural parameters estimation accuracies for different voxel sizes. (a) The vertical resolution of voxel sizes is 0.5 m and the horizontal resolutions range from 1 m to 10 m; (b) the vertical resolution of voxel sizes is 1 m and the horizontal resolution range from 1 m to 10 m; (c) comparison of the difference of rRMSE values (∆rRMSE) of model accuracies between vertical resolution 0.5 m and vertical resolution 1 m. DBH: mean diameter at breast height; hLorey: Lorey’s mean height; N: Stem density; G: Basal area.
Figure 10. Comparison of the forest structural parameters estimation accuracies for different voxel sizes. (a) The vertical resolution of voxel sizes is 0.5 m and the horizontal resolutions range from 1 m to 10 m; (b) the vertical resolution of voxel sizes is 1 m and the horizontal resolution range from 1 m to 10 m; (c) comparison of the difference of rRMSE values (∆rRMSE) of model accuracies between vertical resolution 0.5 m and vertical resolution 1 m. DBH: mean diameter at breast height; hLorey: Lorey’s mean height; N: Stem density; G: Basal area.
Remotesensing 09 00940 g010
Table 1. A summary of plot-level forest structural parameters data.
Table 1. A summary of plot-level forest structural parameters data.
ParametersConiferous Forest (n = 14)Broad-Leaved Forest (n = 14)Mixed Forest (n = 23)
RangeMeanSDRangeMeanSDRangeMeanSD
DBH/cm8.08–19.2212.622.5311.63–20.9915.323.2910.58–19.6913.902.51
hLorey/m4.47–12.979.502.007.70–18.5211.352.757.79–14.1810.791.71
N/(ha−1)656–31671690.64643.15322.00–1833.001126.00428.55689.00–2344.001431.78438.40
G/(m2·ha−1)6.97–34.0723.086.7912.11–28.1021.923.8916.84–35.3723.984.46
V/(m3·ha−1)32.19–178.08116.5334.7590.62–212.45132.7732.3082.78–187.91131.9828.67
AGB/(Mg·ha−1)11.02–127.3969.7427.7632.03–219.6794.2844.9349.65–141.7389.3625.95
Notes: DBH: Mean diameter at breast height; hLorey: Lorey’s mean height; N: Stem density; G: Basal area; V: Volume; AGB: Aboveground biomass.
Table 2. The description of LiDAR canopy metrics.
Table 2. The description of LiDAR canopy metrics.
LiDAR MetricsDescription
Standard metrics
Height-basedPercentile heights (h25, h50, h75 and h95)The percentiles of the canopy height distributions
(25th, 50th, 75th and 95th) of first returns.
Mean height (hmean)Mean height above ground of all first returns.
Coefficient of variation of heights (hcv)Coefficient of variation of heights of all first returns.
Skewness and Kurtosis of heights
(i.e., hskewness and hkurtosis)
The skewness and kurtosis of the heights of all points.
Density-basedCanopy return density (d1, d3, d5, d7 and d9)The proportion of points above the quantiles
(10th, 30th, 50th, 70th and 80th) to total number of points.
Canopy cover above 2 m (CC2m)Percentages of first returns above 2 m.
Canopy metrics
Canopy volumeFilled and Empty zones of CVM
(i.e., Filled and Empty)
The voxels contained point clouds and voxels contained no point clouds within canopy spaces.
Open and Closed gap zones of CVM
(i.e., Open gap (OG) and Closed gap (CG))
The empty voxels located above and below the canopy respectively.
Euphotic and Oligophotic zones of CVM
(i.e., Euphotic (Eu) and Oligophotic (Oligo))
The voxels located within an uppermost percentile (65%) of all filled grid cells of that column, and voxels located below the point in the profile
Weibull-fittedα1 and β1 parameter of Weibull distributionThe scale parameter α and shape parameter β of the Weibull density distribution fitted to CHD.
α2 and β2 parameter of Weibull distributionThe scale parameter α and shape parameter β of the Weibull density distribution fitted to FP.
Table 3. A summary of selected metrics and accuracy assessment results of predictive models.
Table 3. A summary of selected metrics and accuracy assessment results of predictive models.
Forest TypesParametersSM ModelsCM ModelsCombination Models
Standard MetricsAdj-R2RMSErRMSE %Canopy MetricsAdj-R2RMSErRMSE %All MetricsAdj-R2RMSErRMSE %
All plotsDBH/cmh95, d1, d70.60 ***1.7212.33OG, Oligo, Empty, β20.50 ***1.8613.31hcv, h75, d1, Oligo0.61 ***1.6711.97
hLorey/mhcv, h95, d7, d90.75 ***0.979.15Oligo, Filled, Empty, α10.61 ***1.1811.13h50, d1, Empty, β10.77 ***0.908.54
N/(ha−1)hcv, d1, d7, d90.42 ***423.7529.86OG, Eu, Oligo, β10.39 ***415.1729.26d1, Oligo, α1, β10.45 ***410.0228.90
G/(m2·ha−1)h95, d3, d70.44 ***3.9917.23Oligo, Empty, α20.41 ***3.6715.82hkurtosis, h25, h95, Empty0.50 ***3.4714.96
V/(m3·ha−1)hcv, h25, h50, d30.46 ***22.3417.46OG, Eu, Oligo, α10.42 ***22.3617.48h75, Oligo, Empty, β10.58 ***21.0716.47
AGB/(Mg·ha−1)hkurtosis, h95, d3, d90.64 ***19.1722.47OG, Oligo, CG, Empty0.54 ***19.8423.25h95, d3, CC2m, Oligo0.66 ***18.2521.39
Coniferous forestDBH/cmh95, d1, d70.67 **1.209.50OG, Oligo, Empty, β20.541.4011.09hcv, h75, d1, Oligo0.74 **1.088.59
hLorey/mhcv, h95, d7, d90.661.0911.47Oligo, Filled, Empty, α10.641.2112.79h50, d1, Empty, β10.77 **0.9910.43
N/(ha−1)hcv, d1, d7, d90.60315.7818.68OG, Eu, Oligo, β10.58431.6525.53d1, Oligo, α1, β10.64339.2920.07
G/(m2·ha−1)h95, d3, d70.62 **4.5319.63Oligo, Empty, α20.554.7320.48hkurtosis, h25, h95, Empty0.69 **4.2318.32
V/(m3·ha−1)hcv, h25, h50, d30.69 **22.4019.22OG, Eu, Oligo, α10.72 **18.3215.72h75, Oligo, Empty, β10.78 **18.2115.63
AGB/(Mg·ha−1)hkurtosis, h95, d3, d90.72 **16.8624.17OG, Oligo, CG, Empty0.74 **18.5126.55h95, d3, CC2m, Oligo0.81 **14.5320.83
Broad-leaved forestDBH/cmh95, d1, d70.61 **1.7011.12OG, Oligo, Empty, β20.511.8111.79hcv, h75, d1, Oligo0.681.5410.06
hLorey/mhcv, h95, d7, d90.84 ***0.786.91Oligo, Filled, Empty, α10.83 ***0.887.72h50, d1, Empty, β10.88 ***0.726.39
N/(ha−1)hcv, d1, d7, d90.60298.9926.55OG, Eu, Oligo, β10.52299.8226.63d1, Oligo, α1, β10.62273.4924.29
G/(m2·ha−1)h95, d3, d70.542.6211.96Oligo, Empty, α20.502.7412.48hkurtosis, h25, h95, Empty0.632.4911.34
V/(m3·ha−1)hcv, h25, h50, d30.5619.4914.68OG, Eu, Oligo, α10.5818.9714.28h75, Oligo, Empty, β10.6716.6512.54
AGB/(Mg·ha−1)hkurtosis, h95, d3, d90.5726.8028.42OG, Oligo, CG, Empty0.6026.4528.05h95, d3, CC2m, Oligo0.6626.6728.29
Mixed forestDBH/cmh95, d1, d70.48 **1.6611.94OG, Oligo, Empty, β20.481.7812.79hcv, h75, d1, Oligo0.55 **1.5811.34
hLorey/mhcv, h95, d7, d90.81 ***0.605.60Oligo, Filled, Empty, α10.75 ***0.756.94h50, d1, Empty, β10.84 ***0.555.13
N/(ha−1)hcv, d1, d7, d90.48 **336.7328.52OG, Eu, Oligo, β10.44324.7322.68d1, Oligo, α1, β10.50 ***319.0522.28
G/(m2·ha−1)h95, d3, d70.45 **3.0812.86Oligo, Empty, α20.45 ***3.1213.01hkurtosis, h25, h95, Empty0.56 **2.7711.56
V/(m3·ha−1)hcv, h25, h50, d30.60 ***16.7612.70OG, Eu, Oligo, α10.65 ***16.3112.36h75, Oligo, Empty, β10.71 ***15.8712.02
AGB/(Mg·ha−1)hkurtosis, h95, d3, d90.64 ***13.2014.77OG, Oligo, CG, Empty0.71 ***13.0014.55h95, d3, CC2m, Oligo0.79 ***10.8912.19
Notes: Level of significance: NS = not significant (>0.05); ** <0.01; *** <0.001; DBH: mean diameter at breast height; hLorey: Lorey’s mean height; N: Stem density; G: Basal area; V: Volume; AGB: Aboveground biomass. OG: Open gap; Oligo: Oligophotic; Eu: Euphotic; CG: Closed gap.

Share and Cite

MDPI and ACS Style

Zhang, Z.; Cao, L.; She, G. Estimating Forest Structural Parameters Using Canopy Metrics Derived from Airborne LiDAR Data in Subtropical Forests. Remote Sens. 2017, 9, 940. https://doi.org/10.3390/rs9090940

AMA Style

Zhang Z, Cao L, She G. Estimating Forest Structural Parameters Using Canopy Metrics Derived from Airborne LiDAR Data in Subtropical Forests. Remote Sensing. 2017; 9(9):940. https://doi.org/10.3390/rs9090940

Chicago/Turabian Style

Zhang, Zhengnan, Lin Cao, and Guanghui She. 2017. "Estimating Forest Structural Parameters Using Canopy Metrics Derived from Airborne LiDAR Data in Subtropical Forests" Remote Sensing 9, no. 9: 940. https://doi.org/10.3390/rs9090940

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop