Single tree biomass modelling using airborne laser scanning

Accurate forest biomass mapping methods would provide the means for e.g. detecting bioenergy potential, biofuel and forest-bound carbon. The demand for practical biomass mapping methods at all forest levels is growing worldwide, and viable options are being developed. Airborne laser scanning (ALS) is a promising forest biomass mapping technique, due to its capability of measuring the three-dimensional forest vegetation structure. The objective of the study was to develop new methods for tree-level biomass estimation using metrics derived from ALS point clouds and to compare the results with field references collected using destructive sampling and with existing biomass models. The study area was located in Evo, southern Finland. ALS data was collected in 2009 with pulse density equalling approximately 10 pulses/m². Linear models were developed for the following tree biomass components: total, stem wood, living branch and total canopy biomass. ALS-derived geometric and statistical point metrics were used as explanatory variables when creating the models. The total and stem biomass root mean square error per cents equalled 26.3% and 28.4% for Scots pine (Pinus sylvestris L.), and 36.8% and 27.6% for Norway spruce (Picea abies (L.) H. Karst.), respectively. The results showed that higher estimation accuracy for all biomass components can be achieved with models created in this study compared to existing allometric biomass models when ALS-derived height and diameter were used as input parameters. Best results were achieved when adding field-measured diameter and height as inputs in the existing biomass models. The only exceptions to this were the canopy and living branch biomass estimations for spruce. The achieved results are encouraging for the use of ALS-derived metrics in biomass mapping and for further development of the models.     

Using low-density discrete Airborne Laser Scanning data to assess the potential carbon dioxide emission in case of a fire event in a Mediterranean pine forest

The aim of study is to map the carbon dioxide (CO₂) emission of the aboveground tree biomass (AGB) in case of a fire event. The suitability of low point density, discrete, multiple-return, Airborne Laser Scanning (ALS) data and the influence of several characteristics of these data and the study area on the results obtained have been evaluated. A sample of 45 circular plots representative of Pinus halepensis Miller stands were used to fit and validate the model of AGB. The ALS point clouds were processed to obtain the independent variables and a multivariate linear regression analysis between field data and ALS-derived variables allowed estimation of AGB. Then, the influence of several characteristics on the residuals of the model was analyzed. Finally, conversion factors were applied to obtain the CO₂ values. The AGB model presented a R² value of 0.84 with a relative root-mean-square error of 27.35%. This model included ALS variables related to vegetation height variability and to canopy density. Terrain slope, aspect, canopy cover, scan angle and the number of laser returns did not influence AGB estimations at plot level.     

Effect of flying altitude and pulse repetition frequency on laser scanner penetration rate for digital elevation model generation in a tropical forest

In tropical forests, the penetration ability of airborne laser scanning (ALS) may be limited because of highly dense vegetation cover. However, in the typical planning of ALS surveys, the ability of laser pulses to penetrate forests is not considered. Nine round-trip flight lines covering the area of a tropical forest on the northeast side of the Tsengwen Reservoir in Taiwan were designed in this study. Five flight lines flew at altitudes of 1.525, 1.830, 2.135, 2.440, and 2.745 km, and the other four had pulse repetition frequencies (PRFs) of 100, 150, 200, and 250 kHz. The laser penetration index (LPI) is a quantitative index measuring the penetration ability of the ALS and consists of the ratio of the number of laser pulses reaching the forest floor to the total number of laser pulses. The LPI was used to represent the laser penetration rate and investigate the influence of flying altitude and PRF on the LPI. The results showed that as the flying altitude decreased by 1 km, the average LPI increased by 10%, and as the PRF decreased by 50 kHz, the average LPI increased by 2%. The effect of the LPI on digital elevation models (DEMs) was confirmed by visual images obtained by DEMs at five altitudes. The DEM obtained at an altitude of 2.745 km was coarsely textured, whereas that obtained at an altitude of 1.525 km was finely textured. The in-situ height data obtained from the electronic Global Navigation Satellite System (eGNSS) were compared with the data of the ALS-generated DEMs. The results indicated that when the LPI ≥60%, the height difference between the in situ data and DEM data was not prominent. However, when the LPI <60%, the ALS-derived DEM could be higher or lower than the in-situ height; the largest difference between the two was 1.7 m. The LPI of a forest should be considered for ALS survey planning, especially when consistent DEM precision for large tropical forest areas is paramount.     

