# Low-cost precise positioning for automated vehicles

### Carrier-phase differential GNSS produces mass-market centimeter accuracy

**A dense reference network facilitates low-cost**** carrier-phase differential GNSS positioning with rapid integer-ambiguity resolution. This could enable precise lane-keeping for automated vehicles in all weather conditions.**

Strong demand for low-cost precise positioning exists in the mass market. Carrier-phase differential GNSS (CDGNSS) positioning, accurate to within a few centimeters even on a moving platform, would satisfy this demand were its cost significantly reduced. Low-cost CDGNSS would be a key enabler for many demanding consumer applications.

Centimeter-accurate positioning by CDGNSS has been perfected over the past two decades for applications in geodesy, precision agriculture, surveying and machine control. But mass-market adoption of this technology will demand much lower user cost — by a factor of 10 to 100 — yet still require rapid and accurate position fixing. To reduce cost, mass-market CDGNSS-capable receivers will have to make do with inexpensive, low-quality antennas whose multipath rejection and phase center stability are inferior to those of antennas typically used for CDGNSS.

Moreover, there will be a strong incentive to use single-frequency receivers, whereas almost all receivers used for CDGNSS in surveying and similar applications are multi-frequency. Despite these user-side disadvantages, mass-market precise positioning will be expected to demonstrate convergence and accuracy performance rivaling that of the most demanding current precise positioning applications: Users will be dissatisfied with techniques requiring more than a few tens of seconds to converge to a reliable sub-decimeter solution.

Meeting this challenge calls for innovation targeting both the rover (user) equipment and the reference network. Here we examine the challenge from the point of view of the reference network and offer demonstration results for a low-cost end-to-end system.

The recent trend in precise satellite-based positioning has been toward precise point positioning (PPP), whose primary virtue is the sparsity of its reference network. But standard PPP requires several tens of minutes or more to converge to a sub-10-centimeter 95 percent horizontal accuracy. Faster convergence can be achieved by recasting the PPP problem as one of relative positioning, thereby exposing integer ambiguities to the end user.

This technique, known as PPP-RTK or PPP-AR, is mathematically similar to traditional network real-time kinematic (NRTK) positioning. As the network density is increased, sub-minute or even instantaneous convergence is possible with dual-frequency high-quality receivers. Even single-frequency PPP-RTK is possible, with convergence times of approximately 5 minutes for a 40-kilometer network spacing.

For PPP-RTK and NRTK, convergence time is synonymous with the time required to resolve the integer ambiguities that arise in double-difference (DD) carrier-phase measurements, referred to here as time to ambiguity resolution, or TAR. As reference networks become denser, they can better compensate for spatially-correlated variations in signal delay introduced by irregularities in the ionosphere and, to a lesser extent, in the neutral atmosphere. Improvement is manifest as reduced uncertainty in the atmospheric corrections that the network sends to the user. Reduced uncertainty in the atmospheric corrections is key to reducing TAR.

Prior work has established an analytical connection between uncertainty in the ionospheric corrections (denoted *σ _{i }*) and TAR. The existing literature does not, however, offer a satisfactory model for the dependence of σι on network density.

The prevailing model is based on single-baseline CDGNSS, which is inapt for PPP-RTK and NRTK. Moreover, prior work does not address the effect of network-side multipath on the accuracy of the corrections data, which becomes increasingly important as low-cost and poorly-sited reference stations are used to densify the network.

Here, we examine the relationship between ionospheric uncertainty and probability of correct ambiguity resolution, and present the results of an empirical investigation of the relationship between network density and the total uncertainty in network correction data. We developed a simple analytical model relating error variance in network corrections to network density. Our analysis and experiments indicate that for rapid TAR in challenging urban environments with low-cost receivers, network density must be significantly increased. We report on the design and deployment of a dense network in Austin, Texas, and demonstrate a new system that taps into the network to provide reliable vehicle lane-departure warning.

### AMBIGUITY RESOLUTION

