Next Article in Journal
Landsat-TM-Based Discrimination of Lithological Units Associated with the Purtuniq Ophiolite, Quebec, Canada
Next Article in Special Issue
Application of Semi-Automated Filter to Improve Waveform Lidar Sub-Canopy Elevation Model
Previous Article in Journal
Dielectric and Radiative Properties of Sea Foam at Microwave Frequencies: Conceptual Understanding of Foam Emissivity
Previous Article in Special Issue
Comparison of Methods for Estimation of Stem Volume, Stem Number and Basal Area from Airborne Laser Scanning Data in a Hemi-Boreal Forest
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Remote Sensing
Volume 4
Issue 5
10.3390/rs4051190
remotesensing-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Advances in Forest Inventory Using Airborne Laser Scanning

by
Juha Hyyppä
1,*,
Xiaowei Yu
1,
Hannu Hyyppä
2,
Mikko Vastaranta
3,
Markus Holopainen
3,
Antero Kukko
1,
Harri Kaartinen
1,
Anttoni Jaakkola
1,
Matti Vaaja
2,
Jarkko Koskinen
1 and
Petteri Alho
4
1
Department of Remote Sensing and Photogrammetry, Finnish Geodetic Institute, FI-02431 Masala, Finland
2
School of Science and Technology, Aalto University, FI-00076 Aalto, Finland
3
Department of Forest Sciences, University of Helsinki, FI-00014 Helsinki, Finland
4
Department of Geography and Geology, University of Turku, FI-20014 Turku, Finland
*
Author to whom correspondence should be addressed.
Remote Sens. 2012, 4(5), 1190-1207; https://doi.org/10.3390/rs4051190
Submission received: 15 March 2012 / Revised: 23 April 2012 / Accepted: 25 April 2012 / Published: 3 May 2012
(This article belongs to the Special Issue Laser Scanning in Forests)
Browse Figures
Versions Notes

Abstract

:
We present two improvements for laser-based forest inventory. The first improvement is based on using last pulse data for tree detection. When trees overlap, the surface model between the trees corresponding to the first pulse stays high, whereas the corresponding model from the last pulse results in a drop in elevation, due to its better penetration between the trees. This drop in elevation can be used for separating trees. In a test carried out in Evo, Southern Finland, we used 292 forests plots consisting of more than 5,500 trees and airborne laser scanning (ALS) data comprised of 12.7 emitted laser pulses per m2. With last pulse data, an improvement of 6% for individual tree detection was obtained when compared to using first pulse data. The improvement increased with an increasing number of stems per plot and with decreasing diameter breast height (DBH). The results confirm that there is also substantial information for tree detection in last pulse data. The second improvement is based on the use of individual tree-based features in addition to the statistical point height metrics in area-based prediction of forest variables. The commonly-used ALS point height metrics and individual tree-based features were fused into the non-parametric estimation of forest variables. By using only four individual tree-based features, stem volume estimation improved when compared to the use of statistical point height metrics. For DBH estimation, the point height metrics and individual tree-based features complemented each other. Predictions were validated at plot level.

1. Introduction

Trees are important for the carbon balance of the Earth. Forests have great economic and ecological importance. In Finland, about 77% of the country’s land area is forested, which is the highest percentage in Europe. International interest in biomass detection is strongly linked to forest health, photosynthetic activity and other processes related to the carbon cycle [1] and the variability of the climate. There is a growing need and constant shortage of data for improved forest monitoring (e.g., [2]).
Approaches aimed at obtaining forest and forestry data from airborne laser scanning (ALS) data have been divided into two groups [3]: (1) area-based approaches (ABAs) and (2) individual/single-tree detection approaches (ITDs). ABA prediction of forest variables is based on the statistical dependency between the variables measured in the field and the predictor features derived from ALS data [4–12]. The sample unit in the ABA is most often a grid cell, the size of which depends on the size of the field-measured training plot. Stand-level forest inventory results are aggregated by summing and weighting the grid-level predictions inside the stand. When using ITD techniques, individual trees are detected and tree-level variables, such as height and volume, are measured or predicted from the ALS data, i.e., the basic unit is an individual tree. Then the stand-level forest inventory results are aggregated by summing up the tree data. The ABA does not make use of the neighborhood data of laser returns, whereas in ITD [13–19] the neighborhood data are applied using pattern recognition methods to locate individual trees and to derive the physical features of individual trees, such as height, species, crown diameter and crown volume. On the other hand, ABAs are based on the height and density data acquired by ALS, which are highly correlated with the forest variables. Currently, the ABA is operationally applied in the Nordic countries when carrying out standwise forest inventories. Some 5.5 Mha of Finland’s forests have already been inventoried by applying ALS.
Initially, ITD began with the manual interpretation of analogue aerial images [20], followed by attempts to automate this task [21]. In recent years, there have been several attempts to improve both image-based and laser-based ITD [22–71]. An international comparison of the use of ALS for ITD was reported by Kaartinen and Hyyppä [34] and more recently by Vauhkonen et al. [63]. In [34], it was concluded that (1) the ITD method is the main factor impacting on the achieved location, tree height and inventory accuracy, and that (2) the impact of an increase in point density is marginal when compared to the effect of the ITD method. More detailed comparisons of the extraction techniques were carried out in [68], and it was found that methods relying on pixel-based Canopy Height Models (CHMs) are not suitable for finding suppressed trees and that full-waveform or laser-point based techniques should be further investigated (as is also proposed in [42,69,70]).
The ABA has been performed with statistical metrics, calculated from a point cloud, and ITD has been done using the dimensions of the detected trees. However, the methods have been converged for a long time (e.g., [24,26,35,41,56,71]). Already in [16,17], individual tree methodology was proposed and tested at stand level. Hyyppä et al. [24,26] proposed an individual tree-based technique for focusing directly on tree clusters. A similar concept was adapted later on by Breidenbach et al. [44]. Since the work done by Villikka et al. [35], density-related and height-related features, derived from first and last pulse ALS data, have also been used in the estimation of individual tree variables. In the past, the ABA has relied on large, meticulously collected field data from experimental plots, whereas ITD has provided results either without any field date, or with some reference to it [3].
Conventionally, trees are detected from CHMs, which are interpolated from the point data [3,63]. In the present study, we demonstrate that the use of last pulse data is beneficial for tree detection. Firstly, we demonstrate by means of a large, real-life data set that trees can be discriminated better by using last pulse data than by using first pulse data. Secondly, we show by means of examples that the use of individual tree features as predictors increases the accuracy of stem volume in area-based predictions, which suggests that individual tree-based features should be used in addition to point height metrics in the ABA. These two improvements are separately depicted. The use of last pulse is described for tree detection in Sections 2.4 (the method of verification) and 3.1 (results). The improvements in area-based predictions are described in Sections 2.5 (methods) and 3.2 (results). These two improvements are separate cases.

