gamit中的sh_gamit命令批量解算基线

注意： 文件名只能是4个字母或加数字组合

1、只有2017连续站IGS解算

准备工作

mkidr brdc igs rinex

广播导航文件

• 生成N文件

cd brdc
sh_get_nav
sh_get_nav -yr 2017 -doy 010 -ndays 1

精密星历

• 生成SP3文件，而GMIT自动轨道积分生成G文件

cd igs
sh_get_orbits
sh_get_orbits -orbit igsf  -yr 2017 -doy 010  

观测数据

cd ../rinex/
sh_get_rinex -archive sopac cddis  -yr 2017 -doy 010 -ndays 1 -sites urum lhaz shao 

2、链接TABLE表，即station.info制作

sh_setup -yr 2017

• 红色死链接。
• 蓝色可用链接。
• 查看文件

ls -al antmod.dat
ls -al /home/demon/gg/tables/antmod.dat

• 查看解算策略

vim sestbl. 

• 红色改为N
  蓝绿色改为Y
• 制作自己的 station.info
  • 把 station.info 中第一行到第六行的头文件生成到rinenx中去

sed -n '2,6p' station.info
sed -n '2,6p' station.info > ../rinex/station.info

sh_upd_stnfo
sh_upd_stnfo -files *.17o 

• 处理完成后，确认天线类型、天线高、天线高量取方式、接收机类型
• 把station.info 拷贝到tables中
  cp station.info …/tables

3、近似坐标的制作

• 只使用头文件中的tables表中的近似坐标
  cd …/tables/
• 打开处理过程的文件

vim process.defaults
把N改为Y

#Set "Y" to use RINEX header coordinates not in lfile or apr file
 set use_rxc      = "Y"

4 在工程目录下面，用SH批处理命令自动进行GMIT解算,如果出现错误，就进行步骤5,如果依旧错误，那就是表文件没更新。

sh_gamit
批处理
sh_gamit -expt demo -d 2017 010 -orbit IGSF -pres ELEV -noftp -dopt D ao x c > sh_gamit.log

5、更新L文件

在table表目录下执行

grep "APPROX POSITION XYZ"$ ../rinex/*o > lfile.rnx 
# 提取观测文件中的概略坐标

#上面这一步莫似对hkkt和hkks站用不上，暂时不知道啥原因，目前先采用手动添加的方式，打开o文件，找到“APPROX POSITION XYZ”这一行，手动将hkkt和hkks站点近似坐标添加在由上一条命令生成的lfile.rnx内，完成此步骤后继续看下面

rx2apr lfile.rnx 2017 010  # 生成 lfile.rnx.apr 文件

gapr_to_l lfile.rnx.apr lfile.2017 010  # 生成 lfile. 文件

6、更新表文件。

7、解算完成之后，因为解算了两次，会生成

 ls q*
qdemoa.010  qdemop.010
有约束       无约束

 ls o*
odemoa.010  odemop.010

ls hdemoa.17010 
hdemoa.17010