pan11的个人博客分享 http://blog.sciencenet.cn/u/pan11

博文

GAMIT处理第二部分

已有 6636 次阅读 2017-11-18 11:15 |个人分类:GPS学习|系统分类:科研笔记| gps处理, GAMIT处理, GAMIT处理

处理过程中,为了自己的方便使用,需要编写自己的脚本(用红色标出)

在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解




https://wap.sciencenet.cn/blog-651121-1085780.html

上一篇:GAMIT处理第一部分
下一篇:QOCA命令简介以及简单使用
收藏 IP: 101.75.214.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-3-29 01:44

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部