Analysis of Factors Affecting Random Measurement Error in LiDAR Point Cloud Feature Matching Positioning

Measured surface-atmosphere fluxes of energy (sensible heat, H, and latent heat, LE) and CO 2 ( FCO 2) represent the “true” flux plus or minus potential random and systematic measurement errors. Here, we use data from seven sites in the AmeriFlux network, including five forested sites (two of which include “tall tower” instrumentation), one grassland site, and one agricultural site, to conduct a cross-site analysis of random flux error. Quantification of this uncertainty is a prerequisite to model-data synthesis (data assimilation) and for defining confidence intervals on annual sums of net ecosystem exchange or making statistically valid comparisons between measurements and model predictions. We differenced paired observations (separated by exactly 24 h, under similar environmental conditions) to infer the characteristics of the random error in measured fluxes. Random flux error more closely follows a double-exponential (Laplace), rather than a normal (Gaussian), distribution, and increase as a linear function of the magnitude of the flux for all three scalar fluxes. Across sites, variation in the random error follows consistent and robust patterns in relation to environmental variables. For example, seasonal differences in the random error for H are small, in contrast to both LE and FCO 2, for which the random errors are roughly three-fold larger at the peak of the growing season compared to the dormant season. Random errors also generally scale with R n ( H and LE) and PPFD ( FCO 2). For FCO 2 (but not H or LE), the random error decreases with increasing wind speed. Data from two sites suggest that FCO 2 random error may be slightly smaller when a closed-path, rather than open-path, gas analyzer is used.

Light Detection and Ranging (LiDAR) systems have been widely adopted for the acquisition of dense and accurate topographic data over extended areas. Although the utilization of this technology has increased in different applications, the development of standard methodologies for the quality control (QC) of LiDAR data has not followed the same trend. In other words, a lack in reliable, practical, cost-effective, and commonly acceptable QC procedures is evident. A frequently adopted procedure for QC is comparing the LiDAR data to ground control points. Aside from being expensive, this approach is not accurate enough for the verification of horizontal accuracy, unless specifically designed LiDAR targets are used. This paper is dedicated to providing accurate, economical, and convenient internal QC procedures for the evaluation of LiDAR data, which is captured from parallel flight lines. The underlying concept of the proposed methodologies is that, in the absence of systematic and random errors in system parameters and measurements, conjugate surface elements in overlapping strips should perfectly match each other. Consistent incompatibilities and the quality of fit between conjugate surface elements in overlapping strips can be used to detect systematic errors in the system parameters/measurements and to evaluate the noise level in the LiDAR point cloud, respectively. Experimental results from real data demonstrate that all the proposed methods, with one exception, produce compatible estimates of systematic discrepancies between the involved data sets, as well as good quantification of inherent noise.

Accurate vehicle ego-localization is key for autonomous vehicles to complete high-level navigation tasks. The state-of-the-art localization methods adopt visual and light detection and ranging (LiDAR) simultaneous localization and mapping (SLAM) to estimate the position of the vehicle. However, both of them may suffer from error accumulation due to long-term running without loop optimization or prior constraints. Actually, the vehicle cannot always return to the revisited location, which will cause errors to accumulate in Global Navigation Satellite System (GNSS)-challenged environments. To solve this problem, we proposed a novel localization method with prior dense visual point cloud map constraints generated by a stereo camera. Firstly, the semi-global-block-matching (SGBM) algorithm is adopted to estimate the visual point cloud of each frame and stereo visual odometry is used to provide the initial position for the current visual point cloud. Secondly, multiple filtering and adaptive prior map segmentation are performed on the prior dense visual point cloud map for fast matching and localization. Then, the current visual point cloud is matched with the candidate sub-map by normal distribution transformation (NDT). Finally, the matching result is used to update pose prediction based on the last frame for accurate localization. Comprehensive experiments were undertaken to validate the proposed method, showing that the root mean square errors (RMSEs) of translation and rotation are less than 5.59 m and 0.08°, respectively.

