GAMIT学习资料整理(2)

记录各测站的站点、站点名称、开始观测时间、结束观测时间、天线高、 天线量测方法、接收机类型、软硬件版本号、天线类型等信息。严格按照已经给出的跟踪站格式编辑。(注意:要现在当前目录下拷入station.info文件,把内容删掉留下前面的几行,否则后面的步骤没法进行)

具体操作方法有二:

1、手动操作,按照相应文件格式对齐操作;可以从样本格式开始编辑。

2、自动提取观测文件的有关信息,生成 station.info 文件。

在各时段文件目录下操作:

sh_upd_stnfo –files site1doy0.yyo site2doy0.yyo –orbt igsf –u –ref ../tables/station.info

运行上述命令后,就自动从观测文件 0-files中提取重要信息,生成 station.info 文件。

我对station.info文件的准备(武大郭博峰):

1、rinex文件夹下新建文件夹all_station(存放所有点)

2、工程观测文件(o-files)复制到all_station,框架网(o-files)复制到all_station

3、进入Terminal

cd test/rinex/all_station

sh_upd_stnfo –files *.11o -ref station.info(此文件为软链接文件,从~/不常更新表/station.info 链接至指定文件夹下,且只有一个表头)

4、查看station.info文件是否正确

gedit station.info

5、将得到的station.info复制到lfiles文件夹下

2) l-file:

记录各测站的近似坐标: Gamit要求是球面坐标,即l-file。将所有 o 文件整理在一起,然后分别运行如下命令:

grep POSITION *.11o > lfile.rnx

rx2apr lfile.rnx 2011 250 (注:根据所有的o文件中的x、y、z取平均值,文件名为lfile.rnx.apr)

gapr_to_l lfile.rnx.apr lfile. “ “ 2011 250

我对l-file文件的准备(武大郭博峰):

1、在rinex文件夹下新建文件夹all_data(存放所有点)

2、进入该文件夹:cd test/rinex/all_data

3、提取POSITION

grep POSITION *.11o >lfile.xyz

4、查看lfile.xyz文件是否正确

–s ../tables/soltab.2011.J2000 ./soltab.

ln –s ../tables/nutabl.2011 . /nutabl.

ln –s ../tables/pole.usno . /pole.

ln –s ../tables/ut1.usno . /ut1.

ln –s ../tables/pmu.usno . /pmu.

另外,潮汐改正文件也是很容易出错的文件。以10.4版本为例, 默认的潮汐 改正文件名为otl.grid(此文件为损坏文件)(链接前是otl._FES2004.grid )。一定要检查一下 otl.grid 是否是真实链接,出现死链接的情况需要手动链接。命令为:

ln -s ../tables/otl._FES2004.grid ./otl.grid (我在算的时候没有用这个文件,下次试试用)

(注意:链接时,如果该文件夹中已有该文件,无论链接成功与否,要重新链接时都要先删除原文件在进行链接。)

、数据处理阶段

一、批处理

1) 运行makexp 程序,生成输入文件

sh_makexp –expt test –orbt igsf –yr 2011 –doy 250 –sess 99 –srin –nav brdc2500.11n –apr lfile. –sinfo 15 00 00 5760

或者:直接运行 makexp,然后按提示输入

2) 运行sh_sp3fit 脚本,生成轨道初始根数;

(曾由于ut1.链接的不是最新的,导致错误。)

sh_sp3fit –f igs15623.sp3 –o igsf –d 2011 250 –r BERNE –t –u

3) 运行sh_check_sess 脚本,检查卫星一致性;

sh_check_sess –sess 250 –type gfile –file gigsf1.250

4) 运行makej 程序,生成卫星钟差文件;

makej brdc2500.11n jbrdc1.250

5) 运行sh_check_sess,检查卫星一致性;

sh_check_sess –sess 250 –type jfile –file jbrdc1.250

6) 运行makex 程序,生成X 文件;

makex test.makex.batch

7) 运行fixdrv 程序,生成批处理文件;

fixdrv dtest1.250 (进行这步可能发生错误,提示

FATAL :121202:2337:41.0 FIXDRV/armake: Invalid T-file interval = 10 27 2012 10 27 2012 --need either T-file or one X-file to run FIXDR 修改fixdrv/armake.f 文件中的2010为2099或9999)并在此目录下运行 make -f Makefile 进行重新编译。

8) 运行fixdrv 生成的批处理文件

csh btest1.bat 运行此命令前先把GAMIT.FATAL 和 GAMTI.WARNING 删除

二、 gamit 基线处理结果质量评价指标: (检验q文件)

A)标准均方根 nrms解算结果为 0.25左右,若是 nrms大于 0.5 ,表明处理结果是有问题的(例如:周跳没有修复;测站的起算点坐标有问题等)。若是nrms小于 0.5 ,则认为成功解算, 否则需要检查原因,重新处理;

B) 查看gamit基线解算精度,一般要求相对精度达到10-8~10-9

C) 检查所有测站点是否都参与计算

说明:如果有问题,需要查看

1.查看station.info文件: REC#、 ANT#、 HtCod

IGS跟踪站的天线高都是到DHARP的。

外业观测的需要问清楚当时量取天线的位置。

2.查看sestbl文件

1) 高度角 Elevation Cutoff(不能是0°,按规范来为15°);

2) Choice of Observable 长基线选择LC_help、短基线选择L1\L2 Independent

3) 其他小错误该正项

内容版权声明:除非注明,否则皆为本站原创文章。

转载注明出处:https://www.heiqu.com/15502.html