2. Material and Methods

2.1. Test Site and Field Data

The boreal managed forest study area, 5 × 5 km in size, is situated in Evo, Southern Finland. A total of 5,532 trees from 292 plots, each with a radius of 10 m, were used. The plots were measured in 2009. The average stand size in the study area is slightly less than 1 ha. The terrain elevation varies between 125 m and 185 m above sea level. Scots Pine (Pinus sylvestris) and Norway Spruce (Picea abies) are the dominant tree species, contributing respectively 40% and 35% of the total volume. The percentage of deciduous trees is 25% of the total stem volume. The sampling of the field plots in our study was based on pre-stratification of the existing stand inventory data. All the trees in each plot having a DBH greater than 5 cm were recorded and tree height, DBH, the lower limit of the living crown, crown width, and tree species were determined. The plot-level data, summarized in Table 1, were obtained by averaging or summing the tree data. The tree locations were calculated using the coordinates of the plot centers and using the direction and distance of the trees relative to these centers. The plot centers were measured using a Trimble GEOXM 2005 Global Positioning System (GPS) device (Trimble Navigation Ltd., Sunnyvale, CA, USA), and the locations were post-processed with local base-station data. This resulted in an average error of approximately 0.6 m. The tree heights were measured using Vertex clinometers. The DBH was measured using steel calipers.

2.2. Airborne Laser Data Collection

The ALS data were collected in the summer of 2009 using a Leica ALS50-II system operating at a pulse rate of 150 kHz. The data were acquired using a flight altitude of 400 m with FOV of 30 degrees. ALS data comprised of 12.7 emitted laser pulses per m2, including overlaps of strips. A footprint of 8.8 cm in diameter, collected by a beam having a divergence of 0.22 mrad, was achieved. The system was configured to record multiple returns per pulse, i.e., the first of many returns, the last of many returns, an only return and intermediate returns.

2.3. Laser Data Pre-Processing

The ground points were classified using TerraScan and based on the method explained in [72]. Strip adjustment was performed using TerraMatch. A DTM was then calculated using the ground points. Canopy heights were calculated by subtracting the corresponding ground elevation from the elevations of the laser returns.

2.4. Improved Tree Detection

Most of the current approaches for tree detection are based on finding trees from the CHM, which is calculated as a maximum of canopy height values within each raster cell. Thus, the CHM corresponds to the maximum canopy height of the first pulse data. Our approach is based on using the canopy penetration capability of the last pulse returns with overlapping trees. When trees overlap, the surface model corresponding to the first pulse stays high, whereas with last pulse, even a small gap results in a drop in elevation, i.e., the trees can be more readily discriminated. Use of first pulse works, when the whole laser beam penetrates the gap between the crowns, so that the drop is detectable after filtering. With last pulse, the drop in elevation is substantially larger and the drop can be detected even with overlapping trees since the last pulse is more sensitive to lower canopy levels. Since the objective was to demonstrate the usefulness of last pulse data for tree detection using the above-mentioned phenomenon, we used the same tree detection method and verification for both the first and last pulse data. From the last pulse data, we calculate three different products corresponding to the minimum of last returns, the mean of last returns and the maximum of last returns. The tree-detection algorithm in this test was based on finding the local maximum. Finding the local maximum is computationally cost-effective, and this is needed in operational inventories where data processing currently costs about €0.20/ha in Finland. Based on [68], finding the local maximum was one of the most useful methods based on individual tree extraction in the EuroSDR/ISPRS Tree Extraction test.
The comparison consisted of the following steps:
  • Different raster models with 0.5 m × 0.5 m pixel size were created for tree location. The created models were as follows: minimum of last returns (Lmin), maximum of last returns (Lmax), mean of last returns (Lmean), and first returns maximum (Fmax).
  • The raster models were smoothed by means of a Gaussian filter. A 3 m × 3 m window size was selected in order to eliminate minor tree level fluctuations and to avoid the merging of overlapping trees. Given Finland’s forest conditions, larger window sizes lead to the merging of overlapping trees, especially when trying to locate trees from within the suppressed tree storey.
  • Local maxima were sought from the smoothed surface model in a 3 m × 3 m window, and trees were considered to have been detected if the local maxima were greater than 2 m above the ground.
The extracted tree locations were compared with the tree locations measured in the field. The laser-detected trees were automatically matched with trees measured in the field, based on the Hausdorff distance on XY plane. The matching technique is described in detail in [73]. Tree pairs with a distance from each other of less than 2.5 m were considered to be a correct match. The value of 2.5 m was set by taking the distance between the treetop and the roots into consideration.
The methodology and the applied automatic accuracy assessment are further demonstrated in Figure 1, which shows four raster models for one sample plot with the detected trees marked as ‘+’ and field-measured trees as ‘o’. The Fmax surface model corresponds to the commonly accepted way of finding trees. When comparing Lmin and Fmax models, it seems to be easier to discriminate trees from Lmin rather than from Fmax, even though visual processing of laser data in the forest is inferior to the best automated techniques [68].

2.5. Using Point Height Metrics and Individual Tree-Based Features in Area-Based Predictions

