Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Skip to main content

Fast error analysis of continuous GNSS observations with missing data

  • Original Article
  • Published:
Journal of Geodesy Aims and scope Submit manuscript

Abstract

One of the most widely used method for the time-series analysis of continuous Global Navigation Satellite System (GNSS) observations is Maximum Likelihood Estimation (MLE) which in most implementations requires \(\mathcal{O }(n^3)\) operations for \(n\) observations. Previous research by the authors has shown that this amount of operations can be reduced to \(\mathcal{O }(n^2)\) for observations without missing data. In the current research we present a reformulation of the equations that preserves this low amount of operations, even in the common situation of having some missing data.Our reformulation assumes that the noise is stationary to ensure a Toeplitz covariance matrix. However, most GNSS time-series exhibit power-law noise which is weakly non-stationary. To overcome this problem, we present a Toeplitz covariance matrix that provides an approximation for power-law noise that is accurate for most GNSS time-series.Numerical results are given for a set of synthetic data and a set of International GNSS Service (IGS) stations, demonstrating a reduction in computation time of a factor of 10–100 compared to the standard MLE method, depending on the length of the time-series and the amount of missing data.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Subscribe and save

Springer+ Basic
$34.99 /Month
  • Get 10 units per month
  • Download Article/Chapter or eBook
  • 1 Unit = 1 Article or 1 Chapter
  • Cancel anytime
Subscribe now

Buy Now

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6

Similar content being viewed by others

References

  • Agnew DC (1992) The time-domain behaviour of power-law noises. Geophys Res Lett 19(4):333–336

    Article  Google Scholar 

  • Altamimi Z, Collilieux X, Mtivier L (2011) ITRF2008: an improved solution of the international terrestrial reference frame. J Geod 85(8):457–473. doi:10.1007/s00190-011-0444-4

    Article  Google Scholar 

  • Amiri-Simkooei AR (2009) Noise in multivariate GPS position time-series. J Geod 83:175–187. doi:10.1007/s00190-008-0251-8

    Article  Google Scholar 

  • Ammar GS, Gragg WB (1988) Superfast solution of real positive definite Toeplitz systems. SIAM J Matrix Anal Appl 9:61–76

    Article  Google Scholar 

  • Bos MS, Fernandes RMS, Williams SDP, Bastos L (2008) Fast error analysis of continuous GPS observations. J Geod 82:157–166. doi:10.1007/s00190-007-0165-x

    Google Scholar 

  • Brualdi RA, Schneider H (1983) Determinantal identities: Gauss, Schur, Cauchy, Sylvester, Kronecker, Jacobi, Binet, Laplace, Muir, and Cayley. Linear Algebra Appl 52(53):769–791

  • Dow JM, Neilan RE, Rizos C (2009) The International GNSS Service in a changing landscape of Global Navigation Satellite Systems. J Geod 83:191–198. doi:10.1007/s00190-008-0300-3

    Article  Google Scholar 

  • Gohberg IC, Semencul AA (1972) On the inversion of finite Toeplitz matrices and their continuous analogs. Mat Issled 7(2):201–223

    Google Scholar 

  • Hackl M, Malservisi R, Hugentobler U, Wonnacott R (2011) Estimation of velocity uncertainties from GPS time series: examples from the analysis of the South African TrigNet network. J Geophys Res 116(B15):B11404

    Google Scholar 

  • Hosking JRM (1981) Fractional differencing. Biometrika 68:165–176

    Article  Google Scholar 

  • Johnson HO, Agnew DC (1995) Monument motion and measurements of crustal velocities. Geophys Res Lett 22(21):2905–2908

    Article  Google Scholar 

  • Kasdin NJ (1995) Discrete simulation of colored noise and stochastic processes and \(1/f^\alpha \) power-law noise generation. Proc IEEE 83(5):802–827

    Article  Google Scholar 

  • Langbein J (2004) Noise in two-color electronic distance meter measurements revisited. J Geophys Res 109(B04406). doi:10.1029/2003JB002819

  • Langbein J (2010) Computer algorithm for analyzing and processing borehole strainmeter data. Comput Geosci 36(5):611–619. doi:10.1016/j.cageo.2009.08.011

    Article  Google Scholar 

  • Langbein J (2012) Estimating rate uncertainty with maximum likelihood: differences between power-law and flicker-random-walk models. J Geod 86:775–783. doi:10.1007/s00190-012-0556-5

    Article  Google Scholar 

  • Mao A, Harrison CGA, Dixon TH (1999) Noise in GPS coordinate time series. J Geophys Res 104(B2):2797–2816. doi:10.1029/1998JB900033

    Google Scholar 

  • Santamaría-Góme A, Bouin MN, Collilieux X, Wöppelmann G (2011) Correlated errors in GPS position time series: implications for velocity estimates. J Geophys Res 116(B15):B01405

    Google Scholar 

  • Trench WF (1964) An algorithm for the inversion of finite Toeplitz matrices. J Soc Indust Appl Math 12:515–522

    Article  Google Scholar 

  • Webb FH, Zumberge JF (1995) An introduction to GIPSY/OASIS-II. Technical report JPL D-11088, California Institute of Technology, Pasadena, CA

  • Williams SDP (2003) The effect of coloured noise on the uncertainties of rates from geodetic time series. J Geod 76(9–10):483–494. doi:10.1007/s00190-002-0283-4

    Article  Google Scholar 

  • Williams SDP (2008) CATS: GPS coordinate time series analysis software. GPS Solutions 12(2):147–153. doi:10.1007/s10291-007-0086-4

  • Williams SDP, Bock Y, Fang P, Jamason P, Nikolaidis RM, Prawirodirdjo L, Miller M, Johnson DJ (2004) Error analysis of continuous GPS position time series. J Geophys Res 109(B03412). doi:10.1029/2003JB002741

  • Zhang J, Bock Y, Johnson H, Fang P, Williams S, Genrich J, Wdowinski S, Behr J (1997) Southern California Permanent GPS Geodetic Array: error analysis of daily position estimates and site velocities. J Geophys Res 102(B8):18035–18055

    Google Scholar 

  • Zumberge JF, Heflin MB, Jefferson DC, Watkins MM, Webb FH (1997) Precise point positioning for the efficient and robust analysis of GPS data from large networks. J Geophys Res 102(B3):5005–5018. doi:10.1029/96JB03860

    Google Scholar 