Reducing the ionospheric uncertainty σι allows a strong prior constraint to be applied in the ionosphere-weighted model, thereby increasing P( = z), the probability that the estimated and true integer ambiguity vectors are equivalent. It is instructive to consider single-epoch ambiguity resolution (AR), for two reasons.

First, for stationary users with low-cost equipment, multipath errors dominate in the carrier-phase measurement and are strongly correlated over 100 seconds or more. Thus, if single-epoch AR fails then a static user may have to wait an unacceptably long time for multipath errors to decorrelate enough to permit AR. In any case, singe-epoch performance is a strong predictor of multi-epoch performance over an interval short enough (a few tens of seconds) to satisfy impatient mass-market users.

Second, a convenient and accurate analytical model (by Dennis Odijk and PJG Teunissen) for single-epoch AR reveals the dependency of P( = z) on scenario parameters of practical interest: the standard deviation of ionospheric correction errors, the number of visible satellites, the standard deviation of undifferenced carrier- and code-phase measurement errors (including multipath-induced errors), a satellite geometry factor, the number p of free parameters to be estimated (p=3 for negligible tropospheric error, p=4 to estimate a single additional tropospheric parameter), and the number of carrier frequencies broadcast by each of the satellites (1, 2 or 3) along with each carrier’s wavelength.

The model is highly accurate for single-epoch AR, but only approximate for multiple epochs, with accuracy degrading as the data interval lengthens. The model’s inaccuracy results from its assumption that overhead satellites remain static from epoch to epoch, which yields pessimistic results for even fairly short data capture intervals (for example, 30 seconds). Fully accounting for satellite motion in an analytical model for P( = z) is an open problem, which is why studies that wish to account for satellite motion resort to simulation.

Figures 1 and 2 show single-epoch, single-frequency results from the analytical P( = z) model for parameters approximately reflecting the mass-market use case. The most important conclusion to draw from these figures is that for single-epoch, single-frequency AR to be even moderately reliable (PT⩾0.9) over the next few years, the ionospheric uncertainty σι must be held under 2 millimeters. This will relax somewhat as more Galileo and MEO BeiDou satellites come online, but signal blockage in built-up areas will raise the effective elevation mask angle significantly above the 15 degrees assumed here, reducing the number of available satellites. Thus, sub-2-mm ionospheric uncertainty remains desirable for urban environments even as GNSS constellations become fully populated.

Figures 3 and 4 offer results for a dual-frequency (L1-L2) single-epoch scenario. All other scenario parameters are held as for the single-frequency scenario except that, in an attempt to be somewhat more pessimistic, P(m) is based only on GPS satellites. It is assumed that from each satellite the user can extract dual-frequency measurements. As with the single-frequency case, it is evident that dual-frequency PT is strongly dependent on σι. The dual-frequency case is more forgiving, but substantial performance improvement can still be had by reducing σι to under 2 mm.

**Corrections Uncertainty and Network Density.** A key question arises in connection with *σ _{i}*: How is related to reference network density? One expects to decrease with increased network density, but what is the exact relationship?

Dennis Odijk’s work adopts a linear relationship between *σ _{i }*and the distance

*l*between the user and the nearest reference station:

*σ _{i }*=

*βl, 0.3 ≤ β ≤ 3 mm/km*

Parameter *β *depends on ionospheric activity; Odijk recommends determining *β* empirically. Similarly, his other work adopts a linear relation, with *β = *mm/km. But there appears to be no justification for applying this linear model to ionospheric corrections provided to a user by a network of reference receivers. The linear trend corresponds to individual single-baseline solutions involving a single master reference station without network aiding; it is not representative of how *σ _{i }*varies for a rover within a reference network.

Instead of determining how *σ _{i }*varies throughout a reference network, it will be more useful to consider the spatial variation in the variance of aggregate error in network-provided corrections. The aggregate error variance, denoted , can be modeled as the sum of variances associated with (1) residual ionospheric delay error, (2) residual neutral atmospheric (hereafter tropospheric) error, and (3) error due to carrier-phase multipath at the reference network stations:

This model assumes that precise orbital ephemerides are used to eliminate spatially-correlated errors due to satellite ephemeris errors and that the contribution to from reference station carrier-phase thermal noise is negligible compared to reference station carrier-phase multipath error.

Focusing therefore on *σ _{v }*, consider its relationship to reference network density

*γ*, expressed in stations per unit area. This relationship depends on the assumed model for the DD ionospheric and tropospheric delays. Let

*a*denote the master reference station and let S = {

*s*} denote the set of all secondary stations available in the network. Then, for pivot satellite

_{1}, s_{2}, …, s_{N}*i*and alternate satellite

*j*, suppose that the true combined DD atmospheric delay at secondary station

*s*∈

**S**can be accurately modeled as follows, where x

*, y*

_{s}*, and z*

_{s}*represent the secondary station’s east, north, and up displacement from the master:*

_{s }(1)

Dai et al. refer to this model as a linear interpolation model or first-order surface model. The quantities and are the model parameters for the satellite pair

*i, j.*

For the linear model in (1), one can show that if stations are sufficiently uniformly distributed (i.e., no station clumping), then the average value of σ_{v }across a network, denoted * _{ }*, is approximately related to the network density

*γ*by

(3)

where

*q*is a parameter related to the variance of the uncorrelated errors ,

*s*∈

**S**. This approximation becomes highly accurate as

*γ*increases. [See full paper for details.]

It is clear from (3) that, for the linear model (1),* _{ }*can be driven to an arbitrarily small value by increasing the network density

*γ*, and this is true despite the presence of multipath in the reference station carrier-phase measurements. Whether (3) applies in practice depends on whether (1) can be considered an accurate model for , at least over a compact region. The following section examines this question empirically. It further seeks to identify, for an example dense reference network, the density

*γ*beyond which further reduction in

*no longer matters (would no longer improve ) because rover multipath dominates.*

_{}### ANALYSIS OF A DENSE REFERENCE NETWORK

We examined σι as a function of network density using data from several organizations providing GNSS reference station observations: National Geodetic Survey Continuously Operating Reference Stations, UNAVCO, and the California Real Time Network. This combination allowed analysis of a hypothetical reference network of 23 high-quality GNSS receivers with an overall network density of approximately 8 nodes/1,000 km^{2}, or an average inter-station spacing of 14 km. The relative positions of the sites selected to comprise this reference network, located between Los Angeles and Pomona, California, are depicted graphically below.

DD carrier-phase observations from GPS L1 C/A signals spanning GPS weeks 1850 through 1859 were used for the analysis. A minimum satellite elevation mask was enforced at 20 degrees. Any satellite not above the elevation mask and providing carrier-phase observations at both the beginning and end of each processing window was excluded. A step size of 10 minutes was used. The longest available sub-window, meeting the above requirements and providing a minimum of 6 satellite vehicles (1 pivot satellite and 5 others), was selected for processing.

To facilitate batch processing, integer ambiguities were assumed to be resolved correctly when the mean standard deviation of carrier-phase residuals for that solution was less than one quarter wavelength of the GPS L1 frequency. In application, this constraint resulted in rejecting only 0.6 percent of all solutions.

**Network Corrections Estimation. **Estimation of network corrections made use of least-squares estimation applied to carrier-phase residuals measured between master station LONG, denoted *a* hereafter, and secondary reference stations *s*∈**S**, where **S **is now taken to be the set of all stations other than LONG. Consider the following model for the DD carrier-phase measurement, expressed in meters, between master station *a*, secondary station *s*∈**S**, pivot satellite *i*, and alternate satellite *j*:

Here, *λ* is the carrier wavelength; is the DD carrier-phase measurement, in cycles; is the DD range; is the DD integer ambiguity; is the DD combined atmospheric delay, which includes tropospheric and ionospheric delays; and is the DD carrier-phase measurement error, which is dominated by carrier-phase multipath error at* a* and* s*.

Experimental analysis of* _{ }*as a function of network density proceeded as follows. A subset of secondary stations

**S**

*⊂*

_{k}**S**was chosen, together with

*a*, to act as the

*k*th test network. A large number

*K*of subsets

**S**

*of various geographic size and density were analyzed. Let {*

_{k}**S**\

**S**

*} denote the set of secondary stations not in the*

_{k}*k*th test network. For each

**S**

*,*

_{k}*k*= 1, 2,…,

*K,*all secondary stations in {

**S**\

**S**

*} were designated, one at a time, to act as a test station, or rover. Atmospheric delays estimated by the*

_{k}*k*th network for test station

*s*∈{

**S**\

**S**

*} were then differenced from actual delays measured by*

_{k}*s*to evaluate the quality of the atmospheric delay estimates.

Details of the atmospheric delay estimation procedure for the *k*th test network are as follows. For each *s*∈**S*** _{k}*, a DD measurement residual was formed for each pivot satellite

*i*and alternate satellite

*j*as

where was assumed known to sub-millimeter accuracy and was assumed to have been resolved correctly. The true DD atmospheric error contributing to (5) was assumed to vary linearly with geometry over sufficiently short baselines as modeled in (2). The DD multipath error term was assumed to be zero mean, and the component due solely to *s* was assumed to be uncorrelated with all corresponding components .

Under these assumptions, can readily be estimated via least squares. Let be the vector containing the residuals for |**S*** _{k}*|x

**1**. This residuals vector can be modeled as

where *H* is an |**S*** _{k}*|x

**4**matrix whose rows are of the form [x

*y*

_{s }*z*

_{s }

_{s}**1**]. The

**4**x

**1**vector contains the parameters of the hyper-plane to be estimated at each epoch. The |

**S**

*|x*

_{k}**1**vector

*contains DD measurement errors.*

**w**^{ij}An estimate from a least-squares solution of (6) was used to produce a network correction for a test secondary station *s*∈{**S**\**S*** _{k}*}, acting as rover, at location

*x*:

_{s, }y_{s, }z_{s }The subscript *l* on the atmospheric correction indicates that the correction is based on a linear model for DD atmospheric errors; it is used to distinguish the correction from those produced by a quadratic model later on. The correction was applied at test station *s*∈{**S**\**S*** _{k}*} to produce a corrected DD phase measurement

This procedure was repeated at each epoch for each satellite pair *i, j* visible to each test station *s*∈{**S**\**S*** _{k}*} of the

*k*th test network

*, k =*1, 2,…

*K*.

If the assumed models hold, then in the limit as the network density increases, can be modeled as

where is DD phase measurement error due only to multipath at *s*. In other words, as network density increases, application of the network correction eliminates not onlybut also , the component of the DD phase measurement error due to multipath at the master.

**Linear least-squares compared to quadratic-least squares estimation. **To evaluate the assumption that DD tropospheric and ionospheric errors vary proportional to relative position, **c**_{1} was estimated with the full set of secondary stations **S** for single epochs at 300 second intervals. The probability distributions of the contributions of those parameters (e.g., *c*_{xl}x_{s }and not simply *c*_{xl}) are shown below. For comparison, equivalent values are calculated for a quadratic least-squares estimate of the following form:

Here, the subscript *q* of denotes a quadratic model for DD atmospheric delays. The distributions of comparable terms from (9) are also shown in the next two figures. These data represent the collection of all satellites above the elevation mask angle. It is noted that when all satellites are considered together, the expected value of these terms is very near zero.

The next two figures show the same data as the two above, but with each GPS satellite plotted separately. It is noted that the linear parameters, when considering only a particular satellite, are not necessarily zero-mean. This is hypothesized to be a manifestation of the satellite orbit reflected in the tropospheric and ionospheric errors. It is interesting to note that the quadratic terms shown in the second figure below largely exhibit zero mean behavior despite non-zero mean for the associated linear terms.

