Showing papers in "Gps Solutions in 2009"

A simplified yaw-attitude model for eclipsing GPS satellites

Abstract: A simplified yaw-attitude modeling, consistent with Bar-Sever (1996), has been implemented and tested in the NRCan PPP software. For Block IIR GPS satellite it is possible to model yaw-attitude control during eclipsing periods by using the constant hardware yaw rate of 0.20°/s. The Block IIR satellites maintain the nominal yaw attitude even during a shadow crossing (Y. E. Bar-Sever, private communication, 2007), except for the noon and shadow midnight turn maneuvers, both of which can be modeled and last up to 15 min. Thus, for Block IIR satellites it is possible to maintain continuous satellite clock estimation even during eclipsing periods. For the Block II/IIA satellites, it is possible to model satisfactorily the noon turns and also shadow crossing, thanks to the permanent positive yaw bias of 0.5°, implemented in November 1995. However, in order to model the Block II/IIA shadow crossings, satellite specific yaw rates should be used, either solved for or averaged yaw-rate solutions. These yaw rates as estimated by the Jet Propulsion Laboratory (JPL) can differ significantly from the nominal hardware values. The Block II/IIA post-shadow recovery periods, which last about 30 min, should be considered uncertain and cannot be properly modeled. Data from post-shadow recovery periods should, therefore, not be used in precise global GPS analyses (Bar-Sever 1996). For high-precision applications, it is essential that users implement a yaw-attitude model, which is consistent with the generation of the satellite clocks. Initial testing and analyses, based on the IGS and AC Final orbits and clocks have revealed that during eclipsing periods, significant inconsistencies in yaw-attitude modeling still exist amongst the IGS Analyses Centers, which contribute to the errors of the IGS Final clock combinations.

Ambiguity resolution in precise point positioning with hourly data

Abstract: Precise point positioning (PPP) has become a powerful tool for the scientific analysis of Global Positioning System (GPS) measurements. Until recently, ambiguity resolution at a single station in PPP has been considered difficult, due to the receiver- and satellite-dependent uncalibrated hardware delays (UHD). However, recent studies show that if these UHD can be determined accurately in advance within a network of stations, then ambiguity resolution at a single station becomes possible. In this study, the method proposed by Ge et al. J Geod 82(7):389–399, 2007 is adopted with a refinement in which only one single-difference narrow-lane UHD between a pair of satellites is determined within each full pass over a regional network. This study uses the EUREF (European Reference Frame) Permanent Network (EPN) to determine the UHD from Day 245 to 251 in 2007. Then 12 International GNSS Service stations inside the EPN and 15 outside the EPN are used to conduct ambiguity resolution in hourly PPP. It is found that the mean positioning accuracy in all hourly solutions for the stations inside the EPN is improved from (3.8, 1.5, 2.8) centimeters to (0.5, 0.5, 1.4) centimeters for the East, North and Up components, respectively. For the stations outside the EPN, some of which are over 2,000 km away from the nearest EPN stations, the mean positioning accuracy in the East, North and Up directions still achieves (0.6, 0.6, 2.0) centimeters, respectively, when the EPN-based UHD are applied to these stations. These results demonstrate that ambiguity resolution at a single station can significantly improve the positioning accuracy in hourly PPP. Particularly, UHD can be even applied to a station which is up to thousands of kilometers from the UHD-determination network, potentially showing a great advantage over current network-based GPS augmentation systems. Therefore, it is feasible and beneficial for the operators of GPS regional networks and providers of PPP-based online services to provide these UHD estimates as an additional product.

Kalman-filter-based GPS clock estimation for near real-time positioning

Abstract: In this article, an algorithm for clock offset estimation of the GPS satellites is presented. The algorithm is based on a Kalman-filter and processes undifferenced code and carrier-phase measurements of a global tracking network. The clock offset and drift of the satellite clocks are estimated along with tracking station clock offsets, tropospheric zenith path delay and carrier-phase ambiguities. The article provides a brief overview of already existing near-real-time and real-time clock products. The filter algorithm and data processing scheme is presented. Finally, the accuracy of the orbit and clock product is assessed with a precise orbit determination of the MetOp satellite and compared to results gained with other real-time products.

