## Integrating area-based and individual tree detection approaches for estimating tree volume in plantation inventory using aerial image and airborne laser scanning data

iForest - Biogeosciences and Forestry, Volume 10, Issue 1, Pages 296-302 (2016)
doi: https://doi.org/10.3832/ifor1880-009

Research Articles

Remote sensing has been increasingly used to assist forest inventory. Airborne Laser Scanning (ALS) systems can accurately estimate tree height in forests, and are being combined with more traditional optical images that provide further details about the horizontal structure of forests. To predict forest attributes two main techniques are applied to process ALS data: the Area Based Approach (ABA), and the Individual Tree Detection (ITD). The first part of this study was focused on the effectiveness of integrating ALS data and aerial imagery to estimate the wood volume in Eucalyptus urograndis plantations using the ABA approach. To this aim, we analyzed three different approaches: (1) using only ALS points cloud metrics (RMSE = 6.84%); (2) using only the variables derived from aerial images (RMSE = 8.45%); and (3) the integration of both 1 and 2 (RMSE = 5.23%), which underestimated the true volume by 2.98%. To estimate individual tree volumes we first detected individual trees and corrected the density estimate for detecting mean difference, with an error of 0.37 trees per hectare and RMSE of 12.68%. Next, we downscaled the total volume prediction to single tree level. Our approach showed a better result of the overall volume in comparison with the traditional forest inventory. There is a remarkable advantage in using the Individual Tree Detection approach, as it allows for a spatial representation of the number of trees sampled, as well as their volume per unit area - an important metric in the management of forest resources.

# Introduction

Intensively managed forest plantations have an important role, both nationally and internationally, due to their ability to efficiently produce large quantities of biomass. This biomass can be utilized in forms suitable for bioenergy or, after processing into cellulose-derivatives, in materials such as paper. The myriad uses of fibers derived from intensively managed plantation forests can aid in the reduction of anthropogenic pressures on natural forests, while simultaneously generating millions of jobs ([22]).

In Brazil, 4.4 million jobs are generated directly and indirectly by forest plantations, and the total reforested area exceeds 7.6 million hectares ([9]). Of this area, 72% is composed of Eucalyptus spp., a key genus in Brazilian forest plantations. The reasons for the success of Eucalyptus spp. plantations in Brazil include Eucalyptus’ rapid growth, genetic improvement, and its tolerance of a wide range of climatic conditions ([32]).

In order to manage the growth and yield of these plantations, annual forest inventories are required. In general, the traditional approaches to conducting plantation inventories follow the precepts of sampling theory, which is based on ground measurements collected from simple random sampling plots. The forest attributes of interest are then inferred for an entire stand based on observations collected at each sample plot ([21]). However, considering the difficulties in field surveys such as sampling error, bias, and the intensive use of time and labor, remote sensing represents an important additional tool for monitoring large areas accurately and efficiently.

Both active and passive remote sensors have been presented as feasible alternatives for estimating forest variables. High spatial resolution optical data is useful at providing spectral information on species and condition ([35], [18]), while active remote sensing technologies, such as Airborne Laser Scanning (ALS), can provide highly accurate forest height data ([37]). When predicting forest inventory attributes using ALS data or aerial images, two main approaches can be applied: the Area-based Approach (ABA - [19], [36]) and the Individual Tree Detection (ITD) approach ([27], [11]).

To estimate volume in Eucalyptus spp. plantations Packalén et al. ([20]) applied the ABA approach. In the ABA, the response variable (e.g., wood volume per unit area in m3 ha-1) is aggregated over a grid cell of predetermined size or over a plot size ([3]). This response variable is estimated using predictor variables, which are derived from the ALS metrics directly or indirectly, such as the mean height, height percentiles, densities and proportions, or from image, such as the reflectance ([19], [39]).

The ABA approach is most suited to the estimation of a number of crucial forest stand attributes over large areas ([19], [36]). The prediction of stand attribute variables from ALS data, such as biomass, volume, basal area, mean diameter, or average height have been shown to have similar or improved accuracy compared to traditional inventory methods ([36]). However, when tree-level information such as stem number or species is desired, the ABA approach is less suitable ([3]).

In this case, the ITD approach is the most appropriate to estimate individual tree characteristics. However, algorithms normally used in the ITD approach detect, on average, 75% of trees ([26], [18]). Therefore, the use of automatic ITD may lead to the non-detection of some trees, and can interfere with forest attributes prediction ([3]).