The figure above shows the probability distributions of the difference between (7) and (9) (i.e., ) at three representative reference station positions. It can be noticed that despite the increasing baseline distance of LORS and BGIS as compared to CGDM, there is no apparent correlation in these estimation errors. Notice that CGDM and LORS have very similar distributions despite their difference in baselines. BGIS and LORS, with similar baselines, exhibit very different distributions. There is no apparent correlation found between reference station positions and these error terms. Additionally, these distributions are zero-mean for all *s*∈**S **(to within 0.5 mm in each case) with 68.27% boundaries positioned between 1.5-5.5 mm. Because these errors appear indistinguishable from multipath, it is concluded, for this specific network and time period, that linear least-squares estimation is sufficient for estimating tropospheric and ionospheric errors. This is fortunate, because the linear model for atmospheric DD delays provides an averaging effect on multipath present at the reference stations which minimizes the introduction of multipath errors into the estimates produced.

**Uncorrected carrier-phase residuals. **The figure below shows the expected values for DD carrier-phase residual standard deviations for all *s*∈**S **through use of uncorrected observations. These data were produced by averaging the standard deviation of the DD carrier-phase residuals calculated at each epoch across all satellites present in the solution. The fitted curve indicates a linear growth of DD carrier-phase residuals with *β = *0.62 mm/km. Additionally, the mm-level scatter of these data points suggest that position biases of the resolved reference station positions are also mm-level. If the linear fit is shifted down by approximately 4 mm (e.g., taking the minimum data points as those with very little position bias) and extrapolated to 0 km, one can consider this as providing a rough estimate of DD multipath at the reference stations; 4.7 mm (DD) or 3.3 mm (single-difference equivalent).

**Uncorrected Carrier-Phase Residuals.** Figure 5 shows the expected values for DD carrier-phase residual standard deviations for all secondary stations, based on observations that were not corrected for atmospheric delay. These data were produced by averaging the standard deviation of the DD carrier-phase residuals calculated at each epoch across all satellites present in the solution. The fitted curve indicates a linear growth of DD carrier-phase residuals with distance to the master. The mm-level scatter of these data points suggest that biases of the resolved reference station positions are also mm-level.

**Network-Corrected Residuals**. Figure 6 displays similar data to Figure 5, except that the carrier-phase residuals are those that remain after network corrections are applied. Each data point corresponds to a particular subset of secondary stations together with the master, and a particular rover selected at random from the remaining stations. Both the size and specific selection of secondary stations comprising each subset were randomly selected. In all, 70 different network configurations and more than 3.67 million NRTK solutions were analyzed.

Figure 6 shows that carrier-phase residuals after application of network corrections are considerably reduced compared to those original magnitudes seen in Figure 5. With increasing network density, the DD residuals’ deviation asymptotically approaches a minimum value of about 4 mm, which corresponds to an undifferenced deviation of 2 mm. This floor is due to multipath at the rover. Deviations in excess of this floor are caused by residual ionospheric errors and, to a lesser extent, neutral atmospheric errors. Attributing the excess deviation entirely to residual ionospheric errors, and assuming these are uncorrelated with multipath, one can estimate from Figure 6 the undifferenced ionospheric uncertainty. For example, for a 50-km inter-station distance, σ_{ι}=(℘(14^{2} – 4^{2}))/2=6.7mm. To achieve the σ_{ι}<2 mm recommended earlier for fast and reliable AR, station separation should be no more than 22 km, which we round down to a recommended value of 20 km to provide a margin of station redundancy.

### NETWORK DEPLOYMENT

