Geodetics Preprint August 2000 Epoch-by-EpochÔ Positioning and Navigation Paul J. de Jonge, Yehuda Bock, and Michael Bevis Geodetics, Inc. La Jolla, CA BIOGRAPHY applications such as automated airborne landing systems, Paul de Jonge received his Ph.D. in geodesy at the Delft harbor navigation, automated highways, machine control, University of Technology in the Netherlands, and he and GPS/INS integration. Epoch-by-EpochÔ positioning currently works for Geodetics. Dr. Yehuda Bock is the also provides a useful architecture for static positioning, President and CEO of Geodetics. He is also a Research since relatively rare but highly discordant solutions (i.e. Geodesist and Senior Lecturer at the Scripps Institution of solution ‘outliers’ associated with extreme levels of Oceanography where he directs the Scripps Orbit and measurement noise) can be recognized and excluded from Permanent Array Center (SOPAC) and the California any position averaging. When single epoch solutions are Spatial Reference Center (CSRC). Dr. Michael Bevis is the averaged by a median filter, the scatter of the averaged Vice-President of Geodetics and a Professor of solutions is inversely proportional to the logarithm of Geophysics and Ge odesy at the University of Hawaii. averaging time. ABSTRACT INTRODUCTION Geodetics, Inc. has developed a new class of Geodetics, Inc. has developed an innovative instantaneous, real-time GPS positioning algorithms based suite of software modules that implement Epoch-by- on dual frequency phase and pseudorange data, which we EpochÔ analysis of dual-frequency phase and range data call Epoch-by-EpochÔ positioning. The main advantage from two or more GPS receivers. The central feature of this of this approach over traditional real-time methods is that approach is that a high accuracy geodetic solution, based integer-cycle phase ambiguities are independently on instantaneous integer-cycle ambiguity resolution, is estimated for each and every observation epoch, allowing achieved for each measurement epoch using only the precise positions (1-2 cm horizontal and 5-10 cm vertical) observations collected at that epoch. Accordingly, each to be estimated given only a single epoch of data from two solution is independent of the solutions obtained for the or more receivers, in either baseline or network mode. Our previous and following epochs. Because Epoch-by- algorithms are extremely efficient – a single epoch on a EpochÔ analysis routinely resolves the complete set of single baseline can be processed in several milliseconds. double difference ambiguity parameters for each epoch, There is no need for the initialization period associated the technique delivers horizontal accuracies comparable to with conventional RTK algorithms. More importantly there those achieved in conventional kinematic and RTK is no need for reinitialization immediately following loss- processing. Our software can be used to process networks of-lock problems such as occur when a mobile GPS consisting of two or more GPS receivers, and allows receiver passes under a bridge or other obstruction. nearest-neighbor station spacings several times wider Accordingly, our approach minimizes the amount of time a than conventional RTK. receiver spends in a non-productive prepositioning mode Epoch-by-EpochÔ data processing is achieved of operations. Epoch-by-EpochÔ positioning has other by a core library of geodetic routines called the Rapid advantages, including a significantly extended range Network Analysis (RNA) library. All of the application relative to conventional RTK, and a better framework for software developed by Geodetics to date, and several quality control. The Epoch-by-EpochÔ approach can also applications currently under development, are built around be applied to the autonomous real-time determination of the proprietary RNA library. Our algorithms are extremely precise attitude and heading (about 0.05 degrees) of a efficient and allow the simultaneous, real-time network moving platform instrumented with 2-3 receivers and a positioning of at least 10 receivers sampling at a 1 Hz rigid multiple-antenna structure. A combination of this interval on a standard medium-scale PC workstation. A autonomous navigational capability with precise single epoch on a single baseline can be processed in only instantaneous (relative) positioning is ideal for many several milliseconds. Although Epoch-by-EpochÔ Geodetics 2000 Geodetics Preprint August 2000 analysis is instantaneous, in the sense that it produces a CASE STUDY 1: VEHICLE TRACKING solution for each epoch or instant, without direct or Figures 1 and 2 show position and velocity indirect reference to data collected at any other time, it solutions obtained for a vehicle tracking test implemented provides a useful computational architecture for post- by Applanix Corporation on March 3, 2000 in Ontario, processing of multiple-epoch data batches as well as for Canada. The vehicle was a passenger van that traveled up real-time positioning. We illustrate some of these to 6 km away from a GPS base station at speeds of up to 70 applications in the following sections. km/hour. The GPS sampling interval was one second. All results shown in these figures are raw single-epoch KINEMATIC POSITIONING solutions. Figure 1 focuses on a subset of this test in Conventional RTK software partitions the integer which the vehicle made four circuits of the same C-shaped ambiguity problem into an initialization stage that typically section of road. We color code each of the four circuits so lasts 30-40 seconds. During this time phase ambiguities as to provide relative timing information in plots featuring are assessed using multiple epochs of data, and when the only spatial axes. This allows us to confirm that the two integer ambiguities are finally resolved with sufficient tracks apparent in the blow-up map (Fig. 1) correspond to confidence, initialization ends and actual positioning the two sides of the road, with both sides being traversed begins. Should a serious loss of satellite lock occur, as in all four circuits. The 3-D plot shows the first three visits might happen when a mobile GPS receiver passes under a to one end of the track, and one can clearly resolve the bridge, then once the hardware recovers and the GPS data manner in which the vehicle was turned around. Note that stream begins anew, it is necessary for the software to the vehicle arrives well-centered in the right lane, but reinitialize – that is, to return to the problem of resolving leaves the turn less well centered in the adjacent lane (at integer ambiguities as a preliminary to actual positioning. least for the first 40 meters). The total delay associated with loss of lock depends on The velocity solutions (Fig. 2) were obtained by the hardware as well as the software, since the receiver simple differencing of the position solutions, without any also has to reinitialize tracking after a loss of lock. The smoothing or editing. The velocity solutions are produced high performance, dual frequency GPS receivers available at a frequency of 1 Hz, but are offset by 0.5 seconds to civilians can recover from a very brief loss of lock in a relative to the position time series. Because they are second or two, but typically take 10-20 seconds to recover produced by differencing single-epoch positions, velocity from a sustained loss of lock. This hardware down time is time series are strongly sensitive to small errors in followed by additional down time associated with software position. We can detect an artifact in the vertical velocity initialization. Military dual frequency receivers, with direct (Vu) time series about 115 seconds into the test. access to the P code, can relock the satellites much more After completing the four circuits around the C- rapidly than non-classified receivers, with high shaped track, the vehicle stopped for 178 seconds. Both performance military units ‘cold starting’ in less than one the velocity solutions and the positions during this period second. are shown in the lower half of Figure 2. Note that the up Because Epoch-by-EpochÔ algorithms estimate coordinate has an order of magnitude more scatter than the integer ambiguities instantaneously, they completely does the north or east coordinate. In part this reflects the eliminate the software downtime associated with fact that we were estimating atmospheric parameters as initialization and reinitialization. Accordingly Geodetics’ part of each epoch analysis. A shorter static test took software cuts the total hardware plus software place at greater distance from the base station later in the initialization/reinitialization time to something very close to test, and the horizontal scatter associated with these the hardware initialization time. This leads to significantly solutions is shown in the lower left section of Figure 1. improved performance in civilian applications, and dramatically improved performance with military hardware. REAL-TIME ORIENTATION The second architectural advantage of Epoch-by- When three antennas are fixed to a rigid structure EpochÔ positioning is that a mobile receiver can be with spacings of order 1-2 meters, it is possible to positioned against a network of reference stations, and determine the orientation of the antenna array (i.e. its not just a single base station. This makes it easier to pitch, roll and heading angles) without reference to any implement wide-area RTK based on multiple reference external static base station. The solutions for the lengths stations.