Inter-comparison of remote sensing platforms for height estimation of mango and avocado tree crowns

To support the adoption of precision agricultural practices in horticultural tree crops, prior research has investigated the relationship between crop vigour (height, canopy density, health) as measured by remote sensing technologies, to fruit quality, yield and pruning requirements. However, few studies have compared the accuracy of different remote sensing technologies for the estimation of tree height. In this study, we evaluated the accuracy, flexibility, aerial coverage and limitations of five techniques to measure the height of two types of horticultural tree crops, mango and avocado trees. Canopy height estimates from Terrestrial Laser Scanning (TLS) were used as a reference dataset against height estimates from Airborne Laser Scanning (ALS) data, WorldView-3 (WV-3) stereo imagery, Unmanned Aerial Vehicle (UAV) based RGB and multi-spectral imagery, and field measurements. Overall, imagery obtained from the UAV platform were found to provide tree height measurement comparable to that from the TLS (R² = 0.89, RMSE = 0.19 m and rRMSE = 5.37 % for mango trees; R² = 0.81, RMSE = 0.42 m and rRMSE = 4.75 % for avocado trees), although coverage area is limited to 1–10 km² due to battery life and line-of-sight flight regulations. The ALS data also achieved reasonable accuracy for both mango and avocado trees (R² = 0.67, RMSE = 0.24 m and rRMSE = 7.39 % for mango trees; R² = 0.63, RMSE = 0.43 m and rRMSE = 5.04 % for avocado trees), providing both optimal point density and flight altitude, and therefore offers an effective platform for large areas (10 km²–100 km²). However, cost and availability of ALS data is a consideration. WV-3 stereo imagery produced the lowest accuracies for both tree crops (R² = 0.50, RMSE = 0.84 m and rRMSE = 32.64 % for mango trees; R² = 0.45, RMSE = 0.74 m and rRMSE = 8.51 % for avocado trees) when compared to other remote sensing platforms, but may still present a viable option due to cost and commercial availability when large area coverage is required. This research provides industries and growers with valuable information on how to select the most appropriate approach and the optimal parameters for each remote sensing platform to assess canopy height for mango and avocado trees.     

Quantifying individual tree growth and tree competition using bi-temporal airborne laser scanning data: a case study in the Sierra Nevada Mountains,California

Improved monitoring and understanding of tree growth and its responses to controlling factors are important for tree growth modeling. Airborne Laser Scanning (ALS) can be used to enhance the efficiency and accuracy of large-scale forest surveys in delineating three-dimensional forest structures and under-canopy terrains. This study proposed an ALS-based framework to quantify tree growth and competition. Bi-temporal ALS data were used to quantify tree growth in height (ΔH), crown area (ΔA), crown volume (ΔV), and tree competition for 114,000 individual trees in two conifer-dominant Sierra Nevada forests. We analyzed the correlations between tree growth attributes and controlling factors (i.e. tree sizes, competition, forest structure, and topographic parameters) at multiple levels. At the individual tree level, ΔH had no consistent correlations with controlling factors, ΔA and ΔV were positively related to original tree sizes (R?>?0.3) and negatively related to competition indices (R?R|?>?0.7), ΔV was positively related to original tree sizes (|R|?>?0.8). Multivariate regression models were simulated at individual tree level for ΔH, ΔA, and ΔV with the R² ranged from 0.1 to 0.43. The ALS-based tree height estimation and growth analysis results were consistent with field measurements.     

Terrestrial laser scanning to estimate plot-level forest canopy fuel properties

