On the realistic stochastic model of GPS observables: Implementation and Performance

Similar documents
Noise Characteristics in High Precision GPS Positioning

USING THE INTEGER DECORRELATION PROCEDURE TO INCREASE OF THE EFFICIENCY OF THE MAFA METHOD

Approximation of ambiguity covariance matrix for integer de-correlation procedure in single-epoch GNSS positioning

Influence of ambiguity precision on the success rate of GNSS integer ambiguity bootstrapping

It is the purpose of this contribution to present the probability distribution of the 'xed' baseline. This will be done for the case the ambiguities a

ESTIMATING THE RESIDUAL TROPOSPHERIC DELAY FOR AIRBORNE DIFFERENTIAL GPS POSITIONING (A SUMMARY)

SIGMA-F: Variances of GPS Observations Determined by a Fuzzy System

Carrier-phase Ambiguity Success Rates for Integrated GPS-Galileo Satellite Navigation

The effect of an unknown data bias in least-squares adjustment: some concerns for the estimation of geodetic parameters

Least-Squares Variance Component Estimation:

Formulation of L 1 Norm Minimization in Gauss-Markov Models

ARTIFICIAL SATELLITES, Vol. 41, No DOI: /v

Carrier Phase Integer Ambiguity Resolution Recent results and open issues

Design of Adaptive Filtering Algorithm for Relative Navigation

Carrier Phase Techniques

Influences of different factors on temporal correlations of GNSS observations

Impact of Tropospheric Delay Gradients on Total Tropospheric Delay and Precise Point Positioning

ARTIFICIAL SATELLITES, Vol. 41, No DOI: /v THE LAMBDA METHOD FOR THE GNSS COMPASS

Integer least-squares theory for the GNSS compass

GPS cycle slips detection and repair through various signal combinations

Statistical characteristics of L1 carrier phase observations from four low-cost GPS receivers

On the Spectrum of the GPS DD-Ambiguities

Dynamic data processing recursive least-squares

Ionosphere influence on success rate of GPS ambiguity resolution in a satellite formation flying

Least Squares. Ken Kreutz-Delgado (Nuno Vasconcelos) ECE 175A Winter UCSD

STOCHASTIC MODELLING AND ANALYSIS OF IMU SENSOR ERRORS

Autocorrelation Functions in GPS Data Processing: Modeling Aspects

EESC Geodesy with the Global Positioning System. Class 7: Relative Positioning using Carrier-Beat Phase

On the Covariance Matrix of Weighted Total Least-Squares Estimates

v are uncorrelated, zero-mean, white

TIGHT GNSS/INS INTEGRATION AS A CONSTRAINED LEAST-SQUARES PROBLEM. David Bernal, Pau Closas, Eduard Calvo, Juan A. Fernández Rubio

THE EFFECT OF PHYSICAL CORRELATIONS ON THE AMBIGUITY RESOLUTION AND ACCURACY ESTIMATION IN GPS DIFFERENTIAL POSITIONING

Modern Navigation. Thomas Herring

On GNSS Ambiguity Acceptance Tests

SENSITIVITY ANALYSIS OF MULTIPLE FAULT TEST AND RELIABILITY MEASURES IN INTEGRATED GPS/INS SYSTEMS

Maximum Likelihood Diffusive Source Localization Based on Binary Observations

Cramér-Rao Bounds for Estimation of Linear System Noise Covariances

Test Computations

A. Barbu, J. Laurent-Varin, F. Perosanz, F. Mercier and J. Marty. AVENUE project. June, 20

Integer carrier-phase ambiguity resolution is the key to fast. Integer Aperture Estimation. A Framework for GNSS Ambiguity Acceptance Testing

Use of ground-based GNSS measurements in data assimilation. Reima Eresmaa Finnish Meteorological Institute

TOWARDS ROBUST LOCALIZATION OF RTK-GPS TOPOGRAPHIC SURVEYS 23

Huber s M-estimation in relative GPS positioning: computational aspects. 1 Introduction

Expressions for the covariance matrix of covariance data

Influence Analysis of Star Sensors Sampling Frequency on Attitude Determination Accuracy

Iteratively reweighted total least squares: a robust estimation in errors-in-variables models

