记录各测站的站点、站点名称、开始观测时间、结束观测时间、天线高、 天线量测方法、接收机类型、软硬件版本号、天线类型等信息。严格按照已经给出的跟踪站格式编辑。(注意:要现在当前目录下拷入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) 其他小错误该正项