Vastaranta et al. ([34]) applied an alternative method to overcome the shortcomings of ITD detection, involving the aggregation of the tree numbers in a grid cell, according to ABA. By aggregating trees into grid cells it is possible to know how many trees were detected or not in each cell, so long as detection rates are compared with ground truth data and the most appropriate methodology to correct deviations is used. This approach, also used in this study, can lead to a reduction in the number of systematic errors, and to an increased accuracy in the estimation of the number of trees per hectare, a variable extremely important to forest managers ([33], [3]).

According to Packalén et al. ([20]), the advantage of using ITD approaches is that the correct detection of individual trees could reduce the number of plots needed for an accurate and unbiased inventory. Breindenbach & Astrup ([3]) also note that the ITD approach is inherently more intuitive than ABA because the response variables are related to each tree (the smallest unit in a management system), providing coordinates that can be tied to individual tree growth, as well as explicit links to tree growth and yield models. In the case of ALS data, ITD approaches require high density point clouds, which can be impractical because of the cost of data collection ([17]).

The merging of ALS and optical remote sensing data can assist in the accurate estimation of forest attributes, with algorithms based on height or brightness ([8], [19], [38], [35], [29], [18], [20]). Saarela et al. ([30]) integrated ALS and Landsat image data to improve accuracy of timber volume estimation as compared with estimates based on ALS and images variables alone. This comparison is also employed in the present study, though we used high spatial resolution images. For instance, the ABA and ITD approaches can be combined to estimate wood biomass or volume in forest plantations ([19], [2]).

In this study, we use ABA and ITD to estimate stand volume in Eucalyptus spp. plantations using ALS and high spatial resolution multispectral image data. To this goal, we first modeled the stand volume using standard ABA approaches, using different combinations of ALS-based and image-based metrics as independent model variables. Second, we used an ITD-based approach to estimate tree count, applying a correction method to improve single tree detection, which is integrated with the stand volume predicted using ABA. Finally, we merge the two approaches (ITD and ABA), that yields volume predictions at the single tree level. We conclude by discussing the possibilities of using ABA for individual tree volume predictions.

# Materials and methods

## Study area

The study was performed in a forest plantation owned by FIBRIA S.A. company, located in São Luiz do Paraitinga in São Paulo State, Brazil (23° 33′ S, 45° 34′ W and 23° 30′ S, 45° 31′ W - Fig. 1). The dominant species in the study area is Eucalyptus urograndis, which is a hybrid between Eucalyptus urophylla and Eucalyptus grandis ([7]). According to Köppen & Geiger ([12]), the climate is humid subtropical (Cwa) characterized by cold winters and rainy summers, with average annual precipitation of 1210 mm and a mean monthly temperature of 20.2 °C ([5]). The topography of the study area is complex, with elevations varying from 974 to 1266 m a.s.l. and the average slope is 24.3 degrees.

Fig. 1 - The study area near São Luíz do Paraitinga (São Paulo, Brazil), including the 6.5 year-old Eucalyptus urograndis forest plantation analyzed in this investigation.

The forest stands in the study area were planted in December 2006 and harvested between April and May 2013, following a standard 6.5 year rotation period. All harvested wood was taken to a mill. The trees were planted in rows with fixed tree spacing (2 × 3 m), resulting in a density of 1667 stems per hectare.

## Field measurements

Reference data were based on field measurements collected at 48 circular plots (r = 11.28 m, area of 400 m²) between April and May of 2013 prior to harvesting. In each plot, all diameter at breast height (DBH) were recorded, and the height of 10 trees were measured. Heights for all remaining trees in the plot were then estimated by applying the Curtis ([6]) hypsometric model, adjusted based on the 10 measured trees in each plot. Tree volume was estimated with locally-developed volume models based on the model proposed by Schöepfer ([31] - Tab. 1). The geographical location of the plots were obtained using high-precision GPS receivers (ProXRT®, Trimble Inc., Sunnyvale, CA, USA) at the center of each plot.

Tab. 1 - Descriptive statistics of the reference plots (n=48) for the basal area (G), diameter at breast height (DBH), total height (HT), dominant height (HDOM) and cellulose volume without bark (V).

Variable Min Max Mean SD
G (m2 ha-1) 27.08 39.16 33.96 2.75
DBH (cm) 13.80 17.34 16.00 0.65
HT (m) 20.82 27.93 24.89 1.70
HDOM (m) 25.20 33.27 29.79 2.06
V (m3 ha-1) 223.19 391.99 331.16 37.27

After harvesting, all wood arriving at the mill was weighed. The wood density was obtained according to Archimedes’ principle and, therefore, the total volume of wood at the mill could be defined by weight-to-volume (W-V) relationships obtained in the density measurement process. A complete description can be found in Batista ([1]).

## Remote sensing data