Download references

Acknowledgments

The authors would like to thank the assistant editor Jeff Freymueller and three anonymous reviewers for their thorough and constructive review of our manuscript. This work was funded by national funds through FCT in the scope of the Project PesTC/ Mar/LA0015/2011.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to M. S. Bos.

Appendix: Derivation of key equations

Appendix: Derivation of key equations

In this appendix we give the derivations of Eqs. (7) and (8) that were presented in Sect. 2. These derivation make use of Jacobi’s determinant identity and the formula for the inverse of a block matrix (Brualdi and Schneider 1983). To start, assume that our square covariance matrix \(\mathbf C \), size \(n\times n\), can be partitioned into four sub-matrices:

$$\begin{aligned} \mathbf{C } = \begin{pmatrix} \mathbf{C }_{oo}&\mathbf{C }_{om} \\ \mathbf{C }_{mo}&\mathbf{C }_{mm} \end{pmatrix} \end{aligned}$$
(33)

where the subscripts \(m\) and \(o\) represent the rows and columns of missing and observed data respectively. Consequently, \(\mathbf{C }_{mm}\) is a square \(m\times m\) matrix.

Using Gaussian elimination we can reduce \(\mathbf{C }_{mo}\) to a zero matrix:

$$\begin{aligned} \begin{pmatrix} \mathbf{I }^{n-m}&\mathbf{0 } \\ -\mathbf{C }_{mo}\mathbf{C }_{oo}^{-1}&\mathbf{I }^{m} \end{pmatrix} \begin{pmatrix} \mathbf{C }_{oo}&\mathbf{C }_{om} \\ \mathbf{C }_{mo}&\mathbf{C }_{mm} \end{pmatrix} = \begin{pmatrix} \mathbf{C }_{oo}&\mathbf{C }_{om} \\ \mathbf{0 }&\mathbf{C }^{\prime }_{mm} \end{pmatrix} \end{aligned}$$
(34)

where \(\mathbf{C }^{\prime }_{mm}=\mathbf{C }_{mm}-\mathbf{C }_{mo}\mathbf{C }_{oo}^{-1} \mathbf{C }_{om}\) which is called the Schur complement of \(\mathbf{C }_{oo}\) in \(\mathbf C \). When we take the determinant of Eq. (34), we obtain:

$$\begin{aligned} \text{ det}(\mathbf{C })=\text{ det}(\mathbf{C }_{oo})\,\text{ det}(\mathbf{C }^{\prime }_{mm}) \end{aligned}$$
(35)

Using the matrices \(\breve{\mathbf{C }}\) (=\(\mathbf{C }_{oo}\)) and \(\mathbf F \) that were introduced in Sect. 2, we can rewrite Eq. (35) as:

$$\begin{aligned} \text{ det}(\mathbf{C })=\text{ det}(\breve{\mathbf{C }})\, \text{ det}\left((\mathbf{F }^\mathrm{T}\mathbf{C }^{-1}\mathbf{F })^{-1}\right) \end{aligned}$$
(36)

The last term of Eq. (36) can be obtained by realising that the two multiplications with matrix \(\mathbf F \) do nothing more than selecting a submatrix of matrix \(\mathbf{C }^{-1}\). The relation with \(\mathbf{C }^{\prime }_{mm}\) can be found by looking at the formula for the matrix inverse for a block matrix (Brualdi and Schneider 1983):