Antenna phase center calibration for precise positioning of LEO satellites

Abstract: Phase center variations of the receiver and transmitter antenna constitute a remaining uncertainty in the high precision orbit determination (POD) of low Earth orbit (LEO) satellites using GPS measurements. Triggered by the adoption of absolute phase patterns in the IGS processing standards, a calibration of the Sensor Systems S67-1575-14 antenna with GFZ choke ring has been conducted that serves as POD antenna on various geodetic satellites such as CHAMP, GRACE and TerraSAR-X. Nominal phase patterns have been obtained with a robotic measurement system in a field campaign and the results were used to assess the impact of receiver antenna phase patterns on the achievable positioning accuracy. Along with this, phase center distortions in the actual spacecraft environment were characterized based on POD carrier phase residuals for the GRACE and TerraSAR-X missions. It is shown that the combined ground and in-flight calibration can improve the carrier phase modeling accuracy to a level of 4 mm which is close to the pure receiver noise. A 3.5 cm (3D rms) consistency of kinematic and reduced dynamic orbit determination solutions is achieved for TerraSAR-X, which presumably reflects the limitations of presently available GPS ephemeris products. The reduced dynamic solutions themselves match the observations of high grade satellite laser ranging stations to 1.5 cm but are potentially affected by cross-track biases at the cm-level. With respect to the GPS based relative navigation of TerraSAR-X/TanDEM-X formation, the in-flight calibration of the antenna phase patterns is considered essential for an accurate modeling of differential carrier phase measurements and a mm level baseline reconstruction.

A constrained LAMBDA method for GPS attitude determination

Abstract: An improved method to obtain fixed integer ambiguity in GPS attitude determination is presented. Known conditions are utilized as constraints to acquire attitude information when the float solution and its variance–covariance matrix are not accurate enough. The searching ellipsoidal region is first expanded to compensate for errors caused by the inaccurate float solution. Then the constraints are used to shrink the region to a proper size, which maintains the true integer ambiguity. Experimental results demonstrate that this scheme gives a fast search time and a higher success rate in determining the fixed integer ambiguity than the unconstrained method. The accuracy of attitude angles is also improved.

Carrier phase-based integrity monitoring for high-accuracy positioning

Abstract: Pseudorange-based integrity monitoring, for example receiver autonomous integrity monitoring (RAIM), has been investigated for many years and is used in various applications such as non-precision approach phase of flight. However, for high-accuracy applications, carrier phase-based RAIM (CRAIM), an extension of pseudorange-based RAIM (PRAIM) must be used. Existing CRAIM algorithms are a direct extension of PRAIM in which the carrier phase ambiguities are estimated together with the estimation of the position solution. The main issues with the existing algorithms are reliability and robustness, which are dominated by the correctness of the ambiguity resolution, ambiguity validation and error sources such as multipath, cycle slips and noise correlation. This paper proposes a new carrier phase-based integrity monitoring algorithm for high-accuracy positioning, using a Kalman filter. The ambiguities are estimated together with other states in the Kalman filter. The double differenced pseudorange, widelane and carrier phase observations are used as measurements in the Kalman filter. This configuration makes the positioning solution both robust and reliable. The integrity monitoring is based on a number of test statistics and error propagation for the determination of the protection levels. The measurement noise and covariance matrices in the Kalman filter are used to account for the correlation due to differencing of measurements and in the construction of the test statistics. The coefficient used to project the test statistic to the position domain is derived and the synthesis of correlated noise errors is used to determine the protection level. Results from four cases based on limited real data injected with simulated cycle slips show that residual cycle slips have a negative impact on positioning accuracy and that the integrity monitoring algorithm proposed can be effective in detecting and isolating such occurrences if their effects violate the integrity requirements. The CRAIM algorithm proposed is suitable for use within Kalman filter-based integrated navigation systems.

On symmetrical three-dimensional datum conversion