ALS data were acquired by a RIEGL LMS - Q680I® system (Riegl Gmbh, Horn, Austria) simultaneously on 15 January 2012 with aerial photographs under free-cloud conditions. The average ALS point density was 5 pulses per m2. An orthophoto was generated from aerial images with a pixel size of 0.15 cm. The 4 spectral bands were blue (429-514 nm), green (514-600 nm), red (600-676 nm) and near-infrared (695-831 nm).

Initially we separated ALS ground returns from the vegetation returns, to generate a 15 cm resolution DTM and DSM, which is the minimum size resolution based on point density ([15]). The DTM was used to normalize the point cloud, followed by the generation of the most usual metrics (such as tree height, canopy cover, topography, image bands and ratios) from each type of remote sensing data at a 5-meter resolution. This resolution was chosen following the ABA concept that pixel size should be larger than a tree crown ([36]), with Eucalyptus spp. crown diameter varying from 2 to 3 meters in the study area.

## Predictor variables

We extracted LiDAR metrics from the ALS data ([16], [36], [28], [30]) to extract image metrics (bands, ratios between bands and Normalized Difference Vegetation Index), as listed in Tab. 2. We generated the most common metrics (such as tree height, canopy cover, topography, image bands and ratios) from each type of remote sensing data by describing the properties of the point clouds and spectral properties of the images. For the generation of these metrics, only first returns of the ALS point cloud were used and the reflectance value of each pixel in the aerial image was aggregated into a 5 × 5 m grid cell. These metrics were then used in subsequent statistical modeling as explanatory variables. In total, 39 metrics were generated (Tab. 2).

Tab. 2 - List of the variables selected for model inclusion in each of the three approaches adopted in this study.

Data Category ALS metrics Acronym
ALS Height Minimum height hmin
Maximum height hmax
Mean height hmean
Variance of height hvar
Standard deviation of height hsd
Coefficient of variation of height hcv
Mode of height hmode
Kurtosis of height hk
Height percentiles: 1st, 5th, 10th, 20th, 25th, 30th, 40th, 50th, 60th, 70th, 75th, 80th, 90th, 95th, 99th hp01, hp05, hp10, hp20, hp25, hp30, hp40, hp50, hp60, hp70, hp75, hp80, hp90, hp95, hp99
Cover Number of first returns above mean Cab mean
Number of returns above 2 m Cab2m
Percentage of canopy C%
Topography Slope (°) S
Mean altitude A mean
Aerial Image Spectral Properties Band 1 - blue B1
Band 2 - green B2
Band 3 - red B3
Band 4 - near infrared B4
Simple ratio between bands B1/B2, B1/B3, B1/B4, B2/B3, B2/B4, B3/B4 B1/B2, B1/B3, B1/B4, B2/B3, B2/B4, B3/B4
Normalized Difference Vegetation Index NDVI

## Modeling approach to estimate total wood volume

The first step was the estimation of wood volume using the ABA. Three different methods were applied, depending on the input data source. In the first approach, only variables derived from ALS data were used, in the second approach only variables derived from the imagery were used, and in the final one we combined both sets of variables.

One of the most important considerations when working with a large number of predictor variables in a multiple linear regression is variable selection. To that end, we applied the Best Subset Selection (BSS), which matches an appropriate model to each model size ([10], [14], [28]). To select a model among those suggested by BSS, we applied a k-fold cross-validation (k=5). The application of a cross-validation, similar to the leave-one-out method, allows for a better understanding of model behavior, especially when using a database that was not utilized for model adjustment ([10]). After applying the cross-validation we chose the model according to the RMSE. In this study, we used a confidence interval with a p-value equal to 0.05 generated from 1000 simulations.

The total volume of the forest stands was obtained by the mean ABA estimation multiplied by the total area.

### Variable selection and model adjustment

The k-fold estimate of RMSE is computed by averaging the values obtained by eqn. 1:

$$CV_{k}={\frac{1 } {k }}\sum_{i=1}^{k} MSE_{i}$$

where CVk is the k-fold cross validation and MSE is the mean square error. The linear model assumptions were verified according to ([23]), where the assumptions were verified for each of the chosen models with statistical tests and graphical analysis ([24], [28]). The approach was then chosen based on the lowest RMSE and the linear model assumptions analysis, generating a predicted wood volume map for the entire study area. This prediction was validated with stand inventory data and information about tree volume received by the mills, after harvesting.

## Individual tree volume

### Individual Tree Detection based on aerial images

