|||
处理过程中,为了自己的方便使用,需要编写自己的脚本(用红色标出)
在gamit中site.defaults中用自己的站点 localrx,选用的较好的igs站点 用ftprx,若自己的站点没有在itrf08.apr的先验坐标中,最好则要运行sh_gamit两次
第一次为sh_gamit -s 2013024 234 -expt sina -nogifs >sh_gamit.log & rm -r 0??
第二次为sh_gamit -s 2013024 234 -expt sina -nogifs >sh_gamit.log &
第一次为计算大约的坐标用rinex的伪距,这个误差越小,得到的精度会更高,第二部用算好的站点位置估计更好约束站点坐标,当然利用sh_gamit其中有许多命令选项,主要就是可以得到单天解以及保存天空图以及检验单天解的质量,或者命令
sh_gamit -expt expt -d 2012 253 254 255 256 -pres ELEV-orbit IGSF
当然在处理gamit之前或者以后可能用到以下脚本:
sh_nrms ../???/qkasha.??? 得到nrms 计算重复度 重复度并不能代表点位精度
sh_plot_nrms得到单天解的nrms值 cp???/o*a* ooo/
sh_plot_elmean -f../20*/autcln.post.sum -s ALL -R 5/90/-15/15 计算平均相位误差
sh_oneway -sites TASH -autsumautcln.post.sum -prefix DPH -scale 0.5 -page L -grdimage -grdview -elev
进行在doy目录——单天解目录中,进行lc残差检验,生成3D天空图影像,生成1,2,4,6分别生成。
sh_plotatmtide -site GRAS -apr~/gg/tables/itrf00.apr 画超过两天的大气潮汐负载NEU分量
sh_gamit_baseline -o 3??/osmb1a.*-frame 1 -anot 5 可以进行基线的重复性
sh_gamit_baseline -o 3??/osmb1a.* -bCLIF_PIE1 PIE1_VVRS -frame 1 -anot 5 -a CLIF_VVRS # Notethe order
sh_reprep 得到基线重复性的解 将sh_gamit_summary放到ppp文件中然后运行
sh_plot_WL_NL 生成wl与nl的值 在rinex文件下运行plot_station 得到rinex的时间间隔图形
Sh_plot_rinex 与脚本Plot_station 画rinex的时间站点图形
convertc/cortran 转换坐标
在单天解的sh—gamit—summarry中得到评价标准:
1、每个站点的平均的RMS一般小于8
2、一些站点在lfile中得到更新如果调整了1m则进行重新处理。一米先验坐标可以带来1mm的误差。查看q*a文件, 有先验坐标的调整以及误差等
3、post nrms一般小于2左右或者2左右。
4、基线的不连续性应该处于相同量级,坐标约束应用于站点坐标和轨道参数。
5、在sh总结文件中最好的两个站在3-5mm,最差的站点在7-9mm,若是在10-15mm说明是一个高的噪声,
大于15mm代表一个很差的接收机,一个多路径影响很厉害,严重天气问题,或者autcln处理短间隔的问题
应该查看天空图 与高度图,查看噪音来源,若两个最差站点噪音都大于15,则看autcln.post.rms查看其他的站点的rms
grep RMS 027/sh*
6、在antcln.post.sum中看到伪距的rms范围在1-2m,大于5m代表这是坏的站点
7、模糊度结算WL 在p-85为噪音多的天数,大于90%为好的天数 低的WL代表不好的伪距,一个很低的NL百分比相对于WL说明几何学或者大气误差
大于80%的NL也为好的天数 grep WL 027/sh*
8、所有的站点都包括在内 gerp -i stataions 0??/sh* 得到总的x文件也就是rinex的文件个数若x文件少了,则说明tracking的时间段少 舍弃了
可以在process.defaults里设置minxf
9、在q文件的a文件中,adjust大于250px表示,先验坐标需要重新调整,然后重新run sh_gamit
天空图与高度相位图可以评估多路径,水汽,以及天线想为中心模型
在高度与相位图中
常见错误以及操作:
1、没有rinex时段的数据,检查rinex数据文件或者makex.expt.infor
2、太少数据,xfile文件太小没有用到,检查rinex文件大小,改变minxf in process.defaults
3、站点在解后没有数据或者调整,可能先验坐标大于10m,检查autcln.prefit.sum 的范围rmS
sh_rx2apr differentially for several RINEX files
(差的接收机):一般检查C文件或者rinex文件
4、nrms在qa文件大于0.2 解释超出约束,检查GCXvsGLXnrms ,重新运行仅仅一个站点约束时。
经过GAMIT处理后,可以初步首先进行重复性检验
打开tables进行链接
ln -s pmu.bull_a pmu.usno
返回工程目录进行,运用区域站作为框架,进行重复性检验
sh_glred -s 2012 241 2012 266 -expt kash -opt H G E
由于bjfs有几天的数据不能下载,并且数据转换中包含gnss的数据,xian(西安)及其wuhn(武汉)的数据下载基本没有,
lhaz(拉萨)的数据,时间序列NE方向的数据不是很好,然后bish(中亚的一个站点)及其eveb(喜马拉雅的南侧一个站点)没有数据,guao整体的数据质量一般,
最终重新选择hyde,irkt加入 进行gamit处理及其重复性检验,平差后进行,后面然后利用两个cmd进行globk处理,具体详细见以下处理。
将gg的tables里面的编辑globk_combNaNd glorg_combNaNd重新进行序列的运算。
观察时间序列后,删除误差较大的天数数据,重新进行gamit的处理,重新选站点,254天出现异常,怀疑地震的牵扯。
以上仅仅是 gamit处理后,简单的分析时间序列的重复性,检查误差,如果误差不好,则剔除不好的站点或者单天解,重新处理GAMIT解
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-20 01:35
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社