Abstract: A 3-D similarity transformation is frequently used to convert GPS-WGS84-based coordinates to those in a local datum using a set of control points with coordinate values in both systems. In this application, the Gauss-Markov (GM) model is often employed to represent the problem, and a least-squares approach is used to compute the parameters within the mathematical model. However, the Gauss–Markov model considers the source coordinates in the data matrix (A) as fixed or error-free; this is an imprecise assumption since these coordinates are also measured quantities and include random errors. The errors-in-variables (EIV) model assumes that all the variables in the mathematical model are contaminated by random errors. This model may be solved using the relatively new total least-squares (TLS) estimation technique, introduced in 1980 by Golub and Van Loan. In this paper, the similarity transformation problem is analyzed with respect to the EIV model, and a novel algorithm is described to obtain the transformation parameters. It is proved that even with the EIV model, a closed form Procrustes approach can be employed to obtain the rotation matrix and translation parameters. The transformation scale may be calculated by solving the proper quadratic equation. A numerical example and a practical case study are presented to test this new algorithm and compare the EIV and the GM models.

A closed-form formula for GPS GDOP computation

Abstract: Geometric dilution of precision (GDOP) is often used for selecting good satellites to meet the desired positioning precision. An efficient closed-form formula for GDOP has been developed when exactly four satellites are used. It has been proved that increasing the number of satellites for positioning will always reduce the GDOP. Since most GPS receivers today can receive signals from more than four satellites, it is desirable to compute GDOP efficiently for the general case. Previous studies have partially solved this problem with artificial neural network (ANN). Though ANN is a powerful function approximation technique, it needs costly training and the trained model may not be applicable to data deviating too much from the training data. Using Newton’s identities from the theory of symmetric polynomials, this paper presents a simple closed-form formula for computing GDOP with the inputs used in previous studies. These inputs include traces of the measurement matrix and its second and third powers, and the determinant of the matrix.

Ultra-tight GPS/INS/PL integration: a system concept and performance analysis

Abstract: The architecture of the ultra-tight GPS/INS/PL integration is the key to its successful performance; the main feature of this architecture is the Doppler feedback to the GPS receiver tracking loops. This Doppler derived from INS, when integrated with the carrier tracking loops, removes the Doppler due to vehicle dynamics from the GPS/PL signal thereby achieving a significant reduction in the carrier tracking loop bandwidth. The bandwidth reduction provides several advantages such as: improvement in anti-jamming performance, and increase in post correlated signal strength which in turn increases the dynamic range and accuracy of measurements. Therefore, any degradation in the derived Doppler estimates will directly affect the tracking loop bandwidth and hence its performance. The quadrature signals from the receiver correlator, I (in-phase) and Q (quadrature), form the measurements, whereas the inertial sensor errors, position, velocity and attitude errors form the states of the complementary Kalman filter. To specify a reliable measurement model of the filter for this type of integrated system, a good understanding of GPS/PL signal characteristics is essential. It is shown in this paper that phase and frequency errors are the variables that relate the measurements and the states in the Kalman filter. The main focus of this paper is to establish the fundamental mathematical relationships that form the measurement model, and to show explicitly how the system error states are related to the GPS/PL signals. The derived mathematical relationships encapsulated in a Kalman filter, are tested by simulation and shown to be valid.

Geometry-free undifferenced, single and double differenced analysis of single frequency GPS, EGNOS and GIOVE-A/B measurements

Abstract: This paper demonstrates a geometry-free GNSS measurement analysis approach and presents results of single frequency GPS, EGNOS and GIOVE short and zero baseline measurements. The purpose is to separate the different contributions to the measurement noise of pseudo range code and carrier phase observations at the receiver. The influence of multipath on the different combinations of observations is also determined. Quantitative results are presented for the thermal code and phase measurement noise and for the correlation between the observations. Comparison of the results with theoretical approximations confirms the validity of the used approach. Results from field measurements clearly show less thermal noise on the Galileo E1BC observations than on the GPS L1C/A observations due to the new signal modulation. The feasibility of ambiguity resolution with a geometry-free model is also discussed including the significant impact of multipath thereon.

Ionospheric delay corrections for single-frequency GPS receivers over Europe using tomographic mapping