For the second step we used high resolution multispectral images for detecting individual treetops, applying the local maximum filtering method, where the central location of the crown is assigned the highest spectral value ([35], [13]). We applied fixed and varying sizes of moving window, based on the methodology proposed by Kumar ([13]) and Pereira ([25]). The window sizes varied from 0×3 m, 0×5 m, 0×7 m to 0×9 m ([13], [25]). The individual tree detection procedure was performed for each of the 4 spectral bands, and the first principal component (first PCA) of the bands ([35]). According to Wang et al. ([35]), the PCA effectively aggregates a large quantity of information and shows good results for ITD.

To generate the final ITD product, all five intermediate outputs were used. First, the layer with the highest number of detected trees, with their locations visually corrected, was taken as the basis for ITD. From this layer, we removed all the treetops in a 0.7 m buffer around each tree to avoid multiple peaks from single tree. Treetops detected in the second layer, outside this buffer were added creating a new layer of additional ITD points. This layer was then taken as the new basis for ITD. This procedure was repeated until all layers had contributed to the generation of a single ITD layer.

### Merging the ITD approach with ABA

Based on the final ITD layer, stem counts were aggregated into 5 × 5 m grid cells, according to ABA. Following this, the number of trees in each grid cell was identified by local maxima derived from the ITD approach, resulting in a graphical representation of the number of trees per hectare for the entire study area.

### Detection and correction of the mean difference (MD) in total area

The mean difference (MD) detection of the ITD approach was done by comparing the number of detected trees and the number of observed trees in the inventory plots. The difference between the two was considered the response variable, while the ALS and image variables were considered the predictor variables, leading to the modeling of the ITD mean difference, and its prediction for the entire study area. In our case this MD was modeled with 3 variables, maximum height, slope, and ITD detection, as explanatory variables. These variables were selected by applying the Best Subset Selection method and cross-validation, as described earlier. The MD detection model was applied to the entire study area, representing the MD per hectare, to correct the initial estimates of tree count per pixel.

### Estimating individual tree volume

In order to generate the volume per tree in each 5 × 5 m cell, the previously-estimated stand volume for the entire study area (m3 ha-1), obtained by the integration of ALS data and high-resolution imagery, was divided by the number of trees per hectare detected by the ITD approach.

# Results

## Modeling approach to estimate total wood volume

The best fitting models obtained by cross-validation contained 4 parameters, including the intercept, for all the three approaches tested (Tab. 3). In the first approach (using ALS data variables only), the chosen model was based on height percentiles (hp20, hp40, hp50) and elevation (Amean). In the second approach, where only image predictors were used, the selected variables were B3, NDVI, B1/B3 and B1/B4. In the third approach, where both remote sensing data sources were used as predictors, two ALS-based and two image-based variables were selected (hmode, Amean, B4, B2/B3). The RMSE for each approach was 6.84%, 8.45% and 5.23%, respectively. Among the three analyzed approach, we chose the model that combined ALS data and image information, according to the corresponding adjusted R² value (Tab. 3).

Tab. 3 - Fitting statistics of the best regression models used to estimate the total tree volume in the investigated area with three different approaches. (MD): Mean difference; (RMSE): root mean squared error; (Adj): Adjusted. For the acronym of variables included in the models, see Tab. 2.

Approach Regression Model Form Adj-R² MD MD % RMSE RMSE % p-value
(t-test)
(1) ALS data -585.17 + 22.48·(hp20) -165.27·(hp40) + 170.09·(hp50) + 0.22·(Amean) 0.65 0.032 0.001 22.79 6.84 0.993
(2) Aerial images -585.17 + 7.859·(B3) -16167.59·(NDVI) -5101.24·(B1/B3) + 15192.15·(B1/B4) 0.48 0.64 0.19 28.17 8.45 0.883
(3) Both datasets -625.16 + 27.74·(hmode) + 0.32·(Amean) + 1.97·(B4) -377.84·(B2/B3) 0.79 -0.06 -0.018 17.43 5.23 0.983

Fig. 2 displays the relationship between the observed and predicted values of stand volume obtained for each approach. Although the homoscedasticity of residuals was not rejected, their distribution was asymmetric in the first and second approach, suggesting non-normality, while in approach 3 the distribution of residuals showed to match a normal distribution. Therefore, the third approach was selected for wall-to-wall stand volume prediction and subsequent comparisons with the inventory and mill-based volumes as reference values.

Fig. 2 - Scatterplot of the observed and predicted values of stand volume for each of the approaches analyzed in this study. (1) ALS data; (2) aerial images; (3) both datasets.

