含缺值GNSS基准站坐标序列的非插值小波分析与信号提取

受多种因素影响,GNSS基准站坐标序列通常都含有缺值,传统小波分析需要对缺值数据进行内插或补零处理。本文基于小波系数与时间序列观测数据的重构关系,提出了一种非插值的二进小波变换的最小范数解法,导出了相应的计算式,并严格证明了传统的补零处理算法与本文的最小范数解法等价。最后利用中国地壳运动观测网络一期27个基准站实测数据以及模拟数据进行了验证分析。结果表明,本文的非插值算法与插值算法提取的信号差异较小,27个基准站坐标序列的平均残差中误差仅相差2.01%(North),0.54%(East)和1.26%(Up),两种算法提取的信号之差与信号平均方差比仅相差1.16%(North),0.54%(East)和1.62%(Up)。     

TSAnalyzer,a GNSS time series analysis software

GPS Solutions - In geodesy and geophysics, continuous GNSS observations have been used globally. As the number of GNSS observing stations increases, GNSS time series analysis software should be...     

CATS: GPS coordinate time series analysis software

Over the last 10 years, several papers have established that daily estimates of GPS coordinates are temporally correlated and it is therefore incorrect to assume that the observations are independent when estimating parameters from them. A direct consequence of this assumption is the over-optimistic estimation of the parameter uncertainties. Perhaps the perceived computational burden or the lack of suitable software for time series analysis has resulted in many heuristic methods being proposed in the scientific literature for estimating these uncertainties. We present a standalone C program, CATS, developed to study and compare stochastic noise processes in continuous GPS coordinate time series and, as a consequence, assign realistic uncertainties to parameters derived from them. The name originally stood for Create and Analyze Time Series. Although the name has survived, the creation aspect of the software has, after several versions, been abandoned. The implementation of the method is briefly described to aid understanding and an example of typical input, usage, output and the available stochastic noise models are given.     

京津地区GNSS坐标时间序列噪声特性研究

针对地表质量负荷对京津地区GNSS坐标时间序列噪声特性的影响,选取中国大陆构造环境监测网络8个GNSS基准站2012—2014年的坐标时间序列,利用CATS软件计算大气压、非潮汐海洋、积雪和土壤湿度等质量负载改正前后GNSS坐标时间序列的谱指数、最优噪声模型、速度的变化。发现地表质量负载对GNSS坐标时间序列的噪声特性产生了明显影响。结果显示,京津地区GNSS坐标序列包含白噪声和有色噪声,且最优噪声模型具有多样性。扣除质量负载后N、U分量的噪声模型变化明显,主要表现为FN+WN和PL+WN,而N、E分量的谱指数分别趋近于FN和WN。质量负载改正后基准站U方向的线性速度变化较大,且北京地区变化量大于天津地区。研究结果为提高GNSS数据解算精度、精细分析地壳形变提供参考。     

双集合卡尔曼滤波估算时间序列LAI

遥感估算叶面积指数(LAI)时空动态变化对全球气候变化研究具有重要的意义,为了提高遥感估算时间序列叶面积指数的精度,需要耦合遥感观测数据与LAI动态过程模型。本文提出一种基于双集合卡尔曼滤波(Dual EnKF)的时间序列LAI反演方法,同时更新LAI估计值和LAI动态过程模型中的敏感性参数,得到LAI和动态过程模型敏感参数的最优估计值来优化动态过程模型。一方面使得动态过程模型可以更好地描述LAI随时间的变化过程,降低模型预测误差,从而提高LAI动态过程模型的预测能力;另一方面通过耦合动态过程模型和辐射传输模型,集成遥感观测数据与动态过程模型的预测值,进而得到优化的LAI估计值。为检验算法,分别选取作物、草地和林地等典型植被验证站点进行Dual EnKF LAI时间序列估算,并分别与MODIS LAI产品及其SG滤波曲线、集合卡尔曼滤波方法反演LAI、未优化的动态过程模型模拟LAI结果进行比较,并配以一些站点地面实测点数据作为参考。结果表明,采用Dual EnKF方法得到的LAI不但保持了时间上的连续性,而且通过改善动态过程模型的预测能力,即使在缺乏高质量遥感观测数据时,也能够获得符合LAI发展趋势的估算值,没有出现跳跃、波动现象,时间序列曲线较稳定,更符合植被LAI变化规律,表明基于Dual EnKF的时间序列LAI遥感估算方法是提取LAI时间廓线的一种有效途径。     

