Characterizing riparian vegetation and classifying riparian extent using airborne laser scanning data

This dissertation addresses the problem of automated vector extraction from Airborne Laser Scanner (ALS) data in urban areas. The recent popularity of Geographic Information Systems (GIS) has stimulated research on automated object extraction in order to simplify the data acquisition and update process. By automatically generating GIS inputs and updates from a single data source, the cost of data acquisition and processing will be kept to a minimum. Compared with other remote sensing techniques, extraction of objects from ALS data is in its infancy. ALS sensor technology has evolved rapidly and now allows the acquisition of very dense point clouds in a short period of time. ALS data is unique in that it explicitly contains 3D information and is acquired from an active sensor. As such, there are several benefits that can be immediately realised by using an ALS only data approach; data acquisition will not be limited to daylight hours as with other sensors and accurate height information is contained in the data. This means that registration of different data sources is not required and as only one data source is used, acquisition costs are minimised. Apart from these facts, ALS data has some other unique properties that have not been utilised to their full potential. An ALS sensor can record both height and intensity information from multiple returns in a single swath. To date, all the available information has rarely been used. For example, the intensity of a laser return has regularly been dismissed as it was considered under-sampled and noisy. This data still contains usable information and until this information is used, the optimum object recognition results will not be achieved; thus the use of as much of this information as possible is a major focus of this thesis. As ALS is an explicit 3D data source, the early stages of development were primarily focused on topographic mapping of terrain in forested areas in order to generate Digital Terrain Models (DTMs). As sensor technology has improved, so has the achievable resolution of point clouds from ALS data, and methods to extract objects from stand-alone ALS data have emerged. Attempts have been made to create city models and maps from ALS data that included buildings, roads, trees and powerlines. Each of these spatial object types has unique attributes, which means that automatic map creation is not an easy task. For example, buildings in general are easily detected in ALS data but the building outline is not easily delineated. Another difficulty with building detection is the separation of buildings and bridges as they have many of the same properties as observed by an ALS system. Bridges can usually be found in a road network but road extraction techniques typically produce poor detection rates and often require existing data and / or user interaction in semi-automatic techniques. These simple examples highlight the complexity of automatically generating vectorised maps from ALS data or in fact any data source. In this thesis, new methods are presented for the automatic creation of vectorised maps from ALS data. A two-step processing paradigm is adopted for this purpose, namely the classification of the ALS data and the vectorisation of the classification results. A classification strategy is introduced that creates a hierarchy for object detection with respect to the ALS data itself. This approach develops an ontology between the spatial object classes and the ALS data. The hierarchical framework highlights the fact that one object might not be discernable within the data without considering another. New classification algorithms are then presented within this framework. Each algorithm attempts to exploit the attributes of the data that are consistent within the spatial object class being considered as described by the classification framework. New algorithms for the classification of roads, trees and powerlines are all introduced whilst an extension to existing building classification methods is presented. Once classification is complete, a vectorisation process specific to the task at hand can be employed to yield vectorised results. These developed vectorisation processes are new and include an algorithm that has been generalised to allow the vectorisation of thick lines in images by detecting the centreline, direction and width. The primary goal of this thesis is to present a framework of new algorithms that will allow automatic spatial object detection and vectorisation whilst providing results of an acceptable quality. The algorithms presented rely solely on ALS data and require minimal operator knowledge. Each algorithm has been designed to exploit the way in which the object exhibits itself in the data. The new algorithms are integrated into a software package called JTD (Join The Dots) that will facilitate the effective automatic processing of ALS data. The results of the new algorithms have been evaluated over four 2 x 2 km areas that have been sampled with medium resolution ALS data. The results for each area are displayed and analysed to show the applicability of the whole process in an exemplary way.