Across the whole study area, the mean estimated volume per hectare based on the ABA method was 319.9 m3, totaling 44.466.1 m3. When this result was compared with the volume measured at the mill (considered the true volume by the company), we found a difference of 2.98%, while when comparing the forest inventory data with the mill volume, it overestimated by 3.29%. In addition, it was observed that the combination of both ALS and imagery data could provide more information regarding the occurrence of drought and canopy gaps in the forest.

## Semi-ITD approach to estimate wood volume per tree

Overall, using the ABA method 151.048 trees were detected in the study area, with an average of 994 trees ha-1, while the average number of trees per hectare estimated by the forest inventory was 1708 trees ha-1. Therefore, approximately 42% of the trees were not detected using the above approach, reinforcing the need to correct the remotely sensed tree count by modeling the tree count mean difference (MD).

After identifying the predictive mean difference model, we performed a linear adjustment using the least squares method. The model selected used three variables: maximum height (hmax); the slope (S) estimated from ALS data; and the number of trees detected in the aerial images found using the semi-ITD approach (Tab. 3), which had an RMSE of 12.66%.

The average number of trees per hectare detected by the conversion of ITD into ABA (semi-ITD) for the study area after mean difference (MD) correction was 1707.7, while the inventory data provided a figure of 1708.3 trees per hectare. Therefore, the mean volume per tree estimated from study data was 0.188 m3, while the mean volume of an individual tree derived from forest inventory was 0.195 m3. Fig. 3 shows the wall-to-wall tree volume maps generated in this work.

Fig. 3 - Maps of ABA-derived predictions for stand volume (top left), corrected number of trees (top right) and mean tree volume (bottom).

# Discussion

Accurate stand attribute estimates are critical for effective forest plantation management, partly due to their markedly shorter rotation ages and faster growth rates relative to natural forests ([20]). Moreover, for proper management of forest plantations, information on individual trees is of high value ([33]). Remote sensing has proven to be an excellent data source for supporting sustainable forest management ([37]). However, obtaining accurate estimates of forest attributes is often still elusive ([30]). Our results indicate that the combination of different remote sensing data sources to estimate stand volume can lead to improved accuracy of predictions (RMSE = 5.23%).

In this research, we used two approaches to estimate stand volume in Eucalyptus spp. plantations using ALS as well as high spatial resolution multispectral image data. We first modeled the stand volume using ABA with different combinations of ALS- and image-derived metrics as input data, in three different models. Despite ABA slightly underestimated stand volume (2.98%) compared to the reference volume (total wood volumed measured at the mill) the results were close to those obtained by the forest inventory. However this approach could be enhanced by predicting the number of trees using the ITD approach, this introducing a correction method to improve tree count predictions. We then integrated that information with stand volume previously predicted using the ABA method, allowing for the calculation of mean individual tree volume.

The model applied in the volumetric prediction used two variables derived from aerial imagery, the near infrared band (B4) and the green/red (B2/B3) ratio, and two ALS-derived variables of height (hmode) and cover (Amean). Using the ratio between the green and red bands is analogous to using the near infrared and red bands ratio. The reflectance in the near infrared band is linked to scattering in the mesophyll, which occurs due to the interaction of solar radiation with the internal leaf structure ([4]). Similarly, green wavelengths are reflected by mesophilic leaf tissues, while red light is absorbed by chlorophyll ([4]). Therefore, the ratio between these bands is indicative of vegetative vigor and vegetative abundance. Both vigor and abundance information were important complements to the ALS metrics.

Despite the use of high-resolution imagery, the second approach (which used the aerial image variables) showed the poorest performance among the three approaches tested (RMSE = 8.45%).

Our results show that use of information from a single remote sensing technology can meet the predictive needs of plantation managers, so long as expectations of accuracy and precision are attained accordingly. However, our results also highlight that the integration of data from multiple remote sensing technologies, such as the ALS and high-resolution passive imagery data used here, can lead to more accurate results than the use of a single technology.

This improvement in accuracy is achievable even in situations where the difference in time between data acquisition (2012) and field measurement (2013) is up to one year, as is the case in this study. Because ALS data acquisition can be costly, the field measurement-remote sensing hybrid approach can be of remarkable benefit to forest managers. The integration of data from different sources with different acquisition dates can increase the useful life of forest attribute data collected by any means, potentially decreasing the frequency of acquisition and therefore reducing costs.

While ABA provided estimates of total wood volume for the entire study area, individual tree delineation, detected from passive image data, allowed us to downscale the estimates to the individual tree level. Images were used for tree detection instead of ALS data because the trees were planted at a relatively high density (3 × 2 m) resulting in small canopy sizes. Their small size made it challenging for crowns to intercept a sufficient number of LiDAR points for crown extent to be accurately defined. Therefore, the performance of ALS was poorer when detecting individual trees compared to image-based approaches.