Abstract. Utilizing ground control points (GCPs) to georeference photogrammetry-based point cloud data is a common practice in unmanned aerial system (UAS) mapping. Direct georeferencing or integrated sensor orientation (ISO) can be used to obtain georeferenced point clouds from UAS without relying heavily on GCPs. However, the accuracy of the point cloud may be impacted by the accuracy of the trajectory solution obtained by GNSS. To improve point cloud accuracy, post-processing kinematic (PPK) solutions can be applied to the UAS trajectory, which may provide higher accuracy than low-accuracy trajectory solutions and minimize the reliance on GCPs. This study compares the accuracy and precision of two different point clouds generated using different methods. One point cloud was generated using traditional photogrammetric methods with low accuracy Global Navigation Satellite System (GNSS) observations from the UAS and GCPs that have an average accuracy of one to two centimeters, while the other was generated using PPK trajectory solution for the UAS’s trajectory with two software: open-source Emlid Studio and the widely used Inertial Explorer. The use of PPK techniques in UAS mapping may have several potential benefits over traditional methods. By correcting the errors in the UAS's trajectory, a user may only need to depend on fewer ground control points, which can reduce the time and cost associated with fieldwork. This is particularly useful in areas that are difficult to access or have limited ground control point options, such as in urban or forested areas. To evaluate performance, a GNSS receiver is used to obtain measurements on checkpoints, which are used to assess the accuracies of the point clouds. In our experiments, the accuracy of the point clouds generated using PPK trajectory solution with high accuracy GCPs was found to be higher than those generated with low accuracy GNSS observations while aided with high accuracy ground control points. While the use of PPK with GCPs is generally expected to provide more accurate and reliable data than low-accuracy GNSS observations even after adjusting with GCPs, the number and distribution of GCPs can still significantly impact overall accuracy. Therefore, careful consideration of the number of GCPs and their placement is essential to achieve the desired level of efficiency and effectiveness in UAS mapping.

Routine bridge inspections require labor intensive and highly subjective visual interpretation to determine bridge deck surface condition. Light Detection and Ranging (LiDAR) a relatively new class of survey instrument has become a popular and increasingly used technology for providing as-built and inventory data in civil applications. While an increasing number of private and governmental agencies possess terrestrial and mobile LiDAR systems, an understanding of the technology’s capabilities and potential applications continues to evolve. LiDAR is a line-of-sight instrument and as such, care must be taken when establishing scan locations and resolution to allow the capture of data at an adequate resolution for defining features that contribute to the analysis of bridge deck surface condition. Information such as the location, area, and volume of spalling on deck surfaces, undersides, and support columns can be derived from properly collected LiDAR point clouds. The LiDAR point clouds contain information that can provide quantitative surface condition information, resulting in more accurate structural health monitoring. LiDAR scans were collected at three study bridges, each of which displayed a varying degree of degradation. A variety of commercially available analysis tools and an independently developed algorithm written in ArcGIS Python (ArcPy) were used to locate and quantify surface defects such as location, volume, and area of spalls. The results were visual and numerically displayed in a user-friendly web-based decision support tool integrating prior bridge condition metrics for comparison. LiDAR data processing procedures along with strengths and limitations of point clouds for defining features useful for assessing bridge deck condition are discussed. Point cloud density and incidence angle are two attributes that must be managed carefully to ensure data collected are of high quality and useful for bridge condition evaluation. When collected properly to ensure effective evaluation of bridge surface condition, LiDAR data can be analyzed to provide a useful data set from which to derive bridge deck condition information.

A Light Detection and Ranging (LiDAR) sensor mounted on an Unmanned Aerial Vehicle (UAV) can map the overflown environment in point clouds. Mapped canopy heights allow for the estimation of crop biomass in agriculture. The work presented in this paper contributes to sensory UAV setup design for mapping and textual analysis of agricultural fields. LiDAR data are combined with data from Global Navigation Satellite System (GNSS) and Inertial Measurement Unit (IMU) sensors to conduct environment mapping for point clouds. The proposed method facilitates LiDAR recordings in an experimental winter wheat field. Crop height estimates ranging from 0.35–0.58 m are correlated to the applied nitrogen treatments of 0–300 . The LiDAR point clouds are recorded, mapped, and analysed using the functionalities of the Robot Operating System (ROS) and the Point Cloud Library (PCL). Crop volume estimation is based on a voxel grid with a spatial resolution of 0.04 × 0.04 × 0.001 m. Two different flight patterns are evaluated at an altitude of 6 m to determine the impacts of the mapped LiDAR measurements on crop volume estimations.