联合LMD与EMD的GNSS站坐标时间序列去噪方法

针对GNSS站坐标时间序列信噪不易分离的问题,在传统EMD去噪方法的基础上,本文提出了一种联合LMD与EMD的坐标时间序列去噪方法。该方法首先采用LMD分解原始坐标时间序列,基于连续均方误差（CMSE）原则分离高频噪声与低频信号,保持低频分量不变;然后对高频分量进行EMD去噪;最后以2次分解所得低频信号之和作为去噪后时间序列。以仿真数据与8个GNSS基准站实测数据进行试验,通过多种评价指标进行精度评估。结果表明,与传统EMD方法相比,联合LMD与EMD的方法能够更加精确地去除坐标时间序列中的噪声。     

联合EMD-HD和小波分解的GNSS坐标时间序列降噪分析EI北大核心CSCD

针对经验模态分解(empirical mode decomposition,EMD)在GNSS坐标时间序列的降噪过程中存在筛选准则的选取和模态混叠效应等问题,本文引入Hausdorff距离(Hausdorff distance,HD)筛选准则并结合小波分解(wavelet decomposition,WD),提出EMD-HD&WD算法。通过对我国大陆构造环境监测网络149个GNSS测站的垂向坐标时间序列降噪处理,分别利用复合指标T值、测站的速度不确定度和闪烁噪声振幅验证算法的可靠性和普适性。结果显示:HD优于现有的筛选准则;EMD-HD&WD算法对测站的速度不确定度和闪烁噪声振幅的平均改正率均为88.4%。分析表明,本文算法能够有效识别和剔除噪声并且改善EMD的模态混叠效应,提高GNSS垂向坐标时间序列的模型精度。     

LMD和小波阈值的GNSS坐标时间序列降噪应用

局部均值分解方法降噪过于粗糙,将认定为噪声的乘积函数(PF)分量直接剔除,导致有用信息丢失.为了有效提取GNSS站坐标时间序列的有用信息,该文提出一种局部均值分解和小波阈值相结合的降噪方法.通过局部均值分解将坐标时间序列分解为一系列PF分量和余项,依据消除趋势波动分析方法计算各PF分量的Hurst指数,利用小波阈值提取H≤1的PF分量中的有用信息,将提取出的信息与剩余PF分量叠加重构获得最终降噪的坐标时间序列.通过对5个测站的坐标时间序列进行实验,结果表明局部均值分解和小波阈值相结合的方法能够有效提取噪声分量中的有用信息,信噪比提高了27.8％,从而验证了该方法的有效性.     

基于独立分量法的新疆GNSS时间序列共模误差分析CSCD

共模误差(CME)是区域连续全球卫星导航系统(GNSS)网中的主要误差来源之一.针对GNSS时间序列具有非高斯分布特征,基于二阶统计量的主成分分析(PCA)难以准确提取出CME分量问题,采用具有高阶统计量的独立分量分析(ICA)对CME进行提取.以2011—2018年新疆区域GNSS坐标时间序列为例,将PCA滤波效果进行对比验证,分析了CME对GNSS坐标时间序列的影响,并对CME序列进行周期分析.结果表明:前6个独立分量包含CME分量,这可能与卫星轨道、地表质量负荷和时钟误差有关,ICA滤波后东(N)、北(E)、天顶(U)三个方向的均方根(RMS)值分别降低31.83%、32.29%、35.49%,速度不确定度分别降低44.14%、38.49%、43.32%,各测站的周期项振幅较滤波前更一致,有效地剔除了CME,提高了坐标时间序列的精度.     

Bayesian methods for outliers detection in GNSS time series

This article is concerned with the problem of detecting outliers in GNSS time series based on Bayesian statistical theory. Firstly, a new model is proposed to simultaneously detect different types of outliers based on the conception of introducing different types of classification variables corresponding to the different types of outliers; the problem of outlier detection is converted into the computation of the corresponding posterior probabilities, and the algorithm for computing the posterior probabilities based on standard Gibbs sampler is designed. Secondly, we analyze the reasons of masking and swamping about detecting patches of additive outliers intensively; an unmasking Bayesian method for detecting additive outlier patches is proposed based on an adaptive Gibbs sampler. Thirdly, the correctness of the theories and methods proposed above is illustrated by simulated data and then by analyzing real GNSS observations, such as cycle slips detection in carrier phase data. Examples illustrate that the Bayesian methods for outliers detection in GNSS time series proposed by this paper are not only capable of detecting isolated outliers but also capable of detecting additive outlier patches. Furthermore, it can be successfully used to process cycle slips in phase data, which solves the problem of small cycle slips.     

