Crowdsourcing RTK: a new GNSS positioning framework for building spatial high-resolution atmospheric maps based on massive vehicle GNSS data

High-quality spatial atmospheric delay correction information is essential for achieving fast integer ambiguity resolution (AR) in precise positioning. However, traditional real-time precise positioning frameworks (i


Introduction
Precise positioning services using Global Navigation Satellite Systems (GNSS) have become an indispensable part of our daily lives.This technology is widely used in various fields, such as disaster monitoring, autonomous driving, the Internet of Things, and public pedestrian navigation.It provides users with graded positioning services with accuracies ranging from millimeters to meters (Li & Chen, 2022;Liu et al., 2020;Leick et al., 2015;Li et al., 2022;Qin et al., 2022;Yang et al., 2021).
Currently, two methods are predominantly used to provide real-time dynamic centimeter-level GNSS positioning services-traditional network RTK (Dingfa et al.Han, 1997;Hu et al., 2003;Vollath et al., 2000;Huang et al. 2022) and the recently developed PPP-RTK (Li et al., 2023;Teunissen & Khodabandeh, 2015;Wabbena et al., 2005;Zhang et al., 2022aZhang et al., , 2022b) ) technology.Both of these methods rely on Continuous Operational Reference System (CORS) base stations for high-precision positioning, which come with high construction and maintenance costs and are mainly funded by national or local governments.Moreover, the distribution of CORS base stations is relatively sparse, which means that both positioning methods require interpolation to obtain atmospheric delay correction data, which can be affected by the accuracy of the ionosphere during its active period.As a result, the interpolation accuracy can decline sharply, leading to poor positioning performance.These issues collectively lead to high costs and low user experience when using these services.
Autonomous vehicles have become more common and are now being mass-produced.These vehicles are often equipped with GNSS receivers to provide timing services and centimeter-level GNSS positioning services (Cao et al., 2022;Li et al., 2021;Qin et al., 2021;Tang et al., 2022), which provide more potential data sources for precise GNSS positioning.It is noted that the density of the road network is much greater than that of the CORS network.There is a clear imbalance in the effective use of on-board GNSS data on roads.For example, the city of Wuhan has more than 4 million privately owned cars (Security, 2022).During low and peak periods of vehicle travel, approximately 0.8-3.2 million vehicles drive on the road.The density of CORS stations is much less than one ten thousandth of the density of these vehicles.This imbalance shows that we urgently need more intelligent methods to improve the quality and coverage of precise GNSS positioning services while reducing related costs to meet the growing demand for mass precise positioning applications.
Crowdsourcing is a collaborative model that has proven highly effective in various domains, such as autonomous driving (Qin et al., 2021), map creation (Haklay, 2010;Wang et al., 2023), and earthquake warning (Minson et al., 2015).This approach is an innovative way to significantly reduce the expenses associated with gathering professional surveying data by leveraging nonprofessional crowdsourced surveying data.It provides real-time information that is both reliable and cost-effective (Liu et al., 2020).
The purpose of this study is to improve the costeffectiveness and efficiency of positioning services by overcoming the high cost and sparse distribution of traditional CORS base stations and the imbalanced effective utilization of onboard GNSS data on roads.To achieve this goal, we propose a new GNSS positioning framework that uses public GNSS observables to establish a spatial high-resolution atmospheric delay correction map.We also propose the Crowdsourced Atmospheric delay correction Map (CAM), which includes Crowdsourced Ionospheric and Tropospheric delay correction Maps (CIM and CTM).This positioning framework has several unique characteristics.First, it only requires two base stations.Based on this, a set of quality control methods is designed to obtain highprecision atmospheric information at specific times and places.Second, it gathers atmospheric information from the general public and constructs the CAM.Finally, the CAM will further augment the user positioning performance.With the application of this positioning framework, the accuracy of the CAM and the user location performance will mutually improve each other, forming a virtuous cycle that will meet the real-time precise positioning needs of autonomous driving users.