Currently, 3D point clouds are obtained via LiDAR (Light Detection and Ranging) sensors to compute vegetation parameters to enhance agricultural operations. However, such a point cloud is intrinsically dependent on the GNSS (global navigation satellite system) antenna used to have absolute positioning of the sensor within the grove. Therefore, the error associated with the GNSS receiver is propagated to the LiDAR readings and, thus, to the crown or orchard parameters. In this work, we first describe the error propagation of GNSS over the laser scan measurements. Second, we present our proposal to overcome this effect based only on the LiDAR readings. Such a proposal uses a scan matching approach to reduce the error associated with the GNSS receiver. To accomplish such purpose, we fuse the information from the scan matching estimations with the GNSS measurements. In the experiments, we statistically analyze the dependence of the grove parameters extracted from the 3D point cloud -specifically crown surface area, crown volume, and crown porosity- to the localization error. We carried out 150 trials with positioning errors ranging from 0.01 meters (ground truth) to 2 meters. When using only GNSS as a localization system, the results showed that errors associated with the estimation of vegetation parameters increased more than 100 % when positioning error was equal or bigger than 1 meter. On the other hand, when our proposal was used as a localization system, the results showed that for the same case of 1 meter, the estimation of orchard parameters improved in 20 % overall. However, in lower positioning errors of the GNSS, the estimation of orchard parameters were improved up to 50% overall. These results suggest that our work could lead to better decisions in agricultural operations, which are based on foliar parameter measurements, without the use of external hardware.

NRLI-UAV: Non-rigid registration of sequential raw laser scans and images for low-cost UAV LiDAR point cloud quality improvement

A flexible targetless LiDAR–GNSS/INS–camera calibration method for UAV platforms