Despite improvements over ALS data, we still found the accuracy of single tree detection from images to be unsatisfactory (58% correctly detected trees) and so a correction factor was introduced. This correction, based primarily on slope, allowed us to increase the accuracy of single tree detection and achieve unbiased individual tree volume estimates. Slope proved to be a key variable to include in the modeling of tree detection error. With increasing slope more trees were located in the same projected area, leading to an aggregation of the canopies and an increase in the maximum tree height of the inventory plots. Increased maximum tree height, therefore, showed a higher correlation with the standard deviation (i.e., heterogeneity) of a plot, which made tree detection from the images more difficult.

In general, the correction procedure applied to our tree counting process performed well, showing good agreement with existing tree count estimates and leading to individual tree volumes with negligible mean difference and low RMSE (12.66%). Our corrected estimates of tree count were greater than uncorrected values identified from the image-based detection by 713.7 trees per hectare, very close to extant estimates based on the forest inventory. The final estimates of individual tree volume after mean difference correction were 0.188 m3 per tree, 3.5% less than the inventory estimate ([33], [20]).

We demonstrated that a traditional ITD approach can be modified by using ABA, in spite of losing the exact location of each tree, but enabling to obtain the desired tree density information ([3]).

# Conclusions

Accurate information about wood volume at both the stand and individual tree levels is required to support the effective management of forest plantations. In this study we enhanced the traditional area-based approach using individual tree detection. The integration of geometrical and spectral information provided by ALS and aerial imagery, respectively, allows for the reliable predictions of volume at both spatial scales. Moreover, remote sensing data integration results in lower mean difference and RMSE. Our approach allowed for the estimation of individual tree volume with negligible bias.

# Acknowledgements

The authors would like to acknowledge the Brazilian Research Foundation, CNPq, and the National Institute for Space Research (INPE) for supporting this research, FIBRIA S/A cellulose company for supplying data from their plantation, and the Integrated Remote Sensing Studio (IRSS) at the Faculty of Forestry, University of British Columbia, Canada, for its assistance, support and friendship. David Williams provided some editorial assistance.

# References