System overview
As illustrated in Fig. 1, our framework relies on collaboration among users as well as between users and the cloud, primarily comprising the following three parts.The first part involves the initialization of atmospheric information on the user end, where high-precision and reliable ionospheric and tropospheric delay correction is generated.The users solve three baseline parameters and validate their position and ambiguity closure using dual base stations.Subsequently, they upload the validated position along with atmospheric information to the cloud-server end.
In the second part, the cloud gathers atmospheric information uploaded by the general users, creates the CAM, and initializes the on-cloud atmospheric map.
In the third part, the CAM augments user positioning and is subsequently refined.A new positioning request is made by the user, and the cloud-server end provides the atmospheric information within the grid corresponding to the approximate position of the user.Based on this precise atmospheric delay correction information, the user end utilizes an ionosphere-and troposphereweighted model to achieve rapid and precise positioning.After the atmospheric information is validated through the ambiguity closure test, it is sent back to the cloudserver end for real-time updating of the CAM.
This iterative process creates a positive feedback mechanism that continuously improves the accuracy and spatial resolution of the CAM, maintaining the high precision of user positioning.

Single-baseline RTK functional model
In our positioning framework, we add tropospheric constraints to the traditional long-baseline ionosphereweighted single-differenced (SD) RTK function model (Odolinski et al., 2015) to solve the three baselines (two base stations and one user) separately.The specific form is given as: (1) where P s * ru,j and L s * ru,j represent the SD code and phase observables, respectively.s * = 1 * , 2 * , . . ., m * denotes the tracked satellites; * represents the tracked systems, with G for GPS, E for Galileo, and B for BDS; u and r represent the receivers of the baseline; j * = 1 * , . . ., f * denotes the tracked frequencies; κ s * u is the line-of-sight unit vector; and p ru denotes the coordinate increment of receiver r. dt * ru,P j and dt * ru,L j represent receiver clock error terms that absorb pseudorange hardware and phase hardware delays, respectively.I s * ru represents the relative slant ionospheric delay, with a corresponding coefficient µ j * = 2 j / 2 1 ; j is the wavelength for frequency j * ; and T ru and m s * u represent the relative wet Zenith Tropospheric Delay (ZTD) and corresponding mapping function, respectively.N 1s * ru,j * represents double-differenced integer

Single-baseline RTK stochastic models
The variance-covariance (VCV) matrix for the combined GPS/Galileo/BDS ionosphere-weighted and troposphere-weighted model is given as where σ 2 T ru is the variance of (T ru ) prior and Q * I ru represents the single-system VCV matrix of (I s * ru ) prior .Q * yy represents the single-system VCV matrix of the SD code and phase observables and is given as where the code and phase observable a priori variance factors are given by submatrices , respectively, and blkdiag and diag denote the block diag- onal and diagonal matrix, respectively.Here, we assume no cross-correlation between the code and phase or between frequencies.Furthermore, ⊗ is the Kronecker product, D T n is the between-station SD operator and W −1 m * is the elevation-dependent weighting function.This subsection describes the stochastic models for the SD code and phase observables, and later subsections describe the stochastic models for the SD ionospheric and tropospheric delays under different conditions.

Generating atmospheric information with an ambiguity-fixed solution for individual vehicles
During the initialization process, no prior atmospheric information is provided by the cloud server.For individual vehicles, The corresponding stochastic model is given as where σ 2 I s * ru is the SD variance factor of the slant ionospheric delays that can be modeled as a function of the baseline length.σ 2 T ru is a larger constant, which is equivalent to white noise estimation. (2) We utilize (1)-( 5) to obtain the three baseline solutions constructed by the two base stations and the user and attempt to fix the ambiguity.

Atmospheric quality control based on the ambiguity closure test
Once the ambiguities are fixed, the mutual difference between the user's positions calculated by the dual base stations and the ambiguity closure of the three baselines are tested.The expression of the ambiguity closure test is given as where r 1 , r 2 , and u represent the dual base stations and the user, respectively.
If the test is passed, the user's position, along with ionospheric and tropospheric information (including the adjustment value and posterior variance), is uploaded to the cloud server end.

