下面是GNSSer非差事后精密定位测试情况,版本为2018.06之后,对应公开版为 GNSSer 1.4 。
〇、数据说明与测试环境
测试数据来自IGS参考站,采样率为30s,分别采用46个、321个、465个,以及对BJFS站 2744个文件,进行了测试,最后测试了单站8年 2744个文件 的数据,以IGS当日解算为真值,计算结果与之比较,计算偏差。
采用30秒IGS钟差,观测数据采样率30秒,时段为24小时,计算到最后历元(23:59:30);
为了加快计算效率和减少内存使用,可以关闭实时信息输出,关闭历元信息输出;
测试电脑和配置:Windows 10, I7-7700HQ 处理器,内存 16GB 的ThinkPad笔记本,Dell T5810工作站,CPU E5-1620 3.50GHz 28g内存。
一、46个测站计算结果
采用2013.001(GNSSer 1.3 以前默认测试数据),欧洲区域的46个站。
计算结果如下:
参数 | 平均数 | 均方差 |
X | 0.003155 | 0.002615 |
Y | 0.004515 | 0.002587 |
Z | 0.003155 | 0.002615 |
E | 0.003451 | 0.00258 |
N | 0.002263 | 0.004539 |
U | 0.004089 | 0.004097 |
Len | 0.00652 | 0.003376 |
三维绝对偏差平均:0.65 厘米
46个测站的分布图
二、321个GNSSer原始测试站
三维平均偏差 0.72厘米,各分方向平均绝对值偏差: E方向 0.31厘米,N方向0.28厘米,U方向0.59厘米;
三维偏差大于 3 厘米0个,大于 2 厘米 8个(含大于3厘米的);
最大内存耗费378.19MB,平均 6.44 秒一个。
分布图和偏差图如所示
三、收敛时间
将上述321个站按3h时段长度分割,总计321*8个文件,收敛时间的分布如下图所示,收敛的评判标准是当前历元及向后20个历元的E N U均小于1分米,平均收敛时间为50.7个历元,即25.35min。
四、2018001测试情况
下载IGS站2008.001当天的压缩观测数据,共725个文件,850MB,地址:ftp://cddis.gsfc.nasa.gov/pub/gps/data/daily/2018/001/18d/,注:18o文件夹下只有496个文件,数据量1.24GB;
采用工具解压后剩余552个,其它报告为解压缩错误,可能网络传输有误,尚不知具体原因,GNSSer解压和工具解压的数量相同;
采用GNSSer执行格式化,只保留GPS系统,30s采样率,发现有部分文件头部坐标格式错误,但最终都完全格式化成功;
过滤1M和单频数据,余540个文件;
将此540个文件执行PPP计算,取消历元输出,取消Info显示;
成功计算539个(其中MOBK站一直报错,计算结果误差巨大),每个 4.64 秒,最大内存 586MB,计算结果内存454MB;
与 igs18P19821_all.snx 互差比较有371个公共站,超过1米有4个 1 % ,[MAD2,QUI4,QUI3,UNSA], 超过10厘米共8个, 2.1 %, [PTAG,OHI2],超过5厘米共10个, 2.7 %。
去掉互差大于 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 个文件并行计算
解算了bjfs,2010.01.01-2017.12.31,2744个文件,
格式化后只计算GPS系统,1个由于星历下载不全而计算失败
全程自动下载星历钟差(本地产品库有约2年的数据),耗时 1天16小时,平均53s一个,最大耗费内存6.39GB。
弹窗结果和坐标偏差下图
注意:这里53秒一个,包含全自动下载星历和钟差的时间,如果需要获得更快速的计算时间,请提前下载星历和钟差。
剔除10余个粗差,得到以下ENU三方向时序图。从图中,可以明显的发现趋势性和周期性。
注意这里2015年底到2016年初有一段时间高程方向有6cm的偏差,经过分析发现,见下图,IGS结果在这段时间发生了跳跃。
直接将偏差作个分量的计算中误差,结果如下:
RmseOf_dX | 0.0036821307822060968 |
RmseOf_dY | 0.0069533813515795986 |
RmseOf_dZ | 0.0052515676528004018 |
RmseOf_dE | 0.0041634755107285782 |
RmseOf_dN | 0.00232973367856501 |
RmseOf_dU | 0.0081684989311517874 |
RmseOf_length | 0.0094597337240381486 |
下一步将重点,针对误差较大的测站进行分析,并且提供可靠的精度评价和报警方法。
尚无评论!