This paper evaluates the potential of a terrestrial laser scanner (TLS) to characterize forest canopy fuel characteristics at plot level. Several canopy properties, namely canopy height, canopy cover, canopy base height and fuel strata gap were estimated. Different approaches were tested to avoid the effect of canopy shadowing on canopy height estimation caused by deployment of the TLS below the canopy. Estimation of canopy height using a grid approach provided a coefficient of determination of R² = 0.81 and an RMSE of 2.47 m. A similar RMSE was obtained using the 99th percentile of the height distribution of the highest points, representing the 1% of the data, although the coefficient of determination was lower (R² = 0.70). Canopy cover (CC) was estimated as a function of the occupied cells of a grid superimposed upon the TLS point clouds. It was found that CC estimates were dependent on the cell size selected, with 3 cm being the optimum resolution for this study. The effect of the zenith view angle on CC estimates was also analyzed. A simple method was developed to estimate canopy base height from the vegetation vertical profiles derived from an occupied/non-occupied voxels approach. Canopy base height was estimated with an RMSE of 3.09 m and an R² = 0.86. Terrestrial laser scanning also provides a unique opportunity to estimate the fuel strata gap (FSG), which has not been previously derived from remotely sensed data. The FSG was also derived from the vegetation vertical profile with an RMSE of 1.53 m and an R² = 0.87.     

Stochastic gradient boosting classification trees for forest fuel types mapping through airborne laser scanning and IRS LISS-III imagery

This paper presents an application of Airborne Laser Scanning (ALS) data in conjunction with an IRS LISS-III image for mapping forest fuel types. For two study areas of 165 km² and 487 km² in Sicily (Italy), 16,761 plots of size 30-m × 30-m were distributed using a tessellation-based stratified sampling scheme. ALS metrics and spectral signatures from IRS extracted for each plot were used as predictors to classify forest fuel types observed and identified by photointerpretation and fieldwork. Following use of traditional parametric methods that produced unsatisfactory results, three non-parametric classification approaches were tested: (i) classification and regression tree (CART), (ii) the CART bagging method called Random Forests, and (iii) the CART bagging/boosting stochastic gradient boosting (SGB) approach. This contribution summarizes previous experiences using ALS data for estimating forest variables useful for fire management in general and for fuel type mapping, in particular. It summarizes characteristics of classification and regression trees, presents the pre-processing operation, the classification algorithms, and the achieved results. The results demonstrated superiority of the SGB method with overall accuracy of 84%. The most relevant ALS metric was canopy cover, defined as the percent of non-ground returns. Other relevant metrics included the spectral information from IRS and several other ALS metrics such as percentiles of the height distribution, the mean height of all returns, and the number of returns.     

A method for vertical adjustment of digital aerial photogrammetry data by using a high-quality digital terrain model

The accuracy of vertical position information can be degraded by various sources of error in digital aerial photogrammetry (DAP) based point clouds. To address this issue, we propose a relatively straightforward method for automated correction of such point clouds. This method can be used in conjunction with any 3D reconstruction method in which a point cloud is generated from a pair of aerial images. The crux of the method involves separately co-registering each DAP point cloud (formed by the overlap of two or more images) to a common airborne laser scanning (ALS) based digital terrain model. The proposed method has the following essential steps: (1) Ground surface patches are identified in the normalized DAP point clouds by selecting areas in which standard deviation of vertical height is low, (2) height differences between the DAP and ALS point clouds are calculated at these patches, and (3) a correction surface is interpolated from these height differences and is then used to rectify the entire DAP point cloud. The performance of the proposed method is verified using plot data (n = 250) from a forested study area in Eastern Finland. We observed that DAP data from the area corrected using our proposed method resulted in significant increases in prediction accuracy of key forest variables. Specifically, the root mean squared error (RMSE) values for dominant height predictions decreased by up to 23.2%, while the associated model R² values increased by 16.9%. As for stem volume, RMSEs dropped by 20.6%, while the model R² improved by 14.6%, respectively. Hence, prediction accuracies were almost as good as with ALS data. The results suggest that vertically misaligned DAP data, if rectified using an algorithm such as the one presented here, could deliver near ALS data quality at a fraction of the cost.     