Positioning after individual vehicle initialization
After the initialization for individual users, two methods for accessing atmospheric information to augment positioning are available.The first method is to employ the user's own historical atmospheric information weighting via random walk estimation, where t k−1 represents the previous epoch, and the corre- sponding stochastic model is given as where (σ 2 T ru ) t k−1 represent the posterior SD slant ionospheric delays and relative wet ZTD of the previous epoch, prn The second method utilizes atmospheric information downloaded from the cloud-sever end to weight: represent the SD slant ionospheric delays and relative wet ZTD returned from the cloud, respectively.The corresponding random model is: where represent the posterior SD slant ionospheric delays and relative wet ZTD returned from the cloud server end, respectively.

Grid division and interpolation
As shown in Fig. 2, to create an atmospheric map for a certain period, we divide the service area into grids of a certain size, such as 1 km × 1 km, which can be adjusted according to the actual requirements.The atmospheric information provided by the user is comprehensively processed, and the information within each grid is (10) obtained.However, we may not be able to obtain all the atmospheric delay information within the grid initially.Therefore, we use interpolation algorithms such as LIM (Han, 1997) to fill in the gaps and interpolate all the atmospheric information within each grid based on the existing data.This helps us create the initial CAM and complete the initialization process.

Storage and communication of the CAM
Taking Wuhan city as an example, the city spans 132.12 km from east to west and 153.The user end and the cloud-server end communicate bidirectionally to exchange information, enabling the server to provide the atmospheric delay correction information as real-time feedback to the users at the grid's approximate location.

Initialization time
The initialization process can be completed rapidly in practice.First, in terms of spatial distribution, users near the base station can easily obtain an ambiguity-fixed solution (Zhu et al., 2023).This means that atmospheric information in that area can be readily initialized and then extrapolated to provide priors for other newly added users.As the system iteration proceeds, a dense atmospheric map can quickly be constructed.
Second, in terms of temporal distribution, for vehicle users, low-traffic periods typically occur in the morning when ionospheric activity tends to be stable, and the ground truth of the SD ionospheric delay is often near zero, making it easier to perform AR and generate reliable atmospheric delay information.

CAM-augmented user positioning and CAM refinement
New users download atmospheric delay correction information corresponding to their approximate positions within the grid from the cloud-server end.By utilizing Eqs. ( 1)-( 3), ( 9) and ( 10) to solve the long-baseline RTK equations with atmospheric constraints, a single-epoch ambiguity-fixed solution can be obtained.Subsequently, the atmospheric information, validated through ambiguity closure checks, is uploaded back to the cloud.
During system iteration, the crowdsourced ionospheric map and crowdsourced tropospheric map will undergo stepwise refinement, leading to the establishment of a more detailed atmospheric delay correction map.

Experimental results and discussion
We first describe the experimental setup in "Experimental setup" section.Then, "Quality of vehicle GNSS data" and "Assessment of the atmospheric quality generated by long-baseline dynamic users" sections evaluate the various stages of the initialization process, including the data quality, position accuracy, AR, and atmospheric solution accuracy.In "Dynamic baseline (40-300 km) generation of the CAM" section, we extend the baseline to 300 km to explore the spatial applicability of the framework.In the "Simulation experiment based on CAM-augmented positioning" and "Multi-vehicle based augmentation in a real environment" sections, we conduct simulation and multi-vehicle based experiments on CAM-based localization after initialization.