利用改进的Kalman滤波算法进行GNSS单历元姿态测量

传统G PS单历元姿态解算中存在没有利用历元间的相互关系以及不能有效抵制粗差影响等问题,而Kalman滤波则可一定程度上解决这些问题。结合稳健估计理论对GPS单历元姿态解算提出了一种改进的Kalman稳健估计算法,通过实验及与传统的最小二乘解算方法相比较,该算法能够较好地消除历元间计算结果中出现的跳变值,使得解算结果更稳健。     

Analyzing GNSS data in precise point positioning software

This work demonstrates that precise point positioning (PPP) can be used not only for positioning, but for a variety of other tasks, such as signal analysis. The fact that the observation model used for accurate error modeling has to take into consideration the several effects present in GPS signals, and that observations are undifferenced, makes PPP a powerful data analysis tool sensitive to a variety of parameters. The PPP application developed at the University of New Brunswick, which is called GAPS (GPS Analysis and Positioning Software), has been designed and built in order to take advantage of available precise products, resulting in a data analysis tool for determining parameters in addition to position, receiver clock error, and neutral atmosphere delay. These other estimated parameters include ionospheric delays, code biases, satellite clock errors, and code multipath among others. In all cases, the procedures were developed in order to be suitable for real-time as well as post-processing applications. One of the main accomplishments in the development described here is the use of very precise satellite products, coupled with a very complete observation error modeling to make possible a variety of analyses based on GPS data. In this paper, several procedures are described, their innovative aspects are pointed out, and their results are analyzed and compared with other sources. The procedures and software are readily adaptable for using data from other global navigation satellite systems.     

EPN coordinate time series monitoring for reference frame maintenance

The EUREF permanent network (EPN) is a network of continuously operating global navigation satellite system (GNSS) stations, primarily installed for reference frame maintenance. In order to ensure the long term reliability of the EPN products, a thorough station performance monitoring has been initiated and carried out in addition to the routine GNSS data management, processing and analysis. This paper addresses the main factors influencing the quality of the coordinate time series in a permanent GNSS network. Relevant examples, based on the EPN experience are given, the analysis strategy is introduced, the estimated coordinate offsets are published and the importance of this analysis for site velocity estimation is demonstrated. The results are derived from the analysis of the EPN weekly combined solutions covering the period from 1996 to 2003. Our target is the identification, interpretation and elimination of offsets and outliers present in the EPN coordinate time series in order to estimate reliable coordinates and velocities and consequently maintain a high quality kinematic reference network.

Detecting discontinuities in GNSS coordinate time series with STARS: case study,the Bologna and Medicina GPS sites

Global navigation satellite systems (GNSS) data are a fundamental source of information for achieving a better understanding of geophysical and climate-related phenomena. However, discontinuities in the coordinate time series might be a severe limiting factor for the reliable estimate of long-term trends. A methodological approach has been adapted from Rodionov (Geophys Res Lett 31:L09204, 2004; Geophys Res Lett 31:L12707, 2006) and from Rodionov and Overland (J Marine Sci 62:328–332, 2005) to identify both the epoch of occurrence and the magnitude of jumps corrupting GNSS data sets without any a priori information on these quantities. The procedure is based on the Sequential t test Analysis of Regime Shifts (STARS) (Rodionov in Geophys Res Lett 31:L09204, 2004). The method has been tested against a synthetic data set characterized by typical features exhibited by real GNSS time series, such as linear trend, seasonal cycle, jumps, missing epochs and a combination of white and flicker noise. The results show that the offsets identified by the algorithm are split into 48 % of true-positive, 28 % of false-positive and 24 % of false-negative events. The procedure has then been applied to GPS coordinate time series of stations located in the southeastern Po Plain, in Italy. The series span more than 15 years and are affected by offsets of different nature. The methodology proves to be effective, as confirmed by the comparison between the corrected GPS time series and those obtained by other observation techniques.     

