Gamit/GLOBK下载
Gamit/GLOBK数据解算流程
数据准备
在相应的网站下载好观测文件o文件、广播星历n文件、精密星历sp3文件以及电离层文件,可以直接在虚拟机里下载。
- o文件、n文件、sp3文件下载:武汉大学IGS数据中心:http://www.igs.gnsswhu.cn
- 电离层文件下载:https://cddis.nasa.gov
下载后,在虚拟机桌面上新建工程文件夹test(注意不要有数字,最好只有英文),在test文件夹下新建四个文件夹rinex、igs、brdc、ionex。上述操作命令为:
mkdir test
cd test
mkdir rinex igs brdc ionex
然后将o文件移到rinex内,n文件移到brdc内,sp3文件移到igs内,电离层文件移到ionex内。
然后在工程文件夹内使用命令链接tables表文件:
sh_setup -yr **** //****为自己解算的的年份
这样数据就准备好了。
参数设置以及表文件制作
station.info
在工程文件夹下依次输入:
cd tables
sed -n '2,6p' station.info > ../rinex/station.info
cd ../rinex/
ls
查看rinex文件夹下是否有station.info,输入:
sh_upd_stnfo -files *.21o //21根据自己的解算年份修改
cp station.info ../tables/
cd ..
sestbl.
sittbl.
sites.defaults
process.defaults
基线解算
在工程文件夹下,输入:
sh_gamit -expt test -s 2021 124 130 -orbit IGSF -prse ELEV
其中,-expt后为你的工程文件名;-s表示连续解算多天,后面为解算年份、起算日、终止日;-orbit IGSF表示使用最终星历;-pres ELEV表示绘制天空图和相位残差图。
解算完成后,会在工程文件夹下生成以解算年积日命名的文件夹,主要查看o*a.*、h*a.*文件。在figs文件夹下可以查看天空图和相位残差图。
输入如下命令提取所有解算结果的RMS值和模糊度:
grep nrms: ???/o*a.*
grep fixed ???/*.summary
输入如下命令提取重复基线:
grep -h "X N" 1??/o*a.* | sort>bl.inp
cat bl.inp | wbslfilt>bl.lst
网平差
在工程文件夹下输入:
sh_glred -cmd
会在 gsoln文件夹下生成glorg.cmd、globk.cmd两个文件,若没有生成,直接进入gamit/tables下复制这两个文件到gsoln文件夹内。
找到 glorg.cmd 文件, 找到如下内容
# This should match the well-determined sites in the apr_file
stab_site clear
在下面添加
stab_site bjfs chan daej hksl jfng hkws(测站名自己修改)
输入命令进行平差:
sh_glred -expt work -d 2019 001 -opt H G E(一天解算)
sh_glred -expt work -s 2019 001 2019 365 -opt H G E(多天连续解算)
解算完成后在gsoln中生成若干个*.org文件,即可查看解算结果。