Experimental setup
The data were gathered on May 12th, 2023, in Wuhan, covering the urban area from the Second to the Fourth Ring Road.The vehicle test lasted approximately 6 h, and the test trajectory and base station distribution are displayed in Fig. 4. The distance between base stations ranged from 40 to 50 km, and the distance between the base stations and the test trajectory ranged between 0-60 km.We used the traditional network RTK algorithm to calculate the trajectory and obtain the ground truth.During the experiment, we selected dual base stations to create a dynamic triangulation network with onboard users for dynamic long-baseline RTK positioning experiments.Specifically, we chose UBAS and UBHH as the dual base stations for the new positioning framework as the basis for further experiments.
The specific configuration of the test equipment and scenarios is shown in Fig. 5.During the test, we used two GNSS receivers of varying cost-the geodetic receiver Trimble Net-R9 (R9) and the low-cost receiver Septentrio Mosaic X5 (MSX5).The data sampling rate was set to 1 Hz, and we validated the new GNSS positioning framework proposed in this paper using dual-frequency GPS and Galileo data.The test covered multiple typical scenarios, such as elevated roads with broad views, urban roads, and road sections obstructed by trees, tall buildings, and elevated bridges.
The strategies for primary data processing are listed in Table 1.During the initialization process of the system, the study considered only the continuity of ambiguity between epochs, which was estimated as a random walk process.The atmospheric parameters were given a zero mean a priori and were estimated as white noise due to the poor quality of dynamic data compared to static data.This strategy resets the atmospheric information at each epoch.Once the ambiguity closure test, based on dual base stations, is passed, the atmospheric delay at that location is uploaded to the cloud server end to complete the initialization of the CAM.Afterward, users within the service area use the atmospheric delay from the cloud as a prior for localization.If the initialization is not completed and the atmospheric parameters are estimated as

Quality of vehicle GNSS data
This subsection elaborates on satellite visibility, PDOP, and multipath error to assess the feasibility of creating a crowdsourced atmospheric map.
The number of satellites that are available and their corresponding PDOP values have a significant impact on the position and AR performance.During a test conducted with elevation cutoff angles of 15° and 40°, 97.67% and 95% of the periods, respectively, had 10 or more visible satellites, as shown in the first row of Fig. 6.In cases where fewer satellites are visible, there may be severe GNSS electromagnetic signal occlusion under the elevated bridge, but at the same time, there will be a wide field of view on the elevated bridge.The city's ring roads, elevated roads, and other wide roads provide a good spatial distribution of road networks and GNSS observation conditions.This means that in most cases, we will have sufficient satellites for positioning, which creates favorable conditions for establishing a CAM.
Based on the first two rows of Fig. 6, as the elevation cutoff angle increases, the number of available satellites decreases gradually, while the corresponding PDOP increases.Table 2 also shows the HDOP and VDOP values of multiple GNSSs at elevation cutoff angles ranging from 15° to 40°.Notably, as the satellite elevation increases, the DOP value gradually increases.At an elevation cutoff angle of 40°, the HDOP is still only 1.07, whereas the VDOP is 2.85.Thus, at higher elevation cutoff angles, the positioning accuracy of the GNSS in the  horizontal direction is less affected, and the impact in the vertical direction is more significant.In Fig. 6, the third row displays the multipath errors of the R9 and MSX5 receivers and two base stations.The data suggest that the multipath error of static base stations is much lower than that of dynamic users.In dynamic environments, the multipath errors for R9 and MSX5 are approximately the same.This indicates that both R9 and MSX5 can be used to establish a CAM.
It is important to understand that in situations where the availability of GNSS signals is limited, such as in tunnels or in areas with high multipath effects such as urban canyons, even short-baseline RTK positioning can be difficult to achieve.The new GNSS positioning framework that we propose is an exclusive satellite navigation system positioning method, and its highest positioning performance is that of short-baseline RTK positioning.Therefore, the discussion of this study is limited to scenarios where GNSS signals are available and of good quality and where short-baseline RTK positioning can be performed.Situations in which GNSS signals are completely blocked, the signal quality is low, or short-baseline RTK positioning cannot be completed are not within the scope of this study.