局部均值分解和奇异值分解在GNSS站坐标时间序列信号降噪中的应用

为了有效地提取GNSS站坐标时间序列的有用信息，降低噪声干扰，本文提出一种局部均值分解和奇异值分解相结合的信号降噪方法，并利用5个测站的实测坐标时间序列对新方法进行了验证。首先通过局部均值分解将坐标时间序列分解成一系列PF分量和余项，然后利用连续均方误差方法确定高频分量与低频分量的分界点，保持低频分量不变，运用奇异值分解方法对高频分量进行降噪重构，最后将重构的高频分量与低频分量叠加得到最终的降噪坐标时间序列，并对降噪效果进行对比分析。结果表明，与单纯的奇异值分解方法相比，局部均值分解和奇异值分解相结合方法能够自适应地选择合适的奇异值个数进行信号重构，提高了降噪效果。     

基于多通道奇异谱的GNSS坐标序列 粗差探测与数据插值

由于接收机故障、天线更换以及一些未知外界环境因素的影响,导致GPS时间序列中不可避免地存在数据缺失和粗差,数据缺失和粗差会产生各种问题,因此需要鲁棒探测粗差和缺失数据插值来获得连续完整统一的时间序列。传统方法可能需要针对不同类型的数据和不同长度的数据间隙研究不同的插值方法,这使得缺失数据的插值较为困难。针对上述问题,文中采用多通道奇异谱分析(Multichannel Singular Spectrum Analysis,MSSA)对时间序列进行粗差探测和缺失数据插值,重建非均匀采样时间序列的连续可靠模型,且不需要关于时间序列的任何先验信息。在该方法中,粗差探测与数据插值同时进行。模拟数据和实测GNSS时间序列数据测试结果都表明该方法的有效性。     

Choosing the coherent integration time for Kalman filter-based carrier-phase tracking of GNSS signals

Carrier phase–based positioning using Global Navigation Satellite System (GNSS) signals can provide centimeter-level accuracy; however, to do so requires robust, continuous tracking of the phase of the received signal. The phase-locked loop is typically the weakest link in GNSS signal processing, with frequent cycle slips and loss of lock occurring at lower signal-to-noise ratios. One way to improve the signal-to-noise ratio is to increase the coherent integration time; however doing so reduces the loop update rate, thereby degrading performance. This paper investigates this trade-off between sensitivity and loop update rate by investigation of the Kalman filter-based tracking loop. It is shown that it is possible to choose an optimal integration time for a given application. A relatively straightforward procedure is given to determine this optimal value. The results are confirmed through real-time kinematic processing of live satellite signals.     

三种GNSS高程时序降噪方法的效果对比分析

为了探究经验模态分解(EMD)、整体经验模态分解(EEMD)和小波降噪三种方法的降噪性能,以中国区6个国际GNSS服务(IGS)站高程分量的5?a、10?a和20?a时序数据为例,对它们的降噪结果进行比较分析.?首先利用线性拟合分离趋势项,并采用3σ准则剔除异常值,得到满足符合降噪要求的样本序列;然后分别用这三种方法分...     

最小二乘拟合GNSS位置时间序列分析

针对GNSS位置时间序列包含的线性趋势及其变化可能干扰后续分析并掩盖动力学因素信息的问题,该文使用最小二乘方法分析时间序列的线性趋势并减弱时间序列中阶跃的影响。分析了最小二乘方法用于GNSS位置时间序列分析的可行性及利用该方法分别获取趋势变化点前后的线性趋势,据此估计时间序列趋势变化的大小进而可以对时间序列进行修复。该方法用于GNSS位置时间序列的初步分析,可以方便有效地去除线性趋势变化对后续时间序列分析的影响,同时拟合结果本身也能反应出时间序列的变化特征。     

GPS坐标时间序列严密三维噪声模型建立方法

积累多年的GPS基准站坐标时间序列被广泛应用于大地测量学与地球动力学研究。实现坐标时间序列中噪声与信号的有效分离并建立顾及不同分量噪声之间相关性的严密三维噪声模型,有助于确定测站线性以及非线性运动参数估值更加真实的不确定度。这对于精化速度场、合理解释地壳形变特征、建立毫米级的地球参考框架等研究具有重要的科学和现实意义。