We have developed and deployed a low-cost reference network testbed in Austin, Texas, with site hosting courtesy of the Texas Department of Transportation. The Longhorn Reference Network boasts a dozen stations, with plans for 20 (Figure 7). The network’s average inter-station spacing is far shorter than the 20-km spacing recommended earlier. The tighter spacing provides redundancy and flexibility of experimentation. The low-cost reference stations are deployed in environments with greater multipath and signal blockage than those of the high-quality stations studied earlier. Such non-ideal signal environments are to be expected in a dense low-cost reference network, for which choice of station siting is driven largely by opportunity.

The reference station design, pictured in Figure 8 and diagrammed in Figure 9, is novel. Each station is a self-contained, solar-powered node supporting a software-defined dual-frequency, dual-antenna GNSS receiver with an always-on cellular connection to university servers for data collection and software maintenance.

**Live Vehicle Demonstration.** In partnership with Radiosense, an Austin-based precise positioning startup, we have developed and demonstrated a low-cost vehicle lane departure warning system that receives corrections from our dense reference network. The system takes in lane widths from an external database and infers a safe driving corridor within each lane by analyzing the behavior of human drivers on the same road. A vehicle’s proximity to the lane boundary is displayed in real time to the driver and passengers.

For robustness against cycle slips and to provide a baseline against which to compare future improvements, the system currently employs single-epoch CDGNSS positioning without aiding from additional sensors. In choosing a single-epoch approach, the system naively discards information regarding the underlying integer ambiguities at the beginning of each measurement epoch. Still, the system performs well with the typical number of overhead signals in a light urban environment: correct and internally-validated solutions were available in over 92 percent of measurement epochs. When a second rover antenna is included to combat multipath with spatial diversity, this percentage improves to 96. Such good single-epoch performance suggests that, when armed with additional sensor aiding and proper integer ambiguity persistence, reliable and accurate vehicle positioning can be maintained in more challenging environments.

**Demonstration setup.** The live demonstration followed a predetermined route in the vicinity of the University of Texas campus. The 1-mile route (Figure 10) passed through both open-sky and partially-blocked environments.

Prior to the demonstration, the vehicle was driven several times on the same route collecting GNSS measurements to precisely map typical driving trajectories on the route. The ensemble of trajectories was used to build a centimeter-accurate model of the lane center along the route. The sensing equipment employed during this mapping phase is no different than that used during the demonstration, making feasible eventual crowd-sourcing, wherein end-user vehicles generate and update the centerline models.

The demonstration vehicle was outfitted with two dual-frequency GNSS antennas mounted with magnetic bases onto the roof. The first antenna, designated primary, operated as the rover in a single-baseline CDGNSS solution against the master reference station of the Longhorn Reference Network, as illustrated in Figure 11. This baseline provided the geo-referenced, centimeter-accurate vehicle position. The other antenna, designated secondary, was paired with the primary antenna to produce a constrained-baseline CDGNSS solution providing sub-degree-accurate vehicle heading. The secondary antenna also served as a backup when the primary antenna produced a result that did not pass the precise positioning engine’s internal validity testing.

The GNSS antennas were connected to a low-cost, dual-frequency front-end in the trunk of the vehicle (FIGURE 12)which downconverted and digitized the incoming signals and subsequently fed them to a low-cost single-board computer running the precise positioning engine. A cellular modem received real-time measurements from the master reference station, while a WiFi router streamed real-time solutions to several Android devices in the vehicle for real-time visualization of precise within-lane position.

**Demonstration Results.** Figures 13, 14 and 15 show snapshots of the Android application and a still frame of the side of the vehicle in three different scenarios. The large rectangle indicates vehicle position with respect to the modeled lane center, changing color from green, when the vehicle is within the safe driving corridor, to yellow as the vehicle nears the edge of the lane, and finally to red if the vehicle breaks the lane boundary. One could imagine wrapping a control loop around these signals to enable last-moment lane-keeping.

