Next Article in Journal
Variation in Trembling Aspen and White Spruce Wood Quality Grown in Mixed and Single Species Stands in the Boreal Mixedwood Forest
Previous Article in Journal
Extinction Risk of Pseudotsuga Menziesii Populations in the Central Region of Mexico: An AHP Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of Wooded Area Reduction by Airborne Laser Scanning

1
Department of Geodesy and Geoinformation, Vienna University of Technology, Gußhausstraße 27-29, A-1040 Vienna, Austria
2
Faculty of Surveying and Mapping, Hanoi University of Mining and Geology, Hanoi 10000, Vietnam
*
Authors to whom correspondence should be addressed.
Forests 2015, 6(5), 1613-1627; https://doi.org/10.3390/f6051613
Submission received: 20 January 2015 / Revised: 8 April 2015 / Accepted: 27 April 2015 / Published: 7 May 2015

Abstract

:
Airborne Laser Scanning (ALS) data hold a great deal of promise in monitoring the reduction of single trees and forests with high accuracy. In the literature, the canopy height model (CHM) is the main input used frequently for forest change detection. ALS also has the key capability of delivering 3D point clouds, not only from the top canopy surface, but also from the entire canopy profile and also from the terrain. We investigated the use of two additional parameters, which exploit these capabilities for assessing the reduction of wooded area: Slope-adapted echo ratio (sER) and Sigma0. In this study, two ALS point cloud data sets (2005 and 2011) were used to calculate Digital Surface Model (DSM), sER, and Sigma0 in 1.5 km2 forest area in Vorarlberg, Austria. Image differencing was applied to indicate the change in the three difference models individually and in their combinations. Decision trees were used to classify the area of removed trees with the minimum mapping unit of 13 m2. The final results were evaluated by a knowledge-based manual digitization using completeness and correctness measures. The best result is achieved using the combination of sER and DSM, namely a correctness of 92% and a completeness of 85%.

1. Introduction

Forests are an important factor in maintaining the balance in the Earth system. However, the ecological processes are often affected by human activity [1]. In order to control the change of forests under the impacts of deforestation, wind throw, and diseases, it is required for forest managers to apply techniques supporting monitoring and updating forest information regularly. In situ forest inventory and remote sensing technologies are in use for detecting and monitoring these changes. Remote sensing, as one of these techniques, has proven its ability in change detection, automatically, efficiently, and consistently, especially for large areas, this also requiring less manual labor.
Optical remote sensing is a good choice to detect changes of forests as demonstrated in different studies [2,3,4,5]. However, optical acquisition techniques are limited by clouds if they are below the platform. Photograhic techniques additionally provide little information in cast shadow, as well as topographic shadow areas, providing little or no texture and, thus, lead to lower accuracy there. Furthermore, the varying shadow conditions for different acquisition times limit an automatic large area derivation of forest changes.
The image matching technique shows its potential use in forestry [6,7]. Photogrammetric imagery can also be exploited for gaining 3D point clouds [8]. Nevertheless, the limitation of this technique is that information is restricted to the point cloud of the upper canopy and, thus, does not provide ground height [9]. The image-based point cloud quality depends on factors like ground sampling distance, radiometric image resolution, stereo-parameters, viewing geometry, sun-angle, and amount of shadows [9].
With the advantage of penetrating the canopy through small gaps, Airborne Laser Scanner (ALS) is a potential technique for monitoring vegetation changes. Being an active technique, ALS emits its own energy for sensing and is consequently not affected by the ambient illumination (cast shadows, shadows of high clouds). Using 3D point clouds from ALS, the change in both coverage and height can be detected [10,11]. Moreover, not only dense forests, but even single trees can be detected from dense point clouds [12], which can be used to estimate forest biomass [13,14] and generate 3D tree models [15,16].
Although ALS data holds a high promise in vegetation change detection, thus far, research using multi-temporal ALS to detect forest change cover has not yet been fully explored. Yu et al. [17] detected harvested trees using two small footprint, high sampling density ALS acquisitions based on image differencing. Three-dimensional canopy height models (CHM) were calculated for both data sets using raster-based algorithms. The major change of CHM at the same pixel was acquired by a threshold value. They reported that 61 out of 83 field-checked harvested trees were detected automatically. The undetected trees were mainly smaller trees. St-Onge et al. [18] also used the threshold of CHM difference of two medium density LiDAR data acquired in 1998 and 2003 to identify new canopy gaps and assess height growth. With the same data Vepakomma et al. [19] expanded their study in accessing the feasibility of small footprint LiDAR to map the canopy gap expansions and canopy gap closures for the conservation zone of Quebec. Vastaranta et al. [20] developed a ∆CHM method for canopy change detection of snow-damaged trees by applying bi-temporal LiDAR data for the period 2006–2010. Næsset and Gobakken [21] estimated a boreal forest growth over two years by using canopy metrics, i.e., measures of spatial distribution of the acquired point cloud. Nyström et al. [22] employed histogram matching to calibrate the metrics in order to reduce the difference between two ALS datasets and produce change imagery. They controlled the changes by partial and complete tree removal in selected plots. Hollaus et al. [23] assessed the changes in growing stock of 160 km2 mountain forest by two ALS datasets. The model was established with 184 FI (forest inventory) plots.
The above-mentioned studies used the differences found in (two) datasets over a forest area. Alternatively, processes in the forest can be detected using only one dataset by the traces they leave in the site. For example, Mücke et al. [24] used full-waveform ALS data, obtained under leaf-off conditions to detect fallen trees. Here, an echo width model was derived based on the point cloud and normalized echo heights in order to delineate downed stems. Lindberg et al. [25] and Nyström et al. [26] contributed two different methods highlighting the potential of high density ALS data to detect wind-thrown trees under forest canopy. Lindberg et al. [25] used a line template matching method applied directly to the ALS point cloud (69 points/m2), while Nyström et al. [26] used the difference between two elevation models created from the same high density ALS data to detect wind-thrown trees.
In this research, we investigate the ability of forest reduction detection from two different ALS datasets by using image differencing [27]. In comparison to Nyström, who applied histogram matching to account for different sensor characteristics, our aim is to find features of LiDAR point clouds, that are, as much, as possible independent of the sensor characteristics. Unlike other studies mentioned above, our hypothesis is that forest reduction up to individual trees can be observed by the three LiDAR derived models: Digital Surface Model (DSM), Slope-adaptive Echo Ratio (sER), and “Sigma0” (a local roughness measure), as well as their combinations. Our primary interest is, thus, not to demonstrate that ALS can detect forest changes as this was done before [17,22]. Rather, we are interested in finding robust methods in the presence of different ALS mission parameters. The study is done for a mountain forest in Vorarlberg, Austria. Two epochs were acquired with six years difference between the data acquisitions. The DSMs demonstrate the change in height and, thus, indicate that tall objects were removed. sER demonstrates the change in vertical penetrability and indicates that layered objects (e.g., understory and canopy) were removed. Sigma0 demonstrates the change in the vertical dispersion of the points and indicates that objects distributed in height (e.g., trees, bushes) were removed. Using these variables, in different combination, the forest reduction is derived. The results are accessed with their accuracy by using the completeness and correctness measure. All the processes are supported by OPALS [28] and ArcGIS software.

