一、处理前准备
1、在主文件夹内新建test项目文件,项目内新建brdc、igs和rinex三个文件夹,分别存放广播星历,精密星历几观测值文件,所用的命令分别为sh_get_nav、sh_get_orbits和sh_get_rinex (若文件为.Z,用gunzip命令解压,若仍为d,用命令sh_crx2rnx -f 命令解压为o文件)
2、进入test项目文件夹,链接tables,运行
sh_setup -yr 2015 -apr itrf08_com.apr
3、生成station.info文件。
将test/tables下的station.info文件拷贝到rinex文件夹下,打开并编辑,仅保留以#或*开头的前几行,保存并关闭。
打开终端进入rinex文件夹,运行
sh_upd_stnfo -files *.150
4、建立lfile.文件(先验坐标文件)。(非IGS站时,才需此步)
打开终端进入rinex文件夹,用批处理的方式生成lfile.文件,分为三步:
(1) 提取观测值文件的先验坐标:
grep POSITION *.15o > lfile.rnx
(2) 将.rnx文件转化为.apr文件
rx2apr lfile.rnx 2015 001
(3) 由.apr文件生成lfile.文件
gapr_to_l lfile.rnx.apr lfile. “” 2015 001
将station.info和lfile.两个文件复制到test/tables下,覆盖原文件。
5、sestbl.文件的配置。
首先配置 sestbl.这里列出主要需要配置的几项内容:
Choice of Experiment = BASELINE
如果是解算基线的话建议选择 BASELINE,即固定轨道
Choice of Observable = LC_AUTCLN
选择观测值类型,如果是双频接收机并且是大于 5KM 的基线话建
议选取 LC_AUTCLN,如果是小于 5KM 的基线也可以选择 LC_AUTCLN,
其中 LC_HELP 和 LC_AUTCLN 的差别不大,LC_AUTCLN 略好些。
Zenith Delay Estimation = Y
Interval zen = 2
对流层参数的设置,如果是短基线(大概小于 40KM)可以设置
成 N,即不估计对流层天定延迟,采用双差的方法来消除对流层延迟。
如果是大于 40KM 或者 50KM 的基线要设置成 Y。第二项是如果估计对
流层延迟的话,设置多长时间估计一个对流层延迟参数,单位是小时。
Tides applied = 31
Use otl.list = Y
Use otl.grid = Y
设置潮汐改正参数,如果为 31 则考虑海潮改正,为 23 则不加海
潮改正,如果考虑海潮改正的话,需要 tables 下有两个海潮文件,
otl.list 和 otl.grid,其中前者是可有可无的,当两者都采用 Y 时,
程序会根据测站位置首先利用 otl.list(只包括 IGS 站的海潮模型
参数,以及位于该 IGS 站方圆 10KM 范围内都是一套参数)中的海潮
参数,如果该测站没有出现在 otl.list 中,则会根据该测站的位置
和 otl.grid 文件内差出该测站的海潮改正。
以上为 sestbl.中主要的设置项,其他项对基线解算的影响很小,
可以采用默认值。
6、sittbl.文件的配置。
下面设置 sittbl.文件,即设置各个站坐标的约束量,对其算点
坐标采取紧约束,对非起算点的站采用松弛约束。
执行
vi sittbl.
SITE
打开 sittbl.文件。
FIX
–COORD.CONSTR.–
<< default for regional stations >>
ALL
NNN30.
30.
30.
<< IGS core stations >>
BJFS BJFS_GPS NNN 0.050 0.050 0.05
LHAZ LHAZ_GPS NNN 0.050 0.050 0.05
TWTF TWTF_GPS NNN 0.050 0.050 0.05
USUD USUD_GPS NNN 0.050 0.050 0.05
对起算点的 XYZ 坐标分量约束为 5cm,非起算点的站采用松弛约束,
约束量为 30 米。
7、sites.defaults文件的配置。
用来控制需要参与解算的测站,可定义不参与结算的测站。
8、process.default文件的配置。
用来控制处理过程中的细节。
二、理由GAMIT解算基线。
用终端进入test项目文件夹,输入批处理命令:
sh_gamit -expt test -s 2015 001 008 -orbit IGSF -noftp
参数说明:
-d:处理数据的指定日期,-d 2015 001 008 ,处理001和008天的数据
-s:处理数据的指定日期,-s 2015 001 008 ,处理001到008天的数据
结果文件:
qtesta.001:解算记录
otesta.001:解算记录的精简版,一般用此文件
htesta.15001:协方差矩阵、参数平均值
postfit nrms 优于0.3时为最佳,大于1时有问题。
sh_setup -yr -doy -update <links/files/all> -apr -series -localeop
-upd_l -tables