(1)
Batista JLF (2014). Quantificação dos recursos florestais - árvores, arvoredos e florestas [Quantification of forest resources - trees, grove and forests]. Oficina de Textos, São Paulo, Brazil, pp. 384. [in Portuguese]
Gscholar
(2)
Brandtberg T, Warner TA, Landenberger RE, Mcgraw JB (2003). Detection and analysis of individual leaf-off tree crowns in small footprint, high sampling density lidar data from the eastern deciduous forest in North America. Remote Sensing of Environment 85: 290-303.
CrossRef | Gscholar
(3)
Breindenbach J, Astrup R (2014). The semi-individual tree crown approach. In: “Forestry Applications of Airborne Laser Scanninng: Concepts and Case Studies” (Maltamo M ed). Springer Science + Business Media, Dordrecht, Netherlands, pp. 113-133.
Gscholar
(4)
Campbell JB, Wynne RH (2011). Introduction to remote sensing. The Guildford Press, New York, USA, pp. 667.
Online | Gscholar
(5)
CEPAGRI (2014). Clima dos municípios paulistas [Climate of the municipalities of São Paulo]. Web site. [in Portuguese]
Online | Gscholar
(6)
Curtis R (1967). Height-diameter and height-diameter-age equations for second-growth Douglas-fir. Forest science 13 (4): 365-375.
Online | Gscholar
(7)
Ferreira M (1979). Escolha de espécies de eucalipto [Choosing Eucalyptus species]. Circular Técnica IPEF, Piracicaba, Brazil, vol. 47, pp. 1-30. [in Portuguese]
Gscholar
(8)
Hudak AT, Lefsky M, Cohen WB, Berterretche M (2002). Integration of lidar and Landsat ETM+ data for estimating and mapping forest canopy height. Remote Sensing of Environment 82: 397-416.
CrossRef | Gscholar
(9)
IBA (2014). Industria Brasileira de Árvore [Brazilian Industry of trees]. IBA, Brasília, Brazil, pp. 100. [in Portuguese]
Online | Gscholar
(10)
James G, Witten D, Hastie T, Tibshirani R (2014). An introduction to statistical learning with applications in R. Springer Science+Business Media, New York, USA, pp. 426.
Gscholar
(11)
Kaartinen H, Hyppä J, Yu X, Vastaranta M, Hyppä H, Kukko A, Holopainen M, Heipke C, Hirschmugl M, Morsdorf F, Naesset E, Pitkänen J, Popescu S, Solberg BMW, Wu J (2012). An international comparison of individual tree detection and extraction using airborne laser scanning. Remote Sensing 4: 950-974.
CrossRef | Gscholar
(12)
Köppen W, Geiger R (1954). Köppen W, Geiger R (1954): Klima der Erde (Climate of the earth). Wall Map 1:16 Mill. Klett-Perthes, Gotha, Germany.
Online | Gscholar
(13)
Kumar V (2012). Forest inventory parameters and carbon mapping from airborne LiDAR. Dissertation, The University of Twente, Enschede, The Netherlands, pp. 89.
Gscholar
(14)
Lumley T (2014). LEAPS: R package version 2.9 for regression subset selection including exhaustive search. University of Washington, Washington, DC, USA, pp. 8.
Online | Gscholar
(15)
Mccullagh MJ (1988). Terrain and surface modelling systems: theory and practice. Photogrammetric Record 2 (72): 747-749.
CrossRef | Gscholar
(16)
McGaughey RJ (2013). FUSION/LDV: software for LiDAR data analysis and visualization. Pacific Northwest Research Station, USDA Forest Service, WA, USA, pp. 150.
Online | Gscholar
(17)
Mehtätalo L, Nyblom J, Virolainen A (2014). A model-based approach for the recovery of forest attributes using airborne laser scanning data. In: “Forestry Applications of Airborne Laser Scanninng: Concepts and Case Studies” (Maltamo M ed). Springer Science + Business Media, Dordrecht, Netherlands, pp. 193-211.
CrossRef | Gscholar
(18)
Monnet J, Mermim E, Chanussot J, Berger F (2010). Tree top detection using local maxima filtering: a parameter sensitivity analysis. In: Proceedings of the “10th International Conference on LiDAR Application for Assessing Forest Ecossystems”. Freiburg (Germany) Sep 2010. Silvaser 2010, Freiburg, Germany, pp. 9.
Online | Gscholar
(19)
Næsset E (2002). Predicting forest stand characteristics with airborne laser using a practical two-stage procedure field data. Remote Sensing of Environment 80: 88-90.
CrossRef | Gscholar
(20)
Packalén P, Mehtätalo L, Maltamo M (2011). ALS-based estimation of plot volume and site index in a Eucalyptus plantation with a nonlinear mixed-effect model that accounts for the clone effect. Annals of Forest Science 68: 1085.
CrossRef | Gscholar
(21)
Paris C, Bruzzone LA (2015). Three-dimensional model-based approach to the estimation of the tree top height by fusing low-density LiDAR data and very high resolution optical images. IEEE Transactions on Geoscience and Remote Sensing 53 (1): 467-480.
CrossRef | Gscholar
(22)
Payn T, Carnus J, Freer-smith P, Kimberley M, Kollert W, Liu S, Orazio C, Rodriguez L, Silva LN (2015). Changes in planted forests and future global implications. Forest Ecology and Management 352: 57-67.
CrossRef | Gscholar
(23)
Peña EA, Slate EH (2006). Global validation of linear model assumptions. Journal of the American Statistical Association 101 (473): 341-354.
CrossRef | Gscholar
(24)
Peña EA, Slate EH (2014). Gvlma: Global validation of linear models assumptions. R package version 1.0.0.2.
Online | Gscholar
(25)
Pereira JP (2014). Mensuração automática de copas de Araucaria angustifolia (Bertol.) Kuntze a partir de dados LiDAR para estimativa de variáveis dendrométricas [Automatic measurement of Araucaria angustifolia’s crowns (Bertol.) Kuntze from LiDAR data to estimate dendrometric variables]. Dissertation, State University of Santa Catarina, Brazil, pp. 174. [in Portuguese]
Gscholar
(26)
Persson A, Holmgren J, Söderman U (2002). Detecting and measuring individual trees using airborne laser scanner. Photogrammetric Engineering and Remote Sensing 68 (9): 925-932.
Gscholar
(27)
Popescu SC, Wynne RH, Nelson RF (2002). Estimating plot level tree heights with lidar: local filtering with canopy-height based variable window size. Computers and Electronics 37: 71-95.
CrossRef | Gscholar
(28)
R Developing Core Team (2007). R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.
Online | Gscholar
(29)
Roberts SD, Dean TJ, Evans DL, Mcombs JW, Harrington RL, Glass PA (2005). Estimating individual tree leaf area in loblolly pine plantations using LiDAR-derived measurements of height and crown dimensions. Forest Ecology and Management 213: 54-70.
CrossRef | Gscholar
(30)
Saarela S, Grafström A, Stahl G, Kangas A, Holopainen M, Tuominen S, Nordkvist K, Hyyppä J (2015). Model-assisted estimation of growing stock volume using different combinations of LiDAR and Landsat data as auxiliary information. Remote Sensing of Environment 158: 431-440.
CrossRef | Gscholar
(31)
Schöepfer W (1966). Automatisierung des massen, sorten und wertberechnung stenender waldbestände Schriftenreihe Bad [Automation of the mass, varieties and value for calculating forest stock series]. Wurtt-Forstl, Koblenz, Germany.
Gscholar
(32)
Silva CA, Klauberg C, Carvalho SPC, Hudak AT, Rodriguez LCE (2014). Mapping aboveground carbon stocks using LiDAR data in Eucalyptus spp. plantations in the state of São Paulo, Brazil. Scientia Forestalis 42 (104): 591-604.
Online | Gscholar
(33)
Simões D (2008). Avaliação econômica de dois sistemas de colheita florestal mecanizada de eucalipto [Economic evaluation of two systems of mechanized forest harvest of Eucaliptus]. Master Thesis, Faculty of Agronomic Sciences Ciências, UNESP, Botucatu, SP, Brazil, pp. 105. [in Portuguese]
Online | Gscholar
(34)
Vastaranta M, Holopainen M, Yu X, Hyppä J, Mäkinen A, Rasinmäki J, Melkas T, Kaartinen H, Hyyppä H (2011). Effects of individual tree detection error sources on forest management planning calculations. Remote Sensing 3: 1614-1626.
CrossRef | Gscholar
(35)
Wang L, Gong P, Biging GS (2004). Individual tree-crown delineation and treetop detection in high-spatial-resolution aerial imagery. Photogrammetric Engineering and Remote Sensing 70 (3): 351-357.
CrossRef | Gscholar
(36)
White JC, Wulder MA, Varhole A, Vastaranta M, Coops NC, Cook DC, Pitt D, Woods M (2013). A best practices guide for generating forest inventory attributes from airborne laser scanning data using an area-based approach (Version 2.0). Information Report FI-X-010, Natural Resources Canada, Canadian Forest Service and Canadian Wood Fiber Centre, Victoria, BC, Canada, pp. 41.
Online | Gscholar
(37)
Wulder M, Niemann KO, Goodenough DG (2000). Local Maximum Filtering for the extraction of tree locations and basal area from high spatial resolution imagery. Remote Sensing of Environment 73: 103-114.
CrossRef | Gscholar
(38)
Wulder M, Seemann D (2003). Forest inventory height update through the integration of lidar data with segmented Landsat imagery. Canadian Journal of Remote Sensing 29 (5): 536-543.
CrossRef | Gscholar
(39)
Zhu X, Liu D (2015). Improving forest aboveground biomass using seasonal Landsat NDVI time-series. ISPRS Journal of Photogrammetry and Remote Sensing 102: 222-231.
CrossRef | Gscholar