Incorporating the effect of ALS-derived DEM uncertainty for quantifying changes due to the landslide in 2011, Mt. Umyeon,Seoul

ABSTRACT

The landslide, which occurred at Umyeon mountain (Mt. Umyeon) in Seoul, Korea in 2011, was a prime example that raised awareness about the landslide in the highly urbanized area. Although many studies have been done on Umyeon landslide, there is a lack of research that detects the area where the landslide occurred and quantifies the elevation changes through remote sensing data. In this regard, this paper aims to detect and assess topographic changes quantitatively over Mt. Umyeon by using digital elevation models (DEMs) derived from airborne laser scanning (ALS) data. Since Mt. Umyeon was hilly and covered with dense trees during summer, traces of the landslide were detected by estimating the spatially distributed uncertainty of ALS-derived DEMs. The probabilistic analysis with Bayes'? theorem considering the spatially distributed DEM of difference (DoD) uncertainty enabled to detect the landslide traces efficiently and was less affected by the influence of ALS errors. The results indicated that ALS-derived DEMs have the potential to detect landslides with their uncertainty estimation, although the ALS data were acquired in hilly and densely vegetated areas. Moreover, quantifying topographic changes due to landslides with high reliability is considered to be beneficial and practically helpful for disaster recovery.     

基于机载激光雷达点云数据提取林木参数方法研究

本文通过黑河流域遥感—地面观测同步试验,获取林木参数,对机载激光雷达与实地观测获取的林木参数进行对比分析,论证了本文提出的基于机载激光雷达点云数据提取林木参数的算法是可行的。试验通过机载激光雷达点云数据,研究由点云数据生成冠层高度模型(CHM),提出从CHM中提取单株木参数(树高、冠幅等)的关键算法;同时,通过在试验区布设1个100m×100m超级样地和16个25m×25m的子样地,利用DGPS和全站仪对单株木进行精确定位与树木参数测量。     

Deriving airborne laser scanning based computational canopy volume for forest biomass and allometry studies

A computational canopy volume (CCV) based on airborne laser scanning (ALS) data is proposed to improve predictions of forest biomass and other related attributes like stem volume and basal area. An approach to derive the CCV based on computational geometry, topological connectivity and numerical optimization was tested with sparse-density, plot-level ALS data acquired from 40 field sample plots of 500–1000 m² located in a boreal forest in Norway. The CCV had a high correspondence with the biomass attributes considered when derived from optimized filtrations, i.e. ordered sets of simplices belonging to the triangulations based on the point data. Coefficients of determination (R²) between the CCV and total above-ground biomass, canopy biomass, stem volume, and basal area were 0.88–0.89, 0.89, 0.83–0.97, and 0.88–0.92, respectively, depending on the applied filtration. The magnitude of the required filtration was found to increase according to an increasing basal area, which indicated a possibility to predict this magnitude by means of ALS-based height and density metrics. A simple prediction model provided CCVs which had R² of 0.77–0.90 with the aforementioned forest attributes. The derived CCVs always produced complementary information and were mainly able to improve the predictions of forest biomass relative to models based on the height and density metrics, yet only by 0–1.9 percentage points in terms of relative root mean squared error. Possibilities to improve the CCVs by a further analysis of topological persistence are discussed.     

Individual tree crown approach for predicting site index in boreal forests using airborne laser scanning and hyperspectral data