2. Material and Method

2.1. Study Area

Study area is located in the south of Vorarlberg province, Austria, centered on lat. N 47°04′12″, long. E 9°49′12″ (Figure 1). The total covered area is about 1.5 km2 of mountainous region, with the elevations ranges from 1225 m above sea level (a.s.l) in the valleys, to a maximum of 1786 m a.s.l. In general, in this area, approx. 24% is covered by forest and the dominating tree species is Norway spruce (Picea abies). The forests in the study area are managed by the Stand Montafon Forstfonds (http://stand-montafon.at/forst), which operates the local forest inventory. Based on this inventory the tree heights vary between approximately 6 m to 42 m, with the mean height of 27.5 m and the standard deviation is 6.8 m [29]. Due to the topography the majority of the forests in the study area have a protection function against natural hazards, i.e., snow avalanches. Therefore, exploitations of single trees or group of trees are foreseen in the forest management plan, meaning clear cuttings of larger areas are not allowed.
Figure 1. Study area (Orthophoto: office of survey and geoinformation from Vorarlberg, Austria, Political Map of Austria: GEOATLAS.com.).
Figure 1. Study area (Orthophoto: office of survey and geoinformation from Vorarlberg, Austria, Political Map of Austria: GEOATLAS.com.).
Forests 06 01613 g001

2.2. Data

The two ALS data sets were provided by the local office of survey and geo-information from the federal state of Vorarlberg, and are subsets of the official federal state-wide ALS data acquisition campaigns. The ALS data sets were acquired in 2005 and 2011 using an Airborne Laser Terrain Mapper systems (ALTM 1225) and a Trimble Harrier 56 system, respectively. The study area is covered by two flight strips for the first ALS acquisition and three flight strips for the second acquisition. All ALS data sets were acquired under snow-free and leaf-off conditions and were available as georeferenced 3D-point clouds. In Table 1, the relevant sensor characteristics of used ALS systems are summarized.
Table 1. Summary of sensor characteristics of the applied ALS systems.
Table 1. Summary of sensor characteristics of the applied ALS systems.
Sensor CharacteristicsSensors
Acquisition Year 2005Acquisition Year 2011
Optech ALTM 1225Trimble Harrier 56
Beam divergence0.3 mrad<0.5 mrad
Max scan angle (from nadir)20°30°
Wavelength1064 nm1064 nm
Pulse repetition frequency<25 kHz160 kHz
Sensor typeDiscreteFull waveform
Average point density6 echos/m224 echos/m2

2.3. Reference Data

The reference data was derived, based on image interpretation of aerial orthophotos [30] with additional use of 3D point cloud viewing of the raw ALS data. Orthophotos with a spatial resolution of 0.12 m (2012) and 0.5 m (2005) are used to overview the forest cover status. Using additionally various visualizations of the ALS data, such as 3D point cloud visualization in FugroViewer software, nDSM values of the same pixels in two epochs, each single 1 m2 pixel was evaluated and digitized. Due to border effects, a tree crown can be influenced in the DSMs by a slightly different shape. Therefore, small changes of very few m2 were not considered as relevant and the therefore, the minimum mapping area was set to 13 m2. This process took roughly 100 h. For the 223 digitized polygons, the minimum polygon area is 13 m2 and the maximum area is 2351 m2, the average size of a harvested polygon is 152 m2 and the standard deviation is 289 m2. The mean height of harvested polygons is 30.8 m and the standard deviation is 8.2 m. The final digitized output is converted into a binary format with a raster size of 1 × 1 m2 that is used for accuracy assessment.

2.4. Method

In this study 3D ALS point clouds are used as the basis for deriving the following parameters, which are used for delineating harvested trees:
  • DSM,
  • sER, and
  • Sigma0.
For forest change detection the DSMs are an important input because they describe the height changes of the top most canopy surface. This means a decrease of the DSM indicates the loss of trees. To derive the DSMs from the two ALS data sets, the land cover dependent method described, in Hollaus et al. [31], is applied. This method uses the strengths of different algorithms for generating the final DSM by using surface roughness information to combine two DSMs, which are calculated based on the highest echo within a raster cell, and on moving least squares (mls) interpolation with a plane as functional model (i.e., a tilted regression plane fitted through the k-nearest neighbors). For smooth areas (e.g., roof planes, streets, short grassland), noise reduction by moving least squares interpolation is exploited, whereas for rough surfaces (e.g., canopy surface, building edge) the highest point within a raster cell is used (DSM(X, Y)max). The input for moving least squares interpolation is a subset of the point cloud (highest points within 0.5 m rater cells), which ensures that the interpolated surface goes through the tree tops. The derived final DSMs have a spatial resolution of 1 m.
More formally, the DSM calculation runs in the following way (Here we use the trinary operator c ? r1: r2. Its value depends on the condition c. If c is TRUE, the result is r1, otherwise the result is r2.).
z[DSM (Xi, Yi)]2005 = z[σz(Xi, Yi)]2005 < 0.5 or not z[DSM (Xi, Yi)max]2005 ? z[DSM (Xi, Yi)mls]2005:
z[DSM (Xi, Yi)max]2005
z[DSM (Xj, Yj)]2011 = z[σz(Xj, Yj)]2011 < 0.5 or not z[DSM (Xj, Yj)max]2011 ? z[DSM (Xj, Yj)mls]2011:
z[DSM (Xj, Yj)max]2011
ΔDSM (X, Y) = DSM2011 (X, Y) − DSM2005(X, Y)
The sER is a measure that describes the vertical point distribution and thus the penetrability of the surface [32,33]. The echo ratio (not slope-adaptive) is defined as the ratio between the number of neighboring echoes in a fixed search distance, measured in 3D (a sphere, n3D, see appendix for exact definition), and all echoes located within the same search distance in 2D (a vertical cylinder, n2D). To guarantee a correct ER on steep slopes, the search radius of the sphere (r2D) has to be extended considering the slope (r3D = r2D/cos(α)) (i.e., dividing the initial 3D search distance by the cosine of the slope). Thus, the derived ER is the slope-adaptive echo ratio, sER.
s ER = n 3 D ( r / cos α ) n 2 D ( r ) × 100
The sER is computed for each echo in the first processing step and shows for continuous and impenetrable surface (i.e., ground and roof surfaces) values of 100% and for tree canopy points of lower value. For further analyses, the sER is aggregated in 1 m cells using the max value within each cell. The modules opalsEchoRatio and opalsCell were used for this computation.
ΔsER (X, Y) = sER2011(X, Y) − sER2005(X, Y)
For the computation of Sigma0, all echoes are used with a neighborhood size of the ten nearest neighbors. It is derived during the interpolation of a height model using the moving least squares approach. The standard deviation of the residuals in this interpolation is determined at each grid post. This provides a grid congruent with the interpolated heights. In each grid post, this value of Sigma0 indicates how well all of the original points fit to the least squared plane. The grid width was 1 m.
ΔSigma0 (X, Y) = Sigma02011(X, Y) − Sigma02005(X, Y)
As the decreases of the elevation of the canopy surface indicate the loss of trees, the DSM2005 is subtracted from the DSM2011. In addition to the elevation changes, the changes in the vertical echo distribution and penetrability described with the sER model and the changes of the surface roughness, represented by the Sigma0 model, also indicate areas with lost (e.g., harvested) trees. Thus, difference models of the sER and the Sigma0 models are also calculated. To use this information, it is assumed that the sER of point clouds have a lower value for trees than for open ground (e.g., removed trees) and for Sigma0 it receives a higher value.
The thresholds for all change values (i.e., DSM, sER and Sigma0) are assessed empirically. No primary maps of vegetated area were derived for the epochs. In other words, no maps of DSM, sER, and Sigma0 were studied. Instead, we are searching directly for thresholds on the observed differences, thus, in ΔDSM, ΔsER, and ΔSigma0.
All three data sets (ΔDSM, ΔsER, ΔSigma0) indicate lost tree positions. Additionally, combinations of those input variables, for improving the accuracy of the final result, were investigated. Each pair and the triple of variables are used with newly determined empirical thresholds. The change results can be expected to show small errors, localized in single pixels or very small groups of pixels (i.e., along the border of trees or forests). This is caused by the accuracy of the acquired data (i.e., point density, georeferencing), as well as by the interpolation. The results are, therefore, converted into binary format for applying methods of binary mathematical morphology. Closing and then opening morphology with a circular kernel shape with a diameter of 1 and 2 pixels respectively are applied to all output raster datasets to reduce noise and to smooth object outlines. Finally, six change detection outputs were established: DSM only, sER only, Sigma0 only, DSM combined with sER, sER combined with Sigma0 and DSM, and sER and Sigma0. Early in the analysis it became obvious, that the combination of DSM and Sigma0 offers no increase of the achievable accuracy than the other pairs, thus, it was omitted.
Completeness (Comp) and correctness (Corr) [34] are used for the accuracy assessment of the final results.
Comp = TP TP + FN
Corr = TP TP + FP
The forest reduction area in the reference and change detection results are compared, where a true positive (TP) indicate the change in both datasets, false negative (FN) is labeled in the reference data but has no correspondence in the change detection results, and false positive (FP) is labeled in the change detection results and has no corresponding in the reference data.

3. Results and Discussion

The main properties of the primary models (DSM, sER and Sigma0) are summarized in Table 2. Based on the empirical analyses, the thresholds for each of the different image (ΔDSM, ΔsER and ΔSigma0 variables were found and summarized in Table 3. Figure 2 shows the harvested tree detection results of the three variables DSM, sER, and Sigma0 independently, as well as the manual digitized reference map. As can be seen, the downed tree area was detected, more or less, correctly. However, the results are affected by salt-and-pepper type of noise. After applying image morphology operations, the Sigma0 final results are still strongly affected by this and give the worst accuracy (Table 3) compared with the other two variables (Figure 3a–c).
Table 2. Minimum, maximum, mean, standard deviation of DSM, sER, and Sigma0 in the different epochs.
Table 2. Minimum, maximum, mean, standard deviation of DSM, sER, and Sigma0 in the different epochs.
MinMaxMeanStd.dev
DSM2005(m)1226.21817.11516.4124.0
DSM2011(m)1226.02022.1 11517.1124.0
sER2005(%)1.9100.069.929.3
sER2011(%)2.3100.065.529.9
Sigma02005(m)0.028.02.94.2
Sigma02011(m)0.0220.3 22.74.2
1 This value is affected by gross errors. Excluding them leads to a max DSM of 1817.9; 2 Excluding gross error the max Sigma02011 is 27.2.
Figure 2. Forest reduction based on the selected thresholds for the (a) DSM; (b) sER; (c) Sigma0; and (d) Knowledge-based digitization.
Figure 2. Forest reduction based on the selected thresholds for the (a) DSM; (b) sER; (c) Sigma0; and (d) Knowledge-based digitization.
Forests 06 01613 g002
Table 3. Threshold values and accuracy measures.
Table 3. Threshold values and accuracy measures.
Change thresholdDSM (m)sER (%)Sigma0 (m)Corr (%)Comp (%)
ΔDSM<−7.0--------84.690.9
ΔsER---->30----87.587.1
ΔSigma0--------<−7.038.656.8
ΔDSM and sER<−2.0>27----91.985.1
ΔsER and Sigma0---->27<−2.090.980.8
ΔDSM and sER and Sigma0<−2.0>25<−1.092.882.4
ΔDSM and sER and Sigma01<−7.0>30<−7.096.438.2
1 The result of this combination is not shown in the Figure 3.
Figure 3. Final forest reduction after morphological operation (a) DSM; (b) Echo ratio; (c) Sigma0; (d) Echo ratio and DSM; (e) Echo ratio and Sigma0; (f) Echo ratio and DSM and Sigma0.
Figure 3. Final forest reduction after morphological operation (a) DSM; (b) Echo ratio; (c) Sigma0; (d) Echo ratio and DSM; (e) Echo ratio and Sigma0; (f) Echo ratio and DSM and Sigma0.
Forests 06 01613 g003
The limitations of DSM and Sigma0 compared to sER are to detect the reduction of low trees. Using a threshold, which is too low in height change, the lost tree cover is easy to be mixed with the unchanged forest cover. Additionally, Sigma0 may depend much more on flying parameters, such as the flying height, which influences the point density and consequently the Sigma0 values. Regarding to sER, because sER is larger influenced by its neighborhood so it has an increased value for the lost trees in a larger area. This leads to a reduction in the capability of detecting lost trees in dense canopy regions.
In order to overcome the limitation of each variable, all variables are incorporated in combinations in order to achieve improved results. To acquire the threshold of each pair combination between ΔsER versus ΔDSM, ΔsER versus ΔSigma0, and ΔDSM versus ΔSigma0, a feature space (Figure 4) is used to distinguish changed areas. The area of change is in either case in the upper left part of the feature space. It deviates from the distribution of unchanged areas, which is centered on (0, 0). As stated above, the pair ΔDSM versus ΔSigma0 does not show the discrimination in the feature space compared to other pairs, thus it is not used to detect changes. The median value of ΔDSM, ΔsER and ΔSigma0 of the reference data also calculated and plotted into the feature space for delineating the threshold. From the feature space, the threshold for the ΔDSM and ΔsER combination (ΔDSM < −2 AND ΔsER > 27) and the ΔsER and ΔSigma0 combination (ΔsER > 27 AND ΔSigma0 < −2) are found, as a result, the accuracy of change detection is improved (Table 3, Figure 3). The former provides a correctness of 91.9% and completeness of 85.1% higher than the latter with a correctness of 90.9% and a Completeness of 80.8%. The combination of all three variables awards the highest correctness (92.8%) and a lower completeness (82.4%). As could be expected, using the original thresholds of the single variable classifications shows the highest correctness (96%) at the cost of a lower completeness (only 38.2%). The accuracy assessment of seven final change detection results is shown in Table 3.
Figure 4. Feature space of (a) ΔsER versus ΔSigma0; (b) ΔsER versus ΔDSM. The red lines indicate the threshold values for change detection. Brighter colors denote a higher density of points (from blue, to green, to red). Red dots indicated the median value of ΔDSM, ΔsER, and ΔSigma0 of the reference polygons.
Figure 4. Feature space of (a) ΔsER versus ΔSigma0; (b) ΔsER versus ΔDSM. The red lines indicate the threshold values for change detection. Brighter colors denote a higher density of points (from blue, to green, to red). Red dots indicated the median value of ΔDSM, ΔsER, and ΔSigma0 of the reference polygons.
Forests 06 01613 g004
Our final result achieved a accuracy (with a minimum mapping unit of 13 m2) compared to the research of Yu et al. [17], who detected 61 out of 83 harvested trees with an accuracy of 73.5%. St-Onge and Vepakomma [18] used canopy height differences to identify new gaps (especially fallen trees) and the minimum area criterion is 5 m2. Their data has a density of 3 shots/m2 in each epoch. Producer and user accuracies are similar to ours, although a bit higher (95%–98%). However, they apparently confirmed the existence of gaps in the reference (optical images) and the LiDAR result, but not the exact spatial location. Small edge localization errors and gap size differences between our reference and LiDAR result add to lower producer and user accuracy in our case. Additionally, the results of Nyström et al. [22] can be compared to our results. Their overall accuracy in detection is 88%, thus, somewhat lower than our result. However, in the forest-tundra ecotone, in which their study is set, the geometric signal of changes is lower than in ours in the case of taller trees.
Analyzing the Sigma0 values of the different missions indicated that they depend more strongly on the parameters of the data acquisition. Two point clouds of a single tree are displayed and visualized in 3D, and it is realized that the point distributions and point densities of the same tree are different. This will influence the value of the Sigma0 results. Therefore, in this study, surface roughness (Sigma0) alone is not a reasonable measure to detect forest change.
In Figure 4, there are some sER median values for reference polygons located under the selected threshold. It is explained that, in some dense tree positions, the sER values of downed trees are influenced by points of the other surrounding trees.
In the case of deciduous forests with a dense canopy surface and fully developed foliage, there is no penetration into the canopy and to the ground. It leads to high sER and small Sigma0 values for trees. Harvested trees can be better detected with ΔDSM in this case. On the other hand, errors in vertical geo-referencing directly influence the DSM and, therefore, ΔDSM, but it does not affect sER and Sigma0. Therefore, the combination of sER and DSM will provide the highest quality in detecting the reduction of wooded area.
Scan frequency, flying altitude, scan angle, acquisition time (i.e., leaf-on, leaf-off), and applied methods for preprocessing have an influence on data quality [35,36]. Based on the applied method for DSM calculation, the influences of these properties are minimized [31].
For a detailed assessment of forest biomass changes it is important to differentiate between forest growth and exploitation. For the quantification of the exploitation, detailed information about reduced (i.e., harvested trees) forest area is required. The assessment of forest growth is based on changes in the DSM, which requires robust methods to derive DSM from the ALS point clouds that are, as much as possible, independent from sensor characteristics and data acquisition settings [31]. To minimize influences originating from georeferencing issues, advanced georeferencing, including a strip adjustment, is normally required [23]. For the assessment of the biomass change, regression models can be applied to the individual data sets. For the calibration of these regression models [37], reference data, e.g., forest inventory data, are required. Finally, the biomass change can be derived from the differences between the individual biomass maps.

4. Conclusions

Using the image differencing method, a traditional pixel-based change detection method was applied to detect reduction of forest area. In this study, we used the three variables DSM, slope adaptive echo ratio (sER), and Sigma0, derived from two different ALS data sets, to detect downed trees in a forest. While many studies have, thus far, used DSM and its change, we found that sER is a good single predictor for tree cover change. sER is a local measurement, which means that global height differences, e.g., ALS block geo-referencing problems, do not influence the assessment. It is also noted that the threshold value for sER (Table 2) did not change as much as for DSM and Sigma0 in the combined classifications.
The incorporations of two or more variables always improved the quality of detection results, only the combination of DSM and Sigma0 does not provide an improvement. This study opens up a new application of discrete return ALS data in forest change detection and, therefore, in forest management. We conclude that the best results were achieved using sER change with additional consideration of DSM change, namely a correctness of 92% and a completeness of 85%.

Acknowledgments

This work is partly founded from the European Community’s Seventh Framework Programme (FP7/2007–2013) under grant agreement No. 606971, the Advanced_SAR project.

Author Contributions

Giang Tran processed the data, performed the experiment, and drafted the article. Duy Nguyen analyzed the data. Markus Hollaus reviewed and edited the manuscript. Norbert Pfeifer supervised the research, reviewed and edited the manuscript.

Appendix

For each point of the point set P, we computed the number of neighbors n 3 D , i in a 3D spherical neighborhood and the number of neighbors n 2 D , i in a cylindrical neighborhood.
p i P ,     i = 1 , , n ,     p i   R 3 ,     p i =   ( x i ,   y i ,   z i )
n 3 D , i = | { q j P |   p i q j r } |
n 2 D , i = | { q j P |   ( x i x j ) 2 + ( y i y j ) 2 r 2 } | ,   j = 1 , , n
n 3 D , i ( r ) ,   n 2 D , i ( r )
s E R = n 3 D ( r / cos α ) n 2 D ( r ) × 100

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. FAO, State of the World’s Forests; Food Agriculture Organization of the Uniterd States: Rome, Italy, 2012.
  2. Broich, M.; Hansen, M.; Stolle, F.; Potapov, P.; Margono, B.A.; Adusei, B. Remotely sensed forest cover loss shows high spatial and temporal variation across sumatera and kalimantan, indonesia 2000–2008. Environ. Res. Lett. 2011, 6, 014010. [Google Scholar] [CrossRef]
  3. Donoghue, D.N.M.; Watt, P.J.; Cox, N.J.; Dunford, R.W.; Wilson, J.; Stables, S.; Smith, S. An evaluation of the use of satellite data for monitoring early development of young sitka spruce plantation forest growth. Forestry 2004, 77, 383–396. [Google Scholar] [CrossRef]
  4. Huang, X.; Friedl, M.A. Distance metric-based forest cover change detection using modis time series. Int. J. Appl. Earth Obs. Geoinf. 2014, 29, 78–92. [Google Scholar] [CrossRef]
  5. Potapov, P.; Hansen, M.C.; Stehman, S.V.; Loveland, T.R.; Pittman, K. Combining modis and landsat imagery to estimate and map boreal forest cover loss. Remote Sens. Environ. 2008, 112, 3708–3719. [Google Scholar] [CrossRef]
  6. Baltsavias, E.; Gruen, A.; Eisenbeiss, H.; Zhang, L.; Waser, L.T. High-quality image matching and automated generation of 3d tree models. Int. J. Remote Sens. 2008, 29, 1243–1259. [Google Scholar] [CrossRef]
  7. Bohlin, J.; Wallerman, J.; Fransson, J.E.S. Forest variable estimation using photogrammetric matching of digital aerial images in combination with a high-resolution dem. Scand. J. For. Res. 2012, 27, 692–699. [Google Scholar] [CrossRef]
  8. Otepka, J.; Ghuffar, S.; Waldhauser, C.; Hochreiter, R.; Pfeifer, N. Georeferenced point clouds: A survey of features and point cloud management. ISPRS Int. J. Geo-Inf. 2013, 2, 1038–1065. [Google Scholar] [CrossRef]
  9. White, J.; Wulder, M.; Vastaranta, M.; Coops, N.; Pitt, D.; Woods, M. The utility of image-based point clouds for forest inventory: A comparison with airborne laser scanning. Forests 2013, 4, 518–536. [Google Scholar] [CrossRef]
  10. Houldcroft, C.J.; Campbell, C.L.; Davenport, I.J.; Gurney, R.J.; Holden, N. Measurement of canopy geometry characteristics using lidar laser altimetry. IEEE Trans. Geosci. Remote Sens. 2005, 43, 2270–2282. [Google Scholar] [CrossRef]
  11. Yu, X.; Hyyppä, J.; Kukko, A.; Maltamo, M.; Kaartinen, H. Change detection techniques for canopy height growth measurements using airborne laser scanner data. Photogramm. Eng. Remote Sens. 2006, 72, 1339–1348. [Google Scholar] [CrossRef]
  12. Lindberg, E.; Eysn, L.; Hollaus, M.; Holmgren, J.; Pfeifer, N. Delineation of tree crowns and tree species classi fi cation from full-waveform airborne laser scanning data using 3-d ellipsoidal clustering. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 3174–3181. [Google Scholar] [CrossRef]
  13. Jochem, A.; Hollaus, M.; Rutzinger, M.; Hofle, B. Estimation of aboveground biomass in alpine forests: A semi-empirical approach considering canopy transparency derived from airborne lidar data. Sensors 2011, 11, 278–295. [Google Scholar] [CrossRef] [PubMed]
  14. Kim, S.-R.; Kwak, D.-A.; Lee, W.-K.; Son, Y.; Bae, S.-W.; Kim, C.; Yoo, S. Estimation of carbon storage based on individual tree detection in pinus densiflora stands using a fusion of aerial photography and lidar data. Sci. China Life Sci. 2010, 53, 885–897. [Google Scholar] [CrossRef] [PubMed]
  15. Vosselman, G. 3D recontruction of roads and trees for city modelling. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2003, 34, 231–236. [Google Scholar]
  16. Rutzinger, M.; Pratihast, A.K.; Oude Elberink, S.J.; Vosselman, G. Tree modelling from mobile laser scanning data-sets. Photogramm. Rec. 2011, 26, 361–372. [Google Scholar] [CrossRef]
  17. Yu, X.; Hyyppä, J.; Kaartinen, H.; Maltamo, M. Automatic detection of harvested trees and determination of forest growth using airborne laser scanning. Remote Sens. Environ. 2004, 90, 451–462. [Google Scholar] [CrossRef]
  18. St-onge, B.; Vepakomma, U. International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences. In Assessing Forest Gap Dynamics and Growth Using Multi-Temporal Laser-Scanner Data; ISPRS Archives: Freiburg, Germany, 2004; pp. 173–178. [Google Scholar]
  19. Vepakomma, U.; St-Onge, B.; Kneeshaw, D. Spatially explicit characterization of boreal forest gap dynamics using multi-temporal lidar data. Remote Sens. Environ. 2008, 112, 2326–2340. [Google Scholar] [CrossRef]
  20. Vastaranta, M.; Korpela, I.; Uotila, A.; Hovi, A.; Holopainen, M. Mapping of snow-damaged trees based on bitemporal airborne lidar data. Eur. J. For. Res. 2012, 131, 1217–1228. [Google Scholar] [CrossRef]
  21. Nasset, E.; Gobakken, T. Estimating forest growth using canopy metrics derived from airborne laser scanner data. Remote Sens. Environ. 2005, 96, 453–465. [Google Scholar] [CrossRef]
  22. Nyström, M.; Holmgren, J.; Olsson, H. Change detection of mountain birch using multi-temporal als point clouds. Remote Sens. Lett. 2013, 4, 190–199. [Google Scholar] [CrossRef]
  23. Hollaus, M.; Eysn, L.; Karel, W.; Pfeifer, N. Growing stock change estimation using airborne laser scanning data. In Proceedings of the Talk: 13th International Conference on LiDAR Applications for Assessing Forest Ecosystems (SilviLaser 2013), Beijing, China, 9–11 October 2013.
  24. Mücke, W.; Hollaus, M.; Pfeifer, N.; Schroiff, A.; Deák, B. Comparison of discrete and full-waveform als for dead wood detection. ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci. 2013, 199, 199–204. [Google Scholar] [CrossRef]
  25. Lindberg, E.; Hollaus, M.; Mücke, W.; Fransson, J.E.S.; Pfeifer, N. Detection of lying tree stems from airborne laser scanning data using a line template matching algorithm. ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci. 2013, 169, 169–174. [Google Scholar] [CrossRef]
  26. Nyström, M.; Holmgren, J.; Fransson, J.E.S.; Olsson, H. Detection of windthrown trees using airborne laser scanning. Int. J. Appl. Earth Obs. Geoinf. 2014, 30, 21–29. [Google Scholar] [CrossRef]
  27. Hussain, M.; Chen, D.; Cheng, A.; Wei, H.; Stanley, D. Change detection from remotely sensed images: From pixel-based to object-based approaches. ISPRS J. Photogramm. Remote Sens. 2013, 80, 91–106. [Google Scholar] [CrossRef]
  28. Pfeifer, N.; Mandlburger, G.; Otepka, J.; Karel, W. Opals—A framework for airborne laser scanning data analysis. Comput. Environ. Urban Syst. 2014, 45, 125–136. [Google Scholar] [CrossRef]
  29. Hollaus, M.; Wagner, W.; Eberhöfer, C.; Karel, W. Accuracy of large-scale canopy heights derived from lidar data under operational constraints in a complex alpine environment. ISPRS J. Photogramm. Remote Sens. 2006, 60, 323–338. [Google Scholar] [CrossRef]
  30. Lillesand, T.M.; Kiefer, R.W.; Chipman, J.W. Remote Sensing and Image Interpretation; John Wiley & Sons: Hoboken, NJ, USA, 2008. [Google Scholar]
  31. Hollaus, M.; Mandlburger, G.; Pfeifer, N. Land Cover Dependent Derivation of Digital Surface Models from Airborne Laser Scanning Data. IAPRS: Saint-Mandé, France, 2010; pp. 221–226. [Google Scholar]
  32. Höfle, B.; Mücke, W.; Dutter, M.; Rutzinger, M. International Conference on Applied Geoinformatics. In Detection of Building Regions Using Airborne Lidar—A New Combination of Raster and Point Cloud Based Gis Methods Study Area and Datasets; GI-Forum: Washington, DC, USA, 2009; pp. 66–75. [Google Scholar]
  33. Rutzinger, M.; Höfle, B.; Hollaus, M.; Pfeifer, N. Object-based point cloud analysis of full-waveform airborne laser scanning data for urban vegetation classification. Sensors 2008, 8, 4505–4528. [Google Scholar] [CrossRef]
  34. Heipke, C.; Mayer, H.; Wiedemann, C. 3D Recontruction and Modeling of Topographic Objects. In Evaluation of Automatic Road Extraction; IAPRS: Stuttgart, Germany, 1997; pp. 151–160. [Google Scholar]
  35. Hodgson, M.E.; Jensen, J.R.; Schmidt, L.; Schill, S.; Davis, B. An evaluation of lidar- and ifsar-derived elevation models in leaf-on conditions with usgs level 1 and level 2 dems. Remote Sens. Environ. 2003, 84, 295–308. [Google Scholar] [CrossRef]
  36. Friedrich, A. Airborne laser scanning—Present status and future expectations. ISPRS J. Photogramm. Remote Sens. 1999, 54, 64–67. [Google Scholar] [CrossRef]
  37. Hollaus, M.; Wagner, W.; Schadauer, K.; Maier, B.; Gabler, K. Growing stock estimation for alpine forests in austria: A robust lidar-based approach. Can. J. For. Res. 2009, 39, 1387–1400. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Tran, T.H.G.; Hollaus, M.; Nguyen, B.D.; Pfeifer, N. Assessment of Wooded Area Reduction by Airborne Laser Scanning. Forests 2015, 6, 1613-1627. https://doi.org/10.3390/f6051613

AMA Style

Tran THG, Hollaus M, Nguyen BD, Pfeifer N. Assessment of Wooded Area Reduction by Airborne Laser Scanning. Forests. 2015; 6(5):1613-1627. https://doi.org/10.3390/f6051613

Chicago/Turabian Style

Tran, Thi Huong Giang, Markus Hollaus, Ba Duy Nguyen, and Norbert Pfeifer. 2015. "Assessment of Wooded Area Reduction by Airborne Laser Scanning" Forests 6, no. 5: 1613-1627. https://doi.org/10.3390/f6051613

APA Style

Tran, T. H. G., Hollaus, M., Nguyen, B. D., & Pfeifer, N. (2015). Assessment of Wooded Area Reduction by Airborne Laser Scanning. Forests, 6(5), 1613-1627. https://doi.org/10.3390/f6051613

Article Metrics

Back to TopTop