Figure 16 reveals the precision with which the positioning engine was able to locate the vehicle’s driver-side antenna in four repeated passes along the test route. The variation between the four yellow traces is primarily due to driver non-repeatability; actual measurement precision is at the centimeter scale. A small bias in the traces’ registration to the picture is present because Google Earth imagery is only registered to the International Terrestrial Reference Frame with meter-level accuracy.

Figure 17 shows a time history of the vertical deviation from the route mean, in meters. The zoomed view of the vertical deviation shown in Figure 18 allows one to appreciate the precision of the positioning engine: the vertical trajectory is smooth at the centimeter level. Figure 19 shows the DD residuals in carrier phase and pseudorange for GPS PRN 30 during the four loops in Figure 17. One-sigma undifferenced phase and pseudorange deviations are 3.4 mm and 42 cm, respectively.

The figures demonstrate that the precise positioning engine fed by reference data from the Longhorn Reference Network maintained centimeter-accurate knowledge of the vehicle’s position during almost the entire trajectory, despite passing between a large football stadium and parking garage, each of which introduced significant signal blockage and multipath.

For the data shown in Figure 17, 96 percent of the 3.3-Hz measurement epochs resulted in a correct and internally-validated positioning solution. The majority of the remaining solutions were correct but did not pass internal validation. For only 0.6 percent of solutions were the carrier-phase integer ambiguities resolved incorrectly, but all of these incorrect solutions were caught and excluded by the validation algorithm.

Furthermore, the number of overhead signals during the time in which this particular dataset (set A) was taken was average, as seen in the upper plot of Figure 20. 16 signals above 15 degrees elevation were available during this time. In contrast, the number of overhead signals for a second dataset taken 8 days prior (set B) was much worse, with only 12 signals above 15 degrees elevation, as seen in the lower plot.

For insight into the performance of the positioning engine as a function of the number of overhead satellites, Table 1 details the performance of these two datasets (as well as a third dataset) in terms of the percentage of epochs that passed the positioning engine’s internal validation testing, based on a ratio test with a fixed threshold of 2.0. Results are shown for single- and dual-antenna positioning solutions and for dual-antenna vehicle heading solutions.

A large drop-off in positioning performance occurs when the number of overhead signals is reduced below 16, while the constrained-baseline heading determination performance remains good throughout. Fortunately, it will not be long until even more signals are available. Within the next 8 months, the Galileo constellation will add six fully operational satellites. These will bring the number of GPS L1, GPS L2C, Galileo E1, and SBAS signals that are above 15 degrees elevation to 16 or more 95 percent of the time, enabling high-reliability single-epoch CDGNSS positioning.

### CONCLUSIONS

For a sufficiently dense reference network, linear least squares estimation can be applied to the task of reducing uncertainties due to tropospheric and ionospheric delays for the purposes of providing improved positioning accuracy as well as faster time to ambiguity resolution for carrier-phase differential positioning. High network density allows use of a strong linear model for atmospheric delays, which has the virtue of suppressing network-side multipath errors in the provided corrections.

A network of 23 high-quality reference stations in the vicinity of Los Angeles, California, was studied to determine what network density is sufficient to make all network-side error sources negligible compared to rover receiver multipath. A density of three stations per 1,000km2, or an average inter-station spacing of 20 km, was found to drive network-side ionospheric, tropospheric, and multipath errors well below rover receiver multipath.

These findings motivate a significant densification of permanent reference networks, at least in built-up areas where signal blockage and multipath are common, to support mass-market applications for which low user (rover receiver) cost and rapid convergence to a reliable sub-decimeter position are a priority. In a light urban setting, and with the kind of satellite coverage that will soon become the norm, we demonstrated vehicle lane departure warning in a field test that produced highly reliable instantaneous sub-decimeter positioning.

### ACKNOWLEDGMENTS

This work was supported in part by Samsung Research America, by the Data-Supported Transportation Operations and Planning Center (D-STOP), a Tier 1 USDOT University Transportation Center, and by the Texas Department of Transportation under the Connected Vehicle Problems, Challenges and Major Technologies project.

## Follow Us