Abstract: The majority of navigation satellites receivers operate on a single frequency and experience a positioning error due to the ionospheric delay. This can be compensated for using a variety of approaches that are compared in this paper. The study focuses on the last solar maximum. A 4D tomographic imaging technique is used to map the ionospheric electron density over the European region during 2002 and 2003. The electron density maps are then used to calculate the excess propagation delay on the L1 frequency experienced by GPS receivers at selected locations across Europe. The excess delay is applied to correct the pseudo-range single frequency observations at each location and the improvements to the resulting positioning are calculated. The real-time tomographic technique is shown to give navigation solutions that are better than empirical modelling methods and approach the accuracy of the full dual-frequency solution. The improvements in positioning accuracy vary from day to day depending on ionospheric conditions but can be up to 25 m during mid-day during these solar maximum conditions at European mid-latitudes.

Analysis of inversion errors of ionospheric radio occultation

Abstract: The retrieved electron density profile of ionospheric radio occultation (RO) simulation data can be compared with the background model value during the simulation and the inversion error can be obtained exactly. This paper studies the inversion error of ionospheric RO through simulation. The sources of the inversion errors are analyzed. The impacts of measurement errors, such as the errors in phase measurements and satellite orbits, are very small and can be neglected. The approximation of straight-line propagation introduces errors at the height of the F1 layer under solar maximum condition. The spherical symmetry approximation of the electron density distribution is found to be the main source of the inversion error. The statistical results reveal some characteristics of the inversion errors. (1) The relative error increases with enhanced solar activity. (2) It is larger in winter than in equinox season, and it is smallest in summer. (3) For all seasons, it is smaller at middle latitude than at other latitudes. (4) For all seasons and geomagnetic latitudes, it is smaller at daytime than at other times. The NmF2 of the ROs from COSMIC are compared with the measurements of ionosondes, and the relative differences show the same dependencies on season, geomagnetic latitude and local time, as the relative errors of the simulated ionospheric ROs.

Assessment of time-series of troposphere zenith delays derived from the Global Data Assimilation System numerical weather model

Abstract: Troposphere zenith path delays derived from the Global Data Assimilation System (GDAS) numerical weather model (NWM) are compared with those of the International GNSS Service (IGS) solutions over a 1.5-year period at 18 globally distributed IGS stations. Meteorological parameters can be interpolated from the NWM model at any location and at any time after December 2004. The meteorological parameters extracted from the NWM model agree with in situ direct measurements at some IGS stations within 1 mbar for pressure, 3° for temperature and 13% for relative humidity. The hydrostatic and wet components of the zenith path delay (ZPD) are computed using the meteorological parameters extracted from the NWM model. The total ZPDs derived from the GDAS NWM agree with the IGS ZPD solutions at 3.0 cm RMS level with biases of up to 4.5 cm, which can be attributed to the wet ZPDs estimates from the NWM model, considering the less accurate interpolated relative humidity parameter. Based on this study, it is suggested that the availability and the precision of the GDAS NWM ZPD should be sufficient for nearly all GPS navigation solutions.

Optimization of two GPS/MEMS-IMU integration strategies with application to sports

Abstract: The application of low-cost L1 GPS receivers integrated with micro-electro-mechanical system (MEMS) inertial measurement units (IMU) allows the continuous observation of position, velocity and orientation which opens new possibilities for comparison of athletes’ performance throughout a racecourse In this paper, we compare loosely and closely coupled integration strategies under realistic racing scenarios when GPS is partially or completely masked The study reveals that both integration approaches have a similar performance when the satellite constellation is completed or the outages are short However, for less than four satellites, the closely coupled strategy clearly outperforms the loosely coupled approach The second part of the paper is devoted to the important problem of system initialization, because the conventional GPS/IMU alignment methods are no longer applicable when using MEMS-IMU We introduce a modified coarse alignment method and a quaternion estimation method for the computation of the initial orientation Simulations and practical experiments reveal that both methods are numerically stable for any initial orientation of the sensors with the error characteristics of MEMS-IMUs Throughout the paper, our findings are supported by racing experiments with references provided in both, the measurement and the navigation domains

Network-based geometry-free three carrier ambiguity resolution and phase bias calibration

