Search

International audience; Highlights:• A new method for LAD estimation from T-LiDAR in occluded areas is developed.• This method use kriging on LAD estimator (LAD-Kriging) and relies on TLS data only.• LAD-Kriging can be applied to all voxels whatever their reliability.• LAD-Kriging retrieves reliable estimates in not-explored and poorly-sampled voxels.• This method can be fitted to other unbiased estimators extracting metrics in point clouds.Abstract:Terrestrial Laser Scanning (TLS) has been used during the past decade to capture the complexity of 3D forest canopy structures, especially Leaf or Plant Area Density (LAD/PAD). TLS data, i.e. point cloud, can be divided into voxels to estimate the three-dimensional distribution of LAD/PAD. However, the combination effects of vegetation occlusion and shooting pattern of TLS scanners lead to a highly heterogeneous sampling, which limits the reliability of some local estimates, since several voxels are either not explored or poorly-explored by laser beams. In practice, recommendations vary regarding the minimum number of beams crossing voxels or the minimum path lengths required to provide reliable predictions. In addition, assigning a value to non-explored and poorly-explored voxels is still an open question. The present work proposes a new method, called LAD-kriging, to mitigate the impact of non-uniform sampling and to increase the accuracy of LAD estimates in non-explored and poorly-explored voxels. The method takes advantage of i) an unbiased LAD estimator of known variance, which was recently developed; ii) the spatial correlation of the LAD field, which derives from vegetation clumping. LAD-kriging computes kriging weights from mathematical derivations, which takes into account both spatial dependencies in the LAD field and the reliability of the estimate available in each voxel. It was evaluated through numerical experiments, which enabled to validate the algorithm and to evaluate its performance through comparison with true references. An example application to field data shows that such spatial correlations truly exist in the field and that LAD-kriging entails to reduce sampling errors (with respect to full resolution scanning). Although a real validation is impossible due to the lack of precise references at the voxel scale, this example allows to gain confidence for its application to field data. In our realistic numerical experiment, up to 25-30% of voxels can be explored by

The amount and spatial distribution of foliage in a tree canopy have fundamental functions in ecosystems as they affect energy and mass fluxes through photosynthesis and transpiration. They are usually described by the Leaf Area Index (LAI) and the Leaf Area Density (LAD), which can be measured through a variety of methods, including voxel-based methods applied to LiDAR point clouds. A theoretical study recently compared the numerical errors arising from different voxel-based estimation methods for Plant Area Density (PAD) based on Beer’s law-based, contact frequency and Maximum-Likelihood Estimation, showing that the bias-corrected Maximum Likelihood Estimator was theoretically the most efficient. However, this earlier study i) ignored wood volumes; ii) neglected vegetation clumping inside the voxel; iii) ignored instrument characteristics in terms of effective footprint, iv) was limited to a single viewpoint. In practice, retrieving LAD from PAD is not straightforward, vegetation is not randomly distributed in volumes of interest, beams are divergent and forestry plots are usually sampled from more than one viewpoint, to mitigate the effect of occlusion. In the present short communication, we extend the previous efficient formulation to actual field conditions to i) account for the presence of both wood volumes and wood hits, ii) rigorously include correction terms for vegetation and instrument characteristics, iii) integrate multiview data. A numerical comparison with other methods commonly used to combine information from different viewpoints led to error reduction, especially in poorly-explored volumes, which are frequent in actual canopies. Beyond its concision, completeness and efficiency, this new formulation -which can be applied to multiview TLS, but also UAV LiDAR scanning - can help reducing errors in LAD estimation.

International audience; Highlights:• Stand-scale estimates of leaf density from voxelized LiDAR data can be biased.• We identified a negative bias caused by the oversampling of light voxels.• Its magnitude depends on location, vegetation structure, scan design and voxel size• Voxel-scale measurement accuracy is too low in small voxels.• 0.5 m voxel size is the best option because corrections are more straightforward.Abstract:The need for fine scale description of vegetation structure is increasing as Leaf Area Density (LAD, m(2)/m(3)) becomes a critical parameter to understand ecosystem functioning and energy and mass fluxes in heterogeneous ecosystems. Terrestrial Laser Scanning (TLS) has shown great potential for retrieving the foliage area at stand, plant or voxel scales. Several sources of measurement errors have been identified and corrected over the past years. However, measurements remain sensitive to several factors, including, 1) voxel size and vegetation structure within voxels, 2) heterogeneity in sampling from TLS instrument (occlusion and shooting pattern), the consequences of which have been seldom analyzed at the scale of forest plots. In the present paper, we aimed at disentangling biases and errors in plot-scale measurements of LAD with TLS in a simulated vegetation scene. Two negative biases were formerly attributed to (i) the unsampled voxels and to (ii) the subgrid vegetation heterogeneity (i.e. clumping effect), and then quantified, thanks to a the simulation experiment providing known LAD references at voxel scale, vegetation manipulations and unbiased point estimators. We used confidence intervals to evaluate voxel-scale measurement accuracy. We found that the unsampled voxel effect (i) led to underestimations with the "mean layer" method -commonly used to fill unsampled voxels- for small voxels (0.1-0.2 m) and/or low number of scans (

International audience; Estimating leaf and plant area density with Terrestrial Laser Scanners (TLS) continues to be more and more popular, as tridimensional point clouds appear as an appealing measurement technique for heterogeneous environments. Some approaches implement a discretization of the point cloud in a grid (referred to as “voxelbased”) to account for this vegetation heterogeneity and significant work has been done in this recent research field, but no general theoretical analysis is available. Although estimators have been proposed and several causes of biases have been identified, their unbiasedness (zero bias) and efficiency (smallest error) have not been evaluated. Also, confidence intervals are almost never provided. In the present paper, we assumed that the vegetation elements were randomly distributed within voxels and that TLS beams were infinitely thin, in order to focus on the remaining sources of biases and errors. In this simplified context, we both solve the transmittance equation and use the Maximum Likelihood Estimation (MLE), to derive some new estimators for the attenuation coefficient, which is proportional to leaf area density at voxel scale in this idealized context. These estimators include bias corrections and confidence intervals, and account for the number of beams crossing the voxel (beam number), the inequality of path lengths in voxel, the size of vegetation elements, as well as for the variability of element positions between vegetation samples. These theoretical derivations are complemented by numerous numerical simulations for the evaluation of estimator bias and efficiency, as well as the assessment of the coverage probabilities of confidence intervals. Our simulations reveal that the usual estimators are biased and exhibit 95% confidence intervals on the order of±100% of the estimate, when the beam number is smaller than 30. Second, our bias-corrected estimators -especially the bias-corrected MLE- are truly unbiased and efficient in a wider range of validity than the usual ones, even for beam number as low as 5. Third, we found that the confidence intervals can be as high as ≈ ± 50% when the projected area of a single element was on the order of 10% of voxel cross-sectional area and vegetation was dense (optical depth of the voxel equal to 2), even for a beam number larger than 1000. This is explained by the variability of element positions between vegetation samples, which implies that a significant part of residual error is caused by random effects. When LAD estimates are averaged over several small voxels -typically to determine a vertical profile at plot scale or to compute the total leaf area of a single plant-, confidence intervals are typically on the order of±5 to 10% with bias-corrected estimators, which is twice as small as with usual estimators.Our study provides some new ready-to-use estimators and confidence intervals for attenuation coefficients,which are unbiased and efficient within a fairly large range of parameter values. The unbiasedness is achieved for a fairly low beam number, which is promising for application to airborne LiDAR data. They permit to raise the level of understanding and confidence on LAD estimation. Among other applications, their usage should help determine the most suitable voxel size, for given vegetation types and scanning density, whereas existing guidelines are highly variable among studies, probably because of differences in vegetation, scanning design and estimators. The impact of other sources of biases and errors, such as vegetation heterogeneity inside voxels or TLS specifications are not addressed in the present manuscript and would require further investigations.

International audience; Reliable measurements of the 3D distribution of Leaf Area Density (LAD) in forest canopy are crucial for describing and modelling microclimatic and eco-physiological processes involved in forest ecosystems functioning. To overcome the obvious limitations of direct measurements, several indirect methods have been developed, including methods based on Terrestrial LiDAR scanning (TLS). This work focused on various LAD estimators used in voxel-based approaches. LAD estimates were compared to reference measurements at branch scale in laboratory, which offered the opportunity to investigate in controlled conditions the sensitivity of estimations to various factors such as voxel size, distance to scanner, leaf morphology (species), type of scanner and type of estimator. We found that all approaches to retrieve LAD estimates were highly sensitive to voxel size whatever the species or scanner and to distance to the FARO scanner. We provided evidence that these biases were caused by vegetation heterogeneity and variations in the effective footprint of the scanner. We were able to identify calibration functions that could be readily applied when vegetation and scanner are similar to those of the present study. For different vegetation and scanner, we recommend replicating our method, which can be applied at reasonable cost. While acknowledging that the test conditions in the laboratory were very different from those of the measurements taken in the forest (especially in terms of occlusion), this study revealed existence of strong biases, including spatial biases. Because the distance between scanner and vegetation varies in field scanning, these biases should occur in a similar manner in the field and should be accounted for in voxel-based methods but also in gap-fraction methods. This article belongs to the Special Issue Remote Sensing Techniques for Precision Forestry.