Extended Principal Component Analysis for Spatiotemporal Filtering of Incomplete Heterogeneous GNSS Position Time Series

When ordinary principal component analysis (PCA) is employed to analyze the position time series of a regional GNSS station network, the GNSS time series are assumed to be homogeneous, and the missing data in the time series must be restored beforehand. To directly process incomplete and heterogeneous GNSS position time series, we develop the extended PCA (EPCA) and weighted EPCA approaches to solving for the missing values based on the best low-rank approximation in the spatiotemporal domain. The proposed approaches are used to process the real GNSS position time series of 24 stations in North China spanning 2011 to 2019 and successfully extract the common mode errors (CMEs). The proposed approaches are compared with modified PCA (MPCA) and weighted MPCA, in which an additional optimization criterion needs to be introduced in the frequency domain. The results show that EPCA can extract more CMEs than MPCA for both the unweighted and weighted cases. Consequently, EPCA outperforms MPCA in reducing noise and improving the accuracy of site velocity estimates. Repeated simulation experiments show that the CMEs extracted by EPCA are closer to the simulated true values than those extracted by MPCA. When the formal errors of the time series are considered, both weighted EPCA and weighted MPCA outperform their unweighted counterparts, and the former outperforms the latter. In addition, EPCA is computationally more efficient than MPCA since fewer unknowns need to be estimated.