Abstract: Continuously operating reference stations (CORS) are increasingly used to deliver real-time and near-real-time precise positioning services on a regional basis. A CORS network-based data processing system uses either or both of the two types of measurements: (1) ambiguity-resolved double-differenced (DD) phase measurements, and (2) phase bias calibrated zero-differenced (ZD) phase measurements. This paper describes generalized, network-based geometry-free models for three carrier ambiguity resolution (TCAR) and phase bias estimation with DD and ZD code and phase measurements. First, the geometry-free TCAR models are constructed with two Extra-Widelane (EWL)/Widelane (WL) virtual observables to allow for rapid ambiguity resolution (AR) for DD phase measurements without distance constraints. With an ambiguity-resolved WL phase measurement and the ionospheric estimate derived from the two EWL observables, an additional geometry-free equation is formed for the third virtual observable linearly independent of the previous two. AR with the third geometry-free model requires a longer period of observations for averaging than the first two, but is also distance-independent. A more general formulation of the geometry-free model for a baseline or network is also introduced, where all the DD ambiguities can be more rigorously resolved using the LAMBDA method. Second, the geometry-free models for calibration of three carrier phase biases of ZD phase measurements are similarly defined for selected virtual observables. A network adjustment procedure is then used to improve the ZD phase biases with known DD integer constraints. Numerical results from experiments with 24-h dual-frequency GPS data from three US CORS stations baseline lengths of 21, 56 and 74 km confirm the theoretical predictions concerning AR reliability of the network-based geometry-free algorithms.

Accuracy assessment of the National Geodetic Survey's OPUS-RS utility

Abstract: OPUS-RS is a rapid static form of the National Geodetic Survey’s On-line Positioning User Service (OPUS). Like OPUS, OPUS-RS accepts a user’s GPS tracking data and uses corresponding data from the U.S. Continuously Operating Reference Station (CORS) network to compute the 3-D positional coordinates of the user’s data-collection point called the rover. OPUS-RS uses a new processing engine, called RSGPS, which can generate coordinates with an accuracy of a few centimeters for data sets spanning as little as 15 min of time. OPUS-RS achieves such results by interpolating (or extrapolating) the atmospheric delays, measured at several CORS located within 250 km of the rover, to predict the atmospheric delays experienced at the rover. Consequently, standard errors of computed coordinates depend highly on the local geometry of the CORS network and on the distances between the rover and the local CORS. We introduce a unitless parameter called the interpolative dilution of precision (IDOP) to quantify the local geometry of the CORS network relative to the rover, and we quantify the standard errors of the coordinates, obtained via OPUS-RS, by using functions of the form $$ \sigma ({\text{IDOP}},{\text{RMSD}}) = \sqrt {(\alpha \cdot {\text{IDOP}})^{2} + (\beta \cdot {\text{RMSD}})^{2} } $$ here α and β are empirically determined constants, and RMSD is the root-mean-square distance between the rover and the individual CORS involved in the OPUS-RS computations. We found that α = 6.7 ± 0.7 cm and β = 0.15 ± 0.03 ppm in the vertical dimension and α = 1.8 ± 0.2 cm and β = 0.05 ± 0.01 ppm in either the east–west or north–south dimension.

Advantage of velocity measurements on instantaneous RTK positioning

Abstract: Instantaneous real-time kinematic (RTK) techniques are one approach to achieving real-time high-accuracy positioning. The object of this paper is to show the advantage of adding velocity information to instantaneous RTK positioning. In this paper, velocity from Doppler frequency measurements is used to help resolve integer ambiguities in the LAMBDA method. In urban areas, pseudorange measurements can suffer from significant errors due to strong multipath, despite the use of advanced multipath-mitigation techniques. However, Doppler frequency measurements do not deteriorate as much as pseudoranges, because the multipath error on Doppler frequency measurements is only in the range of several centimeters. Our proposed method has been tested by both static and moving users in sub-urban environments. Single-epoch ambiguity fixing performance was improved compared to conventional ambiguity resolution without velocity information.

Total electron content processing from GPS observations to facilitate ionospheric modeling