ABSTRACTWe examine the nearest neighbor (NN) imputation of species-specific logwood volumes using airborne laser scanning (ALS) data and aerial images. We compare different remote sensing (RS) data combinations as predictor variables in an area-based prediction of logwood volumes using separate training and validation data. We include multispectral leaf-on ALS data, bi-temporal leaf-off ALS data and aerial images in the analyses. Two response configurations are used in the NN imputations: (1) simultaneous imputation in which species-specific logwood volumes are response variables, and (2) separate imputation by tree species in which the attributes of one tree species at a time are response variables. Although an unrealistic alternative in practical implementation, the combination of leaf-on and leaf-off ALS metrics as predictors proved to be the most successful RS data combination, according to the RMSE values associated with the predicted species-specific and dominant logwood volumes. The results showed that older leaf-off ALS data perform well in combination with leaf-on ALS data. In general, predictive performance was better with simultaneous imputation than with separate imputation by tree species. Our finding promotes an awareness of how best to utilize various RS data in future forest inventories.

A critical component of the forest ecosystem, the understorey supports the vast majority of wildlife habitat and total ecosystem floristic diversity. Remote sensing data have been developed to provide information at different scales for surveys of forest ecosystems, but obtaining information about the understorey remains a challenge. As rapid and efficient tools for forest structure attribute estimation, Airborne Laser Scanning (ALS) and Terrestrial Laser Scanning (TLS) have attracted much attention. We examine the relationship between ALS and TLS data and detect changes in the forest understorey caused by forest-tending events in the study area. We conducted trials in five plots within a young Khasi pine (Pinus kesiya Royle ex Gord.) plantation in Yunnan province, China, before and after forest tending. We collected bi-temporal ALS data in this area and TLS data from 10 scanning stations. Canopy height profiles were retrieved from ALS and TLS data, and understorey material volume retrieved from filled TLS voxels volume agreed well with the understorey point clouds percentile distribution (PD) obtained from ALS data. The PD value for the understorey from ALS multiplied by the percentage of ALS return points in the overstorey had a stronger correlation (R2 = 0.90) with the TLS-derived understorey material volume than did the ALS PD value for the understorey only (R2 = 0.71). Taking the effect of the overstorey into consideration will improve evaluations of forest understorey parameters from ALS data. This study confirmed the potential of TLS as a validation tool to assess the accuracy of forest understorey material volume estimation at large scales using ALS data.

Improving species-specific plot volume estimates based on airborne laser scanning and image data using alpha shape metrics and balanced field data

Determination of forest fuels characteristics in mortality-affected Pinus forests using integrated hyperspectral and ALS data

Bridging spatio-temporal gaps in ALS data using Landsat time series and forest disturbance-recovery metrics via multi-task neural networks

Tree biomass estimate is essential for carbon accounting, bioenergy feasibility studies, and forest sustainable management. This fact, added to the availability of airborne laser scanning (ALS) information, provided by the Spanish National Plan for Aerial Orthophotography (PNOA), and the existence of little research focusing on the use of ALS technology in Mediterranean Aleppo pine (Pinus halepensis Mill.) forest, determined the main objective of this research. Thus, this study aims to test the suitability of the low point density (0.5 points/m2), discrete, multiple-return, PNOA-ALS data, to estimate and map the total biomass (TB) and its carbon content in Pinus halepensis Mill. forest stands, located in Aragon (north-eastern Spain). TB was calculated in 45 field plots, using allometric equations, and related through a multivariate linear regression analysis with a collection of independent variables extracted from the ALS data. The predictive model was validated using a leave-one-out cross-validation (LOOCV) technique. Then, a regular grid with cell size 25 x 25 m corresponding to the sample plot size was generated by means of GIS, in order to compute TB at stand level and convert biomass to carbon by using the 0.5 conversion factor. The maximum height, kurtosis and the percentage of returns above 1 meter, were the ALS metrics included in the fitted model, which presented a R2 value of 0.89. The implementation of the model in a GIS showed an average of 68633 kg/ha of TB and 34247.95 kg/ha of carbon fixed. The results indicate that despite the low point density of the ALS data, the final model is accurate enough to be used in forestry applications.

