在Linux操作系统下安装与运行GAMIT软件(2)
4、运行实例
本文选取了2006年,年积日为100天的安徽黄山(ahhs)、安徽马鞍山(ahma)、江苏常州(jscz)、江苏高邮(jsgy)、上海奉先(shfx)、上海浦东(shpd)6个国内的GPS观测站的观测资料。
4.1 数据准备
建立rinex、eph、tables目录和以年积日100命名的工作目录。在rienx目录中准备该天的所选用的各个观测站的观测文件o-files。在eph目录中准备导航文件auto1000.06n,GPS卫星星历igu13701.sp3,极移表pole.,TAI-UT1国际时间系统表ut1.,坏卫星信息文件svs_exclude.txt。在tables中建立:测站信息文件station.info,测段信息控制文件sestbl., 测站信息控制文件sittbl.,天线相位中心改正表antmod,dat,周跳的自动探测和修改命令表autcln.cmd,地球形状参数表gdetic.dat,跳秒表leap.sec,月亮表luntab.,章动表nutabl.,太阳表soltab.,星号对照表svnav.dat,接收机及其天线型号对照表rcvan.dat。
上述的文件,可以在网站(ftp://lox.ucsd.edu)下载2006年底100天的卫星星历,以及2005年的太阳表、月亮表和章动表(这些表每年都要更新,因此要下载所处理数据当年的),其余的文件在安装GAMIT时是自带的,复制到对应的文件夹下即可。需要注意的是, station.info, sestbl., sittbl.这三个文件,应该根据实际处理的情况做相应的修改。
编辑station.info文件:首先打开o-files,记录每个选用台站的X、Y、Z方向的高度和接收机型号、天线型号,然后来校正station.info文件中的数据,以求每一个站都有正确的天线高度等信息。开始和结束时间形如:0 0 0 0 24 0 0,则天线信息不仅适用于这一天,而且适用于以后的天,直到有更新的信息。同一站的信息可以出现多次,但不必相邻,但必须是按时间先后顺序排列的。
编辑sestbl.文件:文件选择计算方案和设置与计算方案相应的参数。主要设置下面几项:
A、Choice of Experiment
RELAX:定轨、定位、解ERP;
BASELINE:仅仅定位。
B、Choice of Observable
LC_HELP:用LC观测解模糊度;
LC_RANGE:用LC观测模糊度,但更强调伪距的作用;
LC_ONLY:用LC观测,不解模糊度;
L1_ONLY:仅仅使用L1,解模糊度,对于几公里的小网;
L2_ONLY:仅仅使用L2,解模糊度,对于几公里的小网;
C、Zenith Delay Estimation
YES:解算天顶延迟估计;
NO:不解算天顶延迟估计。
D、Interval Zen
解算天顶延迟参数的时间间隔。
编辑sittbl.文件:首先检查是否包括每个选用测站的信息,缺少的信息要补全,然后对每一个站设置先验约束。这里主要是对台站在X、Y、Z方向摆动的作限制。
4.2 数据处理步骤
(1)用ln命令链接../rinex,../eph,../tables内的文件到work;
(2)执行makexp程序建立所有准备文件的输出及一些模块的输入文件。系统会提示让输入试验名、轨道名、year、doy、session number、概略坐标文件、导航文件名,系统还会提示输入采样间隔、起始时间、历元数。这样,我们依次输入pgga,pgga,2006,100,0,lpgga,auto1000.06n,60 0 0 1440即可。
(3)sh_sp3fit –-f igu13701_18.sp3 –-d 2006 099 100 –-t tpgga6.100
并将生成的*.099文件改名为*.100
(4)makej auto1000.06n jauto6.100
(5)makex pgga.
(6)fixdrv dpgga6.100
(7)csh b*.bat
以上pgga为工程名称,可以根据需要进行替换。
4.3 数据处理结果分析
此次数据处理例子是解算对流层天顶延迟的。结算的精度与可靠性主要看两个方面:
(1)是否使用了足够的资料。所得基线分量的精度是衡量这个标准是否满足的标志。此例中基线边解算精度的统计如表1所示。
由表1可以看到,GAMIT解算得到的基线相对精度均优于对台站坐标和卫星轨道所加的约束10-3。
(2)对资料的拟合模型是否达到它的噪声水平。满足这个条件的判据是解的nrms(normalized rms),即每个自由度的chi-square的平方根。其理论值是1,但实际上这个值在0.25左右才算正常,本例中nrms=0.28。如果这个值大于0.5,则意味着有残余的周跳,或是有严重的其它问题。
表 1
基线边 |
基线长(m) |
精度(m) |
相对精度 |
AHHS- AHMA |
221380.31016 |
0.00386 |
1.74×10-8 |
AHHS- JSCZ |
288514.33161 |
0.00344 |
1.19×10-8 |
AHHS- JSGY |
358731.77533 |
0.00354 |
0.99×10-8 |
AHHS- SHFX |
335520.06794 |
0.00456 |
1.36×10-8 |
AHHS- SHPD |
355131.54150 |
0.00496 |
1.40×10-8 |
AHMA- JSCZ |
139674.63390 |
0.00490 |
3.51×10-8 |
AHMA- JSGY |
149475.83358 |
0.00417 |
2.79×10-8 |
AHMA- SHFX |
294612.54009 |
0.00580 |
1.97×10-8 |
AHMA- SHPD |
293058.38336 |
0.00645 |
2.20×10-8 |
JSCZ- JSGY |
114012.99629 |
0.00343 |
3.01×10-8 |
JSCZ- SHFX |
176465.95374 |
0.00508 |
2.88×10-8 |
JSCZ- SHPD |
165474.09440 |
0.00610 |
3.69×10-8 |
JSGY- SHFX |
281897.78401 |
0.00498 |
1.77×10-8 |
JSGY- SHPD |
264050.84049 |
0.00603 |
2.28×10-8 |
SHFX- SHPD |
32751.45431 |
0.00630 |
1.92×10-7 |
5、结论
读者可以根据以上的步骤进行GAMIT软件的安装与运行,不同的硬件配置,处理结果可能有微小的差异。从以上可以得到结论,GAMIT在LINUX系统上的处理结果是可靠的。然而,GAMIT为一复杂的GPS处理软件,在实际使用时,会遇到各种各样的问题,这就需要在实践中不断地摸索总结,这样,才能熟练掌握它。