In order to demonstrate the usefulness of individual tree features in area-based predictions, which was the second objective of this paper, individual trees were located from the CHM (Fmax model), based on a method employing the minimum curvature object detection [66]. This method was among the best in an extended analysis of the EuroSDR/ISPRS Tree Extraction test [68]. Since the two objectives of this paper were processed in parallel, it was not possible to improve the minimum curvature object detection method with the last pulse data. The DBH was modelled using the tree height and crown area of each tree as the predictors [75]. Tree volume was calculated based on Laasasenaho’s volume equations [74] with laser-derived tree heights (the values of an unfiltered CHM at the tree locations) and DBH as the inputs. Thus, no local data were used to calibrate this conversion.
Plotwise mean height and the mean DBH value were obtained from the arithmetic mean values of the extracted individual tree heights and DBH. The volumes were obtained by summing up the individual tree volumes calculated from Laasasenaho’s equations [74]. Both individual tree-based and point height metrics (Table 2) were used as the inputs for the Random Forest (RF) classifier [66]. RF is a non-parametric regression method in which the prediction is obtained by aggregating regression trees, each constructed using a different random sample of the training data, and choosing splits of the trees from among the subsets of the available features, randomly chosen at each node. The samples that are not used in training are called “out-of-bag” observations. They can be used to estimate the feature’s importance (as applied in Section 3.2) by randomly permutating out-of-bag data across one feature at a time and then estimating the increment in error due to this permutation. The greater the increment, the more important the feature ([66,76]).

3. Results and Discussion

3.1. Tree Detection Accuracy

Figure 2 depicts the percentage of correctly-matched trees using four different raster models for tree location. As mentioned earlier, the conventional way of detecting trees is to use the Fmax surface model (i.e., the CHM). The use of last pulse data gave a higher degree of discrimination between the trees than the use of first pulse data. The use of the raster corresponding to the minimum of last returns resulted in the highest discrimination between trees. An improvement of 6% in individual tree detection is better than that obtained by increasing the pulse density from 2 to 8 pulses per m2 reported in [34]. Figure 3 shows that the improvement in tree detection increases when the density of the forest stand increases. Obviously, the number of commission errors (the percentage of the number of falsely detected trees that would not be matched with field-measured trees) increases, as shown in Figure 4. The method also reveals many small saplings, which are not even recorded in field surveys because of their small diameter (>5 cm DBH is required to qualify for being recorded). For this reason, in future commission errors should be analyzed in detail, taking this into account. The proposed method is prone to creating real commission errors when there are gaps within individual tree crowns. Thus, the filtering has to be specified for the forest type. Figure 5 shows the percentage of correctly matched trees as a function of DBH. The smaller the DBH class, the better the performance of the proposed method (Lmin raster) compared to the use of the CHM (Fmax). With the DBH class 5–10 cm, the last pulse resulted in 10% better detection of trees.
The results confirm that there is also substantial information for tree detection in last pulse data. Currently, in raster-based processing, this information has been largely neglected. The obtained results would even suggest the use of last pulse data for detection, but we assume that a hybrid model utilizing both the first and last pulse data should be developed, even when processing is done at raster level. The advantages of first pulse data obviously include the lower number of commission errors and the high quality of tree separation when the crowns are not overlapping, whereas the advantage of last pulse is in the separation of trees whose crowns overlap. A hybrid model, utilizing the advantages of both pulse types should be developed.
Previously, last pulse data has been demonstrated to be usable for various applications. Liang et al. [77] used a simple technique (the difference of the first and last pulse returns under leaf-off conditions) to discriminate between deciduous and coniferous trees at individual tree level. In Matikainen et al. [78], a comparison between first pulse and last pulse laser scanner data in building detection was carried out. Both visual and numerical quality evaluations showed that the correctness of the results improved when last pulse data were used instead of first pulse data. In building detection using leaf-off data, the clearest difference between the first pulse and last pulse results was the proportion of the area classified as trees, which was considerably larger when the first pulse DSM was used [78].
With leaf-off data, more complex decision rules have to be developed for tree detection. With leaf-off data, the coniferous trees behave similarly to the leaf-on data, but the response of the deciduous trees can vary and needs to be studied in detail. A hybrid technique, utilizing both the first and last pulse data, may provide a working solution for deciduous canopies.
The applied 2.5 m maximum distance for tree matching is a possible error source when matching small trees. A variable distance, based on tree height, could be used in future. Matching using variable distance based on the DBH of the tree is reported in [79].

3.2. The Accuracy of Area-Based Prediction of Forest Variables

We wanted to show that the use of individual tree features as predictors increases the accuracy of stem volume in area-based predictions. Table 3 summarizes the results for mean height, mean DBH and stem volume when using (a) all features (both individual tree-based and point height metrics); (b) point height metrics and (c) individual tree-based features. With regard to volume estimation, there is a substantial improvement when using individual tree-based features as input for area-based prediction when compared to point height metrics. The results were repeated using two point densities and two DBH versus height models, and the combined data set, consisting of individual tree-based and point height metrics, provided the most accurate forest variable predictions in all cases. Since individual tree-based features can easily be added to the ABA method, the results propose a potential improvement in the accuracy of the estimation of forest variables in the ABA.
Individual tree-based features improved the ABA’s accuracy since they had very high correlation, e.g., with the reference stem volume. When calculating the importance of the features, most of the individual tree-based features were among the best features. Figure 6 shows that, when estimating mean height, the best laser-derived feature was the mean height derived by using the individual tree technique. When estimating DBH, the best laser-derived features were (1) mean canopy height and (2) penetration to the ground (as used in [4–6,80,81] and which was confirmed to correlate powerfully with standing volume and biomass in [82]); (3) mean tree height (derived from the extracted individual trees) and (4) mid percentiles. For the estimation of stem volume, the best laser-derived feature was the stem volume derived from extracted individual trees, followed by the basal area derived from extracted individual trees. It is possible to easily derive further laser point height metrics and individual tree-based features. The selected features are those reported quite early on, in [11–13]. Especially for the estimation of DBH, more features should be extracted to improve accuracy. In the field of statistics it is well-known that the use of predictors with high powers of explanation yields better estimates.
Another matter of concern is pulse density. It is probable that individual tree-based features lose some of their explanatory powers when applied at lower pulse densities, and with point densities of about 1 point/m2, the improvement is assumed to be more modest. However, in [34], pulse density did not have a significant impact on the accuracy of individual tree-based results
Several ALS inventory studies have been carried out in the same Evo area as the present study. In [83], the ABA stem volume estimation RMSE (relative stem volume estimation error, Root Mean Squared Error) was 27.1% using 282 field plots for training the k-NN method. In [67], the obtained RMSE for stem volume ranged from 24.8% to 27.2%, which is comparable to the 25.2% obtained in the present study with point height metrics as predictors. An improvement of RMSE value from about 25% to 20.3–20.4% can be considered to be substantial.
Future studies should test the sensitivity of individual tree algorithms in tree finding and the sensitivity of point density for the estimation of forest variables. Since the overall estimation methods based on RF are robust, our preliminary understanding is that estimation is not as sensitive to the tree finding algorithm as has previously been reported in ITD literature. It should be borne in mind that the applied tree extraction method, namely the minimum curvature objector detection method [66], was among the best in the extented EuroSDR/ISPRS comparison [68]. Yet, even that method can be assumed to be improved upon by incorporating last returns, as explained in this paper. We believe that the use of physically-strong features [84], together with statistically proven non-parametric estimation techniques, can also result in savings in practical forest inventory, even when using lower point density data (<1 emitted laser pulses per m2).