More than 50% of the national lands in Japan have been surveyed by airborne laser scanning (ALS) data with different point densities; and developing an effective approach to take full advantage of these ALS data for forest management has thus become an urgent topic of study. This study attempted to assess the utility of ALS data for individual tree detection and species classification in a mixed forest with a high canopy density. For comparison, two types of tree tops and tree crowns in the study area were delineated by the individual tree crown (ITC) approach using the green band of the orthophoto imagery and the digital canopy height model (DCHM) derived from the ALS data, respectively. Then, the two types of tree crowns were classified into four classes—Pinus densiflora (Pd), Chamaecyparis obtusa (Co), Larix kaempferi (Lk), and broadleaved trees (Bl)—by a crown-based classification approach using different combinations of the three orthophoto bands with intensity and slope maps as follows: RGB (red, green and blue); RGB and intensity (RGBI); RGB and slope (RGBS); and RGB, intensity and slope (RGBIS). Finally, the tree tops were annotated with species attributes from the two best-classified tree crown maps, and the number of different tree species in each compartment was counted for comparison with the field data. The results of our study suggest that the combination of RGBIS yielded greater classification accuracy than the other combinations. In the tree crown classifications delineated by the green band and DCHM data, the improvements in the overall accuracy compared to the RGB ranged from 5.7% for the RGBS to 9.0% for the RGBIS and from 8.3% for the RGBS to 11.8% for the RGBIS. The laser intensity and slope derived from the ALS data may be valuable sources of information for tree species classification, and in terms of distinguishing species for the detection of individual trees, the findings of this study demonstrate the advantages of using DCHM instead of optical data to delineate tree crowns. In conclusion, the synthesis of individual tree delineation using DCHM data and species classification using the RGBIS combination is recommended for interpreting forest resources in the study area. However, the usefulness of this approach must be verified in future studies through its application to other forests.

The impact of geolocation uncertainty on GEDI tropical forest canopy height estimation and change monitoring

Classification of tree species or species classes is still a challenge for remote sensing-based forest inventory. Operational use of Airborne Laser Scanning (ALS) data for prediction of forest variables has this far been dominated by area-based methods where laser scanning data have been used for estimation of forest variables within raster cells. Classification of tree species has however not been achieved with sufficient accuracy with area-based methods using only ALS data. Furthermore, analysis of tree species at the level of raster cells with typical size of 15 m × 15 m is not ideal in the case of mixed species stands. Most ALS systems for terrestrial mapping use only one wavelength of light. New multispectral ALS systems for terrestrial mapping have recently become operational, such as the Optech Titan system with wavelengths 1550 nm, 1064 nm, and 532 nm. This study presents an alternative type of area-based method for classification of tree species classes where multispectral ALS data are used in combination with small raster cells. In this “mini raster cell method” features for classification are derived from the intensity of the different wavelengths in small raster cells using a moving window average approach to allow for a heterogeneous tree species composition. The most common tree species in the Nordic countries are Pinus sylvestris and Picea abies, constituting about 80% of the growing stock volume. The remaining 20% consists of several deciduous species, mainly Betula pendula and Betula pubescens, and often grow in mixed forest stands. Classification was done for pine (Pinus sylvestris), spruce (Picea abies), deciduous species and mixed species in middle-aged and mature stands in a study area located in hemi-boreal forest in the southwest of Sweden (N 58°27’, E 13°39’). The results were validated at plot level with the tree species composition defined as proportion of basal area of the tree species classes. The mini raster cell classification method was slightly more accurate (75% overall accuracy) than classification with a plot level area-based method (68% overall accuracy). The explanation is most likely that the mini raster cell method is successful at classifying homogenous patches of tree species classes within a field plot, while classification based on plot level analysis requires one or several heterogeneous classes of mixed species forest. The mini raster cell method also results in a high-resolution tree species map. The small raster cells can be aggregated to estimate tree species composition for arbitrary areas, for example forest stands or area units corresponding to field plots.

Comparing patterns in forest stand structure following variable harvests using airborne laser scanning data

Airborne laser scanner (ALS) data provide an enhanced capability to remotely map two key variables in forestry: leaf area index (LAI) and tree height (H). Nevertheless, the cost, complexity and accessibility of this technology are not yet suited for meeting the broad demands required for estimating and frequently updating forest data. Here we demonstrate the capability of alternative solutions based on the use of low-cost color infrared (CIR) cameras to estimate tree-level parameters, providing a cost-effective solution for forest inventories. ALS data were acquired with a Leica ALS60 laser scanner and digital aerial imagery (DAI) was acquired with a consumer-grade camera modified for color infrared detection and synchronized with a GPS unit. In this paper we evaluate the generation of a DAI-based canopy height model (CHM) from imagery obtained with low-cost CIR cameras using structure from motion (SfM) and spatial interpolation methods in the context of a complex canopy, as in forestry. Metrics were calculated from the DAI-based CHM and the DAI-based Normalized Difference Vegetation Index (NDVI) for the estimation of tree height and LAI, respectively. Results were compared with the models estimated from ALS point cloud metrics. Field measurements of tree height and effective leaf area index (LAIe) were acquired from a total of 200 and 26 trees, respectively. Comparable accuracies were obtained in the tree height and LAI estimations using ALS and DAI data independently. Tree height estimated from DAI-based metrics (Percentile 90 (P90) and minimum height (MinH)) yielded a coefficient of determination (R2) = 0.71 and a root mean square error (RMSE) = 0.71 m while models derived from ALS-based metrics (P90) yielded an R2 = 0.80 and an RMSE = 0.55 m. The estimation of LAI from DAI-based NDVI using Percentile 99 (P99) yielded an R2 = 0.62 and an RMSE = 0.17 m2/m−2. A comparative analysis of LAI estimation using ALS-based metrics (laser penetration index (LPI), interquartile distance (IQ), and Percentile 30 (P30)) yielded an R2 = 0.75 and an RMSE = 0.14 m2/m−2. The results provide insight on the appropriateness of using cost-effective 3D photo-reconstruction methods for targeting single trees with irregular and heterogeneous tree crowns in complex open-canopy forests. It quantitatively demonstrates that low-cost CIR cameras can be used to estimate both single-tree height and LAI in forest inventories.

Modelling top height growth and site index using repeated laser scanning data

This paper investigates the possibilities of improving aboveground forest biomass and basal area weighted mean height estimates from optical multispectral satellite data by adding Canopy Height Models (CHMs) obtained from matching multiple view-angle satellite data. The analysis was carried out using data collected over the Remningstorp test site in southern Sweden from 2008–2011 and used training and validation data from airborne laser scanning and field data. CHMs were produced by subtracting a Digital Elevation Model (based on airborne laser scanning data) from the Digital Surface Models created by matching multiple view-angle SPOT-5 HRS and ALOS PRISM data. By modeling biomass and height using regression analysis on multispectral data from SPOT-5 HRG in combination with height metrics from the CHMs, an improved root mean squared error (RMSE) was attained, compared with using the individual satellite data sources alone. A comparison between SPOT-5 HRS and ALOS PRISM CHMs in combination with multispectral data was made at stand level using biomass and height estimates from laser scanning data as reference data. For biomass, the relative RMSE improved from 32.9% when using only multispectral data, to 29.2% and 22.4% when adding the CHM from SPOT-5 HRS and ALOS PRISM, respectively. The corresponding improvements for height were from 16.1% to 15.3% with the SPOT-5 HRS CHM and to 12.9% with the ALOS PRISM CHM. A further analysis of combining the ALOS PRISM CHM and multispectral data was made at sub-stand level with field measurements as reference data. This combination gave a relative RMSE of 20.6% for biomass and 10.5% for height. In conclusion, the estimation accuracy for aboveground biomass and basal area weighted mean height was improved by adding CHM data to multispectral data from optical satellites.

This study aimed to map the stem biomass of an even-aged eucalyptus plantation in southeastern Brazil based on canopy height profile (CHPs) statistics using wall-to-wall discrete return airborne laser scanning (ALS), and compare the results with alternative maps generated by ordinary kriging interpolation from field-derived measurements. The assessment of stem biomass with ALS data was carried out using regression analysis methods. Initially, CHPs were determined to express the distribution of laser point heights in the ALS cloud for each sample plot. The probability density function (pdf) used was the Weibull distribution, with two parameters that in a secondary task, were used as explanatory variables to model stem biomass. ALS metrics such as height percentiles, dispersion of heights, and proportion of points were also investigated. A simple linear regression model of stem biomass as a function of the Weibull scale parameter showed high correlation (adj.R2 = 0.89). The alternative model considering the 30th percentile and the Weibull shape parameter slightly improved the quality of the estimation (adj.R2 = 0.93). Stem biomass maps based on the Weibull scale parameter doubled the accuracy of the ordinary kriging approach (relative root mean square error = 6 % and 13 %, respectively).