Performance analysis of dynamic long-baseline positioning
In this section, we first focus on position accuracy, which is of utmost concern to the general public.We plot two types of receivers to show the difference in position between epochs 442,500 and 444,500, with a sequence of 2000s.Figure 7 shows the epochs that match well; the position difference calculated near epoch 443,100 is relatively large and is not shown.We evaluated the positions of two receivers under long-baseline dynamic positioning and found that both types of receivers, which have different costs, exhibited slight errors compared with the ground truth.Taking the dynamic long baselines constructed by R9, MSX5, and UBAS as examples, the positioning accuracies (STDs) in the three ENU directions are (1.1, 1.18, 4.84) cm and (0.91, 1.06, 4.74) cm, respectively, as shown in the first row of Fig. 7.The third row of  Fig. 7 shows that for a single user, the position difference of the baseline solutions of the dual base stations is relatively small.The mutual difference accuracies of the two types of receivers are (0.89, 0.83, 2.18) cm and (1.02, 1.07, 3.19) cm, respectively.The difference in the single-baseline solutions between R9 and MSX5 is also small, indicating that both types of receivers perform consistently in long-baseline dynamic positioning.
In addition, by comparing the first and third lines of Fig. 7, we observe that the dynamic long-baseline solution provides better consistency in the vertical direction than does the traditional network RTK solution.Although the network RTK user-end model has fewer parameters to estimate, its redundancy is the same as that of the atmosphere-weighted long baseline model.Moreover, its user-end model loses accuracy information when utilizing atmospheric information provided by the server, and the interpolation error of the atmospheric delay is absorbed by the position, resulting in slightly poorer consistency compared to the long baseline atmosphere-weighted model.
This experiment verified the reliability of the dynamic long-baseline positioning of dual base stations.The results further support the potential of using massive vehicle user data to generate atmospheric information.It is reasonable for us to use the mutual difference in the same user's position calculated by dual base stations as the basis for establishing a dynamic reference network.This process will aid in AR and the creation of high-quality atmospheric maps.

Performance of dynamic long-baseline AR
Correct AR is a prerequisite for achieving high-precision positioning at the centimeter level for the GNSS.We plotted the position difference, ADOP, and ratio sequences of the two types of receivers between epochs 442,500 and 442,800 with a sequence of 300 s in Fig. 8.The graph shows that there are three significant fluctuations in the position differences.The ADOP and ratio indicators show significant changes in the first two areas, while the fluctuations are insignificant in the third area.This makes it challenging for dynamic positioning users to determine whether the ambiguity is correctly fixed based on the theoretical indicators of the ADOP and ratio.For this reason, we introduced dual base stations.This allows us to test the user positioning accuracy directly through mutual position differences.Additionally, the ambiguity parameter in the equation of the phase observables is closely correlated with the ionospheric parameters.Incorrect AR will result in corresponding errors in the ionospheric parameters.Therefore, we examine the ambiguity closure of the dynamic triangulation network to determine whether the AR is correct and then identify and eliminate erroneous ionospheric parameter estimations.
In traditional network RTK or PPP-RTK, the server side relies on the correlation between ambiguity and ionospheric delay to refine ionospheric information for initialization.Additionally, the user end uses interpolated atmospheric delay correction to achieve fast integer ambiguity resolution without contributing to the server's atmospheric delay refinement.In contrast, the new GNSS positioning framework enables the user to use the server's atmospheric delay information to achieve instantaneous integer ambiguity resolution while also contributing to the refinement of the server's atmospheric delay.This approach allows us to provide users with atmospheric correction services that offer higher spatial resolution and greater accuracy than interpolated atmospheric correction by CORS stations.