4. Conclusions

This paper reports two improvements for laser-based forest inventory. The first improvement is based on using the last returns for discriminating between overlapping trees. Using last pulse data and the test, which included 292 forests plots and more than 5,500 trees, an improvement of 6% for individual tree detection was obtained when compared to using first pulse data. In the 5–10 cm diameter breast height class, the use of last pulse resulted in a 10% better detection of trees than when using first pulse data. The results confirm that there is substantial information for tree detection in last pulse data, which should not be neglected even when using raster-based processing.
The second improvement is based on the use of individual tree-based features, in addition to statistical point height metrics, in the area-based prediction of forest variables. By using individual tree-based features as the input in non-parametric estimation, the root mean squared error, when compared to solely point height metrics, were reduced from about 25% to 20% at plot level. Point height metrics and individual tree-based features complemented each other in basal area estimation. The results confirmed the high usability of individual tree level features in the area-based estimation of forest variables.

Acknowledgments

The Academy of Finland is acknowledged for its financial support in the form of the projects “Science and Technology Towards Precision Forestry” and “Towards Improved Characterization of Map Objects”.

References

  1. Sexton, J.O.; Bax, T.; Siqueira, P.; Swenson, J.J.; Hensley, S. A comparison of lidar, radar, and field measurements of canopy height in pine and hardwood forests of south-eastern North America. For. Ecol. Manage 2009, 257, 1136–1147. [Google Scholar]
  2. Houghton, R.A.; Hall, F.; Goetz, S.J. Importance of biomass in the global carbon cycle. J. Geophys. Res 2009, 114. [Google Scholar] [CrossRef]
  3. Hyyppä, J.; Hyyppä, H.; Leckie, D.; Gougeon, F.; Yu, X.; Maltamo, M. Review of methods of small-footprint airborne laser scanning for extracting forest inventory data in boreal forests. Int. J. Remote Sens 2008, 29, 1339–1366. [Google Scholar]
  4. Maclean, G.; Krabill, W. Gross-merchantable timber volume estimation using and airborne LiDAR system. Can. J. Remote Sens 1986, 12, 7–18. [Google Scholar]
  5. Nelson, R.; Krabill, W.; Maclean, G. Determining forest canopy characteristics using airborne laser data. Remote Sens. Environ 1984, 15, 201–212. [Google Scholar]
  6. Nelson, R.; Krabill, W.; Tonelli, J. Estimating forest biomass and volume using airborne laser data. Remote Sens. Environ 1988, 24, 247–267. [Google Scholar]
  7. Næsset, E. Determination of mean tree height of forest stands using airborne laser scanner data. ISPRS J. Photogramm 1997, 52, 49–56. [Google Scholar]
  8. Means, J.; Acker, S.; Harding, D.; Blair, J.; Lefsky, M.; Cohen, W.; Harmon, W.; McKee, A. Use of large-footprint scanning airborne lidar to estimate forest stand characteristics in the western cascades of oregon. Remote Sens. Environ 1999, 67, 298–308. [Google Scholar]
  9. Lim, K.; Treitz, P.; Baldwin, K.; Morrison, I.; Green, J. Lidar remote sensing of biophysical properties of northern tolerant hardwood forests. Can. J. Remote Sens 2002, 29, 658–678. [Google Scholar]
  10. Lim, K.; Treitz, P.; Wulder, M.; St.-Onge, B.; Flood, M. LIDAR remote sensing of forest structure. Prog. Phys. Geogr 2003, 27, 88–106. [Google Scholar]
  11. 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]
  12. Maltamo, M.; Eerikäinen, K.; Packalen, P.; Hyyppä, J. Estimation of stem volume using laser scanning-based canopy height metrics. Forestry 2006, 79, 217–229. [Google Scholar]
  13. Hyyppä, J.; Inkinen, M. Detecting and estimating attributes for single trees using laser scanner. Photogramm. J. Fin 1999, 16, 27–42. [Google Scholar]
  14. Brandtberg, T. Automatic individual tree-based analysis of high spatial resolution remotely sensed data. Acta Univ. Agric. Sueciae Silvestria 1999, 118, 47. [Google Scholar]
  15. Ziegler, M.; Konrad, H.; Hofrichter, J.; Wimmer, A.; Ruppert, G.; Schardt, M.; Hyyppä, J. Assessment of forest attributes and single-tree segmentation by means of laser scanning. Laser Radar Technol. Appl. V 2000, 4035, 73–84. [Google Scholar]
  16. Hyyppä, J.; Kelle, O.; Lehikoinen, M.; Inkinen, M. A segmentation-based method to retrieve stem volume estimates from 3-dimensional tree height models produced by laser scanner. IEEE Trans. Geosci. Remote Sens 2001, 39, 969–975. [Google Scholar]
  17. Hyyppä, J.; Schardt, M.; Haggrén, H.; Koch, B.; Lohr, U.; Scherrer, H.U.; Paananen, R.; Luukkonen, H.; Ziegler, M.; Hyyppä, H.; et al. HIGH-SCAN: The first European-wide attempt to derive single-tree information from laserscanner data. Photogramm. J. Fin 2001, 17, 58–68. [Google Scholar]
  18. Hyyppä, J. Method for Determination of Stand Attributes and a Computer Program to. Perform the Method. US Patent 6792684 B1 2004. [Google Scholar]
  19. Persson, Å.; Holmgren, J.; Söderman, U. Detecting and measuring individual trees using an airborne laser scanner. Photogramm. Eng. Remote Sensing 2002, 68, 925–932. [Google Scholar]
  20. Gougeon, F.; Moore, T. Classification Individuelle des Arbres à Partir d’Images à Haute Résolution Spatiale. Proceedings of 6e Congrès de L'association Québécoise de Télédétection, Sherbrooke, QC, Canada, 4–6 May 1988; VI, pp. 185–196.
  21. Erikson, M.; Olofsson, K. Comparison of three individual tree crown detection methods. Mach. Vis. Appl 2005, 16, 258–265. [Google Scholar]
  22. Pitkänen, J.; Maltamo, M.; Hyyppä, J.; Yu, X. Adaptive Methods for Individual Tree Detection on Airborne Laser Based Canopy Height Model. Proceedings of ISPRS WG VIII/2 Laser-Scanners for Forest and Landscape Assessment, Freiburg, Germany, 3–6 October 2004; 36, pp. 187–191.
  23. Bortolot, Z.; Wynne, R. Estimating forest biomass using small footprint LiDAR data: An individual tree-based approach that incorporates training data. ISPRS J. Photogramm 2005, 59, 342–360. [Google Scholar]
  24. Hyyppä, J.; Mielonen, T.; Hyyppä, H.; Maltamo, M.; Yu, X.; Honkavaara, E.; Kaartinen, H. Using Individual Tree Crown Approach for Forest Volume Extraction with Aerial Images and Laser Point Clouds. Proceedings of ISPRS WG III/3, III/4, V/3 Workshop “Laser scanning 2005”, Enschede, The Netherlands, 12–14 September 2005; 36, pp. 144–149.
  25. Suarez, J.C.; Ontiveros, C.; Smith, S.; Snape, S. Use of airborne LiDAR and aerial photography in the estimation of individual tree heights in forestry. Comput. Geosci 2005, 31, 253–262. [Google Scholar]
  26. Hyyppä, J.; Yu, X.; Hyyppä, H.; Maltamo, M. Methods of Airborne Laser Scanning for Forest Information Extraction. Proceedings of the Workshop on 3D Remote Sensing in Forestry, Vienna, Austria, 14–15 February 2006; p. 16.
  27. Koch, B.; Heyder, U.; Weinacker, H. Detection of individual tree crowns in airborne lidar data. Photogramm. Eng. Remote Sensing 2006, 72, 357–363. [Google Scholar]
  28. Korpela, I.; Tokola, T. Potential of aerial image-based monoscopic and multiview single-tree forest inventory: A simulation approach. For. Sci 2006, 52, 136–147. [Google Scholar]
  29. Barilotti, A.; Sepic, F.; Abramo, E.; Crosilla, F. Assessing the 3D Structure of the Single Crowns in Mixed Alpine Forests. Proceedings of Photogrammetric Image Analysis PIA’07, Munich, Germany, 19–21 September 2007; XXXVI, pp. 129–134.
  30. Hirschmugl, M.; Ofner, M.; Raggam, J.; Schardt, M. Single tree detection in very high resolution remote sensing data. Remote Sens. Environ 2007, 110, 533–544. [Google Scholar]
  31. Popescu, S.C. Estimating biomass of individual pine trees using airborne lidar. Biomass Bioenergy 2007, 31, 646–655. [Google Scholar]
  32. Kwak, D.-A.; Lee, W.-K.; Lee, J.-H.; Biging, G.S.; Gong, P. Detection of individual trees and estimation of tree height using LiDAR data. J. For. Res 2007, 12, 425–434. [Google Scholar]
  33. Heurich, M. Automatic recognition and measurement of single trees based on data from airborne laser scanning over the richly structured natural forests of the Bavarian Forest National Park. For. Ecol. Manage 2008, 255, 2416–2433. [Google Scholar]
  34. Kaartinen, H.; Hyyppä, J. EuroSDR/ISPRS Project, Commission II, “Tree Extraction”, Final Report; European Spatial Data Research Official Publication: Amsterdam, The Netherlands, 2008; No. 53. [Google Scholar]
  35. Villikka, M.; Maltamo, M.; Packalén, P.; Vehmas, M.; Hyyppä, J. Alternatives for predicting tree-stem volume of Norway Spruce using airborne laser scanning. Photogramm. J. Fin 2008, 20, 33–42. [Google Scholar]
  36. Wang, Y.; Weinacker, H.; Koch, B. A lidar point cloud based procedure for vertical canopy structure analysis and 3D single tree modelling in forest. Sensors 2008, 8, 3938–3951. [Google Scholar]
  37. Balic, N.; Koch, B. Canscan—An algorithm for automatic extraction of canyons. Remote Sens 2009, 1, 197–209. [Google Scholar]
  38. Forzieri, G.; Guarnieri, L.; Vivoni, E.R.; Castelli, F.; Preti, F. Multiple attribute decision making for individual tree detection using high-resolution laser scanning. For. Ecol. Manage 2009, 258, 2501–2510. [Google Scholar]
  39. Hudak, A.; Evans, J.; Smith, A. LiDAR utility for natural resource managers. Remote Sens 2009, 1, 934–951. [Google Scholar]
  40. Kato, A.; Moskal, L.M.; Schiess, P.; Swanson, M.E.; Calhoun, D.; Stuetzle, W. Capturing tree crown formation through implicit surface reconstruction using airborne lidar data. Remote Sens. Environ 2009, 113, 1148–1162. [Google Scholar]
  41. Maltamo, M.; Peuhkurinen, J.; Malinen, J.; Vauhkonen, J.; Packalén, P.; Tokola, T. Predicting tree attributes and quality characteristics of scots pine using airborne laser scanning data. Silva Fenn 2009, 43, 507–521. [Google Scholar]
  42. Reitberger, J.; Schnörr, C.; Krzystek, P.; Stilla, U. 3D segmentation of single trees exploiting full waveform LIDAR data. ISPRS J. Photogramm 2009, 64, 561–574. [Google Scholar]
  43. Angelo, J.; Duncan, B.; Weishampel, J. Using lidar-derived vegetation profiles to predict time since fire in an Oak scrub landscape in East-Central Florida. Remote Sens 2010, 2, 514–525. [Google Scholar]
  44. Breidenbach, J.; Næsset, E.; Lien, V.; Gobakken, T.; Solberg, S. Prediction of species specific forest inventory attributes using a nonparametric semi-individual tree crown approach based on fused airborne laser scanning and multispectral data. Remote Sens. Environ 2010, 114, 911–924. [Google Scholar]
  45. Holopainen, M.; Mäkinen, A.; Rasinmäki, J.; Hyyppä, J.; Hyyppä, H.; Kaartinen, H.; Viitala, R.; Vastaranta, M.; Kangas, A. Effect of tree level airborne laser scanning accuracy on the timing and expected value of harvest decisions. Eur. J. For. Res 2010, 129, 899–910. [Google Scholar]
  46. Holopainen, M.; Vastaranta, M.; Rasinmäki, J.; Kalliovirta, J.; Mäkinen, A.; Haapanen, R.; Melkas, T.; Yu, X.; Hyyppä, J. Uncertainty in timber assortment estimates predicted from forest inventory data. Eur. J. For. Res 2010, 129, 1131–1142. [Google Scholar]
  47. Holopainen, M.; Vastaranta, M.; Haapanen, R.; Yu, X.; Hyyppä, J.; Kaartinen, H.; Viitala, R.; Hyyppä, H. Site-type estimation using airborne laser scanning and stand register data. Photogramm. J. Fin 2010, 22, 16–32. [Google Scholar]
  48. Gupta, S.; Weinacker, H.; Koch, B. Comparative analysis of clustering based approaches for 3-D single tree detection using airborne fullwave lidar data. Remote Sens 2010, 2, 968–989. [Google Scholar]
  49. Kantola, T.; Vastaranta, M.; Yu, X.; Lyytikäinen-Saarenmaa, P.; Holopainen, M.; Talvitie, M.; Kaasalainen, S.; Solberg, S.; Hyyppä, J. Classification of defoliated trees using tree-level airborne laser scanning data combined with aerial images. Remote Sens 2010, 2, 2665–2679. [Google Scholar]
  50. Koch, B. Status and future of laser scanning, synthetic aperture radar and hyperspectral remote sensing data for forest biomass assessment. ISPRS J. Photogramm 2010, 65, 581–590. [Google Scholar]
  51. Korpela, I.; Orka, H.; Hyyppa, J.; Heikkinen, V.; Tokola, T. Range and AGC normalization in airborne discrete-return LiDAR intensity data for forest canopies. ISPRS J. Photogramm 2010, 65, 369–379. [Google Scholar]
  52. Korpela, I.; Orka, H.; Maltamo, M.; Tokola, T.; Hyyppä, J. Tree species classification using airborne lidar-effects of stand and tree parameters, downsizing of training set, intensity normalization, and sensor type. Silva Fenn 2010, 44, 319–339. [Google Scholar]
  53. Kwak, D.-A.; Lee, W.-K.; Cho, H.-K.; Lee, S.-H.; Kafatos, M.; Kim, S.-R. Estimating stem volume and biomass of Pinus koraiensis using LiDAR data. J. Plant Res 2010, 123, 421–432. [Google Scholar]
  54. Puttonen, E.; Litkey, P.; Hyyppä, J. Individual tree species classification by illuminated—Shaded area separation. Remote Sens 2010, 2, 19–35. [Google Scholar]
  55. Vauhkonen, J.; Korpela, I.; Maltamo, M.; Tokola, T. Imputation of single-tree attributes using airborne laser scanning-based height, intensity, and alpha shape metrics. Remote Sens. Environ 2010, 114, 1263–1276. [Google Scholar]
  56. Yu, X.; Hyyppä, J.; Holopainen, M.; Vastaranta, M. Comparison of area-based and individual tree-based methods for predicting plot-level forest attributes. Remote Sens 2010, 2, 1481–1495. [Google Scholar]
  57. Edson, C.; Wing, M. Airborne light detection and ranging (LiDAR) for individual tree stem location, height, and biomass measurements. Remote Sens 2011, 3, 2494–2528. [Google Scholar]
  58. Heinzel, J.; Weinacker, H.; Koch, B. Prior-knowledge-based single-tree extraction. Int. J. Remote Sens 2011, 32, 4999–5020. [Google Scholar]
  59. Jung, S.-E.; Kwak, D.-A.; Park, T.; Lee, W.-K.; Yoo, S. Estimating crown variables of individual trees using airborne and terrestrial laser scanners. Remote Sens 2011, 3, 2346–2363. [Google Scholar]
  60. Listopad, C.; Drake, J.; Masters, R.; Weishampel, J. Portable and airborne small footprint LiDAR: Forest canopy structure estimation of fire managed plots. Remote Sens 2011, 3, 1284–1307. [Google Scholar]
  61. Tuominen, S.; Haapanen, R. Comparison of grid-based and segment-based estimation of forest attributes using airborne laser scanning and digital aerial imagery. Remote Sens 2011, 3, 945–961. [Google Scholar]
  62. Ussyshkin, V.; Theriault, L. Airborne lidar: Advances in discrete return technology for 3D vegetation mapping. Remote Sens 2011, 3, 416–434. [Google Scholar]
  63. Vauhkonen, J.; Ene, L.; Gupta, S.; Heinzel, J.; Holmgren, J.; Pitkänen, J.; Solberg, S.; Wang, Y.; Weinacker, H.; Hauglin, K.; et al. Comparative testing of single-tree detection algorithms under different types of forest. Forestry 2011. [Google Scholar] [CrossRef]
  64. Vastaranta, M.; Holopainen, M.; Yu, X.; Hyyppä, J.; Mäkinen, A.; Rasinmäki, J.; Melkas, T.; Kaartinen, H.; Hyyppä, H. Effects of individual tree detection error sources on forest management planning calculations. Remote Sens 2011, 3, 1614–1626. [Google Scholar]
  65. Vehmas, M.; Eerikäinen, K.; Peuhkurinen, J.; Packalén, P.; Maltamo, M. Airborne laser scanning for the site type identification of mature boreal forest stands. Remote Sens 2011, 3, 100–116. [Google Scholar]
  66. Yu, X.; Hyyppä, J.; Vastaranta, M.; Holopainen, M.; Viitala, R. Predicting individual tree attributes from airborne laser point clouds based on the random forests technique. ISPRS J. Photogramm 2011, 66, 28–37. [Google Scholar]
  67. Vastaranta, M.; Kankare, V.; Holopainen, M.; Yu, X.; Hyyppä, J.; Hyyppä, H. Combination of individual tree detection and area-based approach in imputation of forest variables using airborne laser data. ISPRS J. Photogramm 2012, 67, 73–79. [Google Scholar]
  68. Kaartinen, H.; Hyyppä, J.; Yu, X.; Vastaranta, M.; Hyyppä, H.; Kukko, A.; Holopainen, M.; Heipke, C.; Hirschmugl, M.; Morsdorf, F.; et al. An international comparison of individual tree detection and extraction using airborne laser scanning. Remote Sens 2012, 4, 950–974. [Google Scholar]
  69. Wang, Y.; Weinacker, H.; Koch, B.; Sterenczak, K. Lidar Point Cloud Based Fully Automatic 3D Single Tree Modelling in Forest and Evaluations of the Procedure. Proceedings of the 21st ISPRS Congress, Beijing, China, 3–11 July 2008; 37, pp. 45–51.
  70. Gupta, S.; Koch, B.; Weinacker, H. Tree Species Detection Using Full Waveform Lidar Data in a Complex Forest. In Technical Commission VII Symposium—100 Years of ISPRS; Wagner, W., Székely, B., Eds.; Vienna University of Technology: Vienna, Austria, 2010; pp. 249–254. [Google Scholar]
  71. Holmgren, J.; Wallerman, J. Estimation of Tree Size Distributions by Combining Vertical and Horizontal Distribution of Laser Measurements with Extraction of Individual Trees. Proceedings of the Workshop on 3D Remote Sensing in Forestry, Vienna, Austria, 14–15 February 2006; pp. 168–173.
  72. Axelsson, P. DEM Generation From Laser Scanner Data Using Adaptive TIN Models. Proceedings of XIXth ISPRS Congress, Amsterdam, The Netherlands, 16–22 July 2000; 33, pp. 110–117.
  73. Yu, X.; Hyyppä, J.; Kukko, A.; Maltamo, M.; Kaartinen, H. Change detection techniques for canopy height growth measurements using airborne laser scanning data. Photogramm. Eng. Remote Sensing 2006, 72, 1339–1348. [Google Scholar]
  74. Laasasenaho, J. Taper Curve and Volume Functions for Pine, Spruce and Birch; The Finnish Forest Research Institute: Helsinki, Finland, 1982; Volume 108, pp. 1–74. [Google Scholar]
  75. Kalliovirta, J.; Tokola, T. Functions for estimating stem diameter and tree age using tree height, crown width and existing stand database information. Silva Fenn 2005, 39, 227–248. [Google Scholar]
  76. Breiman, L. Random forests. Mach. Learn 2001, 45, 5–32. [Google Scholar]
  77. Liang, X.; Hyyppä, J.; Matikainen, L. Deciduous-Coniferous Tree Classification Using Difference between First and Last Pulse Laser Signatures. ISPRS Workshop on Laser Scanning 2007 and SilviLaser 2007, Espoo, Finland, 12–14 September 2007; XXXVI, pp. 253–257.
  78. Matikainen, L.; Hyyppä, J.; Kaartinen, K. Comparison between first pulse and last pulse laser scanner data in the automatic detection of buildings. Photogramm. Eng. Remote Sensing 2009, 75, 133–146. [Google Scholar]
  79. Olofsson, K.; Lindberg, E.; Holmgren, J. A Method for Linking Field-Surveyed and Aerial-Detected Single Trees Using Cross Correlation of Position Images and the Optimization of Weighted Tree List Graphs. Proceedings of the SilviLaser 2008: 8th International Conference on LiDAR Applications in Forest Assessment and Inventory, Edinburgh, UK, 17–19 September 2008; pp. 95–104.
  80. Hyyppä, J.; Pulliainen, J.; Hallikainen, M.; Saatsi, A. Radar-derived standwise forest inventory. IEEE Trans. Geosci. Remote Sens 1997, 35, 392–404. [Google Scholar]
  81. Hyyppä, J.; Hallikainen, M. Applicability of airborne profiling radar to forest inventory. Remote Sens. Environ 1996, 57, 39–57. [Google Scholar]
  82. Ahokas, E.; Hyyppä, J.; Yu, X.; Holopainen, M. Transmittance of airborne laser scanning pulses for boreal forest elevation modeling. Remote Sens 2011, 3, 1365–1379. [Google Scholar]
  83. Holopainen, M.; Haapanen, R.; Tuominen, S.; Viitala, R. Performance of Airborne Laser Scanning- and Aerial Photograph-Based Statistical and Textural Features in Forest Variable Estimation. Proceedings of the SilviLaser 2008: 8th International Conference on LiDAR Applications in Forest Assessment and Inventory, Edinburgh, UK, 17–19 September 2008; pp. 105–112.
  84. Hollaus, M. Large Scale Applications of Airborne Laser Scanning for a Complex Mountainous Environment. 2006; p. 111. [Google Scholar]