Abstract: With the increasing global distribution of high rate dual-frequency global positioning system (GPS) receivers, the production of a real-time atmospheric constituent definition, total electron content (TEC), has become a beneficial contributor to the modeling applications used in the assessment of GPS position accuracy and the composition of the ionosphere, plasmasphere, and troposphere. Historically, TEC measurements have been obtained through post processing techniques to produce the quality of data necessary for modeling applications with rigorous error estimate requirements. These procedures necessitated the collection of large volumes of data to address the various abnormalities in the computation of TEC associated with the use of greater data quality controls and source selection while real-time modeling environments must rely on autonomous controls and filtration techniques to prevent the production of erroneous model results. In this paper we present methods for processing TEC in real time, which utilize several procedures including the application of an ionospheric model to automatically perform quality control on the TEC output and the computational techniques used to address receiver multipath, faulty receiver observations, cycle-slips, segmented processing, and receiver calibrations. The resulting TEC measurements are provided with rigorous error estimates validated using the vertical TEC from the Jason satellite mission.

Carrier phase wind-up in GPS reflectometry

Abstract: Changes in GPS transmitter and receiver antenna orientations induce variations in observed carrier phase values. An analytic formula for this well-known carrier phase wind-up correction is derived which generalizes a previous result. In addition, it is shown that in GPS reflectometry the wind-up values of direct and coherently reflected rays may differ by up to several centimeters. The results are discussed on the basis of simulated measurements.

MATLAB Algorithm Availability Simulation Tool

