Review Article  Open Access
LiDAR Forest Inventory with SingleTree, Double, and SinglePhase Procedures
Abstract
Light Detection and Ranging (LiDAR) data at 0.5–2 m postings were used with doublesample, stratified procedures involving singletree relationships in mixed, and single species stands to yield sampling errors ranging from % to %. LiDAR samples were selected with focal filter procedures and heights computed from interpolated canopy and DEM surfaces. Tree dbh and height data were obtained at various ratios of LiDAR, ground samples for DGPS located ground plots. Dbhheight and groundLiDAR height models were used to predict dbh and compute Phase 2 estimates of basal area and volume. Phase 1 estimates were computed using the species probability distribution from ground plots in each strata. Phase 2 estimates were computed by randomly assigning LiDAR heights to species groups using a Monte Carlo simulation for each ground plot. There was no statistical difference between volume estimates from 0.5 m and 1 m LiDAR densities. Volume estimates from singlephase LiDAR procedures utilizing existing tree attributes and height bias relationships were obtained with sampling errors of 1.8% to 5.5%.
1. Introduction
Light detection and ranging (LiDAR) is a relatively new remote sensing tool that has the potential for use in the acquisition of measurement data for inventories of standing timber. LiDAR systems have been used in a variety of forestry applications [1–3] for the quantification of biomass [4], basal area, and tree and stand height estimates. Stand level, LiDAR inventory procedures involving average values of tree attributes such as dominant height, mean diameter, basal area, and volume have been applied to obtain unbiased stand level predictions [5–7]. Since LiDAR has the capability to detect individual trees and measure tree height with predictable bias when correlated with ground measurements [8, 9], stratalevel inventory estimates involving individual tree, doublesample inventory procedure have been used by researchers from Mississippi State University in conifer and mixed hardwood stands in the Northwest and Southeast [10–16]. The individual tree approach to stand inventory when combined with doublesample, ground procedures permits relatively precise estimates of volume with a simple prediction function for groundLiDAR height bias, and groundbased attribute relationship functions for tree diameter and total height which can be used with any standard, standing tree volume function. Stand level approaches involving average tree attribute values for sampling units require more sophisticated prediction models than an individual tree approach and procedures that differ radically from traditional groundbased inventory methods. The objective of this paper is to summarize and discuss the procedures, models, and advantages/disadvantages of the singletree approach to using LiDAR data in double and singlephase forest inventory methods.
2. Methods
2.1. Flight Planning for an LiDAR Inventory
Smallfootprint, multireturn LiDAR data have been acquired with various sensors to attain nominal posting spacings of 0.5–2.0 m, 0.25–4 points/m^{2}, and footprint sizes of 0.122–0.330 m for two returns per pulse (Table 1). Aircraft altitudes of 600–1000 m and swath widths of 189–609 m were used. The minimum required density of LiDAR hits is a function of the crown size, average height, and spatial density of the sample trees in the primary canopy. Acceptable sampling statistics were attained for sparse densities of large crown conifers in Idaho (11.5% sampling error at the 95% confidence level with a standard error of 5 m^{3} [11] with 0.25 points/m^{2}; however, 1 point/m^{2} was required to achieve acceptable inventory results in natural pine and mixed pinehardwood stands (7.6% sampling error, [12]) and 2 points/m^{2} in young (6+ years) pine plantations (2.2% sampling error, [14]) in the Southeast. Increasing LiDAR density from 2 to 4 points/m^{2} did not statistically improve the volume estimation precision, and the increased statistical “noise” in the highdensity LiDAR data translated into additional sampling error about the volume estimate.

Target aircraft altitude is a function of the desired swath width and scan angle for the LiDAR pulse generator and sensor and the technical ability of the sensor to achieve the desired posting density. The swath width diminishes as the desired posting density increases, but swath widths of 150–300 m can be achieved with 2–4 points/m^{2} and are of sufficient width to accommodate traditional size, sample ground plots within the swath interior, and minimizing “edge effects” of LiDAR. An important factor influencing desired swath width was the size of the groundbased sample plots used in the inventory procedure. The swath should be sufficiently wide to encompass the sample ground and LiDAR plots within the center onethird of the swath so as to minimize the “edge effects” of the LiDAR data. Tree attribute measurements are severely compromised at the extreme edges of the swath and scan angle because the density of pulse hits is reduced, and the pulse pattern is quite irregular.
Percent LiDAR coverage is a function of economics and inventory design. LiDAR data is relatively expensive to obtain and complete area coverage that is normally not required for most timber inventory designs. In some instances, the costs of complete LiDAR coverage to produce an accurate uptodate Digital Elevation Model (DEM) may be more justifiable than the expense for a timber inventory. The use of a current Geographical Information System (GIS) to locate flight lines that cross the desired sampling strata can minimize the percent coverage of the LiDAR area. Most forested areas to be inventoried can be flown with ten percent or less LiDAR area coverage by orienting flight lines so as to cross the target inventory strata at desired flight line intervals.
2.2. Field and LiDAR Plot Design and Procedures
Inventory design for the singletree LiDAR applications involved the use of circular [12, 14] or rectangular plots [11] with all plots being Phase 1 LiDAR plots and every th plot as a Phase 2 ground plot. Field designs varied from a ratio of LiDAR to ground circular plots in a nested arrangement to a ratio with rectangular or circular plots along a flight line (Figure 1). UTM coordinates were established at the center of each circular Phase 2 plot or at the endpoints of rectangular plots for navigation with a realtime Differential Global Positioning System (DGPS). Differential corrections from either the U.S. Government WAAS or private enterprise OmniStar geostationery satellite were obtained satisfactorily under tree canopies by using a large dome antenna. Based on informal field tests on surveyed bench marks, field locations were obtained with approximately 1 m accuracies with both systems.
2.3. LiDAR Surfacing for Tree Location and Height Determination
The LiDAR data were processed to produce a ground surface or digital terrain model (DTM) and a tree surface for determination of sample tree locations and tree heights within the sample field and LiDAR plot areas. LiDAR datasets were surfaced to produce 1st return canopy and last return DTM with 0.2 m cell sizes using a linear interpolation technique. Tree locations and heights were determined with algorithms and focal filter procedures developed by McCombs et al. [16] that used a variable search window radius based on relative tree density. These procedures used moving 2.5, 4.0, or 5.5 ft radius search windows to identify each tree peak as the point that is higher than 85% of the surrounding maxima from one of the three search window, radius files. Tree height was interpreted as the difference between canopy and DTM values at each tree peak location. Tree heights were converted to point coverages and clipped to sample area boundaries using UTM coordinates to describe sample plot locations and sizes.
A spatial filtering technique derived from image analysis called smoothing was used to reduce commission errors by minimizing the abrupt elevation changes in the initial canopy surface. The Focal analysis option in ERDAS’ Imagine software performed smoothing based on userdefined inputs for window size and preferred statistical procedure. A 5by5 pixel window was used to create a 1 m^{2} filter that would avoid removal of small peaks in the canopy surface (small trees), while maximizing the smoothing function. The filter moved across the LiDAR canopy surface, pixel by pixel, averaged the values within the window, and placed the result in the center pixel.
Smoothing heights on LiDAR surfaces improved the relationship between LiDAR and ground tree heights in terms of , reduced height biases for hardwoods, increased height biases for pines, and improved target recognition in terms of trees/ac estimates. There were, however, no statistical differences () between doublesample regression volume estimates with smoothed versus unsmoothed LiDAR surfaces from low or highdensity LiDAR. Standard errors and sampling errors of the regression estimates were lower for all unsmoothed LiDAR data models than with smoothed data models. Thus, smoothing heights on LiDAR surfaces did not produce a statistical gain for volume estimation using doublesample procedures.
3. Results
3.1. DoubleSample, Regression Estimator Procedures
The doublesample model widely used with groundbased point sampling [17] and aerial photogrammetric inventories and adapted for these studies was
With traditional aerial photogrammetric inventories, the and variables are photographic volume/unit area and ground volume/unit area from Phase 1 and Phase 2 plots, respectively, and is the regression slope coefficient for (ground volume) over (photo volume) on ground plots. Thus, a Phase 1 (large sample) variable such as remotely sensed (i.e., photographic or LiDAR derived) volume has a strong, identifiable relationship with a Phase 2 (small sample) variable such as ground volume.
In applications of the doublesample model with singletree LiDAR data, Phase 2 sample tree measures of dbh and height were used to derive heightdbh and dbhheight equations of the model type:
where was ground measured height, dbh was ground measured dbh, and age was average stand age (years) from GIS data. The age variable in Models (2) and (3) was removed when age did not contribute significantly to the relationship. Models (2) and (3) were derived from the groundmeasured sample trees, but one is not a back transformation of the other. Model (2) was applied to groundplot trees, where dbh was measured on all trees and heights on a subsample. Model (3) was applied to LiDAR derived tree heights to obtain a dbh for singletree volume computation.
Generally only 2 trees per ground plot were measured for height; dbh was measured on all trees. The heightdbh Model (2) was applied to trees on the ground plots for which height was not measured to obtain a height for singletree volume computation. Sample tree heights from the Phase 2 ground plots were used to predict ground height of target trees identified on LiDAR surfaces. The dbhheight Model (3) was applied to the biasadjusted, singletree LiDAR height from the groundLiDAR height bias Model:
where was measured ground height of trees on Phase 2 plots, and was interpolated height of the same trees from the LiDAR surface.
Derived dbh on LiDAR plots and derived height on ground plots permitted the use of a standard, standingtree volume equation with dbh and height as variables to predict volume. Thus, the doublesample models used in this study involved LiDAR mean estimates of basal area (LiBA from Phase 1 and liba from Phase 2 with matching ground plot), and volume (LiVOL from Phase 1 and livol from Phase 2 with matching ground plot) for the variables as with variance:
where bar was Phase 2 mean ground volume, was the regression slope coefficient for (ground volume/unit area) over (LiDAR volume/unit area or basal area/unit area on ground plot), and was volume or basal area on the LiDAR plot. Data were fitted to Models (5) for all data combined (i.e., nonstratified), each agesclass strata, and combined strata. Combined strata, linear regression estimates of volume, and associated standard error of each doublesample model were obtained by
where and were Phase 1 and 2 sample sizes, respectively, for stratum , to strata.
All doublesample volume computations were performed with the Windows software program LIDARDS (LiDAR DoubleSample) developed by Parker [18]. The software allowed the user to specify dbh limits for speciesproduct classes, regression coefficients for the dbhheight and ground heightLiDAR height models, stratum definitions of beginning and ending plot numbers and average age, and to enter comma delimited data files of Phase 1 LiDAR heights and Phase 2 groundplot trees (species, product, dbh, and height of sample trees). A speciesproduct class is a userdefined combination of a species (e.g., fir, spruce, etc.) and a merchantable tree product (e.g., pulpwood, sawtimber, veneer, etc.). LiDAR heights in the Phase 1 data were allocated in a Monte Carlo simulation to speciesproduct classes on each matching Phase 2 ground plot on the basis of percent distribution by numbers on the ground plot. Since species and dbh of the LiDAR trees are unknown, the Monte Carlo simulation (50 iterations) would randomly allocate the LiDAR derived trees (dbh predicted from adjusted LiDARtoground height) to speciesproduct classes and obtain a mean basal area and volume estimate for the speciesproduct class. Thus, basal area and volume estimates from Phase 1 LiDAR plots that had a matching Phase 2 ground plot became Phase 2 LiDAR plots. Phase 1 LiDAR heights that did not have a matching Phase 2 ground plot were randomly allocated to encountered species classes in each stratum in a single iteration and used to compute mean estimates of numbers of trees, basal area, and volume. Phase 2 tree measures of dbh and height were used to compute LiDAR estimates of mean basal area and volume by using fieldderived dbhheight equations to predict dbh from LiDAR height and volume. Predicted dbh and height were used in a singletree volume function to predict individual/single tree volume. Doublesample volume estimates and associated precision statistics were computed with Models (5) for each stratum and with Models (7) for combined strata.
3.2. SinglePhase Inventory Procedures
A recent study by Williams [15] investigated the minimal data inputs for an LiDARbased timber inventory with singlephase procedures. Since the ground phase of a doublesample field procedure is both expensive and time consuming, Williams investigated the feasibility of using LiDAR data in singletree approach to obtain volume estimates by stratum with a singlephase inventory procedure. Previous studies have shown that LiDAR can provide precise, but biased estimates of tree numbers and heights. If the assumptions are made that (1) the LiDAR height bias is known and relatively constant for a given speciesorigin class (e.g., pine plantations) and (2) previously established tree attribute relationships are also known, inventory estimates of volume can be obtained with a singletree approach and singlephase procedures from LiDAR data only.
The tree attribute relationship from Model (3) was developed from ground measurements of dbh and height within Continuous Forest Inventory (CFI) plots and from Phase 2 ground plots in the doublesample approach by Parker and Evans [14]. Sample trees were randomly selected from the original datasets in groups of 75 and fitted to tree attribute Model (3) under the assumption that ground data were available from previous studies. The groundLiDAR height bias equation obtained from Model (4) [14] was assumed to be constant and known. LiDAR derived heights from Phase 1 plots were adjusted for bias with Model (4) then used with Model (3) to obtain singletree dbh estimates for use with a singletree volume function in a conventional, singlephase inventory processor. The stratum volume estimates and precision statistics were compared to estimates obtained from Phase 2 ground plots by Parker and Evans [14].
The singletree, singlephase volume estimates from LiDAR data compared favorably with ground plot estimates (Table 2). At the tract level for 20 age class strata on 10 443 ha, there was no statistical difference between the singlephase LiDAR estimates and the ground plot estimate. Singlephase volume estimates were obtained for the full dataset of Phase 1 LiDAR plots, the Phase 1 LiDAR plots that had a matching Phase 2 ground plot, and 5 iterations of reducing the LiDAR dataset to a ratio with ground plots within each stratum. The Williams [15] study found no statistical difference () between tractlevel, singlephase volume estimates, where the singletree relationship model was developed with 1539 sample trees from the Phase 2 ground plots by Parker and Evans [14], 1509 trees from the regional CFI plots, or 5 iterations of 75 randomly selected trees from the Phase 2 dataset. The study concluded that precise singlephase LiDAR inventory estimates is feasible with minimal inputs of ground data for establishing tree attribute relationships. A potential application of the singletree, singlephase inventory procedures would be the rapid posthinning inventory of pine plantations and periodic inventories of forested holdings.

4. Discussion about SingleTree LiDAR Inventory Procedures
LiDAR provides precise , , and coordinate data that can be used to extract tree heights and locations; however, there are several sources of bias that can impact the accuracy of a per unit area volume estimate. Height bias is primarily caused by the failure of the laser pulse to hit the terminal leader, but this bias can be predicted with acceptable success in conifers but not in hardwoods with broad rounded crowns. Height bias can also be introduced by the interpolation of tree heights with mixed linear and nonlinear procedures within the same dataset. Tree count bias has at least two sources of origin; trees in the mid and lower canopy layers are hidden from the laser pulse by a dominant canopy, and tree maxima locations may not be interpreted correctly during the LiDAR surfacing and height extraction process.
The precision of volume estimates with singletree LiDAR procedures in a doublesample process is not affected by the height or tree count bias inherent in LiDAR data. These biases are effectively adjusted during the doublesample inventory procedures. The height bias can be adjusted prior to singletree volume computations with LiDAR derived heights or afterwards during the doublesample volume adjustment process. Height bias correction prior to volume computation improves the accuracy of the resulting per unit area volume estimate. The tree count bias caused by canopy coverage or LiDAR surfacing/processing is also effectively adjusted through the doublesample volume computations. Tree count bias, however, has a major impact on the accuracy of volume estimates in a singletree procedure computed with singlephase inventory methods. Thus, if there is any doubt about the validity of the tree counts or height bias during LiDAR processing, a doublesample volume computation process should be used.
Ground sample tree measurements are needed to establish the groundLiDAR height bias and the relationship between standing tree height and dbh. The number of sample trees is dependent upon the number of parameters in the regression models and variation in the data. A reasonable rule of thumb is 25 samples per parameter estimated in a regression model. Since the LiDAR height bias is relatively constant for a species in a local area and the dbhheight relationship for a given set of speciessite conditions is also relatively stable, the sample trees can be obtained from either the LiDAR inventory or surrounding forested areas.
Establishing the LiDAR to ground tree height bias requires the matching of trees on the ground and on the LiDAR surface. Past experience has shown that the location of a plot center must be done with a real time, DGPS, and the distance and direction to the sample trees from the plot center should be obtained with a laser so that the , coordinates of the sample trees can be located on the LiDAR surface.
Sample plot size and shape on the ground and on the LiDAR surface should be a function of tree density on the ground and the LiDAR processing procedures employed. Experience has shown that the ground plot size should be adjusted such that a minimum of 6 and a maximum of approximately 15 trees should be selected. The minimum number is associated with the within and between plot variation, and the maximum is a logistical consideration for minimizing omission/commission errors in tallying trees. Rectangular plots are easier to handle during the LiDAR processing, but more difficult to establish in the field and to use in establishing distance and direction to sample trees from a DGPS location. The cost of an LiDAR inventory can be minimized by flying only a portion of the desired inventory area. As long as the LiDAR swaths cover the desired strata, ground plots could be located randomly or systematically within strata within swaths. DGPS permits the location of sample plots and trees with relative ease and precision.
References
 S. Magnussen and P. Boudewyn, “Derivations of stand heights from airborne laser scanner data with canopybased quantile estimators,” Canadian Journal of Forest Research, vol. 28, no. 7, pp. 1016–1031, 1998. View at: Publisher Site  Google Scholar
 M. A. Lefsky, W. B. Cohen, S. A. Acker, G. G. Parker, T. A. Spies, and D. Harding, “LiDAR remote sensing of the canopy structure and biophysical properties of Douglasfir western hemlock forests,” Remote Sensing of Environment, vol. 70, no. 3, pp. 339–361, 1999. View at: Publisher Site  Google Scholar
 J. E. Means, S. A. Acker, B. J. Fitt, M. Renslow, L. Emerson, and C. J. Hendrix, “Predicting forest stand characteristics with airborne scanning LiDAR,” Photogrammetric Engineering & Remote Sensing, vol. 66, no. 11, pp. 1367–1371, 2000. View at: Google Scholar
 R. Nelson, M. A. Valenti, A. Short, and C. Keller, “A multiple resource inventory of Delaware using airborne laser data,” BioScience, vol. 53, no. 10, pp. 981–992, 2003. View at: Publisher Site  Google Scholar
 E. Næsset, “Predicting forest stand characteristics with airborne scanning laser using a practical twostage procedure and field data,” Remote Sensing of Environment, vol. 80, no. 1, pp. 88–99, 2002. View at: Publisher Site  Google Scholar
 E. Næsset, “Accuracy of forest inventory using airborne laser scanning: evaluating the first nordic fullscale operational project,” Scandinavian Journal of Forest Research, vol. 19, no. 6, pp. 554–557, 2004. View at: Publisher Site  Google Scholar
 S. C. Popescu, R. H. Wynne, and R. F. Nelson, “Estimating plotlevel tree heights with LiDAR: local filtering with a canopyheight based variable window size,” Computers and Electronics in Agriculture, vol. 37, no. 1–3, pp. 71–95, 2003. View at: Publisher Site  Google Scholar
 Å. Persson, J. Holmgren, and U. Söderman, “Detecting and measuring individual trees using an airborne laser scanner,” Photogrammetric Engineering & Remote Sensing, vol. 68, no. 9, pp. 925–932, 2002. View at: Google Scholar
 J. Holmgren, “Prediction of tree height, basal area and stem volume in forest stands using airborne laser scanning,” Scandinavian Journal of Forest Research, vol. 19, no. 6, pp. 543–553, 2004. View at: Publisher Site  Google Scholar
 C. A. Collins, Comparing integrated LiDAR and multispectral data with field measurements in hardwood stands, M.S. thesis, Mississippi State University, Starkville, Miss, USA, 2003.
 R. C. Parker and D. L. Evans, “An application of LiDAR in a doublesample forest inventory,” Western Journal of Applied Forestry, vol. 19, no. 2, pp. 95–101, 2004. View at: Google Scholar
 R. C. Parker and P. A. Glass, “High versus lowdensity LiDAR in a doublesample forest inventory,” Southern Journal of Applied Forestry, vol. 28, no. 4, pp. 205–210, 2004. View at: Google Scholar
 R. C. Parker and A. L. Mitchel, “Smoothed versus unsmoothed LiDAR in a doublesample forest inventory,” Southern Journal of Applied Forestry, vol. 29, no. 1, pp. 40–47, 2005. View at: Google Scholar
 R. C. Parker and D. L. Evans, “Industrial application of a LiDAR doublesample forest inventory,” Research Bulletin, Forest and Wildlife Research Center, Mississippi State University, Starkville, Miss, USA, 2006. View at: Google Scholar
 E. C. Williams, Minimum data input for a LiDAR based timber inventory, M.S. thesis, Mississippi State University, Starkville, Miss, USA, 2006.
 J. W. McCombs, S. D. Roberts, and D. L. Evans, “Influence of fusing LiDAR and multispectral imagery on remotely sensed estimates of stand density and mean tree height in a managed loblolly pine plantation,” Forest Science, vol. 49, no. 3, pp. 457–466, 2003. View at: Google Scholar
 T. E. Avery and H. E. Burkhart, Forest Measurements, McGrawHill, New York, NY, USA, 5th edition, 2002.
 R. C. Parker, “Computer automation of a LiDAR doublesample forest inventory,” Research Bulletin, Forest and Wildlife Research Center, Mississippi State University, Starkville, Miss, USA, 2005. View at: Google Scholar
Copyright
Copyright © 2009 Robert C. Parker and David L. Evans. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.