Assessment of the atmospheric quality generated by long-baseline dynamic users
We chose UBAS as the base station, while R9 and MSX5 were selected as rover stations.We then calculated ambiguity-fixed baseline solutions for UBAS-R9 and UBAS-MSX5 separately.Afterward, we differenced the baseline solutions for UBAS-R9 and UBAS-MSX5, which included the position, single-differenced ionospheric delay, and single-differenced tropospheric delay.The results are shown in Fig. 9.
First, the position estimation of dynamic long baselines can converge steadily with centimeter-level accuracy.Second, the ionospheric mutual differences corresponding to most satellites are within 2 cm, but there are still a few cases that exceed this range.This is because the current criterion used to determine the fixed solution of the ionospheric delay relies only on the theoretical index (ratio), which is not entirely reliable.We will further provide more accurate criteria to determine whether the ionospheric delay was correctly estimated.Finally, the mutual differences in the tropospheric delay show a relatively stable trend, which is maintained at the millimeter level.This is due to the correlation between the troposphere and each observed satellite, which has good observability and high redundancy.
In addition, in the static scenario of nearly 200 s between epochs 43,500 and 43,700, there is a smaller difference between the position and atmospheric parameters, which can be effectively attributed to typical vehicle scenarios: users are waiting at traffic lights or experiencing traffic congestion, which provides favorable Figure 10 further illustrates the agreement in the single-differenced ionosphere delay between R9 and MSX5 under different judgment criteria.It can be observed that position constraints do not have a significant effect on the localization and removal of ionospheric gross errors, but ambiguity closure can effectively eliminate satellites with ionospheric solution errors.Taking satellite G08 as an example, combined with the data in Table 3, the percentage of epochs that passed the ambiguity closure test (i.e., the success rate of ambiguity resolution, calculated as the number of epochs that passed the ambiguity closure test divided by the fixed number of epochs) between MSX5 and the dual base stations is 66.42%; in contrast, the percentage for R9 is only 20.83%, which is almost not closed around 443,500 epochs.Therefore, there is a significant difference in the estimated ionospheric delay values around this period.After screening for ambiguity closure, the standard deviation (STD) of To demonstrate how the single-differenced ionosphere delay affects ambiguity resolution, additional simulation experiments are included.An ultrashort baseline comprising CUTA and CUT0, which was collected at Curtin University on January 2, 2023 (with an ionospheric ground truth of 0), was chosen as an example.Figure 11 illustrates the impact of different levels of single-differenced ionospheric noise on the ADOP.Based on prior experience, an ADOP of less than 0.12 cycles corresponds to a success rate of ambiguity resolution greater than 99.9% (Odijk & Teunissen, 2008).A single-differenced ionospheric error of 5 cm is required for ambiguity resolution, while the ionospheric accuracy generated by our dynamic long baseline reaches 0.74 cm.Such high-precision ionospheric information can assist in instantaneous ambiguity resolution for almost any long baseline.
Thus far, we have generated a centimeter-level betweenstation single-differenced ionospheric delay and a  millimeter-level between-station single-differenced tropospheric delay using vehicle dynamic data.

Dynamic baseline (40-300 km) generation of the CAM
This section investigates the feasibility of generating a CAM using dynamic long baseline data of over 100 km.We achieved this by increasing the distance between base stations and dynamic users.We selected base station data that were collected at distances of 200 km and 300 km from the motion trajectory and used precise ephemeris data for longbaseline solutions.We evaluated the most crucial factor of correct AR, as illustrated Fig. 12 and Table 4.The data presented in Fig. 12 reveal two clear trends.First, as the satellite elevation increases, the visibility of the satellite improves significantly, resulting in an increase in the success rate of ambiguity resolution.Second, as the baseline length increases, the success rate generally decreases, although there are some exceptions for the G19 and E33 satellites at low elevations.These two trends suggest that while fixing the dynamic long baseline ambiguities, it is essential to prioritize selecting data with high satellite elevations.This would increase the accuracy and stability of atmospheric parameter solutions.
When the baseline distances are set at 40 km and 200 km and the elevation is 25 degrees or greater, the success rate is greater than 50%.In contrast, for a 300 km baseline, within the same elevation range, the success rate varies from 25% to 50%.Importantly, the success rates for  the low-cost receiver MSX5 and the geodetic receiver R9 are similar, which provides valuable reference for popularizing device selection for the general public in practical applications.
The core of the new framework is that it accurately calculates the atmospheric parameters of a low-cost user dynamic long baseline to complete system initialization.By utilizing the correlation between ambiguity and the ionospheric delay, we can transform this problem into that of answering two key questions: whether the ambiguity can be resolved and whether it is resolved correctly.Figure 12 shows that a dynamic long baseline of less than 300 km has the possibility of correct AR, and the success rate exceeds 25% at a cutoff elevation of 25°.Although this percentage may seem relatively low for individual dynamic users, there is enough redundant information for overall mass users.By processing and merging data effectively, we can enhance the overall quality and coverage of the CAM.