Abstract: This paper describes a set of MATLAB functions currently being developed for Space Based Augmentation System (SBAS) availability analysis. This toolset includes simulation algorithms that are constantly being developed and updated by various working groups. This set of functions is intended for use as a fast, accurate, and highly customizable experimental test bed for algorithm development. A user-friendly interface has also been developed for the tool. It is open source and can be downloaded from the Stanford GPS Research Laboratory web site (http://waas.stanford.edu/~wwu/maast/maast.html). Therefore, it provides a common ground for different working groups to compare their results. This paper demonstrates the utility of this toolset by analyzing the SBAS service volume models for the Conterminous United States. The paper describes the functionality provided within the software, and shows an example set of contour plots which are the means for investigating how different algorithms and parameters impact availability.

Tropospheric delay gradients from numerical weather prediction models: effects on GPS estimated parameters

Abstract: Several numerical weather prediction (NWP) models provide information on the 3D state of the neutral atmosphere which has enabled GNSS researchers to have improved a priori information of the delay induced in the GNSS signals. However, the quality of weather models on the one hand and computational difficulties on the other, are motivations to develop an algorithm based partly on NWP models, while still estimating the remaining residual delay through GNSS processing strategies. An algorithm has been developed to estimate horizontal delay gradients from Meteorological Service of Canada NWP models. The GNSS software “Bernese” has also been modified to handle these gradients, as well as zenith delay and mapping functions based on NWP models in phase and code observation equations. Month-long precise point positioning results show strong correlation between north–south hydrostatic gradients and latitude differences, with significant but less strong correlation with the height and zenith total delay parameters. The longitude components were not sensitive to the implementation of gradients. High precision GNSS applications such as long term geodynamics studies, realization of terrestrial reference frames and climatology and consequential interpretations may be affected by ignoring the asymmetry of the neutral atmosphere. In addition to estimating the gradients, implementing a priori information on gradients in the processing software may have an impact on estimated results and consequential interpretations.

A MATLAB toolbox for attitude determination with GPS multi-antenna systems

Abstract: In this paper a MATLAB toolbox for determining the attitude of a rigid platform by means of multiple non-dedicated antennas using global positioning system is presented. The programs embedded in this toolbox cover the RINEX data analysis, single point positioning, differential positioning, coordinate conversion, attitude determination, and other auxiliary functions. After forming the baselines through double-differenced (carrier phase smoothed) code observables, the attitude parameters are obtained by applying the direct attitude computation and the least squares attitude estimation. The theoretical background is summarized, and some hints regarding the software implementation are given in the paper. Moreover, improvements yielding an expanded functionality are proposed.

Improving strain rate estimation from velocity data of non-permanent GPS stations: the Central Apennine study case (Italy)

Abstract: An efficient procedure is proposed to define realistic lower limits of velocity errors of a non-permanent GPS station (NPS), i.e., a station where the antenna is installed and operates for short time periods, typically 10–20 days per year. Moreover, the proposed method is aimed at being independent of standard GPS data processing. The key is to subsample appropriately the coordinate time series of several continuous GPS stations situated nearby or inside the considered NPS network, in order to simulate the NPS behavior and to estimate the velocity errors associated with the subsampling procedure. The obtained data are used as lower limits to accept or correct the error estimates provided by standard data processing. The proposed approach is applied to data from the dense, non-permanent network in the Central Apennine of Italy based on a sequence of solutions for the overlapping time spans 1999–2003, 1999–2004, 1999–2005 and 1999–2007. Both the original and error-corrected velocity patterns are used to compute the strain rate fields. The comparison between the corresponding results reveals large differences that could lead to divergent interpretations about the kinematics of the study area.

The total optimal search criterion in solving the mixed integer linear model with GNSS carrier phase observations

Abstract: Existing algorithms for GPS ambiguity determination can be classified into three categories, i.e. ambiguity resolution in the measurement domain, the coordinate domain and the ambiguity domain. There are many techniques available for searching the ambiguity domain, such as FARA (Frei and Beutler in Manuscr Geod 15(4):325–356, 1990), LSAST (Hatch in Proceedings of KIS’90, Banff, Canada, pp 299–308, 1990), the modified Cholesky decomposition method (Euler and Landau in Proceedings of the sixth international geodetic symposium on satellite positioning, Columbus, Ohio, pp 650–659, 1992), LAMBDA (Teunissen in Invited lecture, section IV theory and methodology, IAG general meeting, Beijing, China, 1993), FASF (Chen and Lachapelle in J Inst Navig 42(2):371–390, 1995) and modified LLL Algorithm (Grafarend in GPS Solut 4(2):31–44, 2000; Lou and Grafarend in Zeitschrift fur Vermessungswesen 3:203–210, 2003). The widely applied LAMBDA method is based on the Least Squares Ambiguity Search (LSAS) criterion and employs an effective decorrelation technique in addition. G. Xu (J Glob Position Syst 1(2):121–131, 2002) proposed also a new general criterion together with its equivalent objective function for ambiguity searching that can be carried out in the coordinate domain, the ambiguity domain or both. Xu’s objective function differs from the LSAS function, leading to different numerical results. The cause of this difference is identified in this contribution and corrected. After correction, the Xu’s approach and the one implied in LAMBDA are identical. We have developed a total optimal search criterion for the mixed integer linear model resolving integer ambiguities in both coordinate and ambiguity domain, and derived the orthogonal decomposition of the objective function and the related minimum expressions algebraically and geometrically. This criterion is verified with real GPS phase data. The theoretical and numerical results show that (1) the LSAS objective function can be derived from the total optimal search criterion with the constraint on the fixed integer ambiguity parameters, and (2) Xu’s derivation of the equivalent objective function was incorrect, leading to an incorrect search procedure. The effects of the total optimal criterion on GPS carrier phase data processing are discussed and its practical implementation is also proposed.

Simplified equivalent multiple baseline solutions with elevation-dependent weights

Abstract: Since the assumption of all stations tracking the same satellites with identical weights was previously employed by Shen and Xu (GPS Solut 12:99–108, 2008) to derive the simplified GNSS single- and double-differenced equivalent equations, this supplementary paper expands these simplified equations in the case of each station tracking different satellites with elevation-dependent weights. Numerical experiments are performed to demonstrate the computational efficiency of the simplified equivalent algorithm relative to the traditional method in various scenarios of multi-baseline solutions with tracking different satellites. The fast computational speed of the simplified equivalent algorithm will potentially benefit the local, regional and even global GNSS multi-baseline solutions as well as the combined GNSS application.

Correction of apparent position shifts caused by GNSS antenna changes

Abstract: Antenna changes at GNSS reference stations frequently produce discontinuities in the coordinate time series. These apparent position shifts are mainly caused by changes of carrier-phase multipath effects and different errors in the antenna phase center corrections. A monitoring method was developed and successfully tested, which requires additional GNSS observations from a local, temporary reference station. Changes of carrier-phase measurement errors due to the antenna change are determined and stored in L1 and L2 phase maps. These phase maps provide corrections to be applied either to the observation data obtained before the antenna change or to the observation data obtained after the antenna change. The observation corrections are able to remove coordinate discontinuities independent of the selected coordinate estimation algorithm.

A data archive of GPS navigation messages

Abstract: Since 18 June 2007 navigation data messages transmitted by the GPS constellation are recorded by five receivers within GeoForschungsZentrum’s global groundstation network. We describe the recording, processing, validation, analysis and archiving of the navigation data. During the 197 days between 18 June 2007 and 31 December 2007 a total of 125,723,666 subframes were collected. By taking into consideration that the same data set frequently is observed by two or more receivers concurrently, 65,153,955 unique subframes could be extracted from the observations. With an estimated 88,099,200 subframes transmitted by the constellation during this time period a data yield of about 74% was achieved. Simulation studies suggest that with two additional GPS receivers, which are scheduled for addition to the network in 2008, about 95% of the transmitted subframes will be retrieved. The message data archive is open to the scientific community for non-commercial purposes and may be accessed through GFZ’s Information System and Data Center (http://isdc.gfz-potsdam.de).

Pitch and roll attitude estimation of a small-scaled helicopter using single antenna GPS with gyroscopes

Abstract: A method is presented for estimating the roll and pitch attitude of a small-scaled unmanned helicopter based on the velocity measurements of the global positioning system (GPS). The small-scaled helicopter is a radio controlled (RC) model which is readily available and affordable for academic laboratories as a research platform. Only one single antenna GPS receiver is equipped on the RC helicopter to acquire the velocity measurements needed for the attitude estimation. The velocity information is recorded by the onboard computer for post-processing. An attitude and heading reference system (AHRS) is used to provide the reference attitudes. The required angular rates and heading for this study are also given by the gyroscopes and compass of the AHRS for the sake of system’s simplification. The Kalman filter is applied to estimate the helicopter’s accelerations by using the GPS velocity measurements. The estimated accelerations form the fundamental elements of synthesizing the pseudo-roll and the pseudo-pitch. With some legitimate simplifications and assumptions, the relation between the helicopter’s attitudes and the accelerations estimated from the GPS velocity measurements can be developed. Furthermore, to enhance the accuracy of the pseudo-attitudes, the angular rates acquired from the gyroscopes are incorporated into the estimation algorithm of pseudo-attitudes by using a complementary filter.

A new multi-C/A code acquisition method for GPS

Abstract: An acquisition method is proposed which saves processing time and rapidly finds the Global Positioning System (GPS) satellite with high receiving power. The idea is to combine two or more Coarse/Acquisition (C/A) codes to process acquisition, thus called “multi-C/A code acquisition method.” The proposed method will change the procedure of conventional Fast Fourier Transform acquisition slightly, and it has the capabilities to combine with other acquisition methods to perform signal acquisition simultaneously. In order to show how this multi-C/A code acquisition method may improve the Time To First Fix at the cold start, the real raw Intermediate Frequency data from a GPS software receiver are used to validate this multi-C/A code acquisition method. The results show that this method can save at least 23% processing time and it is able to detect the satellite more rapidly.

On exact solutions of the multi-constellation GNSS navigation problem

Abstract: In Global Navigation Satellite System (GNSS) positioning, the receiver measures the pseudorange with respect to each observable navigation satellite and determines the position and clock bias. In addition to the GPS, several other navigation satellite constellations including Glonass, Galileo and Compass can/will also be used to provide positioning, navigation, and timing information. The paper is concerned with the solvability of the navigation problem when the receiver attempts to process measurements from different constellations. As two different constellations may not be time-synchronized, the navigation problem involves the determination of position of the receiver and clock bias with respect to each constellation. The paper describes an analytic approach to account for the two-constellation navigation problem with three measurements from one constellation and two measurements from another constellation. It is shown that the two-constellation GNSS navigation problem becomes the solving of a set of two simultaneous quadratic equations or, equivalently, a quartic equation. Furthermore, the zero-crossover of the leading coefficient and the sign of the discriminant of the quartic equation are shown to play a significant role in governing the solvability, i.e., the existence and uniqueness of the navigation solutions.