The variation in point cloud density is driven by many different factors. This variation is expected to affect the quality of the information extracted from the point clouds, however, the extent to which these variations impact the ability to accurately extract and assess geometric features of highways from point cloud data are unknown. This paper investigates the impacts of point density reduction on the extraction and assessment of four critical geometric features. The density of light detection and ranging (LiDAR) data was first reduced and the different features were extracted at varying levels of point density and on a selection of different highway segments in Alberta, Canada. The information obtained at lower point density was then compared to what was obtained at 100% point density. It was found that clearance assessments and sight distance assessments had low sensitivity to reductions in point density (i.e., reducing the point density to as low as 10% of the original data (30 ppm2 on the pavement surface) yielded results comparable to what was obtained at 100% density (300 ppm2) In contrast, for cross section slope estimation and curve attribute estimation higher sensitivity to point density was observed. These findings are critical for transportation agencies considering the adoption of LiDAR technology to manage elements of their infrastructure and for researchers developing data processing tools and algorithms for the semantic segmentation of transportation features from remotely sensed point clouds.

During the past few decades, the presence of global navigation satellite systems (GNSSs) such as GPS, GLONASS, Beidou and Galileo has facilitated positioning, navigation and timing (PNT) for various outdoor applications. With the rapid increase in the number of orbiting satellites per GNSS, enhancements in the satellite-based augmentation systems (SBASs) such as EGNOS and WAAS, as well as commissioning new GNSS constellations, the PNT capabilities are maximized to reach new frontiers. Additionally, the recent developments in precise point positioning (PPP) and real time kinematic (RTK) algorithms have provided more feasibility to carrier-phase precision positioning solutions up to the third-dimensional localization. With the rapid growth of internet of things (IoT) applications, seamless navigation becomes very crucial for numerous PNT dependent applications especially in sensitive fields such as safety and industrial applications. Throughout the years, GNSSs have maintained sufficiently acceptable performance in PNT, in RTK and PPP applications however GNSS experienced major challenges in some complicated signal environments. In many scenarios, GNSS signal suffers deterioration due to multipath fading and attenuation in densely obscured environments that comprise stout obstructions. Recently, there has been a growing demand e.g. in the autonomous-things domain in adopting reliable systems that accurately estimate position, velocity and time (PVT) observables. Such demand in many applications also facilitates the retrieval of information about the six degrees of freedom (6-DOF - x, y, z, roll, pitch, and heading) movements of the target anchors. Numerous modern applications are regarded as beneficiaries of precise PNT solutions such as the unmanned aerial vehicles (UAV), the automatic guided vehicles (AGV) and the intelligent transportation system (ITS). Hence, multi-sensor fusion technology has become very vital in seamless navigation systems owing to its complementary capabilities to GNSSs. Fusion-based positioning in multi-sensor technology comprises the use of multiple sensors measurements for further refinement in addition to the primary GNSS, which results in high precision and less erroneous localization. Inertial navigation systems (INSs) and their inertial measurement units (IMUs) are the most commonly used technologies for augmenting GNSS in multi-sensor integrated systems. In this article, we survey the most recent literature on multi-sensor GNSS technology for seamless navigation. We provide an overall perspective for the advantages, the challenges and the recent developments of the fusion-based GNSS navigation realm as well as analyze the gap between scientific advances and commercial offerings. INS/GNSS and IMU/GNSS systems have proven to be very reliable in GNSS-denied environments where satellite signal degradation is at its peak, that is why both integrated systems are very abundant in the relevant literature. In addition, the light detection and ranging (LiDAR) systems are widely adopted in the literature for its capability to provide 6-DOF to mobile vehicles and autonomous robots. LiDARs are very accurate systems however they are not suitable for low-cost positioning due to the expensive initial costs. Moreover, several other techniques from the radio frequency (RF) spectrum are utilized as multi-sensor systems such as cellular networks, WiFi, ultra-wideband (UWB) and Bluetooth. The cellular-based systems are very suitable for outdoor navigation applications while WiFi-based, UWB-based and Bluetooth-based systems are efficient in indoor positioning systems (IPS). However, to achieve reliable PVT estimations in multi-sensor GNSS navigation, optimal algorithms should be developed to mitigate the estimation errors resulting from non-line-of-sight (NLOS) GNSS situations. Examples of the most commonly used algorithms for trilateration-based positioning are Kalman filters, weighted least square (WLS), particle filters (PF) and many other hybrid algorithms by mixing one or more algorithms together. In this paper, the reviewed articles under study and comparison are presented by highlighting their motivation, the methodology of implementation, the modelling utilized and the performed experiments. Then they are assessed with respect to the published results focusing on achieved accuracy, robustness and overall implementation cost-benefits as performance metrics. Our summarizing survey assesses the most promising, highly ranked and recent articles that comprise insights into the future of GNSS technology with multi-sensor fusion technique.

Differential LiDAR (Light Detection and Ranging) from repeated surveys has recently emerged as an effective tool to measure three-dimensional (3D) change for applications, such as quantifying slip and spatially distributed warping associated with earthquake ruptures, and examining the spatial distribution of beach erosion after hurricane impact. Currently, the primary method for determining 3D change from LiDAR is through the use of the iterative closest point (ICP) algorithm and its variants. However, all current studies using ICP have assumed that all LiDAR points in the compared point clouds have uniform accuracy. This assumption is simplistic given that the error for each LiDAR point is variable, and dependent upon time varying factors such as target range, angle of incidence, and aircraft trajectory accuracy. Therefore, to rigorously determine spatial change, it would be ideal to model the random error for every LiDAR observation in the differential point cloud, and use these error estimates as apriori weights in the ICP algorithm. To test this approach, we implemented a rigorous LiDAR observation error propagation method to generate estimated random error for each point in a LiDAR point cloud, and then determine 3D displacements between two point clouds using an anisotropic weighted ICP (A-ICP) algorithm. The algorithm was evaluated by qualitatively and quantitatively comparing point clouds with synthetic fault ruptures between a uniform weight and anistropically weighted ICP algorithm. Then post-earthquake slip is estimated for the 2010 El Mayor-Cucapah Earthquake (EMC), using pre- and post-event LiDAR. Based on the analysis, Moving Window A-ICP is able to better estimate the synthetic surface ruptures, and provides a smoother estimate of actual displacement for the EMC earthquake.

Abstract. For the azimuth observation to be made at its magnetic observatories routinely, Japan Meteorological Agency (JMA) has adopted a traditional method based on Polaris sighting. Due to its difficulty to implement under overcast weather conditions and to its demand on observers, for overtime work into the evening, we are motivated to seek for an alternative method based on the Global Navigation Satellite Systems (GNSS) observation that might potentially relieve those two disadvantages. An experiment is made at Kakioka to assess the eligibility and effectiveness of the GNSS method for JMA's unmanned observatories, Memambetsu and Kanoya. The GNSS observations themselves achieve as high a precision as approximately 1 arcsec, as far as they are analyzed with Static mode. Derived from the results of GNSS observation and some supplementary horizontal angle measurements, the azimuth of the azimuth mark for the absolute measurement is determined with a precision of a few arcsecond, which is comparable to the azimuth precision achieved by the Polaris sighting. However, we end up with their significant difference by about 10 arcsec. We discuss this discrepancy to be possibly due to a local geoid gradient. The Polaris observation is made with a theodolite tilted in the gravitational direction, also known as the vertical line deviation, whereas the GNSS observations are based on the azimuth of the compliant ellipsoid plane.

The study of climate change is an important field of research. Monitoring of atmospheric variability especially the tropospheric precipitable water vapor (PWV) is a powerful way to investigate climate change. Global navigation satellite systems (GNSS) provide a good tool for studying atmospheric parameters as GNSS signals along its path from the satellites to the ground based receivers suffer a significant delay due to the refractivity of earth’s atmosphere. GNSS signals are not only delayed but also refracted in the neutral atmosphere. The zenith wet delays (ZWD) caused by the troposphere can be estimated during the geodetic processing of GNSS signals. Since the ZWD is tightly correlated to the PWV, GNSS observations can be used to study the changes of PWV. In this study GNSS data from the Egyptian permanent GNSS network (EPGN), International GNSS Service (IGS), EUREF Permanent Network (EPN) and Scripps Orbit and Permanent Array Center (SOPAC), for the years 2013 and 2014, were used. These GNSS stations provided a good geometrical coverage over the respective region. All GNSS data were processed using Bernese V 5.2 software. The needed product from GNSS data processing is the tropospheric zenith total delay (ZTD) over each GNSS site. The ZTD is divided based on physical parameters into zenith hydrostatic delay (ZHD) and zenith wet delay (ZWD). The ZWD is the basic observable used to calculate PWV. The values of the PWV were calculated and its variation over the study area was investigated. The quality of calculated PWV values from GNSS data evaluated against traditional Radiosonde (RS) measurements. The results of GNSS PWV showed good agreement with RS PWV when the PWV values were low, but this agreement became worse at high PWV values as the differences between the two techniques increased. The correlation coefficient between RS PWV and GNSS PWV varied from 0.31 to 0.84. Standard deviation of the differences between RS PWV and GNSS PWV ranged from 2.369 to 5.973 mm. The PWV estimated from GNSS observations had annual cycle, and its cycle in 2013 was different from that in 2014. The PWV differences between the 2 years clarified that the water vapor content over the east Mediterranean in most of 2014 days was higher than that in 2013. Furthermore, PWV variations were noted on both temporal and spatial scales. The highest temporal variation value was 25.41 mm whereas the maximum value of the spatial variation was 19.67 mm. The present study illustrated the importance of using geodetic networks to provide atmospheric information.

This paper presents a new algorithm for building an indoor map by integrating point clouds of 2D light detection and ranging (LIDAR) scanners in indoor environments. Iterative closest point (ICP) algorithm is one of the well-known methods for such purpose and often used for mobile robot SLAM. However, the algorithm is designed based on dense (or continuous) measurement of the same space with known relative positions of measurement points and angles, and it does not often work efficiently if the measurement is sparse and/or LIDAR locations are unknown. Such situations are seen when some installed LIDARs in a room are used to build a background indoor map for object tracking, or mobile LIDARs are used by technicians to build a digital indoor map, where each space is captured only at a few locations with different angles. To tackle this issue, our method extracts line segments and edge points as features from given LIDAR point clouds and finds shape coincidences commonly contained in a pair of given point clouds to identify positional relationships between LIDARs. By this information, these point clouds can be integrated into common 2D coordinates. Indoor map generation is realized by sequentially applying this integration procedure to every pair of point clouds. Also, the experiments on real data show that our method can identify the relative positions of LIDARs with 10cm-order errors in average, and by sequentially applying the point-cloud integration, the generated maps have only 3% errors.