Coordinate time sequence period and noise analysis of Shandong CORS
-
摘要: 针对不同区域的周期和噪声特性存在差异,以27个山东连续运行参考站(CORS)2015—2018年的坐标时间序列为研究对象,进行周期和噪声特征分析. 结果表明:坐标时间序列并没有严格的周期信号,U方向的周期运动较为规律;坐标分量的谱指数显示大部分站点可用白噪声(WN)+闪烁噪声(FN)的噪声模型进行描述;采用最大似然估计法(MLE)确定山东CORS站坐标时间序列N方向的最佳噪声模型为WN+FN,E、U方向的最佳噪声模型为WN+FN+随机漫步噪声(RWN).
-
关键词:
- 坐标时间序列 /
- CORS /
- 噪声特征 /
- 最大似然估计(MLE) /
- 谱指数
Abstract: In view of the difference of period and noise characteristics in different regions, this paper takes the coordinate time series of 27 CORS in Shandong Province from 2015 to 2018 as the research object, and analyzed the period and noise of coordinate components of each station on the basis of data preprocessing. The results show that, there is no strict whole periodic signal in coordinate time series, and the periodic motion in U direction is relatively regular. The spectrum index of coordinate component shows that the noise model of most stations is white noise (WN) + flicker noise (FN). Using the maximum likelihood estimation criterion, the best noise model of coordinate time series in Shandong CORS station is WN + FN, and the best noise model in E and U direction is WN + FN. The best noise model is WN + FN + random walk noise (RWN). -
0. 引 言
GPS坐标时间序列不仅反映坐标分量的整体趋势,还反映出基准站的非线性变化. 相关研究表明,由于全球卫星导航系统(GNSS)基准站坐标时间序列受到地球物理因素及系统误差的影响,不仅存在线性变化,还存在周期变化,尤其在高程方向表现得最为明显[1-6]. 分析坐标时间序列的非线性变化特征,有利于研究基准站的稳定性和基于非线性运动模型获得基准站的瞬时位置. 坐标时间序列的噪声信息会影响坐标解算精度,造成测站非线性运动模型有一定的偏差. GNSS坐标时间序列不仅存在白噪声(WN)还存在有色噪声(CN),目前最优随机模型广泛地被描述为WN+闪烁噪声(FN)[7-14]. 就山东部分连续运行参考站(CORS)而言,唐江森等[15]研究了两年坐标时间序列的噪声特征,表明其最佳噪声模型也为WN+FN. 通常认为时间序列跨度大于2.5年,才能准确估计时间序列的线性项、周期项及其精度[16].
本文在已有研究的基础上,以山东CORS 2015-01-01—2018-12-31连续四年的坐标时间序列为研究对象,深入分析CORS站坐标时间序列的周期与噪声特征,以确定山东CORS站点的最佳噪声模型.
1. 数据处理策略
选取了较为均匀分布在山东地区的CORS站,如图1所示.
受到外界条件、多路径效应或传输信号干扰等因素的影响,观测时间序列存在粗差,不仅影响数据质量,而且影响站点速度的计算及噪声特性的分析,有必要将其剔除. 采用拉依达准则法(
$3\sigma $ 准则)探测并剔除观测序列的粗差. 外界观测条件、站点受损或网络故障等都会影响时间序列的连续性和完整性,加之剔除粗差的历元,坐标时间序列会出现在某天或某几天的观测数据缺失的现象,经过对多种插值方法反复实验对比,采用三次多项式插值补全数据.坐标时间序列在较长时间内受某些因素的影响,形成一种整体变动的趋势,该趋势使得坐标时间序列为非平稳序列. 对坐标时间序列进行零均值化处理,实际上就是平稳化的过程,其拟合模型如下:
$${y _t} = {a} + {bt} + {{x} _t}.$$ (1) 式中:
${y _t}$ 为坐标时间序列;${a} ,\;b$ 分别为坐标时间序列拟合的常数值和一次项系数值;${{x} _t}$ 为去除趋势项的残差序列.2. CORS站坐标时间序列分析
2.1 周期特征分析
小波分析方法是一种同时在时、频两域内分析时间序列的方法,实现坐标时间序列周期项的探测与提取. 其基本原理如下:
$$\left\{ {\begin{array}{*{20}{c}} {{{{\mathit{\boldsymbol{A}}}}} _{{K} + {\rm{1}}}} = {{\mathit{\boldsymbol{H}}}}{{{{\mathit{\boldsymbol{A}}}}} _{K} } \;\;\;{{K} = {\rm{1}},{\rm{2}}, \cdots ,{N} } \\ {{{{{\mathit{\boldsymbol{D}}}}} _{{K} + {\rm{1}}}} = {{\mathit{\boldsymbol{G}}}}{{{{\mathit{\boldsymbol{D}}}}} _{K} } \;\;\;{{K} = {\rm{1}},{\rm{2}}, \cdots ,{N} } } \end{array}} \right..$$ (2) 式(2)为小波分解公式. 式中:
${{{{\mathit{\boldsymbol{A}}}}} _{K} }$ ,${{{{\mathit{\boldsymbol{D}}}}} _{K} }$ 分别表示分解后第${K} $ 层的低频和高频部分;H、G分别表示小波低通和高通滤波器;${N} $ 表示信号长度.$${{{{\mathit{\boldsymbol{A}}}}} _{K} } = {{{\mathit{\boldsymbol{H}}}}^{*}}{{{{\mathit{\boldsymbol{A}}}}} _{{K} + {\rm{1}}}} + {{{{\mathit{\boldsymbol{G}}}}^{*}}{{{{\mathit{\boldsymbol{D}}}}} _{{K} + {\rm{1}}}}} \;\;{{K} = {N} , \cdots ,{\rm{2}},{\rm{1}}}. $$ (3) 式(3)为小波重构公式. 式中:H*、G*为H、G的共轭转置.
通过小波分析提取出QUFU站的年、半年和季节周期项,如图2所示.
由图2可知,QUFU站N、E、U三个坐标方向均表现出一定的周期特征,具有明显的年周期项和半年周期项变化,整体周期性变化较为稳定. 对27个CORS站的周期特征进行分析可得:1)U方向的半年周期变化特征最为明显,其次是年周期特征,最后是3~4个月的季节周期特征;2)N、E方向具有明显的年周期或半年周期,不同站点,差异较大;3)坐标时间序列的周期信号并不是严格的年周期、半年周期等周期项,而是接近于年周期、半年周期等周期项. 不同测站所表现出来的周期特征存在着差异,同一测站的不同方向坐标分量的周期规律也有所不同,这是因为测站所处的地理位置不同,环境因素成为影响它们周期规律不同的主要因素.
2.2 噪声特征分析
2.2.1 噪声序列的获取
GPS基准站的拟合模型[17]可用下式表示
$$ \begin{aligned} y \left( {{{t} _i}} \right) = & {a} + {{bt} _i} + {c} {\rm{sin}}\left( {2{\text{π}}{t_i}} \right) + {d} {\rm{cos}}\left( {2{\text{π}}{t_i}} \right) +\\ & {e} {\rm{sin}} \left( {4{\text{π}}{t_i}} \right) + {f} {\rm{cos}}\left( {4{\text{π}}{t_i}} \right) +\\ & \displaystyle \sum\limits_{{j} = 1}^{{{n} _j}} {{{g} _i}{{H}}\left( {{{t} _i} - {t_{g_j} }} \right)} {\rm{ + }}{{v} _i} \end{aligned} .$$ (4) 式中:
$y \left( {{{t} _i}} \right)$ 为CORS站原始坐标时间序列;${a} $ 为观测序列的初始位置;${b} $ 为线性速度;${t_i}$ 为时间;c、d、e、f分别为年、半年周期项系数;${{g} _i}$ 为${{t} _{g_{j}}}$ 时刻的阶跃;H为阶梯函数;${{v} _i}$ 为噪声序列. 若顾及更复杂的情况,如某时刻测站运动速率的改变或某一事件发生后测站运动速率呈指数衰减的情况等,可采用文献[8]提供的坐标分量每日解观测序列参数模型. 在数据处理策略中可实现坐标时间序列线性趋势项、阶跃的去除,再利用小波分析提取坐标时间序列的周期项,由此获得坐标分量的噪声序列.2.2.2 谱指数与最大似然估计
噪声的功率谱密度
${P} \left( {f} \right)$ 与噪声频率${f} $ 之间存在着某种幂次关系,即$${P} \left( {f} \right) \propto {{f} ^\alpha },$$ (5) 式中,
$\alpha $ 为谱指数. 式(5)两边取对数为$${\rm{ln}}{P} \left( {f} \right) \propto \alpha {\rm{ln}} {f} .$$ (6) 明显地,
$\alpha $ 在双对数直角坐标系中为拟合直线的斜率. 不同的谱指数对应不同的噪声类型,$\alpha = {\rm{0}}$ 对应WN,$\alpha = - 1$ 对应FN,$\alpha = - 2$ 对应RWN.最大似然估计(MLE)是确定不同噪声模型下噪声分量的方法,使得在相应噪声模型下噪声序列与其协方差的概率密度最大[18],即
$$\begin{array}{l} {L} \left( {{{\mathit{\boldsymbol{X}}}},{{{\mathit{\boldsymbol{\sigma }}}}_{\text{w}}},{{{\mathit{\boldsymbol{\sigma }}}}_{{\text{rw}}}},{{{\mathit{\boldsymbol{\sigma }}}}_{\text{f}}}} \right) = {\left( {2{\text{π}}} \right)^{ - {\rm{0}}{\rm{.5}}{N} }}{\rm{*}}{\left| {{{{\mathit{\boldsymbol{Q}}}}_{{{XX}}}}} \right|^{ - {\rm{0}}{\rm{.5}}}}{\rm{*}} \\ \qquad\qquad\qquad \qquad \exp \left[ { - 0.5{{{\mathit{\boldsymbol{X}}}}^{\rm{T}}}{{\mathit{\boldsymbol{Q}}}}_{{{XX}}}^{{\rm{ - 1}}}{{\mathit{\boldsymbol{X}}}}} \right] \\ \end{array} .$$ (7) 式中:
${{\mathit{\boldsymbol{X}}}}$ 为时间序列列向量;${{{\mathit{\boldsymbol{\sigma }}}}_{\text{w}}},{{{\mathit{\boldsymbol{\sigma }}}}_{{\text{rw}}}},{{{\mathit{\boldsymbol{\sigma }}}}_{\text{f}}}$ 分别为WN、RWN和FN;${N} $ 为序列长度;${{{\mathit{\boldsymbol{Q}}}}_{{XX} }}$ 为${{\mathit{\boldsymbol{X}}}}$ 的协方差.根据式(6)谱指数的定义,计算求得CORS站各坐标分量的谱指数,如表1所示.
由表1可知,CORS站水平方向的谱指数均在−1~0之间,垂直方向除ZAZH站外,谱指数均在−1~0之间,说明CORS站各坐标分量噪声均不具有纯WN的特性,且同一测站不同方向的噪声模型有所不同. 为确定CORS站坐标分量的最佳噪声模型,根据最大似然估计准则,以ZAZH站为例,计算了五种噪声模型:WN、WN+FN、WN+RWN、FN+RWN和WN+FN+RWN的MLE值,结果如表2所示.
表 1 CORS站各方向谱指数测站 N E U 测站 N E U ANQI −0.637 −0.713 −0.835 QUFU −0.698 −0.704 −0.823 CAOX −0.681 −0.670 −0.836 RZRS −0.943 −0.965 −0.890 DEZH −0.620 −0.698 −0.936 SDTC −0.909 −0.854 −0.946 DOMI −0.752 −0.651 −0.792 SDZB −0.627 −0.696 −0.827 FEIX −0.592 −0.564 −0.869 SHID −0.739 −0.671 −0.845 HEKO −0.515 −0.592 −0.673 SHXI −0.737 −0.715 −0.935 JUCH −0.833 −0.644 −0.952 TERZ −0.670 −0.723 −0.869 JUNA −0.871 −0.873 −0.898 WEIH −0.594 −0.692 −0.770 JUXI −0.657 −0.631 −0.856 WEND −0.669 −0.797 −0.662 JUYE −0.781 −0.679 −0.832 YAXI −0.602 −0.683 −0.826 LAIW −0.624 −0.650 −0.847 YINA −0.715 −0.718 −0.789 LINQ −0.612 −0.700 −0.826 YISH −0.698 −0.678 −0.969 LISH −0.792 −0.720 −0.916 ZAZH −0.746 −0.798 −1.011 NIJI −0.655 −0.606 −0.769 由表2可知,不同噪声模型对应的MLE值不同,WN模型MLE值最小,说明坐标分量中不仅含有WN,还含有CN. 水平方向上,WN+FN模型和WN+FN+RWN模型的MLE值最大且相等,说明在没有RWN模型存在的情况下,WN+FN模型与WN+FN+RWN模型所估计出来的噪声分量是相同的. 因此,ZAZH站水平方向的最佳噪声模型为WN+FN或WN+FN+RWN;垂直方向上,MLE值为−4202.15所对应的最佳噪声模型为WN+FN+RWN. 由此可以看出同一CORS站的不同方向可以有不同的最佳噪声模型. 假设CORS站的噪声模型为WN+FN+RWN,在此模型下计算了山东地区27个CORS站的噪声分量,如表3所示.
表 2 ZAZH噪声模型MLE值方向 WN WN+FN WN+RWN FN+RWN WN+FN+RWN N −3 175.58 −2 901.89 −2 957.20 −2 904.16 −2 901.89 E −3 642.48 −3 249.34 −3 299.46 −3 251.23 −3 249.34 U −5 365.97 −4 223.56 −4 218.40 −4 208.10 −4 202.15 表 3 WN+FN+RWN噪声模型下的噪声分量估计值测站 N E U WN FN RWN WN FN RWN WN FN RWN ANQI 0.853 6.939 0.000 1.279 8.070 0.000 1.664 19.089 0.000 CAOX 1.099 7.309 0.000 1.771 7.294 4.108 2.318 22.512 9.085 DEZH 1.048 6.274 0.000 1.353 7.770 0.000 0.996 22.032 6.664 DOMI 0.850 7.513 0.000 1.340 7.959 0.000 1.976 21.472 0.000 FEIX 1.204 6.736 0.000 1.718 6.800 0.000 2.124 20.632 6.112 HEKO 1.171 6.675 0.000 1.556 7.339 0.000 2.795 16.744 0.000 JUCH 0.947 7.743 2.991 1.548 7.269 3.405 2.174 21.862 14.566 JUNA 0.400 8.019 0.000 0.301 9.533 0.000 0.000 18.214 0.000 JUXI 1.175 6.952 0.000 1.654 6.895 3.277 1.868 20.698 7.005 JUYE 1.094 7.671 2.345 1.586 7.216 3.062 1.029 23.330 3.978 LAIW 1.180 6.449 0.000 1.524 7.461 2.166 1.916 22.091 7.233 LINQ 1.022 5.783 0.000 1.276 7.324 0.000 1.055 22.669 0.000 LISH 0.541 7.523 0.000 0.991 8.469 0.000 0.000 19.475 0.000 NIJI 1.201 6.882 0.000 1.609 8.112 0.000 1.489 21.411 0.000 QUFU 0.839 7.271 0.000 11.572 0.126 0.000 0.000 20.898 0.000 RZRS 0.000 7.337 0.000 0.000 9.628 0.000 0.000 15.363 0.000 SDTC 0.000 9.738 0.000 1.048 10.140 3.640 0.000 20.185 0.000 SDZB 0.906 6.047 0.000 0.983 8.176 0.000 0.000 20.053 0.000 SHID 0.618 7.696 0.000 1.170 8.300 0.000 2.345 17.315 8.187 SHXI 0.929 6.903 0.000 1.421 7.667 4.859 0.000 22.410 10.722 TERZ 0.773 6.843 0.000 1.114 8.574 0.000 0.000 21.085 0.000 WEIH 0.870 6.622 0.000 1.084 8.658 0.000 2.319 18.461 0.000 WEND 0.660 7.751 0.000 0.817 9.334 0.000 2.506 15.586 0.000 YAXI 1.180 6.625 0.000 1.412 8.755 0.000 1.681 20.653 0.000 YINA 0.659 6.926 0.000 1.208 8.563 0.000 1.548 17.635 0.000 YISH 1.005 6.568 0.000 1.305 7.965 0.000 0.000 18.823 7.103 ZAZH 0.660 7.029 0.000 0.785 9.028 0.000 2.013 14.735 20.745 由表3可知,观测时间序列坐标分量具有不同的噪声特征. N方向上,有25个站可用模型WN+FN描述,2个站可用模型WN+FN+RWN描述;E方向上,20个站可用模型WN+FN描述,7个站可用模型WN+FN+RWN描述;U方向上,16个站可用模型WN+FN描述,11个站可用模型WN+FN+RWN描述. 因此山东CORS站的最佳噪声模型为:N方向上采用WN+FN模型,E和U方向上采用WN+FN+RWN模型. 同时,U方向的噪声分量估计值要明显高于其它两个方向,这与高程方向精度低于水平方向精度的结论是一致的.
3. 结束语
本文以山东27个CORS站的坐标时间序列为研究对象,通过小波分析方法提取周期项进行分析,得出了坐标时间序列中并不严格地存在年周期项、半年周期项等周期信号,而是与之接近的周年、半周年项等周期信号;另外,在观测序列的坐标分量中,U方向的周期信号较水平方向明显.
CORS站的谱指数表明坐标时间序列中不仅存在WN,还存在CN. 根据MLE方法,估计了CORS站各坐标分量在WN+FN+RWN模型下的噪声分量,得出了CORS站各坐标分量可以有不同的噪声模型,N方向的最佳噪声模型为WN+FN,E、U方向的最佳噪声模型为WN+FN+RWN的结论.
-
表 1 CORS站各方向谱指数
测站 N E U 测站 N E U ANQI −0.637 −0.713 −0.835 QUFU −0.698 −0.704 −0.823 CAOX −0.681 −0.670 −0.836 RZRS −0.943 −0.965 −0.890 DEZH −0.620 −0.698 −0.936 SDTC −0.909 −0.854 −0.946 DOMI −0.752 −0.651 −0.792 SDZB −0.627 −0.696 −0.827 FEIX −0.592 −0.564 −0.869 SHID −0.739 −0.671 −0.845 HEKO −0.515 −0.592 −0.673 SHXI −0.737 −0.715 −0.935 JUCH −0.833 −0.644 −0.952 TERZ −0.670 −0.723 −0.869 JUNA −0.871 −0.873 −0.898 WEIH −0.594 −0.692 −0.770 JUXI −0.657 −0.631 −0.856 WEND −0.669 −0.797 −0.662 JUYE −0.781 −0.679 −0.832 YAXI −0.602 −0.683 −0.826 LAIW −0.624 −0.650 −0.847 YINA −0.715 −0.718 −0.789 LINQ −0.612 −0.700 −0.826 YISH −0.698 −0.678 −0.969 LISH −0.792 −0.720 −0.916 ZAZH −0.746 −0.798 −1.011 NIJI −0.655 −0.606 −0.769 表 2 ZAZH噪声模型MLE值
方向 WN WN+FN WN+RWN FN+RWN WN+FN+RWN N −3 175.58 −2 901.89 −2 957.20 −2 904.16 −2 901.89 E −3 642.48 −3 249.34 −3 299.46 −3 251.23 −3 249.34 U −5 365.97 −4 223.56 −4 218.40 −4 208.10 −4 202.15 表 3 WN+FN+RWN噪声模型下的噪声分量估计值
测站 N E U WN FN RWN WN FN RWN WN FN RWN ANQI 0.853 6.939 0.000 1.279 8.070 0.000 1.664 19.089 0.000 CAOX 1.099 7.309 0.000 1.771 7.294 4.108 2.318 22.512 9.085 DEZH 1.048 6.274 0.000 1.353 7.770 0.000 0.996 22.032 6.664 DOMI 0.850 7.513 0.000 1.340 7.959 0.000 1.976 21.472 0.000 FEIX 1.204 6.736 0.000 1.718 6.800 0.000 2.124 20.632 6.112 HEKO 1.171 6.675 0.000 1.556 7.339 0.000 2.795 16.744 0.000 JUCH 0.947 7.743 2.991 1.548 7.269 3.405 2.174 21.862 14.566 JUNA 0.400 8.019 0.000 0.301 9.533 0.000 0.000 18.214 0.000 JUXI 1.175 6.952 0.000 1.654 6.895 3.277 1.868 20.698 7.005 JUYE 1.094 7.671 2.345 1.586 7.216 3.062 1.029 23.330 3.978 LAIW 1.180 6.449 0.000 1.524 7.461 2.166 1.916 22.091 7.233 LINQ 1.022 5.783 0.000 1.276 7.324 0.000 1.055 22.669 0.000 LISH 0.541 7.523 0.000 0.991 8.469 0.000 0.000 19.475 0.000 NIJI 1.201 6.882 0.000 1.609 8.112 0.000 1.489 21.411 0.000 QUFU 0.839 7.271 0.000 11.572 0.126 0.000 0.000 20.898 0.000 RZRS 0.000 7.337 0.000 0.000 9.628 0.000 0.000 15.363 0.000 SDTC 0.000 9.738 0.000 1.048 10.140 3.640 0.000 20.185 0.000 SDZB 0.906 6.047 0.000 0.983 8.176 0.000 0.000 20.053 0.000 SHID 0.618 7.696 0.000 1.170 8.300 0.000 2.345 17.315 8.187 SHXI 0.929 6.903 0.000 1.421 7.667 4.859 0.000 22.410 10.722 TERZ 0.773 6.843 0.000 1.114 8.574 0.000 0.000 21.085 0.000 WEIH 0.870 6.622 0.000 1.084 8.658 0.000 2.319 18.461 0.000 WEND 0.660 7.751 0.000 0.817 9.334 0.000 2.506 15.586 0.000 YAXI 1.180 6.625 0.000 1.412 8.755 0.000 1.681 20.653 0.000 YINA 0.659 6.926 0.000 1.208 8.563 0.000 1.548 17.635 0.000 YISH 1.005 6.568 0.000 1.305 7.965 0.000 0.000 18.823 7.103 ZAZH 0.660 7.029 0.000 0.785 9.028 0.000 2.013 14.735 20.745 -
[1] 于亚杰, 唐江森. CORS基准站坐标时间序列周期性分析[J]. 城市勘测, 2019(4): 132-135. DOI: 10.3969/j.issn.1672-8262.2019.04.031 [2] 魏冠军, 陈晨, 黄逸宇, 等. 基于小波变换的CORS高程坐标时序周期性分析与降噪方法研究[J]. 全球定位系统, 2019, 44(4): 61-67. [3] 戴海亮, 孙付平, 朱新慧. 中国区域内IGS站时间序列的非线性变化分析[J]. 全球定位系统, 2018, 43(6): 76-81. [4] 姜卫平, 李昭, 刘鸿飞, 等. 中国区域IGS基准站坐标时间序列非线性变化的成因分析[J]. 地球物理学报, 2013, 56(7): 2228-2237. DOI: 10.6038/cjg20130710 [5] 蔡富, 孙付平, 戴海亮, 等. 小波和傅里叶变换在坐标时间序列分析中的应用[J]. 全球定位系统, 2019, 44(4): 40-46. [6] RAY J, ALTAMIMI Z, COLLILIEUX X, et al. Anomalous harmonics in the spectra of GPS position estimates[J]. GPS solutions, 2008, 12(1): 55-64. DOI: 10.1007/s10291-007-0067-7.
[7] 明锋. GPS坐标时间序列分析研究[J]. 测绘学报, 2019, 48(10): 1340-1340. DOI: 10.11947/j.AGCS.2019.20190006 [8] LANGBEIN J. Noise in GPS displacement measurements from southern California and southern Nevada[J]. Journal of geophysical research atmosphere, 2008, 113(B5): B05405. DOI: 10.1029/2007JB005247.
[9] 袁林果, 丁晓利, 陈武, 等. 香港GPS基准站坐标序列特征分析[J]. 地球物理学报, 2008, 51(5): 1372-1384. DOI: 10.3321/j.issn:0001-5733.2008.05.011 [10] 蒋志浩, 张鹏, 秘金钟, 等. 顾及有色噪声影响的CGCS2000下我国CORS站速度估计[J]. 测绘学报, 2010, 39(4): 355-363. [11] 李昭, 姜卫平, 刘鸿飞, 等. 中国区域IGS基准站坐标时间序列噪声模型建立与分析[J]. 测绘学报, 2012, 41(4): 496-503. [12] 曹奇, 张迎燕, 岳东杰, 等. 城市CORS坐标时间序列噪声模型建立与分析[J]. 大地测量与地球动力学, 2015, 35(6): 1001-1006. [13] 金波文, 王慧, 刘玉龙, 等. 中国沿海海洋站GNSS坐标时间序列噪声模型的建立与分析[J]. 大地测量与地球动力学, 2020, 40(5): 476-481. [14] 许家琪. GNSS坐标时间序列噪声模型建立影响因素研究[D]. 上海: 东华理工大学, 2019. [15] 唐江森, 曲国庆, 苏晓庆. 山东CORS基准站坐标时间序列噪声分析[J]. 测绘科学, 2016, 41(10): 63-68. [16] 王健, 许安安, 周伯烨. 顾及共模误差的大区域GPS网坐标时间序列噪声分析[J]. 测绘通报, 2018(4): 6-9, 56. [17] NIKOLAIDIS R. Observation of geodetic and seismic deformation with the global position system[D]. California: University of California, 2002.
[18] AlVARO S-G, BOUIN M-N, COLLIlLIEUX X, et al. Correlated errors in GPS position time series: implications for velocity estimates[J]. Journal of geophysical research, 2011, 116(B1): B01405. DOI: 10.1029/2010JB007701.
-
期刊类型引用(1)
1. 朱文强,李建章,隋哲民,高志钰. 华中华东陆态网坐标时序噪声模型及环境负载分析. 导航定位学报. 2023(02): 147-152 . 百度学术
其他类型引用(1)