Site productivity is essential information for sustainable forest management and site index (SI) is the most common quantitative measure of it. The SI is usually determined for individual tree species based on tree height and the age of the 100 largest trees per hectare according to stem diameter. The present study aimed to demonstrate and validate a methodology for the determination of SI using remotely sensed data, in particular fused airborne laser scanning (ALS) and airborne hyperspectral data in a forest site in Norway. The applied approach was based on individual tree crown (ITC) delineation: tree species, tree height, diameter at breast height (DBH), and age were modelled and predicted at ITC level using 10-fold cross validation. Four dominant ITCs per 400 m² plot were selected as input to predict SI at plot level for Norway spruce (Picea abies (L.) Karst.) and Scots pine (Pinus sylvestris L.). We applied an experimental setup with different subsets of dominant ITCs with different combinations of attributes (predicted or field-derived) for SI predictions. The results revealed that the selection of the dominant ITCs based on the largest DBH independent of tree species, predicted the SI with similar accuracy as ITCs matched with field-derived dominant trees (RMSE: 27.6% vs 23.3%). The SI accuracies were at the same level when dominant species were determined from the remotely sensed or field data (RMSE: 27.6% vs 27.8%). However, when the predicted tree age was used the SI accuracy decreased compared to field-derived age (RMSE: 27.6% vs 7.6%). In general, SI was overpredicted for both tree species in the mature forest, while there was an underprediction in the young forest. In conclusion, the proposed approach for SI determination based on ITC delineation and a combination of ALS and hyperspectral data is an efficient and stable procedure, which has the potential to predict SI in forest areas at various spatial scales and additionally to improve existing SI maps in Norway.     

Integration of Landsat ETM+ Data with Field Measurements for Mapping Leaf Area Index in the Grasslands of Central Kazakhstan

A fine-resolution leaf area index (LAI) data set over a 150 km × 150 km region in central Kazakhstan is retrieved using Landsat ETM+ imagery and ground-based LAI inferred from hemispherical photography and direct measurements. Regression analysis and geostatistics are applied for developing empirical models of LAI from Landsat ETM+ data. The best accuracy is achieved using a model employing a canonical index that combines all the contributions of individual Landsat ETM+ bands into a single index (R ² = 0.67; RMSE = 0.21). This model is then applied for mapping LAI at a regional scale.     

Effect of flying altitude,scanning angle and scanning mode on the accuracy of ALS based forest inventory

Airborne laser scanning (ALS) is a widely used technology in the mapping of environment and forests. Data acquisition costs and the accuracy of the forest inventory are closely dependent on some extrinsic parameters of the ALS survey. These parameters have been assessed in numerous studies about a decade ago, but since then ALS devices have developed and it is possible that previous findings do not hold true with newer technology. That is why, the effect of flying altitudes (2000, 2500 or 3000 m), scanning angles (±15° and ±20° off nadir) and scanning modes (single- and multiple pulses in air) with the area-based approach using a Leica ALS70HA-laser scanner was studied here. The study was conducted in a managed pine-dominated forest area in Finland, where eight separate discrete-return ALS data were acquired. The comparison of datasets was based on the bootstrap approach with 5-fold cross validation. Results indicated that the narrower scanning angle (±15° i.e. 30°) led to slightly more accurate estimates of plot volume (RMSE%: 21–24 vs. 22.5–25) and mean height (RMSE%: 8.5–11 vs. 9–12). We also tested the use case where the models are constructed using one data and then applied to other data gathered with different parameters. The most accurate models were identified using the bootstrap approach and applied to different datasets with and without refitting. The bias increased without refitting the models (bias%: volume 0 ± 10, mean height 0 ± 3), but in most cases the results did not differ much in terms of RMSE%. This confirms previous observations that models should only be used for datasets collected under similar data acquisition conditions. We also calculated the proportions of echoes as a function of height for different echo categories. This indicated that the accuracy of the inventory is affected more by the height distribution than the proportions of echo categories.     

Modeling biophysical properties of broad-leaved stands in the hyrcanian forests of Iran using fused airborne laser scanner data and ultraCam-D images

