Forest mapping, one of the main components of performing a forest inventory, is an important driving force in the development of laser scanning. Mobile laser scanning (MLS), in which laser scanners are installed on moving platforms, has been studied as a convenient measurement method for forest mapping in the past several years. Positioning and attitude accuracies are important for forest mapping using MLS systems. Inertial Navigation Systems (INSs) and Global Navigation Satellite Systems (GNSSs) are typical and popular positioning and attitude sensors used in MLS systems. In forest environments, because of the loss of signal due to occlusion and severe multipath effects, the positioning accuracy of GNSS is severely degraded, and even that of GNSS/INS decreases considerably. Light Detection and Ranging (LiDAR)-based Simultaneous Localization and Mapping (SLAM) can achieve higher positioning accuracy in environments containing many features and is commonly implemented in GNSS-denied indoor environments. Forests are different from an indoor environment in that the GNSS signal is available to some extent in a forest. Although the positioning accuracy of GNSS/INS is reduced, estimates of heading angle and velocity can maintain high accurate even with fewer satellites. GNSS/INS and the LiDAR-based SLAM technique can be effectively integrated to form a sustainable, highly accurate positioning and mapping solution for use in forests without additional hardware costs. In this study, information such as heading angles and velocities extracted from a GNSS/INS is utilized to improve the positioning accuracy of the SLAM solution, and two information-aided SLAM methods are proposed. First, a heading angle-aided SLAM (H-aided SLAM) method is proposed that supplies the heading angle from GNSS/INS to SLAM. Field test results show that the horizontal positioning accuracy of an entire trajectory of 800 m is 0.13 m and is significantly improved (by 70%) compared to that of a traditional GNSS/INS; second, a more complex information added SLAM solution that utilizes both heading angle and velocity information simultaneously (HV-aided SLAM) is investigated. Experimental results show that the horizontal positioning accuracy can reach a level of six centimetres with the HV-aided SLAM, which is a significant improvement (by 86%). Thus, a more accurate forest map is obtained by the proposed integrated method.
The detailed structure information under the forest canopy is important for forestry surveying. As a high-precision environmental sensing and measurement method, terrestrial laser scanning (TLS) is widely used in high-precision forestry surveying. In TLS-based forestry surveys, stem-mapping, which is focused on detecting and extracting trunks, is one of the core data processing tasks and the basis for the subsequent calculation of tree attributes; one of the most basic attributes is the diameter at breast height (DBH). This article explores and improves the methods for stem mapping and DBH estimation from TLS data. Firstly, an improved 3D stem mapping algorithm considering the growth direction in random sample consistency (RANSAC) cylinder fitting is proposed to extract and fit the individual tree point cloud section. It constructs the hierarchical optimum cylinder of the trunk and introduces the growth direction into the establishment of the backbone buffer in the next layer. Experimental results show that it can effectively remove most of the branches and reduce the interference of the branches to the discrimination of trunks and improve the integrity of stem extraction by about 36%. Secondly, a robust least squares ellipse fitting method based on the elliptic hypothesis is proposed for DBH estimation. Experimental results show that the DBH estimation accuracy of the proposed estimation method is improved compared with other methods. The mean root mean squared error (RMSE) of the proposed estimation method is 1.14 cm, compared with other methods with a mean RMSE of 1.70, 2.03, and 2.14 cm. The mean relative accuracy of the proposed estimation method is 95.2%, compared with other methods with a mean relative accuracy of 92.9%, 91.9%, and 90.9%.
The rapid movement of low Earth orbit (LEO) satellite can improve geometric diversity, which contributes to the rapid convergence of Global Navigation Satellite System (GNSS) precise point positioning (PPP). However, the LEO onboard receiver clock cannot be used directly by PPP users as the LEO satellite clock because the LEO onboard receiver clock and LEO satellite clock absorb different code delays when receiving and transmitting signals. In this study, a real-time estimation approach for the LEO satellite clock based on ground tracking stations was proposed for the first time. The feasibility for the rapid convergence of the LEO satellite clock was analyzed using the satellite time dilution of precision (TDOP) that one satellite is relative to multiple ground tracking stations. The LEO constellation of 168 satellites and observations for 15 ground tracking stations were simulated to verify the proposed method. The experiment results showed that the average convergence time was 31.21 min for the Global Positioning System (GPS) satellite clock, whereas the value for the LEO satellite clock was only 2.86 min. The average root mean square (RMS) and standard deviation (STD) values after convergence were 0.71 and 0.39 ns for the LEO satellite clock, whereas the values were 0.31 and 0.13 ns for the GPS satellite clock. The average weekly satellite TDOP for the LEO satellite was much smaller than that for the GPS satellite. The average satellite TDOPs for all LEO and GPS satellites were 19.13 and 1294.70, respectively. However, the average delta TDOPs caused by satellite motion for all LEO and GPS satellites were both 0.10. Therefore, the rapid convergence of the LEO satellite clock resulted from the better geometric distribution of the LEO satellite relative to ground stations. Despite errors and the convergence time of the LEO satellite clock, the convergence time and positioning accuracy for LEO-augmented GPS and BeiDou Navigation Satellite System (BDS) PPP with the real-time estimated LEO satellite clock can still reach 10.63 min, 1.94 cm, 1.44 cm, and 4.18 cm in the east, north, and up components, respectively. The improvements caused by LEO satellite for GPS/BDS PPP were 59%, 30%, 31%, and 33%, respectively.
Simultaneous localization and mapping (SLAM) has been investigated in the field of robotics for two decades, as it is considered to be an effective method for solving the positioning and mapping problem in a single framework. In the SLAM community, the Extended Kalman Filter (EKF) based SLAM and particle filter SLAM are the most mature technologies. After years of development, graph-based SLAM is becoming the most promising technology and a lot of progress has been made recently with respect to accuracy and efficiency. No matter which SLAM method is used, loop closure is a vital part for overcoming the accumulated errors. However, in 2D Light Detection and Ranging (LiDAR) SLAM, on one hand, it is relatively difficult to extract distinctive features in LiDAR scans for loop closure detection, as 2D LiDAR scans encode much less information than images; on the other hand, there is also some special mapping scenery, where no loop closure exists. Thereby, in this paper, instead of loop closure detection, we first propose the method to introduce extra control network constraint (CNC) to the back-end optimization of graph-based SLAM, by aligning the LiDAR scan center with the control vertex of the presurveyed control network to optimize all the poses of scans and submaps. Field tests were carried out in a typical urban Global Navigation Satellite System (GNSS) weak outdoor area. The results prove that the position Root Mean Square (RMS) error of the selected key points is 0.3614 m, evaluated with a reference map produced by Terrestrial Laser Scanner (TLS). Mapping accuracy is significantly improved, compared to the mapping RMS of 1.6462 m without control network constraint. Adding distance constraints of the control network to the back-end optimization is an effective and practical method to solve the drift accumulation of LiDAR front-end scan matching.
For GPS medium-long baseline real-time kinematic (RTK) positioning, the troposphere parameter is introduced along with coordinates, and the model is ill-conditioned due to its strong correlation with the height parameter. For BeiDou Navigation Satellite System (BDS), additional difficulties occur due to its special satellite constellation. In fact, relative zenith troposphere delay (RZTD) derived from high-precision empirical zenith troposphere models can be introduced. Thus, the model strength can be improved, which is also called the RZTD-constrained RTK model. In this contribution, we first analyze the factors affecting the precision of BDS medium-long baseline RTK; thereafter, 15 baselines ranging from 38 km to 167 km in different troposphere conditions are processed to assess the performance of RZTD-constrained RTK. Results show that the troposphere parameter is difficult to distinguish from the height component, even with long time filtering for BDS-only RTK. Due to the lack of variation in geometry for the BDS geostationary Earth orbit satellite, the long convergence time of ambiguity parameters may reduce the height precision of GPS/BDS-combined RTK in the initial period. When the RZTD-constrained model was used in BDS and GPS/BDS-combined situations compared with the traditional RTK, the standard deviation of the height component for the fixed solution was reduced by 52.4% and 34.0%, respectively.
The precision of railway map is becoming a significant issue for autonomous train scheduling, monitoring and maintenance, related location-based service (LBS), and further ensuring travel safety. Mobile 3D laser scanning is an efficient method for making relative high-precision railway track maps, particularly during the night period of railway maintenance, for light detection and ranging (LiDAR) can work without ambient light. In this paper, we propose an efficient and accurate railway track vectorization method based on the LiDAR point clouds from the self-built train Mobile Laser Scanning (MLS) system. Our method takes full use of railway track geometry and reflection intensity feature of LiDAR, without any trajectory prior information. Firstly, clear track points are filtered by intensity; then, a K-means clustering fused Region-Grow Fitting algorithm is applied. It can not only extract the line vector of railway track, but also can tell the track branches apart, especially on bends and turnout. Experiments were carried on using point clouds with an average density of 490 points per square meter. The experimental results show that the method not only can quickly extract linear objects such as railway track and catenary, but also can detect the railways even in complex real-world topologies such as at bends and turnouts. The precision of the detection area in bends and turnouts are 90.32% and 81.31% respectively, the sensitivity is 83.27% and 83.33%, respectively. Moreover, it can identify the track networks.