estimator is a random vector. Conceptually we may write a = F (y), where a denotes the integer ambiguity vector, y the vector of observables, and F (:

SPATIAL-TEMPORAL TECHNIQUES FOR PREDICTION AND COMPRESSION OF SOIL FERTILITY DATA

OPTIMAL TIME TRANSFER

FERMENTATION BATCH PROCESS MONITORING BY STEP-BY-STEP ADAPTIVE MPCA. Ning He, Lei Xie, Shu-qing Wang, Jian-ming Zhang

EESC Geodesy with the Global Positioning System. Class 4: The pseudorange and phase observables

Recursive On-orbit Calibration of Star Sensors

Bias Correction in the Balanced-half-sample Method if the Number of Sampled Units in Some Strata Is Odd

Reliability of the Conventional Deformation Analysis Methods for Vertical Networks

Generalized Attitude Determination with One Dominant Vector Observation

Velocity (Kalman Filter) Velocity (knots) Time (sec) Kalman Filter Accelerations 2.5. Acceleration (m/s2)

Interferometric GPS ambiguity resolution

Analysis of incremental RLS adaptive networks with noisy links

at least 50 and preferably 100 observations should be available to build a proper model

Publication IV Springer-Verlag Berlin Heidelberg. Reprinted with permission.

Estimation and Prediction Scenarios

2938 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 46, NO. 11, NOVEMBER Integer Parameter Estimation in Linear Models with Applications to GPS

Constructing high-resolution, absolute maps of atmospheric water vapor by combining InSAR and GNSS observations

Modified Leaky LMS Algorithms Applied to Satellite Positioning

CONTROL SYSTEMS, ROBOTICS, AND AUTOMATION - Vol. V - Recursive Algorithms - Han-Fu Chen

Statistics 910, #5 1. Regression Methods

On the reliability of errors-in-variables models

Simplified Filtering Estimator for Spacecraft Attitude Determination from Phase Information of GPS Signals

Deformation measurement using SAR interferometry: quantitative aspects

Geodetic Network Adjustment Examples

DATA ASSIMILATION FOR FLOOD FORECASTING

A New Kalman Filtering Method to Resist Motion Model Error In Dynamic Positioning Chengsong Zhoua, Jing Peng, Wenxiang Liu, Feixue Wang

Kinematic Orbit Determination of Low Earth Orbiting Satellites using GPS Observables

Quadratic Extended Filtering in Nonlinear Systems with Uncertain Observations

A STATE-SPACE APPROACH FOR THE ANALYSIS OF WAVE AND DIFFUSION FIELDS

Satellite Navigation error sources and position estimation

VISION TRACKING PREDICTION

A NOVEL OPTIMAL PROBABILITY DENSITY FUNCTION TRACKING FILTER DESIGN 1

IMPRECISION IN GEODETIC OBSERVATIONS CASE STUDY GPS MONITORING NETWORKS

Estimating and Testing the US Model 8.1 Introduction

Estimation of Troposphere Decorrelation Using the Combined Zenith-dependent Parameter

Least-Squares Spectral Analysis Theory Summary

Musical noise reduction in time-frequency-binary-masking-based blind source separation systems

GPS 10, GPS , GPS. (J PL ) Zumbeger. mm, cm ; ( Global RT K) [3 ]

GPS Multipath Detection Based on Sequence of Successive-Time Double-Differences

1 EM algorithm: updating the mixing proportions {π k } ik are the posterior probabilities at the qth iteration of EM.

ADAPTIVE INVERSE CONTROL BASED ON NONLINEAR ADAPTIVE FILTERING. Information Systems Lab., EE Dep., Stanford University

S-system theory applied to array-based GNSS ionospheric sensing

1 Introduction to Generalized Least Squares

Least squares: introduction to the network adjustment

8.2 Harmonic Regression and the Periodogram

Nonlinear Parameter Estimation for State-Space ARCH Models with Missing Observations

Improving Adaptive Kalman Estimation in GPS/INS Integration

3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE

Performance of a Dynamic Algorithm For Processing Uncorrelated Tracks

Evaluation of Some Techniques for Forecasting of Electricity Demand in Sri Lanka

Evaluating the Consistency of Estimation

Econ 424 Time Series Concepts

Prediction-based adaptive control of a class of discrete-time nonlinear systems with nonlinear growth rate

Prof. Dr. Roland Füss Lecture Series in Applied Econometrics Summer Term Introduction to Time Series Analysis

Transcription:

he International Archives of the Photogrammetry, Remote Sensing Spatial Information Sciences, Volume XL-/W5, 05 International Conference on Sensors & Models in Remote Sensing & Photogrammetry, 3 5 Nov 05, Kish Isl, Iran On the realistic stochastic model of GPS observables: Implementation Performance F. Zangeneh-Nead a, *, A. R. Amiri-Simooei b, M. A. Sharifi a, J. Asgari b a School of Surveying Geosptial Engineering, Research Institute of Geoinformation echnology (RIG), College of Engineering, University of ehran, Iran- (f.zangenehnead, sharifi@ut.ac.ir) c Department of Geomatics Engineering, Faculty of Engineering, University of Isfahan, 8746-7344 Isfahan, Iran- (amiri, asgari@eng.ui.ac.ir) KEY WORDS: GPS stochastic model, Least-squares variance component estimation (LS-VCE), noise assessment, integer ambiguity resolution (IAR), Success rate, GPS observables. ABAC: High-precision GPS positioning requires a realistic stochastic model of observables. A realistic GPS stochastic model of observables should tae into account different variances for different observation types, correlations among different observables, the satellite elevation dependence of observables precision, the temporal correlation of observables. Least-squares variance component estimation (LS-VCE) is applied to GPS observables using the geometry-based observation model (GBOM). o model the satellite elevation dependent of GPS observables precision, an exponential model depending on the elevation angles of the satellites are also employed. emporal correlation of the GPS observables is modelled by using a first-order autoregressive noise model. An important step in the high-precision GPS positioning is double difference integer ambiguity resolution (IAR). he fraction or percentage of success among a number of integer ambiguity fixing is called the success rate. A realistic estimation of the GNSS observables covariance matrix plays an important role in the IAR. We consider the ambiguity resolution success rate for two cases, namely a nominal a realistic stochastic model of the GPS observables using two GPS s collected by the rimble R8 receiver. he results confirm that applying a more realistic stochastic model can significantly improve the IAR success rate on individual frequencies, either on L or on L. An improvement of 0% was achieved to the empirical success rate results. he results also indicate that introducing the realistic stochastic model leads to a larger stard deviation for the baseline components by a factor of about.6 on the s considered.. INRODCION GPS data processing is usually performed by the least squares adustment method for which both the functional stochastic models must be correctly specified. he functional model, describing the mathematical relation between the observations the unnown parameters, is usually well nown either in the relative positioning or in the single/precise point positioning (Seeber, 003; Hofmann-Wellenhof et al., 008; Leic, 004; eunissen Kleusberg, 998; Rizos, 997). However, the stochastic model, expressing the statistical characteristics of the GPS observations by means of a covariance matrix, is still a challenging problem. Misspecification in the stochastic model leads to unreliable suboptimal estimates. he realistic stochastic model should thus be utilized aiming at obtaining reliable least squares estimates. here are a few common stochastic models for GPS observables: () Equal-weight stochastic model in which identical variances for each measurement type are chosen (i.e. the same variance for the code the same variance for the phase observations). his structure ignores the correlations among different observables, () Satellite elevation-angle dependent model in which the observations weighting procedure is performed using trigonometric or exponential functions. he observations at lower elevation angles have larger variances than those at higher elevation angles, (3) Signal-to-noise ratio (SNR) based model in which the observations weighting procedure depends on SNR values, a quality indicator for the GPS observations. hese models are in fact simple rudimentary stochastic models. Satirapod Wang (000) found that in some cases both the use of SNR satellite elevation angle information failed to reflect the reality of data quality. herefore, it still remains necessary to investigate rigorous methods for constructing a more reliable covariance matrix of the GPS observables. An appropriate stochastic model for GPS observables may include different variances for each observation type, the correlation among different observables, the satellite elevation dependence of observables precision, the possible temporal correlation of the GPS observables (Amiri-Simooei et al., 009, 03, 05). herefore, we first define a general form of the covariance matrix of the GPS observables, which includes the abovementioned variance covariance components into the stochastic model. his general form of the covariance matrix is expressed as an unnown linear combination of nown cofactor matrices. he estimation of the unnown (co)variance parameters is referred to as variance component estimation (VCE). here are many different methods for VCE among which we mae use of the least-squares variance component estimation (LS-VCE), in a straightforward manner, to assess the noise characteristics of the GPS observables. LS-VCE was originally developed by eunissen (988). Further elaboration development of the LS-VCE are provided by eunissen * Corresponding author doi:0.594/isprsarchives-xl--w5-755-05 755

he International Archives of the Photogrammetry, Remote Sensing Spatial Information Sciences, Volume XL-/W5, 05 International Conference on Sensors & Models in Remote Sensing & Photogrammetry, 3 5 Nov 05, Kish Isl, Iran Amiri-Simooei (008) Amiri-Simooei (007). he geometry-based observation model (GBOM), employed as the functional model, is the widely used mathematical model for high precision GPS positioning from code phase observations using a relative GPS receiver setup (Odi, 008). P ( t ) ( t ) e ( t ) () s, s, s, r,, L i r,, L i r,, L i An important step in the high-precision GPS positioning is double difference (DD) integer ambiguity resolution (IAR). he fraction or percentage of success among a number of integer ambiguity fixing is called the success rate. A realistic estimation of the GNSS observables covariance matrix plays an important role in the IAR. When compared to a realistic stochastic model, a simple stochastic model results in poorer precision of the estimated parameters hence a lower integer ambiguity success rate. In this paper, the reliability of the resolved ambiguities is investigated through the computation of the IAR success rate for two cases, namely a nominal a realistic stochastic model of the GPS observables. For the nominal case, the same stard deviations for the phase observables on L L the same stard deviations for the code observables on L L is assumed as 3 mm P 30 cm P P, respectively. his structure ignores the correlation among different observables, the satellite elevation dependence of observables precision, the possible temporal correlation of the GPS observables. where ( t ) ( t ) a ( t ) () s, s, s, s, r,, L i r,, L i L r,, L r,, L i (.) s r is an abbreviation for s s s (.) (.) (.) ((.) (.) ), L corresponds to either L r r r or L frequency, P is the DD pseudorange observation on the L or L frequency, denotes the DD receiver-satellite range, indicates the pseudorange measurement errors on the L or e L frequency t i refers to the measurement time instant. In Eq. (), indicates the DD carrier phase observation on the L or L frequency, is the corresponding wavelength, a denotes the DD integer ambiguities on L or L, is the phase measurement errors on L or L. he DD pseudorange carrier-phase observation equations on the L or L frequency can be written in the following matrix form: he baseline components uncertainties expressed by the estimated covariance matrix of the unnowns depend on the choice made for the weight matrix of the observables. We thus assess the effect of the stochastic model on the efficiency of the estimated baseline uncertainties. o this end, the baseline uncertainties are also computed for two cases, nominal a realistic stochastic model of the GPS observables. s, P ( t ) r,, i s, g P ( t ) r,, i s, E A a s, r,, ( t ) r,, i s, a s, r,, ( t ) r,, i (3) his paper is organized as follows. Section provides a brief review on a few theoretical concepts followed in this research. his section reviews the functional part of the GPS GBOM, a realistic stochastic model of GPS observables. Section 3 presents the theory of the LS-VCE method its application to determine an appropriate stochastic model. In section 4, to investigate the performance of the realistic stochastic model obtained by the LS-VCE approach, two GPS baseline s of short baselines are employed. he reliability of the resolved ambiguities is then investigated through the computation of the IAR success rate for two cases, namely a nominal a realistic stochastic model of the GPS observables. We also assess the effect of the stochastic model on a realistic stard deviation for the baseline components. We further highlight the effect of an optimal stochastic model in other fields of GPS applications such as the precise point positioning method. Finally, we mae conclusions in section 6.. GPS GEOMERY-BASED OBSERVAION MODEL (GBOM). Functional model Consider two receivers r simultaneously trac the same satellites s. he functional part of GBOM is based on the DD GPS observables. Ignoring the DD atmospheric delays, which is a valid assumption for zero very short baselines, the observation equation for the DD GPS pseudorange (code) carrier phase observations are respectively defined as (eunissen Kleusberg, 998) relating to two satellites which were observed by two receivers. In Eq. (3), E sts for the mathematical expectation operator. he expectation of pseudorange carrier-phase measurement errors is assumed to be zero. In Eq. (3), the design matrix A is of the form G G A G G 0 0 0 0 0 0 s in which the geometry matrix G is G ( u u ), where u ( r r ) / r r is the unit direction vector between receiver satellite (line-of-sight vector). he vectors r r (4) are the geocentric vector of the satellite receiver, respectively. In Eq. (3), vector g consists of the unnown baseline components between the reference rover receivers, which is defined as g X Y Z. A Realistic Stochastic Model. r r r A realistic stochastic model for DD GPS observables is of the form (Amiri-Simooei et al., 009, 03, 05) doi:0.594/isprsarchives-xl--w5-755-05 756

he International Archives of the Photogrammetry, Remote Sensing Spatial Information Sciences, Volume XL-/W5, 05 International Conference on Sensors & Models in Remote Sensing & Photogrammetry, 3 5 Nov 05, Kish Isl, Iran D Q (5) y C E where indicates the Kronecer product of two matrices. For four observation types C covariance components of the form C includes 0 unnown variance P P, P P, P, P, P P P, P, P, P,, P P,,, (6) Assume the linear model of observation equations as E( y) Ax; D( y) Q Q (9) p y where D denotes the mathematical dispersion operator, y is the observation vector of order m, is the n-vector of unnown parameters, A is the design matrix describing the functional relation between the observations vector the unnown parameters,,..., p are the unnown m n (co)variance components. Also, x Q y is the covariance matrix of the observables expressed as an unnown linear combination of nown cofactor matrices ( Q,,..., p ). he diagonal off-diagonal elements of C are the variances the covariance between the GPS observables, respectively. Also, a oeplitz matrix represents temporal correlation of the GPS observables. It is a K K matrix, where K denotes the number of the observation epochs. may consists of K number of unnown parameters as follows Finally, E is a (0) () ( K ) () (0) ( K ) ( K) ( K) (0) (7) matrix representing the dependence of GPS observables precision on satellite elevation, where is the number of satellites. he satellite number [] is assumed as the reference satellite. E E is of the form [] [ ] [] [] [] [] [3] [] [] [] [] [ ] he stochastic model presented in Eq. (5) considers all items of a realistic covariance matrix of GPS observables pointed out in the introduction. In this paper, we estimate the unnown (co)variance components using LS-VCE which is briefly presented in the next section. 3. SOCHASIC MODEL ESIMAION 3. Least-Squares Variance Component Estimation (LS- VCE) he LS-VCE was originally developed by eunissen (988). Further elaboration development of the method is provided by eunissen Amiri-Simooei (008) Amiri-Simooei (007). (8) Based on Eq. (9), the covariance matrix of the observables is partly unnown as it is expressed as an unnown linear combination of nown cofactor matrices. o estimate the unnown (co)variance components, different VCE methods have been invented. Among them the LS-VCE method is employed in the present contribution. he unnown (co)variance components can be estimated using LS-VCE as p -vector l n ˆ N l, where the p p are respectively obtained as i tr( Q Q P Q Q P i y A y A ) matrix N the (0) ˆ e Q Q Q e i y i y ˆ l () where eˆ P y is the m-vector of the least squares residuals, in which proector, A A m y y P I A( A Q A) A Q is an orthogonal tr denotes the trace of a matrix. N o compute the matrix, the vector the residuals vector ê the covariance matrix of the observables is required. LS-VCE is thus an iterative method that needs the approximate values of the (co)variance components. he iteration proceed until the converge criterion of ˆ l ˆ is met, i.e., the difference of two consecutive solutions is less than the tolerance. he covariance matrix of the estimated (co)variance components can automatically be obtained as Qˆ N. When the number of observations is large, we can divide the entire observations into a few (r) groups. One can show that the unnown (co)variance parameters ( ; : p) can be separately estimated for each group the final estimates are obtained by averaging the groupwise estimates over all (r) groups, i.e., r ( i ) ˆ ˆ,,..., ( ) p where ˆ i is the r i doi:0.594/isprsarchives-xl--w5-755-05 757

he International Archives of the Photogrammetry, Remote Sensing Spatial Information Sciences, Volume XL-/W5, 05 International Conference on Sensors & Models in Remote Sensing & Photogrammetry, 3 5 Nov 05, Kish Isl, Iran th co(variance) components of the i th group (Amiri-Simooei et al., 009). 3.. Estimation of components of E : After estimating the 3. Determination of a realistic stochastic model of GPS observables We now investigate the noise characteristics of GPS observables pointed out in the introduction. For this purpose, we apply LS-VCE to the GPS observables. he LS-VCE method may be applied to GPS GBOM as the functional model in three ways, namely Float ambiguities, Fixed ambiguities Fixed ambiguities fixed baseline (Amiri-Simooei et al. 03, 05). In this paper, we carried out the Fixed ambiguities method. After fixing the DD integer ambiguities via the leastsquares ambiguity decorrelation adustment (LAMBDA) method (eunissen, 993, 995) or the lattice theory (Jazaeri et al., 0), the fixed DD integer ambiguities were moved to the left h side of Eq. (3). In this case, A simplifies to A [ G G G G ]. Equation (5) presents a realistic stochastic model of GPS observables. Estimation of the (co)variance components of this model may be performed in there steps () Estimation of the (co)variance components of. components E C, () Estimation of the, (3) Estimation of the components of For more details on the application of LS-VCE to noise assessment of the GPS observables, the reader is referred to Amiri-Simooei et al. (009, 03, 05). 3.. Estimation of (co)variance components of first stage, the (co)variance components C C : In the including the different variances of observation types the possible correlation among the GPS observables are estimated. For this purpose, the temporal correlation the satellites elevation dependence of GPS observables precision are ignored. On the other h, is assumed as an identity matrix of size K, i.e.,, I K E E simplifies to 4 4 4 (5) meaning that all satellites with different elevation angles have equal variance factors (i.e., ). [] [ ] [ ] By using the structure introduced in Eq. (5) for I, the (co)variance components of K C E can be estimated by LS-VCE. he approximate values for the covariance components was assumed to be zero, un-correlated observables, the approximate values for the stard deviation of GPS code (C P) phase observations (L L) were taen as 30 cm 3 mm, respectively. Also, the tolerance for the convergence of the LS-VCE solution is set to = 0-6. (co)variance components of C, we may consider the satellites elevation dependence of GPS observables precision by estimating the components of matrix C E using LS-VCE. o this end,, which is already estimated, is introduced into the stochastic model of Eq. (5). he temporal correlation is still assumed to be absent at this stage, i.e., components of E I K. he can then be estimated using LS-VCE. Also, the approximate values for the variance factors of assumed identical for all satellites as. [] [ ] [ ] o consider the dependence of GPS observables precision on the satellites elevation, an exponential elevation-dependent model may also be used. his model is similar to the Euler Goad elevation-dependent model (Euler Goad, 99), is expressed as follows where E 0 a a e (6) is the variance of the observables, two unnown parameters to be estimated, elevation angle of the satellites 0 a a are are denotes the is an unnown reference elevation angle. he unnown parameters of the nonlinear model given in (6) can be obtained by a nonlinear least squares fit to the estimated variance components of E using LS-VCE. herefore, the exponential elevation-dependent model was implemented to express the dependence of the GPS observable precision on the satellite elevation angle (Amiri-Simooei et al., 05). 3..3 Estimation of components of : In the last step, time correlation of the GPS observables is investigated. At this stage, the components of purpose, are to be estimated. For this matrices, which were estimated in the two previous steps, are assumed nown hence introduced into Eq. (5). he components of can then be estimated by using the empirical autocovariance function (ACF), which is shown to be a by product of the LS-VCE method, if the weight matrix is chosen as the identity matrix. he unbiased sample autocovariance function of a zero-mean stationary time-series can be obtained as (Amiri-Simooei 007, p. 4) K ee ˆˆ i i i ˆ 0,,..., K K (7) where is the time lag (the time interval between two samples) in seconds, is the covariance among the observables with doi:0.594/isprsarchives-xl--w5-755-05 758

he International Archives of the Photogrammetry, Remote Sensing Spatial Information Sciences, Volume XL-/W5, 05 International Conference on Sensors & Models in Remote Sensing & Photogrammetry, 3 5 Nov 05, Kish Isl, Iran lag 0 K is the number of epochs. he coefficient at lag equals the diagonal elements (variance factors) of As a predefined model, the temporal correlation of the GPS observables is approximated by an exponential auto covariance function (first-order auto regressive, AR()) as e, 0,,..., K (8) where are the unnown parameters of the model. he unnown parameters of the nonlinear model given in (8) can be obtained by a nonlinear least squares fit to the estimated ˆ obtained by Eq. (7). herefore, the exponential auto covariance function was implemented to express the temporal correlation of the GPS observables. 4. NUMERICAL RESULS AND DISCUSSIONS o investigate the effect of the realistic stochastic model of the GPS observables on the IAR success rate on the estimated baseline uncertainties, two real GPS s collected by a geodetic GPS receiver, namely rimble R8, are employed. he receiver used is a geodetic dual frequency receiver which provides the code phase observations on both the L L frequencies (namely, C-P-L-L). he first GPS data experiment was collected on November 8, 00 with one second interval over a short baseline of.589 m from 3::9 to 3:5:8 UC consisting of 50 epochs. he second was also measured over a short baseline of 0.56 m with one second interval on November 8th, 00, from 3:7:8 to 3::7 UC consisting of 00 epochs. his geodetic GPS receiver has been recently used in Amiri- Simooei et al. (05) for two different GPS s. We now use other datasets using the same receiver under a different situation, for example different days different satellite configuration. However, we use the stard deviation of the phase code observations as well as the correlation coefficient among different observations as those estimated by Amiri-Simooei et al. (05). According to Amiri-Simooei et al. (05), the estimated stard deviation of the four GPS observation types are ˆ 39. mm, ˆ 65.4 mm, ˆ 5. mm C P L ˆ 4.94 mm. able provides the covariances L correlation coefficients between the GPS observables. he results indicate a significant correlation between observations. A positive correlation of 0.4 is observed between the phase observations on L L for rimble R8, which is considered to be significant (see Amiri-Simooei et al., 03). he estimated values of the parameters given in (6) are a 0.6, a 8.90 the reference angle is estimated as 7.6, 0 for the rimble R8 receiver. able. Estimated covariance ˆ ( mm ) i, correlation coefficient ˆi, among different observations for rimble R8 Between ˆ ( mm ) i, ˆi,. C-P 4500.39 0.8 C-L 385.59 0.9 C-L 63.99 0.4 P-L 99.83 0.06 P-L 50.35 0.08 L-L 0.69 0.4 After introducing the estimated components of C the mathematical model of the GPS baseline (ignoring time correlation of GPS observables in this stage), the integer ambiguity success rate () is then computed compared with its nominal counterpart for each of the GPS s for three following cases: Case (): L only; Case (): L only Case (3): L + L. he empirical success rate computed for three cases are given in able (). he results indicate that applying the estimated components of C E E to to the mathematical model of the GPS baselines can significantly improve the success rate; an improvement of 0% is clear when using only L or L. able. Estimated empirical success rate () for each of GPS s in three cases. Data set L only L only L+L First Second Nominal Realistic Nominal Realistic 7.8 % 6.0 % 00 % 8.6 % 86.4 % 00 % 8.0 % 76.4 % 00 % 83.6 % 87.6 % 00 % Using different GPS s, Amiri-Simooei et al. (05) also showed that applying a more realistic stochastic model can significantly improve the IAR success rate on individual frequencies, either on L or on L. his wor can be considered as a follow-up to the wor carried out by Amiri-Simooei et al. (05) in which some supplementary results, for example considering the effect of introducing the realistic stochastic model on the estimation of the two baseline components uncertainties, are also presented in this contribution. We now aim to consider the effect of introducing the realistic stochastic model on the estimation of the baseline components uncertainties. he baseline uncertainties can be regarded as the square root of the diagonal components of the unnowns covariance matrix approximated by Q ( A Q A) ˆ in which Q y can be considered in three different cases: Case I: Nominal stochastic model of the GPS observables, Case II: A realistic stochastic model of the GPS observables by introducing the matrices C E estimated by the LS-VCE method, x y doi:0.594/isprsarchives-xl--w5-755-05 759

he International Archives of the Photogrammetry, Remote Sensing Spatial Information Sciences, Volume XL-/W5, 05 International Conference on Sensors & Models in Remote Sensing & Photogrammetry, 3 5 Nov 05, Kish Isl, Iran Case III: Considers the C E as estimated in Case II. In addition, the realistic stochastic model of the GPS observables considers time correlation of GPS observables based on the empirical ACF fitting an exponential auto covariance function to the ACF. able (3) provides the estimated baseline uncertainties in the three cases for each of the GPS s. he results indicate that the estimated baseline uncertainties in the case of the nominal stochastic model lead to the optimistic uncertainties of the estimated baseline components when compared to a realistic stochastic model (cases II III). In other words, introducing the realistic stochastic model results in a larger stard deviation for the baseline components by a factor of about.6 on the s considered. able 3. Estimated baseline uncertainties in three cases for two s: Case I: Nominal, Case II: considering Case III: considering Data set First Second Case C X, E (mm) Y. (mm) C Z I 0.7 0.34 0.36 II 0.35 0.50 0.43 III 0.66 0.85 0.90 I 0.7 0.36 0.38 II 0.34 0.49 0.45 III 0.76.03.09 E (mm) In this paper, we investigated the effect of an optimal stochastic model in the case of relative positioning. he choice of a realistic stochastic model of the GPS observables is also important in other fields of GPS applications such as the precise point positioning (PPP) technique. 5. CONCLUSIONS A realistic stochastic model for the GPS observables should include a few important issues lie different variances for each GPS observation type, the correlation between different observations, the satellite elevation dependence of the GPS observables weights, the temporal correlation of the GPS observables. he use of a proper (correct) stochastic model of the GPS observables leads to the best linear unbiased estimators (BLUE) in the high precision GNSS positioning. he LS-VCE method was applied to GPS observables using the GBOM as the functional model. he impact of a realistic GPS stochastic model on IAR success rate was discussed using two GPS short baseline s. he empirical success rate was computed for two cases: () using the nominal covariance matrix, () using the estimated covariance matrix by LS-VCE. For the latter case, the success rate is determined by taing into consideration the variances of the GPS observables, the covariance/correlations of GPS observables, the satellites elevation dependence of the GPS observable precision. he results showed that applying a more realistic stochastic model can significantly improve the IAR success rate on individual frequencies, either on L or on L. An improvement of 0% was achieved to the success rate results. his idea can thus be implemented in attitude determination using single-frequency single-epoch of GPS observations in order to improve /or speed up the ambiguity fixing procedure. Finally, the effect of the realistic stochastic model on the estimated baseline uncertainties was investigated. he results indicated that introducing the realistic stochastic model leads to a larger stard deviation for the baseline components by a factor of about.6 on the s considered. ACKNOWLEDGEMENS We would lie to acnowledge Dr. Shahram Jazaeri for providing us with the GPS s we used in this paper. REFERENCES Amiri-Simooei, AR., 007. Least-squares variance component estimation: heory GPS applications. Ph.D. thesis, Delft Univ. of echnology, Delft, he Netherls. Amiri-Simooei, AR., eunissen, PJG., iberius, CCJM., 009. Application of least-squares varinace component estimation to GPS observables. J Surv Eng., 35(4):49-60. Amiri-Simooei, AR., Zangeneh-Nead, F., Asgari, J., 03. Least-squares variance component estimation applied to GPS geometry-based observation model. J Surv Eng., 39(4):76-87. Amiri-Simooei, AR., Zangeneh-Nead, F., Jazaeri, S., Asgari, J., 05. Role of stochastic model on GPS integer ambiguity resolution success rate, GPS Solutions, DOI 0.007/s09-05-0445-5. Euler, HJ., Goad, C., 99. On optimal filtering of GPS dual frequency observations without using orbit information. Bull. Geod, 65: 30 43. Hofmann-Wellenhof, B., Lichtenegger, H., Wasle, E., 008. GNSS Global Navigation Satellite Systems GPS, GLONASS, Galileo & more. Springer, Wien New Yor. Jazaeri, S., Amiri-Simooei, AR., Sharifi, MA., 0. Fast integer least-squares estimation for GNSS high-dimensional ambiguity resolution using lattice theory. J Geodesy, 86():3 36. Leic, A., 004. GPS Satellite Surveying. 3rd edition, John Wiley, New Yor. Odi, D., 008. GNSS Solutions: What does geometry-based geometry-free mean in the context of GNSS? Inside GNSS, 3(), 4. Rizos, C., 997. Principles practice of GPS surveying. Monograph 7, School of Geomatic Engineering, he University of New South Walse, Sydney. Satirapod, C., Wang, J., 000. Comparing the Quality Indicators of GPS Carrier Phase Observations, Geomatics Research Australasia. Seeber, G., 003. Satellite Geodesy. De Gruyter, Berlin. doi:0.594/isprsarchives-xl--w5-755-05 760

he International Archives of the Photogrammetry, Remote Sensing Spatial Information Sciences, Volume XL-/W5, 05 International Conference on Sensors & Models in Remote Sensing & Photogrammetry, 3 5 Nov 05, Kish Isl, Iran eunissen, PJG., 993. Least squares estimation of the integer GPS ambiguities. Proc., IAG General Meeting, Series No. 6, Delft Geodetic Computing Centre. eunissen, PJG., 995. he least-squares ambiguity decorrelation adustment: A method for fast GPS integer ambiguity estimation. J Geodesy, 70( ):65 83. eunissen, PJG., 988. owards a least-squares framewor for adusting testing of both functional stochastic model. Internal research memo, Geodetic Computing Centre, Delft. A reprint of original 988 report is also available in 004, No. 6 eunissen, PJG., Kleusberg, A., 998. Section 7: Quality control GPS. GPS for Geodesy, PJG, eunissen A, Kleusberg, eds., Springer, Berlin. eunissen, PJG., Amiri-Simooei, AR., 008. Least-squares variance component estimation. J. Geodesy, 8(), 65 8. doi:0.594/isprsarchives-xl--w5-755-05 76