Simulation experiment based on CAM-augmented positioning
The experiments described above verify the initialization of the system.The results show that 12 satellites continuously observed by dynamic users can be combined with dual base stations ambiguity closure tests to generate high-precision atmospheric delay correction information.This section assesses whether the generated atmospheric information can be used to improve the positioning performance and assesses the ambiguity resolution of users for a single-epoch solution.The primary experimental settings are shown in Table 5.Other parameters, such as the elevation and SNR masks, remain the same as those in Table 1.
Figure 13 shows the position, ADOP, and ratio sequences of UBHH-R9 under different atmospheric constraints.The results indicate that simultaneously constraining ionospheric and tropospheric delays is the most effective approach.This suggests that our proposed framework can utilize low-cost receivers for generating atmospheric products, enabling fast integer ambiguity resolution in precise positioning.Constraining only the ionosphere also yields good positioning augmentation, while there is almost no significant gain when only the tropospheric delay is constrained.We observed noteworthy fluctuations in the position and ADOP sequences, particularly at approximately 443,800 s, coinciding with the transition of the vehicle's dynamic conditions from static to dynamic.This phenomenon was particularly pronounced when employing only the ionospheric constraint strategy.The observed instability can be attributed to the smaller number of prior parameters in this scheme compared to that of scheme one, leading to challenges when transitioning between epochs.Hence, in dynamic environments, it becomes imperative to meticulously consider and optimize the configuration of prior parameters to ensure the stability of the position solution.

Multi-vehicle based augmentation in a real environment
We further conduct multi-vehicle experiments to evaluate the enhanced positioning effects of atmospheric information generated by different vehicles.The experimental setup is consistent with real-world applications.
The data were collected on March 13, 2024.The vehicle trajectory was within 5 km of the reference station JFNG, and the ground-truth trajectory was provided by shortbaseline RTK/INS.Dual base stations HB01 and HB05, located 30 km and 40 km away from the vehicle, respectively, were selected.The experimental configuration is shown in Fig. 14, and the strategy settings are detailed in Table 6.We continue to compare the four atmospheric augmentation schemes from "Simulation experiment based on CAM-augmented positioning" section.

Fig. 1
Fig. 1 Illustration of the system structure.The system consists of three parts.The first part is the atmosphere initialization for individual vehicles.The second part is the atmosphere initialization for the cloud-sever end.The third part is the CAM augmenting user positioning and CAM refinement and t is the time interval between epochs.
97 km from north to south.It can be roughly divided into a grid of 133 × 154 (1 km × 1 km) squares, each of which has four float attributes: the values and variances of the single-differenced ionospheric delay and single-differenced tropospheric delay.The storage requirement for individual grid data is 16 bytes.At the same time, approximately 20 observable satellites contribute to the creation of the CAM, with a data size of approximately 12 MB.These data are stored on the cloud-server end.
Figure 3 illustrates this process.

Fig. 3
Fig. 3 CAM refinement and user positioning with CAM augmentation

Fig. 4
Fig. 4 Vehicle experiment test trajectory and base station distribution

Fig. 6
Fig.6The first and second lines show the numbers of visible satellites and PDOP at an elevation cutoff angle of 15° to 40° for the rover station, respectively.The third line is the multipath error of the base station and the rover station

Fig. 7
Fig. 7 Differences in dynamic long-baseline positioning.The first two lines show differences between the long-baseline solutions and network RTK of two rover stations and dual base stations.The third line calculates the mutual differences between the same rover stations from the baselines of the dual base stations

Fig. 8
Fig. 8 Theoretical indicators of the AR and position differences

Fig. 9
Fig. 9 Consistency between the single-differenced ionosphere and troposphere under various dynamic conditions

Fig. 11
Fig. 11 Effect of different single-differenced ionospheric noise on ambiguity resolution

Fig. 12
Fig. 12 Statistics of ambiguity resolution success rates for dynamic long baselines from 40 to 300 km

Table 1
Data processing strategies used in this study

Table 2
Percentage of visible satellites (NSat ≥ 10) and average DOP at elevation cutoff angles of 15° to 40° for rover stations

Table 3
Percentage passed the ambiguity closure test within 2000 epochs

Table 4
Ambiguity resolution success rates for dynamic long baselines from 40 to 300 km

Table 5
Main data processing strategies Fig. 13Single-epoch ambiguity-fixed position differences, ratio, and ADOP sequences of UBHH-R9 for four schemes