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.

As the field of high-precision applications based on carriers continues to expand, the development of low-cost, small, modular receivers and their application in diverse scenarios and situations with complex data quality has increased the requirements of carrier-phase data preprocessing. A new geometry-based cycle slip detection and repair method based on Global Position System (GPS) + BeiDou Navigation Satellite System (BDS) is proposed. The method uses a Time-differenced Carrier Phase (TDCP) model, which eliminates the Inner-System Bias (ISB) between GPS and BDS, and it is conducive to the effective combination of GPS and BDS. It avoids the interference of the noise of the pseudo-range with cycle slip detection, while the cycle slips are preserved as integers. This method does not limit the receiver frequency number, and it is applicable to single-frequency data. The process is divided into two steps to detect and repair cycle slip. The first step is cycle slip detection, using the Improved Local Analysis Method (ILAM) to find satellites that have cycle slips; The second step is to repair the cycle slips, including estimating the float solution of changes in ambiguities at the satellites that have cycle slips with the least squares method and the integer solution of the cycle slips by rounding. In the process of rounding, in addition to the success probability, a decimal test is carried out to validate the result. Finally, experiments with filed test data are carried out to prove the effectiveness of this method. The results show that the detectable cycle slips number with GPS + BDS is much greater than that with GPS. The method can also detect the non-integer outliers while fixing the cycle slip. The maximum decimal bias in repair is less than that with GPS. It implies that this method takes full advantages of multi-system.