#### Authors’ Affiliation

(1)
Emily T Shinzato
Yosio E Shimabukuro
Department of Remote Sensing, National Institute for Space Research - INPE, São Paulo (Brazil)
(2)
Nicholas C Coops
Piotr Tompalski
Department of Forest Resources Management, University of British Columbia, Victoria, BC (Canada)
(3)
Esthevan AG Gasparoto
Department of Forest Sciences, Luiz de Queiroz College of Agriculture, University of São Paulo , Piricicaba (Brazil)

#### Citation

Shinzato ET, Shimabukuro YE, Coops NC, Tompalski P, Gasparoto EAG (2016). Integrating area-based and individual tree detection approaches for estimating tree volume in plantation inventory using aerial image and airborne laser scanning data. iForest 10: 296-302. - doi: 10.3832/ifor1880-009

#### Paper history

Accepted: Aug 26, 2016

First online: Dec 15, 2016
Publication Date: Feb 28, 2017
Publication Time: 3.70 months

© SISEF - The Italian Society of Silviculture and Forest Ecology 2016

#### Breakdown by View Type

(Waiting for server response...)

#### Article Usage

Total Article Views: 12298
(from publication date up to now)

Breakdown by View Type
HTML Page Views: 9205
Abstract Page Views: 547

Web Metrics
Days since publication: 1298
Overall contacts: 12298
Avg. contacts per week: 66.32

Article citations are based on data periodically collected from the Clarivate Web of Science web site
(last update: Jun 2020)

Total number of cites (since 2017): 3
Average cites per year: 0.75

#### iForest Database Search

Search By Author

Search By Keyword

Citing Articles

Search By Author

Search By Keywords

#### PubMed Search

Search By Author

Search By Keyword