$$\begin{aligned}&\mathbf{C }&^{-1}=\begin{pmatrix} \mathbf{C }_{oo}^{-1} \!+\! \mathbf{C }_{oo}^{-1} \mathbf{C }_{om}(\mathbf{C }^{\prime }_{mm})^{-1}\mathbf{C }_{mo}\mathbf{C }_{oo}^{-1}&\!-\!\mathbf{C }_{oo}^{-1}\mathbf{C }_{om}(\mathbf{C }^{\prime }_{mm})^{-1}\nonumber \\ \!-\!(\mathbf{C }^{\prime }_{mm})^{-1}\mathbf{C }_{mo}\mathbf{C }_{oo}^{-1}&(\mathbf{C }^{\prime }_{mm})^{-1} \end{pmatrix}\\ \end{aligned}$$
(37)

one will note the equivalence of this submatrix with the inverse of the Schur complement. Thus, we have:

$$\begin{aligned} \mathbf{C }^{\prime }_{mm} = \left(\mathbf{F }^\mathrm{T}\mathbf{C }^{-1}\mathbf{F }\right)^{-1} \end{aligned}$$
(38)

By taking the logarithm of Eq. (36), and using the relation \(\text{ det}(\mathbf{C })=1/\text{ det}(\mathbf{C }^{-1})\), Eq. (8) of Sect. 2 is obtained.

Of course the missing data will normally not occur after all observations have been made. Determinants have the property that rows and columns of the matrix can be interchanged without changing the value except for a change in sign. However, when rows \(i\) and \(j\) are swapped of the covariance matrix, the columns \(i\) and \(j\) are swapped at the same time which means that the sign does not change. Therefore, Eq. (8) is valid for any sequence of missing data.

Next, using Eq. (38) we have:

$$\begin{aligned} \mathbf{F }(\mathbf{F }^\mathrm{T}\mathbf{C }^{-1}\mathbf{F })^{-1}\mathbf{F }^\mathrm{T} = \begin{pmatrix} \mathbf{0 }&\mathbf{0 } \\ \mathbf{0 }&\mathbf{C }^{\prime }_{mm} \end{pmatrix} \end{aligned}$$
(39)

Using Eq. (37), this leads to the following result:

$$\begin{aligned}&\mathbf{C }^{-1}\mathbf{F }(\mathbf{F }^\mathrm{T}\mathbf{C }^{-1}\mathbf{F })^{-1}\mathbf{F }^\mathrm{T}\mathbf{C }^{-1} \nonumber \\&\quad = \begin{pmatrix} \mathbf{C }_{oo}^{-1}\mathbf{C }_{om}(\mathbf{C }^{\prime }_{mm})^{-1} \mathbf{C }_{mo}\mathbf{C }_{oo}^{-1}&-\mathbf{C }_{oo}^{-1}\mathbf{C }_{om}(\mathbf{C }^{\prime }_{mm})^{-1}\\ -(\mathbf{C }^{\prime }_{mm})^{-1}\mathbf{C }_{mo}\mathbf{C }_{oo}^{-1}&(\mathbf{C }^{\prime }_{mm})^{-1} \end{pmatrix}\nonumber \\ \end{aligned}$$
(40)

Changing the sign of Eq. (40) and adding \(\mathbf{C }^{-1}\) finally gives us:

$$\begin{aligned} \mathbf{C }^{-1} - \mathbf{C }^{-1}\mathbf{F }(\mathbf{F }^\mathrm{T}\mathbf{C }^{-1}\mathbf{F })^{-1} \mathbf{F }^\mathrm{T}\mathbf{C }^{-1} = \begin{pmatrix} \mathbf{C }_{oo}^{-1}&\mathbf{0 } \\ \mathbf{0 }&\mathbf{0 } \end{pmatrix} \end{aligned}$$
(41)

This provides us the proof for Eq. (7) in Sect. 2:

$$\begin{aligned} {\breve{\mathbf{r }}}^\mathrm{T} {\breve{\mathbf{C }}}^{-1}{\breve{\mathbf{r }}} = \mathbf{r }_o^\mathrm{T} \left( {\mathbf{C }^{-1} - \mathbf{C }^{-1} \mathbf{F } (\mathbf{F }^\mathrm{T}\mathbf{C }^{-1} \mathbf{F })^{-1} \mathbf{F }^\mathrm{T} \mathbf{C }^{-1} }\right) \mathbf{r }_o \end{aligned}$$
(42)

Swapping rows \(i\) and \(j\) of the covariance matrix \(\mathbf C \) will cause a swap of columns \(i\) and \(j\) of matrix \(\mathbf{C }^{-1}\). Swapping the columns of \(\mathbf C \) will cause a swap of the rows of \(\mathbf{C }^{-1}\) in a similar way. Thus, we can reshuffle our set of observations and missing data, while at the same time adjusting matrix \(\mathbf F \), to obtain the form of Eq. (34) which proves the general validity of Eq. (42) for any sequence of missing data.

Rights and permissions

Reprints and permissions

About this article

Cite this article

Bos, M.S., Fernandes, R.M.S., Williams, S.D.P. et al. Fast error analysis of continuous GNSS observations with missing data. J Geod 87, 351–360 (2013). https://doi.org/10.1007/s00190-012-0605-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s00190-012-0605-0

Keywords