Inventories of mixed broad-leaved forests of Iran mainly rely on terrestrial measurements. Due to rapid changes and disturbances and great complexity of the silvicultural systems of these multilayer forests, frequent repetition of conventional ground-based plot surveys is often cost prohibitive. Airborne laser scanning (ALS) and multispectral data offer an alternative or supplement to conventional inventories in the Hyrcanian forests of Iran. In this study, the capability of a combination of ALS and UltraCam-D data to model stand volume, tree density, and basal area using random forest (RF) algorithm was evaluated. Systematic sampling was applied to collect field plot data on a 150 m × 200 m sampling grid within a 1100 ha study area located at 36°38′- 36°42′N and 54°24′–54°25′E. A total of 308 circular plots (0.1 ha) were measured for calculation of stand volume, tree density, and basal area per hectare. For each plot, a set of variables was extracted from both ALS and multispectral data. The RF algorithm was used for modeling of the biophysical properties using ALS and UltraCam-D data separately and combined. The results showed that combining the ALS data and UltraCam-D images provided a slight increase in prediction accuracy compared to separate modeling. The RMSE as percentage of the mean, the mean difference between observed and predicted values, and standard deviation of the differences using a combination of ALS data and UltraCam-D images in an independent validation at 0.1-ha plot level were 31.7%, 1.1%, and 84 m³ ha⁻¹ for stand volume; 27.2%, 0.86%, and 6.5 m² ha⁻¹ for basal area, and 35.8%, −4.6%, and 77.9 n ha⁻¹ for tree density, respectively. Based on the results, we conclude that fusion of ALS and UltraCam-D data may be useful for modeling of stand volume, basal area, and tree density and thus gain insights into structural characteristics in the complex Hyrcanian forests.     

Combining UAV-based plant height from crop surface models,visible, and near infrared vegetation indices for biomass monitoring in barley

In this study we combined selected vegetation indices (VIs) and plant height information to estimate biomass in a summer barley experiment. The VIs were calculated from ground-based hyperspectral data and unmanned aerial vehicle (UAV)-based red green blue (RGB) imaging. In addition, the plant height information was obtained from UAV-based multi-temporal crop surface models (CSMs). The test site is a summer barley experiment comprising 18 cultivars and two nitrogen treatments located in Western Germany. We calculated five VIs from hyperspectral data. The normalised ratio index (NRI)-based index GnyLi (Gnyp et al., 2014) showed the highest correlation (R² = 0.83) with dry biomass. In addition, we calculated three visible band VIs: the green red vegetation index (GRVI), the modified GRVI (MGRVI) and the red green blue VI (RGBVI), where the MGRVI and the RGBVI are newly developed VI. We found that the visible band VIs have potential for biomass prediction prior to heading stage. A robust estimate for biomass was obtained from the plant height models (R² = 0.80–0.82). In a cross validation test, we compared plant height, selected VIs and their combination with plant height information. Combining VIs and plant height information by using multiple linear regression or multiple non-linear regression models performed better than the VIs alone. The visible band GRVI and the newly developed RGBVI are promising but need further investigation. However, the relationship between plant height and biomass produced the most robust results. In summary, the results indicate that plant height is competitive with VIs for biomass estimation in summer barley. Moreover, visible band VIs might be a useful addition to biomass estimation. The main limitation is that the visible band VIs work for early growing stages only.     

Two algorithms for extracting building models from raw laser altimetry data

Two new techniques for the determination of building models from laser altimetry data are presented. Both techniques work on the original laser scanner data points without the requirement of an interpolation to a regular grid. Available ground plan information may be used, but is not required. Closed solutions for the determination of the parameters of a standard gable roof type building model based on invariant moments of 2 1/2-D point clouds are shown. In addition, the analysis of deviations between point cloud and model does allow for modelling asymmetries such as dorms on a gable roof. By intersecting planar faces nonparametric buildings with more complex roof types can also be modelled. The techniques were applied to a FLI-MAP laser scanner dataset covering an area of 500×250 m² with a density of more than 5 points/m². Within this region, all but one building could be modelled. An analysis of the variance of the parameters within a group of buildings indicates a precision in the range of 0.1–0.2 m.     

Mapping individual trees with airborne laser scanning data in an European lowland forest using a self-calibration algorithm

