gamit/globk解算数据质量精度指标解释及常见问题解释
补充解算策略
本文进行Gamit/Globk数据处理时,下载邻近测区周边9个IGS的事后精密星历、广播星历文件,采用以下策略:
1.process.defaults中修改“set sint”采样间隔为5s,“set nepc”观测历元为2160s,set.aprf=control.apr平差时的控制点;
2.sestbl.中修改Elevation Cutoff高度截止角改为15°,基线解算类型Choice of Experiment=BASELINE,基线处理类型Choice of Observable=L1,L2_INDEPENDENT,由于测站不邻近海洋取消海洋负荷潮Use met.grid=N,Use otl.grid=N,Use atml.grid=N;
3.station.info手动更新控制点与待解算点测站信息,天线高量取方式DHARP;
4.sittbl.中修改所有测站约束平面0.02,高程0.05;
5.lfile.中配置所有站点的近似坐标,对于IGS测站的坐标从itrf08.apr框架中提取,对于待解算测站的初始坐标利用双差定位法提取;
6.globk_comb.cmd中修改禁用地震文件eq_file,增加输入卫星记录文件make_svs @.svs Z,改变控制点文件apr_file ../tables/controlpoint.apr路径,改变prt_opt GDLF CMDS自由网平差输出文件,改变org_opt PSUM CMDS GDLF约束网平差输出文件,启用所有控制点参与平差use_site all;
7.glorg_comb.cmd修改stab_it网平差迭代次数为5次,设置已知点stab_site bjfs chan chum irkj lhaz shao tash ulab urum,apr_file ../tables/controlpoint.apr设定网平差控制文件路径,禁用source。
1. Q文件分析
qhn01a.094
在该平差结果文件中有对解算精度和可靠性进行评估的指标,其中均方根残差nrms是衡量单天解质量的重要指标之一,通常比较理想的值应小于0.2,如果该值大于0.3,就意味着处理过程中未除去在大的周跳或某一参数的解算存在很大偏差,或者解算模型设定有误,具体原因可进一步在autcln.post.sum中查找。
2.autcln.sum文件
尽管nrms是合限的,还是建议仔细查看autcln.sum(autcln.post.sum/ autcln.out)文件,进一步确定观测站数据和卫星数据有无异常情况,其中最重要的内容之一是对接收机时钟的稳定性统计情况.
检查autcln.post.sum中的观测数据有效性,右边表格最后一列可用数据个数,下面表格最后good列代表有效地数据个数。当有效数/可用数较高,并且有效数绝对量大于18000以上结果基本可用。
3. 天空图精度指标
同一个地点在各个的时间段有高残差预示着存在多路径效应
同一个地点在特定的时间段有高残差预示着存在水波折射
整体趋势线在有的高角度离开蓝色线,表明错误引用了天线相位模型,或卫星信号受到了电磁干扰
正常:红色整体趋势线平顺,在蓝色线上下均匀波动
qexpta.doy-Adjust
正常 Adjust 不等于0
若等于0,则该测点未参与计算
qexpta.doy-Receiver/Antanna/Ht to APR
qexpta.doy
正常 N E < 1cm ,U<2cm
如果 N E U > 3cm ,则数据质量有问题
qexpta.doy-nrms
一般要求 Postfit nrms < 0.3
若Postfit nrms > 1,则表明计算的结果不可靠
autcln.post.sum
sh_gamit_doy.summary
1) Number of stations used:表示参与计算的测站数
2) Total xfiles:创建X文件的总数
注:通常1、2相等,若不等,则有的测站观测时间太短,没有参与计算
也是可能process.defaults中设置了某些测站没有参与计算
1) 最好测站RMS[3 5]mm
2) 最差测站RMS[7 9]mm
3) 有较大的噪声[10 15]mm,大于15mm则说明测站数据有问题
注:如果最好的测站的RMS为0,则该测站已经被移除
1) WL(宽巷模糊度) > 90%
2) NL(窄巷模糊度) > 80%
4 常见问题分析和处理
本部分来自:GAMIT软件GPS测量数据处理流程及常见问题分析和处理_陈永祥 一文
在利用GAMIT软件进行基线解算过程中,由于观测数据、共用表及测站相关文件和过程控制文件,准备中可能存在一些错误,将会导致软件运行非正常中断,出错的原因可以通过分析Gami t. fatal文件中的错误提示来判断,并做出相应的改正。
Gami t. fatal文件中的错误:
1 Cannot find site code 168P on L-file vg_ in
2 Neither T- nor G-file available (Name trism9. 105)
3 Error opening navigation file: erism9. 105 ERROR 2
4 Error opening file: risen. makex. batch ERROR 2
5 Receiver code TPPLEG not found in rcvant.dat
6 input antenna type TPPC3D with alias TPPC3D not in rcant.dat
7 Antenna code (TPSC3D DHPAA) for:168P 2009 105 1 not in hi.dat
8 Error reading OCLEANLOD values ERROR 501
9 Globk平差后无结果或中误差超限