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

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

【技术实现步骤摘要】
一种GNSS区域网中异常点探测方法及系统
本专利技术属于GNSS数据精密处理
，具体涉及一种基于GNSS多测站坐标时间序列的异常点探测技术方案。
技术介绍
GNSS多测站坐标时间序列反映了测站位置随时间变化的规律及在空间上的分布特征。利用GNSS多测站坐标时间序列进行地学现象空间分布特征已成为现实(Hammondetal.,2016)，可构建兼顾时间和空间分辨率的地壳垂向形变影像图。但是，在空间上的异常点会严重影响地壳形变影像图的效果及其地学解释。其中，GNSS为全球导航卫星系统简称。本专利技术提供了一种基于GNSS影像的异常点探测方法，以每个测站最临近点的速度及其不确定性作为空间滤波的基础，结合欧拉矢量反演结果，以欧拉矢量和应变量作为关键参数，利用聚类方法(SavageandSimpson,2013)实现对异常点的探测。
技术实现思路
针对现有技术存在的问题，本专利技术提供了一种以每个测站最临近点的速度及其不确定性作为空间滤波的基础，结合欧拉矢量反演结果，以欧拉矢量聚类实现异常点探测的技术方案。为解决本文档来自技高网...

【技术保护点】
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