软件教程(tutorial)

GNSSer 非差定位(PPP)测试情况 2018.6

2018-5-3 20:09:07 来源: GNSSer 发布人: czsgeo


下面是GNSSer非差事后精密定位测试情况,版本为2018.06之后,对应公开版为 GNSSer 1.4 。


〇、数据说明与测试环境

  1. 测试数据来自IGS参考站,采样率为30s,分别采用46个、321个、465个,以及对BJFS站 2744个文件,进行了测试,最后测试了单站8年 2744个文件 的数据,以IGS当日解算为真值,计算结果与之比较,计算偏差。

  2. 采用30秒IGS钟差,观测数据采样率30秒,时段为24小时,计算到最后历元(23:59:30);

  3. 为了加快计算效率和减少内存使用,可以关闭实时信息输出,关闭历元信息输出;

  4. 测试电脑和配置:Windows 10, I7-7700HQ 处理器,内存 16GB 的ThinkPad笔记本,Dell T5810工作站,CPU E5-1620 3.50GHz 28g内存。


一、46个测站计算结果


采用2013.001(GNSSer 1.3 以前默认测试数据),欧洲区域的46个站。

计算结果如下:

参数 平均数 均方差
X0.0031550.002615
Y0.0045150.002587
Z0.0031550.002615
E0.0034510.00258
N0.0022630.004539
U0.0040890.004097
Len0.006520.003376

三维绝对偏差平均:0.65 厘米

201806-46ENU.jpg

46个测站的分布图




二、321个GNSSer原始测试站

  1. 三维平均偏差 0.72厘米,各分方向平均绝对值偏差: E方向 0.31厘米,N方向0.28厘米,U方向0.59厘米;

  2. 三维偏差大于 3 厘米0个,大于 2 厘米 8个(含大于3厘米的);

  3. 最大内存耗费378.19MB,平均 6.44 秒一个

分布图和偏差图如所示










三、收敛时间


将上述321个站按3h时段长度分割,总计321*8个文件,收敛时间的分布如下图所示,收敛的评判标准是当前历元及向后20个历元的E N U均小于1分米,平均收敛时间为50.7个历元,即25.35min。




四、2018001测试情况

  1. 下载IGS站2008.001当天的压缩观测数据,共725个文件,850MB,地址:ftp://cddis.gsfc.nasa.gov/pub/gps/data/daily/2018/001/18d/,注:18o文件夹下只有496个文件,数据量1.24GB;

  2. 采用工具解压后剩余552个,其它报告为解压缩错误,可能网络传输有误,尚不知具体原因,GNSSer解压和工具解压的数量相同;

  3. 采用GNSSer执行格式化,只保留GPS系统,30s采样率,发现有部分文件头部坐标格式错误,但最终都完全格式化成功;

  4. 过滤1M和单频数据,余540个文件;

  5. 将此540个文件执行PPP计算,取消历元输出,取消Info显示;

  6. 成功计算539个(其中MOBK站一直报错,计算结果误差巨大),每个 4.64 秒,最大内存 586MB,计算结果内存454MB;

  7. 与       igs18P19821_all.snx  互差比较有371个公共站,超过1米有4个 1 % ,[MAD2,QUI4,QUI3,UNSA], 超过10厘米共8个,  2.1 %, [PTAG,OHI2],超过5厘米共10个, 2.7 %。

  8. 去掉互差大于 5 厘米的测站,各个分量的绝对平均值如下表所示,单位:厘米。

    Site

    X

    Y

    Z

    E

    N

    U

    Len

    绝对值平均

    0.031

    0.055

    0.031

    -0.010

    -0.133

    0.315

    0.952

    以下是三方向的偏差和三维绝对偏差图



五、单站 8 年 2744 个文件并行计算


  1. 解算了bjfs,2010.01.01-2017.12.31,2744个文件,

  2. 格式化后只计算GPS系统,1个由于星历下载不全而计算失败

  3. 全程自动下载星历钟差(本地产品库有约2年的数据),耗时 1天16小时,平均53s一个,最大耗费内存6.39GB。  

  4. 弹窗结果和坐标偏差下图

结果统计.jpg

注意:这里53秒一个,包含全自动下载星历和钟差的时间,如果需要获得更快速的计算时间,请提前下载星历和钟差。

剔除10余个粗差,得到以下ENU三方向时序图。从图中,可以明显的发现趋势性和周期性。


注意这里2015年底到2016年初有一段时间高程方向有6cm的偏差,经过分析发现,见下图,IGS结果在这段时间发生了跳跃。





直接将偏差作个分量的计算中误差,结果如下:

RmseOf_dX0.0036821307822060968
RmseOf_dY0.0069533813515795986
RmseOf_dZ0.0052515676528004018
RmseOf_dE0.0041634755107285782
RmseOf_dN0.00232973367856501
RmseOf_dU0.0081684989311517874
RmseOf_length0.0094597337240381486

下一步将重点,针对误差较大的测站进行分析,并且提供可靠的精度评价和报警方法。





阅读次数: 2354

下一篇: GNSSer 桌面公开版 1.3 更新说明
上一篇: 在Windows上安装GNSSer网络计算服务

评论次数: 1
序号 用户 标题 时间
1 楼 lly 计算精度更新 2018-6-4 15:01:21
6.4,lly发现gpt2w服务出现错误,key错了,所以所有站用的同一个对流层改正。 改正前 2013001 327个站的精度(E N U):0.003224241 0.002755141 0.013339781. 改正后 2013001 327个站的精度(E N U):0.003357 0.002708 0.006325。 U方向精度提升显著(1倍多)。
返回上一页面