Remotesensing 04 01190f1a 1024Remotesensing 04 01190f1b 1024
Figure 1. Four raster models for one example plot (radius 10 m) with the detected trees marked as ‘+’ and the field-measured trees marked as ‘o’. The trees designated by A, B, C and D were detectable on the Lmin image but not on the Fmax image. Lmin refers to the minimum of last returns, Lmax to the maximum of last returns, Lmean to the mean of last returns and Fmax to the first returns maximum. The return height is color coded.
Figure 1. Four raster models for one example plot (radius 10 m) with the detected trees marked as ‘+’ and the field-measured trees marked as ‘o’. The trees designated by A, B, C and D were detectable on the Lmin image but not on the Fmax image. Lmin refers to the minimum of last returns, Lmax to the maximum of last returns, Lmean to the mean of last returns and Fmax to the first returns maximum. The return height is color coded.
Remotesensing 04 01190f1aRemotesensing 04 01190f1b
Remotesensing 04 01190f2 1024
Figure 2. The percentage of correctly matched trees when tested with 5,532 trees surveyed in the field.
Figure 2. The percentage of correctly matched trees when tested with 5,532 trees surveyed in the field.
Remotesensing 04 01190f2
Remotesensing 04 01190f3 1024
Figure 3. The percentage of correctly matched trees as a function of plot density when tested with 5,532 trees surveyed in the field. The number of trees refers to a plot with a radius of 10 m.
Figure 3. The percentage of correctly matched trees as a function of plot density when tested with 5,532 trees surveyed in the field. The number of trees refers to a plot with a radius of 10 m.
Remotesensing 04 01190f3
Remotesensing 04 01190f4 1024
Figure 4. The percentage of commission errors for different surface models when using the local maximum finding as the tree detection algorithm.
Figure 4. The percentage of commission errors for different surface models when using the local maximum finding as the tree detection algorithm.
Remotesensing 04 01190f4
Remotesensing 04 01190f5 1024
Figure 5. The percentage of correctly matched trees as a function of diameter breast height (DBH).
Figure 5. The percentage of correctly matched trees as a function of diameter breast height (DBH).
Remotesensing 04 01190f5
Remotesensing 04 01190f6 1024
Figure 6. The feature importance (in the left column) and a scatter plot of the predicted versus observed values (in the right column) for mean height, mean DBH and volume estimates based on all features. The order of the feature index is the same as in Table 2 (24–27 are the individual tree-based features).
Figure 6. The feature importance (in the left column) and a scatter plot of the predicted versus observed values (in the right column) for mean height, mean DBH and volume estimates based on all features. The order of the feature index is the same as in Table 2 (24–27 are the individual tree-based features).
Remotesensing 04 01190f6
Table
Table 1. Descriptive statistics of the field plots accessed in the study.
Table 1. Descriptive statistics of the field plots accessed in the study.
Mean Height (m)Mean Diameter (cm)Stem Volume m3/ha
Minimum3.97.60.4
Maximum31.750.8586.2
Mean18.018.3148.2
Standard deviation6.16.9110.7
Table
Table 2. The features used in predicting the forest attributes.
Table 2. The features used in predicting the forest attributes.
No.FeatureExplanation
Point Height Metrics
1meanHMean canopy height calculated as the arithmetic mean of the heights from the point cloud
2stdHStandard deviations of heights from the point cloud
3PPenetration calculated as a proportion of ground returns to total returns
4COVCoefficient of variation
5H1010th percentile of canopy height distribution
6H2020th percentile of canopy height distribution
7H3030th percentile of canopy height distribution
8H4040th percentile of canopy height distribution
9H5050th percentile of canopy height distribution
10H6060th percentile of canopy height distribution
11H7070th percentile of canopy height distribution
12H8080th percentile of canopy height distribution
13H9090th percentile of canopy height distribution
14maxHMaximum height
15D1010th canopy cover percentile computed as the proportion of returns below 10% of the total height
16D2020th canopy cover percentile computed as the proportion of returns below 20% of the total height
17D3030th canopy cover percentile computed as the proportion of returns below 30% of the total height
18D4040th canopy cover percentile computed as the proportion of returns below 40% of the total height
19D5050th canopy cover percentile computed as the proportion of returns below 50% of the total height
20D6060th canopy cover percentile computed as the proportion of returns below 60% of the total height
21D7070th canopy cover percentile computed as the proportion of returns below 70% of the total height
22D8080th canopy cover percentile computed as the proportion of returns below 80% of the total height
23D9090th canopy cover percentile computed as the proportion of returns below 90% of the total height
Individual Tree Based Features
24LHMean height of all extracted trees
25LDMean DBH of all detected trees, derived from the extracted heights and crown areas
26LBBasal area of the plot, derived from the extracted DBH
27LVVolume of the plot, derived from the extracted DBH and height
Table
Table 3. The bias, RMSE, and correlation coefficient (R) between the predicted and observed values were calculated as measures of the accuracy of area-based inventory based on three different feature sets.
Table 3. The bias, RMSE, and correlation coefficient (R) between the predicted and observed values were calculated as measures of the accuracy of area-based inventory based on three different feature sets.
BiasRMSERMSE (%)R
With all features
Mean height (m)−0.001.106.150.98
Mean DBH (cm)0.002.9116.070.89
Volume (m3/ha)0.2430.0520.320.96
With point height metrics
Mean height (m)−0.031.256.990.98
Mean DBH (cm)0.023.0216.650.88
Volume (m3/ha)0.1337.5625.410.93
Individual tree-based features
Mean height (m)−0.001.246.970.98
Mean DBH (cm)−0.063. 5419.540.83
Volume (m3/ha)−1.0630.1620.400.96

Share and Cite

MDPI and ACS Style

Hyyppä, J.; Yu, X.; Hyyppä, H.; Vastaranta, M.; Holopainen, M.; Kukko, A.; Kaartinen, H.; Jaakkola, A.; Vaaja, M.; Koskinen, J.; et al. Advances in Forest Inventory Using Airborne Laser Scanning. Remote Sens. 2012, 4, 1190-1207. https://doi.org/10.3390/rs4051190

AMA Style

Hyyppä J, Yu X, Hyyppä H, Vastaranta M, Holopainen M, Kukko A, Kaartinen H, Jaakkola A, Vaaja M, Koskinen J, et al. Advances in Forest Inventory Using Airborne Laser Scanning. Remote Sensing. 2012; 4(5):1190-1207. https://doi.org/10.3390/rs4051190

Chicago/Turabian Style

Hyyppä, Juha, Xiaowei Yu, Hannu Hyyppä, Mikko Vastaranta, Markus Holopainen, Antero Kukko, Harri Kaartinen, Anttoni Jaakkola, Matti Vaaja, Jarkko Koskinen, and et al. 2012. "Advances in Forest Inventory Using Airborne Laser Scanning" Remote Sensing 4, no. 5: 1190-1207. https://doi.org/10.3390/rs4051190

APA Style

Hyyppä, J., Yu, X., Hyyppä, H., Vastaranta, M., Holopainen, M., Kukko, A., Kaartinen, H., Jaakkola, A., Vaaja, M., Koskinen, J., & Alho, P. (2012). Advances in Forest Inventory Using Airborne Laser Scanning. Remote Sensing, 4(5), 1190-1207. https://doi.org/10.3390/rs4051190

Article Metrics

Back to TopTop