一种GNSS区域网中异常点探测方法及系统技术方案

本发明专利技术提供一种GNSS区域网中异常点探测方法及系统，包括获取GNSS测站坐标时间序列北方向与东方向的观测值，并获取GNSS测站坐标；建立各测站GNSS坐标时间序列北方向与东方向观测值的速度模型，获取速度估值及其不确定性；利用GNSS成像方法，分别生成研究区域的北方向与东方向的GNSS影像；依据GNSS影像结果，使用聚类分析方法对研究区域进行块体划分；对块体任选3个测站计算欧拉矢量与应变参数，分别形成聚类空间，并分别计算各聚类点至聚类中心的距离；利用异常值探测方法，获取欧拉矢量和应变率参数的异常聚类点，分别记录异常欧拉矢量与异常应变参数对应的测站，并统计各测站出现的次数，大于阈值则认为异常点，对相应异常点集合取交集则获得异常速率点集合。

【技术实现步骤摘要】
一种GNSS区域网中异常点探测方法及系统
本专利技术属于GNSS数据精密处理
，具体涉及一种基于GNSS多测站坐标时间序列的异常点探测技术方案。
技术介绍
GNSS多测站坐标时间序列反映了测站位置随时间变化的规律及在空间上的分布特征。利用GNSS多测站坐标时间序列进行地学现象空间分布特征已成为现实(Hammondetal.,2016)，可构建兼顾时间和空间分辨率的地壳垂向形变影像图。但是，在空间上的异常点会严重影响地壳形变影像图的效果及其地学解释。其中，GNSS为全球导航卫星系统简称。本专利技术提供了一种基于GNSS影像的异常点探测方法，以每个测站最临近点的速度及其不确定性作为空间滤波的基础，结合欧拉矢量反演结果，以欧拉矢量和应变量作为关键参数，利用聚类方法(SavageandSimpson,2013)实现对异常点的探测。
技术实现思路
针对现有技术存在的问题，本专利技术提供了一种以每个测站最临近点的速度及其不确定性作为空间滤波的基础，结合欧拉矢量反演结果，以欧拉矢量聚类实现异常点探测的技术方案。为解决上述技术问题，本专利技术提供的技术方案为一种GNSS区域网中异常点探测方法，包括以下步骤，步骤1，获取GNSS测站坐标时间序列北方向与东方向的观测值，并获取每个GNSS测站的坐标(Bi,Li)，其中，i＝1,2,…,n，n为总测站数；Bi，Li为第i个测站在大地坐标系下的纬度与经度；步骤2，建立各测站GNSS坐标时间序列北方向与东方向观测值的速度模型，获取各测站北方向与东方向的速度估值vni,vei及其不确定性σni,σei，i＝1,2,…,n；步骤3，利用GNSS成像方法，分别生成研究区域的北方向与东方向的GNSS影像；步骤4，依据步骤3所得北方向与东方向的GNSS影像结果，使用聚类分析方法对研究区域进行块体划分，设得到l个块体，记各块体内测站数量为numo，o＝1,2,…,l；步骤5，对第o个块体的numo个测站中，任选3个测站计算欧拉矢量EVo,p与应变参数Straino,p，o＝1,2,…,l；步骤6，基于聚类分析方法，分别形成欧拉矢量与应变参数的聚类空间，设获得第o个块体的欧拉矢量聚类中心CEVo与应变参数聚类中心CStraino，并分别计算聚类空间中各聚类点EVo,p至聚类中心CEVo的距离，计算各聚类点Straino,p至聚类中心CStraino的距离；步骤7，基于各聚类点至其聚类中心的距离，利用异常值探测方法，获取欧拉矢量和应变率参数的异常聚类点，设获得对应的异常欧拉矢量EVo,q与异常应变参数Straino,r，其中，q＝1,…,abno,1；r＝1,…,abno,2，abno,1与abno,2为第o个块体中异常欧拉矢量与异常应变参数个数；步骤8，分别记录异常欧拉矢量与异常应变参数对应的测站，并统计各测站出现的次数，若大于设置的阈值，则认为该测站为可能的异常点，分别获得异常点集合{OutlierEV}和{Outlierstrain}；步骤9，对异常点集合{OutlierEV}与{OutlierStrain}取交集，则获得异常速率点集合{Outlier}。而且，步骤2中所述建立各测站GNSS坐标时间序列北方向与东方向观测值的速度模型，通过最小二乘估计或极大似然估计方法实现。而且，步骤3中，利用基于站间相关性空间结构函数构建的GNSS成像方法，分别生成研究区域的北方向与东方向的GNSS影像。而且，步骤4中，采用k均值聚类方法对研究区域进行块体划分。而且，步骤6中，基于k均值聚类分别形成欧拉矢量与应变参数的聚类空间。而且，步骤7中，所述利用异常值探测方法获取欧拉矢量和应变率参数的异常聚类点，是基于拉依达准则或巴尔达准则实现。而且，步骤8中，阈值优选设置为3。本专利技术还相应提供一种GNSS区域网中异常点探测系统，用于如上所述GNSS区域网中异常点探测方法。与现有技术相比，本专利技术具有特点：本专利技术基于临近点的速度及其不确定性作为空间滤波的基础，形成的GNSS影像是基于实测数据的结果，具有较高的可靠性；而以GNSS影像结果对GNSS测站覆盖区域进行块体划分，很大程度上保证了划分块体的整体变形特征，尽可能减少了由于块体划分不准确所引起的欧拉矢量聚类分析的出错率；与另一重要物理量——应变结合的异常点探测，进一步保证了结果的准确性。附图说明图1为本专利技术实施例方法的具体流程示意图。具体实施方式为使本专利技术目的、技术方案及有益效果更加清楚明白，下面将结合附图及实施例，进一步说明本专利技术。参见图1，本专利技术实施例提供一种GNSS多测站坐标时间序列异常点探测方法，包括步骤：步骤1，获取GNSS测站坐标时间序列北方向与东方向的观测值，并获取每个GNSS测站的坐标(Bi,Li)(i＝1,2,…,n)，其中，n为总测站数；Bi，Li为第i个测站在大地坐标系下的纬度与经度；GNSS测站坐标时间序列反映了测站位置随时间变化的规律特征。具体实施时，可以通过双差定位软件工具、精密单点定位软件工具或IGS分析中心获取积累的GPS测站坐标时间序列观测值。步骤2，建立各测站GNSS坐标时间序列北方向与东方向观测值的速度模型，获取各测站北方向与东方向的速度估值vni,vei及其不确定性σni,ei(i＝1,2,…,n)；具体实施时，可基于现有技术(如最小二乘估计、极大似然估计等)，建立各测站GNSS坐标时间序列的速度模型，获取速度估值vni,vei及其不确定性σni,σei。步骤3，利用GNSS成像方法，分别生成研究区域的北方向与东方向的GNSS影像；本步骤实现可参见现有技术，例如申请人的专利技术《基于站间相关性空间结构函数构建的GNSS成像方法》，其中提出一种基于测站对相关系数的空间结构函数的构建方案，实现GNSS成像。该方法包括输入GNSS测站坐标时间序列观测值，以及每个GNSS测站的速度及不确定性；结合研究区域的地质学与大地测量学结果(如地质断层位置信息、基于长期GPS观测的活动块体划分模型等)，对GNSS网内的测站进行聚类划分，得到聚类区域；计算各聚类区域内任意两测站所构成测站对之间的相关系数，根据测站间相关系数进行数据池划分，获取每个聚类区域中若干数据池及在各数据池内的GNSS测站对；在每个聚类区域内，计算每个数据池中所有测站对相关系数的中值和绝对中位差，构建各聚类区域的空间结构函数，标准化形成整个GNSS网最终的空间结构函数；依据速度不确定性和空间结构函数，确定研究范围内所有测站的权，利用空间插值方法进行空间插值，形成影像。因此，本专利技术优选利用该技术，实现构建GNSS的形变影像图。步骤4，依据步骤3中北方向与东方向的GNSS影像结果，使用聚类分析方法对研究区域进行块体划分，设得到l个块体，记各块体内测站数量为numo(o＝1,2,…,l)；具体实施时，可对北方向与东方向的GNSS影本文档来自技高网...

【技术保护点】
1.一种GNSS区域网中异常点探测方法，其特征在于：包括以下步骤，/n步骤1，获取GNSS测站坐标时间序列北方向与东方向的观测值，并获取每个GNSS测站的坐标(B

【技术特征摘要】
1.一种GNSS区域网中异常点探测方法，其特征在于：包括以下步骤，
步骤1，获取GNSS测站坐标时间序列北方向与东方向的观测值，并获取每个GNSS测站的坐标(Bi，Li)，其中，i＝1，2，…，n，n为总测站数；Bi，Li为第i个测站在大地坐标系下的纬度与经度；
步骤2，建立各测站GNSS坐标时间序列北方向与东方向观测值的速度模型，获取各测站北方向与东方向的速度估值vni，vei及其不确定性σni，σei，i＝1，2，…，n；
步骤3，利用GNSS成像方法，分别生成研究区域的北方向与东方向的GNSS影像；
步骤4，依据步骤3所得北方向与东方向的GNSS影像结果，使用聚类分析方法对研究区域进行块体划分，设得到l个块体，记各块体内测站数量为numo，o＝1，2，…，l；
步骤5，对第o个块体的numo个测站中，任选3个测站计算欧拉矢量EVo，p与应变参数Straino，p，o＝1，2，…，l；
步骤6，基于聚类分析方法，分别形成欧拉矢量与应变参数的聚类空间，设获得第o个块体的欧拉矢量聚类中心CEVo与应变参数聚类中心CStraino，并分别计算聚类空间中各聚类点EVo，p至聚类中心CEVo的距离，计算各聚类点Straino，p至聚类中心CStraino的距离；
步骤7，基于各聚类点至其聚类中心的距离，利用异常值探测方法，获取欧拉矢量和应变率参数的异常聚类点，设获得对应的异常欧拉矢量EVo，q与异常应变参数Straino，r，其中，q＝1，…，abno，1；r＝1，…，abno，2，abno，1与abno，2为第o个块体中异常欧拉矢量与异常应变参数个数；
步骤8，分别记录异常欧拉矢量与异常应...

【专利技术属性】
技术研发人员：周晓慧，姜卫平，杨艺林，陈华，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北;42