Traditional field-based forest inventories tend to be expensive, time-consuming, and cover only a limited area of a forested region. Remote sensing (RS), especially airborne laser scanning (ALS) has opened new possibilities for operational forest inventories, particularly at the single-tree level, and in the prediction of single-tree characteristics. Throughout the world, forests have varying characteristics that necessitate the development of modern, effective, and versatile tools for ALS data processing. To address this need, we aimed to develop a tool for individual tree detection (ITD) utilising a self-calibrating algorithm procedure and to verify its accuracy using the complicated forest structure of near natural forests in the temperate zone.This study was carried out in the Polish part of the Białowieża Forest (BF). The airborne laser scanner (ALS) and color-infrared (CIR) datasets were acquired for more than 60 000 ha. Field-based measurements were performed to provide reference data at the single tree level. We introduced a novel ITD method that is self-calibrated and uses a hierarchical analyses of the canopy height model.There were more than 20 000 000 of trees in first layer in BF above 7 m height. Trees visible from above were divided into coniferous, deciduous and mixed trees that were then matched with an accuracy of 85 %, 85 % and 75 %, respectively. Compared to existing methods, the proposed method is more flexible and achieves better results, especially for deciduous species. Before application of the presented method to other regions, the calibration based on the developed optimisation procedure is needed.     

Modeling Mediterranean forest structure using airborne laser scanning data

The conservation of biological diversity is recognized as a fundamental component of sustainable development, and forests contribute greatly to its preservation. Structural complexity increases the potential biological diversity of a forest by creating multiple niches that can host a wide variety of species. To facilitate greater understanding of the contributions of forest structure to forest biological diversity, we modeled relationships between 14 forest structure variables and airborne laser scanning (ALS) data for two Italian study areas representing two common Mediterranean forests, conifer plantations and coppice oaks subjected to irregular intervals of unplanned and non-standard silvicultural interventions. The objectives were twofold: (i) to compare model prediction accuracies when using two types of ALS metrics, echo-based metrics and canopy height model (CHM)-based metrics, and (ii) to construct inferences in the form of confidence intervals for large area structural complexity parameters.Our results showed that the effects of the two study areas on accuracies were greater than the effects of the two types of ALS metrics. In particular, accuracies were less for the more complex study area in terms of species composition and forest structure. However, accuracies achieved using the echo-based metrics were only slightly greater than when using the CHM-based metrics, thus demonstrating that both options yield reliable and comparable results. Accuracies were greatest for dominant height (Hd) (R² = 0.91; RMSE% = 8.2%) and mean height weighted by basal area (R² = 0.83; RMSE% = 10.5%) when using the echo-based metrics, 99th percentile of the echo height distribution and interquantile distance. For the forested area, the generalized regression (GREG) estimate of mean Hd was similar to the simple random sampling (SRS) estimate, 15.5 m for GREG and 16.2 m SRS. Further, the GREG estimator with standard error of 0.10 m was considerable more precise than the SRS estimator with standard error of 0.69 m.     

A bias-free geodetic boundary value problem approach to height datum unification

A geodetic boundary value problem (GBVP) approach has been formulated which can be used for solving the problem of height datum unification. The developed technique is applied to a test area in Southwest Finland with approximate size of 1.5° × 3° and the bias of the corresponding local height datum (local geoid) with respect to the geoid is computed. For this purpose the bias-free potential difference and gravity difference observations of the test area are used and the offset (bias) of the height datum, i.e., Finnish Height Datum 2000 (N2000) fixed to Normaal Amsterdams Peil (NAP) as origin point, with respect to the geoid is computed. The results of this computation show that potential of the origin point of N2000, i.e., NAP, is (62636857.68 ± 0.5) (m²/s²) and as such is (0.191 ± 0.003) (m) under the geoid defined by W ₀ = 62636855.8 (m²/s²). As the validity test of our methodology, the test area is divided into two parts and the corresponding potential difference and gravity difference observations are introduced into our GBVP separately and the bias of height datums of the two parts are computed with respect to the geoid. Obtaining approximately the same bias values for the height datums of the two parts being part of one height datum with one origin point proves the validity of our approach. Besides, the latter test shows the capability of our methodology for patch-wise application.