CN114114332B - 一种有效探测gnss基准站坐标时间序列不连续点的方法 - Google Patents
一种有效探测gnss基准站坐标时间序列不连续点的方法 Download PDFInfo
- Publication number
- CN114114332B CN114114332B CN202111294457.4A CN202111294457A CN114114332B CN 114114332 B CN114114332 B CN 114114332B CN 202111294457 A CN202111294457 A CN 202111294457A CN 114114332 B CN114114332 B CN 114114332B
- Authority
- CN
- China
- Prior art keywords
- reference station
- gnss reference
- sequence
- station coordinate
- time sequence
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000001514 detection method Methods 0.000 claims abstract description 13
- 238000005457 optimization Methods 0.000 claims abstract description 9
- 238000012360 testing method Methods 0.000 claims abstract description 8
- 238000000528 statistical test Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 7
- 239000006185 dispersion Substances 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 5
- 238000012300 Sequence Analysis Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000005295 random walk Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000012731 temporal analysis Methods 0.000 description 1
- 238000000700 time series analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/23—Testing, monitoring, correcting or calibrating of receiver elements
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Complex Calculations (AREA)
Abstract
本发明提出了一种有效探测GNSS基准站坐标时间序列不连续点的方法。本发明以GNSS基准站坐标时间序列作为输入数据;对坐标时间序列中相邻时刻的观测值作差得到坐标差值序列,使用四分位距方法对差值序列进行探测得到初始不连续点;构造坐标时间序列的参数估计模型,通过优化求解得到未知参数的估计值及残差;构建检验统计量,对初始不连续点进行统计检验,求解获取新的未知参数估值;检查新的未知参数的估值与初始未知参数的估计值是否相等;输出参数估计值,即为GNSS基准站坐标时间序列中不连续点。本发明无需先验信息,即可实现不连续点的自动探测,可靠性和适应性强,可为基准站坐标时间序列的精细建模和分析提供准确的基础信息。
Description
技术领域
本发明属于时间序列分析领域,并具体涉及一种有效探测GNSS基准站坐标 时间序列不连续点的方法。
背景技术
GNSS基准站坐标时间序列中的阶跃(Offsets)现象,在时间序列分析中也 被称为不连续点(Discontinuities),通常是由于接收机天线和天线罩更换、大地 震造成的同震位移、GNSS软件处理策略变化以及一些未知的人为因素和环境因 素引起的。一方面,GNSS时间序列中不连续点的存在会降低GNSS坐标时间序 列的精度,如果未知的不连续点没有被探测并正确修正,则会使速度产生偏差,进而影响地球物理学估计量,降低GNSS成果的可靠性和其建立及维持地球参考 框架的精度。另一方面,由于时间序列中多个不连续点类似于一种随机游走过程, 不连续点的存在会改变坐标时间序列的噪声特性。因此,不连续点的探测分析已 经成为GNSS成果分析与应用中的重要问题。
GNSS基准站中的不连续点处理,通常依靠基准站事件记录的方式来设定阶 跃参数,例如:地震事件、接收机天线更换或人为事件等,实际处理时在事件发 生的时间点上增加一个待估参数。但事件记录往往存在疏漏,一些未知的偶然事 件也无法一一记录,这时只能从时间序列本身获取缺失的阶跃信息,因此有必要研究基于数据本身而无需先验信息的阶跃探测方法。随着GNSS基准站网规模的 不断增加、观测数据的不断积累,实现基准站坐标时间序列中不连续性的快速可 靠探测势在必行。
发明内容
本发明主要是针对目前GNSS基准站坐标时间序列分析预处理时,难以准确 自动探测和识别序列中不连续点这一问题,提供一种无需先验信息实现自动探测 与识别序列中的不连续的方法,为后续坐标时间序列的精细建模和分析提供准确 的基础信息,进而获取更为可靠真实的基准站运动特征。
针对上述问题,本发明的主要技术解决方案为:一种有效探测GNSS基准站 坐标时间序列不连续点的方法。所述方法按以下步骤进行:
步骤1,输入GNSS基准站坐标时间序列;
步骤2,将GNSS基准站坐标时间序列中相邻时刻的观测值作差得到GNSS 基准站坐标差值序列,将GNSS基准站坐标差值序列使用四分位距方法进行探测 得到多个初始不连续点;
步骤3,构造GNSS基准站坐标时间序列的参数估计模型,以步骤1所述的 GNSS基准站坐标时间序列作为输入数据,以所述GNSS基准站坐标时间序列的 参数估计模型的截断项、所述GNSS基准站坐标时间序列的参数估计模型的线性 速度、所述GNSS基准站坐标时间序列的参数估计模型的年周期项正弦系数、所 述GNSS基准站坐标时间序列的参数估计模型的年周期项余弦系数、所述GNSS 基准站坐标时间序列的参数估计模型的半年周期项的正弦系数、所述GNSS基准站坐标时间序列的参数估计模型的半年周期项的余弦系数、步骤2所述的多个初 始不连续点作为未知参数,通过优化求解得到未知参数的最优线性无偏估计值及 所述GNSS基准站坐标时间序列的残差序列;
步骤4,根据步骤3所述的残差序列构建检验统计量,对步骤2所述的多个 初始不连续点进行统计检验,每次从步骤3所述构造的GNSS基准站坐标时间序 列的参数估计模型中剔除不显著的初始不连续点,直到所有的初始不连续点经过 统计量检验都显著为止,以所述GNSS基准站坐标时间序列的参数估计模型的截 断项、所述GNSS基准站坐标时间序列的参数估计模型的线性速度、所述GNSS 基准站坐标时间序列的参数估计模型的年周期项正弦系数、所述GNSS基准站坐 标时间序列的参数估计模型的年周期项余弦系数、所述GNSS基准站坐标时间序列的参数估计模型的半年周期项的正弦系数、所述GNSS基准站坐标时间序列的 参数估计模型的半年周期项的余弦系数、检验显著的不连续点作为新的未知参数,通过步骤3所述优化求解得到新的未知参数的估计值及残差序列;
步骤5,检查步骤4获取的新的未知参数的最优线性无偏估计值Xnew与步骤 3获取的未知参数的最优无偏估计值Xinitial是否相等,若相等则进入步骤6,否则 用Xnew替换Xinitial重复步骤3~步骤4。
步骤6,输出新的未知参数的最优线性无偏估计值Xnew,即为GNSS基准站 坐标时间序列中不连续点。
步骤1所述GNSS基准站坐标时间序列为:
Si
i∈[1,m0]
其中,m0表示观测值的数量,Si表示第i个时刻观测值;
步骤2所述GNSS基准站坐标差值序列为:
i∈[1,m0],且
其中,为第i个时刻观测值与第i-1个时刻观测值之差,m0表示观测值的数 量,即时刻的数量,Si表示第i个时刻观测值;
步骤2所述使用四分位距方法进行探测得到多个初始不连续点为:
将GNSS基准站坐标差值序列按从小到大顺序排列,得到排列后GNSS基准 站坐标差值序列;
将排列后GNSS基准站坐标差值序列中位于数据1/2位置处的数据为中位值, 所述中位值用于对排列后GNSS基准站坐标差值序列的集中程度进行度量;
排列后GNSS基准站坐标差值序列位于低四分位点的数值为Q1;
排列后GNSS基准站坐标差值序列位于高四分位点的数值为Q3;
IQR=Q3-Q1,所述IQR用于对排列后GNSS基准站坐标差值序列的分散程 度进行度量;
NIQR=IQR×0.7413,表示排列后GNSS基准站坐标差值序列的一倍中误差:
其中,median表示返回给定序列的中值,w为时间窗口;相应的探测条件 为:
满足以上条件的观测值所在序列中的位置即为不连续点出现的位置,也即完 成了初始不连续点的探测。
步骤2中多个初始不连续点为:
Dj
j∈[1,n0]
其中,n0表示探测出的不连续点总数,Dj表示第j个初始不连续点;
步骤3所述GNSS基准站坐标时间序列的参数估计模型如下:
式中,Si为第i个时刻的观测值,ti为第i个时刻的观测值对应的观测时间, 以年为单位,a为截断项,b为线性速度,c为年周期项正弦系数,d为年周期 项余弦系数,e为半年周期项的正弦系数,f为半年周期项的余弦系数;n0表示 不连续点个数,Dj为第j个初始不连续点,Tj为第j个不连续点出现的观测时刻, H为海维西特阶梯函数,vi为第i个时刻的模型残差,i∈[1,T0], j=1,2,…,n0,T0表示观测值的数量,n0表示不连续点个数。
步骤3所述通过优化求解得到未知参数的估计值及残差序列,具体如下:
步骤3所述的未知参数为:
步骤3所述的GNSS基准站坐标时间序列的参数估计模型进一步转换为:
式中,A为系数矩阵,为残差序列,vi为第i个时刻的所述 GNSS基准站坐标时间序列对应的残差值,T0表示观测值的数量,n0表示不连续 点个数,i∈[1,T0];
所述A的定义如下:
其中,Tj,j=1,2,…,n0为第j个不连续点在所述GNSS基准站坐标时间序列 Si中的位置,n0表示不连续点个数,有n0<m0,m0表示GNSS基准站坐标时间 序列的数量;
步骤3所述未知参数的最优线性无偏估计值为:
X*=(ATA)-1ATS
其中,X*表示未知参数的最优线性无偏估计值,A表示系数矩阵,S为所 述GNSS基准站坐标时间序列,S=Si,i∈[1,m0]
步骤3所述残差序列为:
V=S-AX*
其中,V表示残差序列,X*表示未知参数的最优线性无偏估计值,A表示 系数矩阵,S为所述GNSS基准站坐标时间序列;
作为优选,步骤4构建检验统计量如下:
其中,T为统计量,V为步骤3获得的残差序列,Urem为设置的阈值,一般 Urem=0.01,当T小于指定的阈值Urem时,说明该参数不显著,参数有无对结果的 影响不大,可予以剔除;
本发明具有如下优点:
不依赖准确的先验信息输入,通过模型构建和参数显著性检验,能够自动探 测和剔除不显著参数,识别指定阈值标准以上的未知阶跃,具有较强的可靠性;
遵循从整体到局部,从大到小的原则,实现不连续性的自动探测,增强模型 的适应性。
附图说明
图1:方法流程图
图2:WUHN站坐标时间序列及探测结果;
具体实施方式
下面结合附图和相关实施案例,对本发明的技术方案做进一步的叙述。
实施案例
本发明的实施案例提出了一种有效探测GNSS基准站坐标时间序列不连续 点的方法。所述方法按以下步骤进行:
步骤1,输入GNSS基准站坐标时间序列;
步骤1所述GNSS基准站坐标时间序列为:
Si
i∈[1,m0]
其中,m0=6075表示观测值的数量,Si表示第i个时刻观测值;
步骤2,将GNSS基准站坐标时间序列中相邻时刻的观测值作差得到GNSS 基准站坐标差值序列,将GNSS基准站坐标差值序列使用四分位距方法进行探测 得到多个初始不连续点;
步骤2所述GNSS基准站坐标差值序列为:
i∈[1,m0],且
其中,为第i个时刻观测值与第i-1个时刻观测值之差,m0表示观测值的数 量,即时刻的数量,m0=6075表示观测值的数量,Si表示第i个时刻观测值;
步骤2所述使用四分位距方法进行探测得到多个初始不连续点为:
将GNSS基准站坐标差值序列按从小到大顺序排列,得到排列后GNSS基准 站坐标差值序列;
将排列后GNSS基准站坐标差值序列中位于数据1/2位置处的数据为中位值, 所述中位值用于对排列后GNSS基准站坐标差值序列的集中程度进行度量;
排列后GNSS基准站坐标差值序列位于低四分位点的数值为Q1;
排列后GNSS基准站坐标差值序列位于高四分位点的数值为Q3;
IQR=Q3-Q1,所述IQR用于对排列后GNSS基准站坐标差值序列的分散程 度进行度量;
NIQR=IQR×0.7413,表示排列后GNSS基准站坐标差值序列的一倍中误差:
其中,median表示返回给定序列的中值,w为时间窗口;相应的探测条件 为:
满足以上条件的观测值所在序列中的位置即为不连续点出现的位置,也即完 成了初始不连续点的探测。
步骤2中多个初始不连续点为:
Dj
j∈[1,n0]
其中,n0表示探测出的不连续点总数,n0=6,Dj表示第j个初始不连续点;
步骤3,构造GNSS基准站坐标时间序列的参数估计模型,以步骤1所述的 GNSS基准站坐标时间序列作为输入数据,以所述GNSS基准站坐标时间序列的 参数估计模型的截断项、所述GNSS基准站坐标时间序列的参数估计模型的线性 速度、所述GNSS基准站坐标时间序列的参数估计模型的年周期项正弦系数、所 述GNSS基准站坐标时间序列的参数估计模型的年周期项余弦系数、所述GNSS 基准站坐标时间序列的参数估计模型的半年周期项的正弦系数、所述GNSS基准站坐标时间序列的参数估计模型的半年周期项的余弦系数、步骤2所述的多个初 始不连续点作为未知参数,通过优化求解得到未知参数的最优线性无偏估计值及 所述GNSS基准站坐标时间序列的残差序列;
步骤3所述GNSS基准站坐标时间序列的参数估计模型如下:
式中,Si为第i个时刻的观测值,ti为第i个时刻的观测值对应的观测时间, 以年为单位,a为截断项,b为线性速度,c为年周期项正弦系数,d为年周期 项余弦系数,e为半年周期项的正弦系数,f为半年周期项的余弦系数;n0表示 不连续点个数,Dj为第j个初始不连续点,Tj为第j个不连续点出现的观测时刻, H为海维西特阶梯函数,vi为第i个时刻的模型残差,i∈[1,T0], j=1,2,…,n0,T0表示观测值的数量,n0表示不连续点个数。
步骤3所述通过优化求解得到未知参数的估计值及残差序列,具体如下:
步骤3所述的未知参数为:
步骤3所述的GNSS基准站坐标时间序列的参数估计模型进一步转换为:
式中,A为系数矩阵,为残差序列,vi为第i个时刻的所述 GNSS基准站坐标时间序列对应的残差值,T0表示观测值的数量,n0表示不连续 点个数,i∈[1,T0];
所述A的定义如下:
其中,Tj,j=1,2,…,n0为第j个不连续点在所述GNSS基准站坐标时间序列 Si中的位置,n0表示不连续点个数,有n0<m0,m0表示GNSS基准站坐标时间 序列的数量;
步骤3所述未知参数的最优线性无偏估计值为:
X*=(ATA)-1ATS
其中,X*表示未知参数的最优线性无偏估计值,A表示系数矩阵,S为所 述GNSS基准站坐标时间序列,S=Si,i∈[1,m0]
步骤3所述残差序列为:
V=S-AX*
其中,V表示残差序列,X*表示未知参数的最优线性无偏估计值,A表示 系数矩阵,S为所述GNSS基准站坐标时间序列;
步骤4,根据步骤3所述的残差序列构建检验统计量,对步骤2所述的多个 初始不连续点进行统计检验,每次从步骤3所述构造的GNSS基准站坐标时间序 列的参数估计模型中剔除不显著的初始不连续点,直到所有的初始不连续点经过 统计量检验都显著为止,以所述GNSS基准站坐标时间序列的参数估计模型的截 断项、所述GNSS基准站坐标时间序列的参数估计模型的线性速度、所述GNSS 基准站坐标时间序列的参数估计模型的年周期项正弦系数、所述GNSS基准站坐 标时间序列的参数估计模型的年周期项余弦系数、所述GNSS基准站坐标时间序列的参数估计模型的半年周期项的正弦系数、所述GNSS基准站坐标时间序列的 参数估计模型的半年周期项的余弦系数、检验显著的不连续点作为新的未知参数,通过步骤3所述优化求解得到新的未知参数的估计值及残差序列,新的参数估计 值记为Xnew;
步骤4构建检验统计量如下:
其中,T为统计量,V为步骤3获得的残差序列,Urem为设置的阈值,一般 Urem=0.01,当T小于指定的阈值Urem时,说明该参数不显著,参数有无对结果的 影响不大,可予以剔除;
步骤5,检查步骤4获取的新的未知参数的最优线性无偏估计值Xnew与步骤 3获取的未知参数的最优无偏估计值Xinitial是否相等,若相等则进入步骤6,否则 用Xnew替换Xinitial重复步骤3~步骤4。
步骤6,输出新的未知参数的最优线性无偏估计值Xnew,即为GNSS基准站 坐标时间序列中不连续点。
WUHN站的坐标时间序列来自SOPAC分析中心,参数显著性水平保持控制 Urem=0.01,序列水平方向设置探测的最小阶跃为10mm,最小速率变化为 5mm/yr,序列垂直方向能够探测的最小阶跃为10mm,最小速率变化为10mm/yr, 周期信号的振幅最小值为1mm。利用本发明算法探测不连续点的结果如图2所 示,圆点表示原始时间序列,曲线为拟合曲线,垂直短虚线表示粗差点,垂直长 虚线表示阶跃及阶跃搜索函数曲线。探测的结果如表1所示,N、E方向序列不 存在明显的阶跃和速率变化,仅存在少部分粗差点。U方向存在两处阶跃,其中 历元2002.0699(即2002年1月26)发生明显的阶跃,该时间恰好更换了天线 罩,且天线高垂向改正数比原来大0.041m,分析认为可能是由于当时天线信息 未及时更新,分析中心解算时仍沿用之前的天线高计算,从而造成发生阶跃现象。而另一处阶跃发生在历元2008.8511(即2008年11月7日),设备更换表在当天 未有设备更换事件发生,此阶跃产生的原因尚未可知,但从图2中可见,该阶跃 是显著存在的。
表1:为本发明算法探测的不连续点信息及其估值。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属领 域的技术人员可以对所描述的具体实施例替换成其他基准站序列,做各种各样的 修改或补充,或采用相似方式替代,但并不会偏离本发明的精神或者超越所附权 利要求书所定义的范围。
Claims (3)
1.一种有效探测GNSS基准站坐标时间序列不连续点的方法,其特征在于,包括以下步骤:
步骤1,输入GNSS基准站坐标时间序列;
步骤2,将GNSS基准站坐标时间序列中相邻时刻的观测值作差得到GNSS基准站坐标差值序列,将GNSS基准站坐标差值序列使用四分位距方法进行探测得到多个初始不连续点;
步骤3,构造GNSS基准站坐标时间序列的参数估计模型,以步骤1所述的GNSS基准站坐标时间序列作为输入数据,以所述GNSS基准站坐标时间序列的参数估计模型的截断项、所述GNSS基准站坐标时间序列的参数估计模型的线性速度、所述GNSS基准站坐标时间序列的参数估计模型的年周期项正弦系数、所述GNSS基准站坐标时间序列的参数估计模型的年周期项余弦系数、所述GNSS基准站坐标时间序列的参数估计模型的半年周期项的正弦系数、所述GNSS基准站坐标时间序列的参数估计模型的半年周期项的余弦系数、步骤2所述的多个初始不连续点作为未知参数,通过优化求解得到未知参数的最优线性无偏估计值及所述GNSS基准站坐标时间序列的残差序列;
步骤4,根据步骤3所述的残差序列构建检验统计量,对步骤2所述的多个初始不连续点进行统计检验,每次从步骤3所述构造的GNSS基准站坐标时间序列的参数估计模型中剔除不显著的初始不连续点,直到所有的初始不连续点经过统计量检验都显著为止,以所述GNSS基准站坐标时间序列的参数估计模型的截断项、所述GNSS基准站坐标时间序列的参数估计模型的线性速度、所述GNSS基准站坐标时间序列的参数估计模型的年周期项正弦系数、所述GNSS基准站坐标时间序列的参数估计模型的年周期项余弦系数、所述GNSS基准站坐标时间序列的参数估计模型的半年周期项的正弦系数、所述GNSS基准站坐标时间序列的参数估计模型的半年周期项的余弦系数、检验显著的不连续点作为新的未知参数,通过步骤3所述优化求解得到新的未知参数的估计值及残差序列;
步骤5,检查步骤4获取的新的未知参数的最优线性无偏估计值Xnew与步骤3获取的未知参数的最优无偏估计值Xinitial是否相等,若相等则进入步骤6,否则用Xnew替换Xinitial重复步骤3~步骤4;
步骤6,输出新的未知参数的最优线性无偏估计值Xnew,即为GNSS基准站坐标时间序列中不连续点;
步骤3所述GNSS基准站坐标时间序列的参数估计模型如下:
式中,Si为第i个时刻的观测值,ti为第i个时刻的观测值对应的观测时间,以年为单位,a为截断项,b为线性速度,c为年周期项正弦系数,d为年周期项余弦系数,e为半年周期项的正弦系数,f为半年周期项的余弦系数;n0表示不连续点个数,Dj为第j个初始不连续点,Tj为第j个不连续点出现的观测时刻,H为海维西特阶梯函数,vi为第i个时刻的所述GNSS基准站坐标时间序列对应的残差值,i∈[1,m0],j=1,2,…,n0,m0表示观测值的数量;
步骤3所述通过优化求解得到未知参数的估计值及残差序列,具体如下:
步骤3所述的未知参数为:
步骤3所述的GNSS基准站坐标时间序列的参数估计模型进一步转换为:
式中,A为系数矩阵,为残差序列;
所述A的定义如下:
其中,Tj*,j=1,2,…,n0为第j个不连续点在所述GNSS基准站坐标时间序列Si中的位置,有n0<m0,m0表示GNSS基准站坐标时间序列的数量;
步骤3所述未知参数的最优线性无偏估计值为:
X*=(ATA)-1ATS
其中,X*表示未知参数的最优线性无偏估计值,S=Si,i∈[1,m0];
步骤3所述残差序列为:
V=S-AX*。
2.根据权利要求1所述的有效探测GNSS基准站坐标时间序列不连续点的方法,其特征在于,
步骤2所述GNSS基准站坐标差值序列为:
其中,为第i个时刻观测值与第i-1个时刻观测值之差;
步骤2所述使用四分位距方法进行探测得到多个初始不连续点为:
将GNSS基准站坐标差值序列按从小到大顺序排列,得到排列后GNSS基准站坐标差值序列;
将排列后GNSS基准站坐标差值序列中位于数据1/2位置处的数据为中位值,所述中位值用于对排列后GNSS基准站坐标差值序列的集中程度进行度量;
排列后GNSS基准站坐标差值序列位于低四分位点的数值为Q1;
排列后GNSS基准站坐标差值序列位于高四分位点的数值为Q3;
IQR=Q3-Q1,所述IQR用于对排列后GNSS基准站坐标差值序列的分散程度进行度量;
NIQR=IQR×0.7413,表示排列后GNSS基准站坐标差值序列的一倍中误差:
其中,median表示返回给定序列的中值,w为时间窗口;相应的探测条件为:
满足以上条件的观测值所在序列中的位置即为不连续点出现的位置,也即完成了初始不连续点的探测。
3.根据权利要求2所述的有效探测GNSS基准站坐标时间序列不连续点的方法,其特征在于,步骤4构建检验统计量如下:
其中,T为统计量,V为步骤3获得的残差序列,Urem为设置的阈值,Urem=0.01,当T小于指定的阈值Urem时,说明该参数不显著,参数有无对结果的影响不大,予以剔除。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111294457.4A CN114114332B (zh) | 2021-11-03 | 2021-11-03 | 一种有效探测gnss基准站坐标时间序列不连续点的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111294457.4A CN114114332B (zh) | 2021-11-03 | 2021-11-03 | 一种有效探测gnss基准站坐标时间序列不连续点的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114114332A CN114114332A (zh) | 2022-03-01 |
CN114114332B true CN114114332B (zh) | 2024-06-14 |
Family
ID=80380763
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111294457.4A Active CN114114332B (zh) | 2021-11-03 | 2021-11-03 | 一种有效探测gnss基准站坐标时间序列不连续点的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114114332B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117388872B (zh) * | 2023-09-05 | 2024-03-19 | 武汉大学 | 一种北斗地基增强***参考站坐标框架维持方法和*** |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104537034B (zh) * | 2014-12-22 | 2017-11-10 | 国家电网公司 | 基于时间序列分析的输变电设备的状态监测数据清洗方法 |
CN106066901B (zh) * | 2016-04-22 | 2019-05-07 | 中南大学 | 一种gnss自动化变形监测的基准点稳定性分析方法 |
CN106772498B (zh) * | 2016-11-21 | 2019-05-10 | 华东交通大学 | 一种gps位置时间序列噪声模型建立方法 |
CN106814378B (zh) * | 2017-01-17 | 2019-05-10 | 华东交通大学 | 一种gnss位置时间序列周期特性挖掘方法 |
CN107102342B (zh) * | 2017-04-28 | 2018-09-21 | 武汉大学 | 基于共模误差的gps坐标时间序列不连续性的补足方法 |
CN109188466A (zh) * | 2018-09-29 | 2019-01-11 | 华东交通大学 | 一种顾及非线性变化的gnss基准站地壳运动速度场估计方法 |
CN110058236B (zh) * | 2019-05-21 | 2023-04-07 | 中南大学 | 一种面向三维地表形变估计的InSAR和GNSS定权方法 |
-
2021
- 2021-11-03 CN CN202111294457.4A patent/CN114114332B/zh active Active
Non-Patent Citations (2)
Title |
---|
L1范数与IQR统计量组合的GNSS坐标序列粗差探测算法;明锋;曾安敏;景一帆;;测绘科学技术学报;20160701(第02期);127-132 * |
徐克科.《GNSS监测地壳形变理论与方法》.北京:测绘出版社,2019,(第978-7-5030-4265-2版),80-84. * |
Also Published As
Publication number | Publication date |
---|---|
CN114114332A (zh) | 2022-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106021710B (zh) | 基于大气电离层参数的震前卫星轨道异常识别方法 | |
CN108307301B (zh) | 基于rssi测距和轨迹相似性的室内定位方法 | |
CN106871776A (zh) | 一种基于gnss的实时变形监测*** | |
US20080147327A1 (en) | Method of Seismic Signal Processing | |
CN114114332B (zh) | 一种有效探测gnss基准站坐标时间序列不连续点的方法 | |
CN109029368A (zh) | 像方补偿的遥感影像/sar影像高精度几何定位后处理方法 | |
CN107831516B (zh) | 融合gnss和地面监测网获取大坝实时高精度位移的方法 | |
CN114936201A (zh) | 基于自适应分块神经网络模型的卫星降水数据校正方法 | |
CN117538503A (zh) | 一种实时智能的土壤污染监测***及方法 | |
CN113742927B (zh) | 一种气象预报数据质量检测的方法 | |
CN113157684B (zh) | 一种水利海量数据的查错方法 | |
CN114019139A (zh) | 一种农用地土壤重金属异常数据的探测方法 | |
Ebeling | Ground-based deformation monitoring | |
CN115545112B (zh) | 一种大量地下水实时自动监测数据自动识别和处理的方法 | |
CN106897683B (zh) | 一种遥感图像的地物检测方法及*** | |
CN117970380A (zh) | 一种铁路机车的卫星定位完好性实时判别方法及装置 | |
Aguilar et al. | A theoretical approach to modeling the accuracy assessment of digital elevation models | |
CN116738375A (zh) | 基于单条带测深数据的诱导升沉误差探测消除方法及*** | |
CN112732773B (zh) | 一种继电保护缺陷数据的唯一性校核方法及*** | |
CN115201864A (zh) | 一种检测卫星钟差跳变的方法、装置、存储介质及设备 | |
CN114563771A (zh) | 基于聚类分析的双阈值激光雷达云层检测算法 | |
CN116448062B (zh) | 一种桥梁沉降变形检测方法、装置、计算机及存储介质 | |
CN114020725B (zh) | 一种考虑空间分布的窗口滑动gpm数据订正方法 | |
CN108303713A (zh) | 一种参考站稳定性的自动检验方法 | |
CN109782324A (zh) | 一种铁路巡检定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |