GAMIT-GLOBK数据处理报告
GAMIT-GLOBK数据处理报告
一．处理任务
利用GAMIT-GLOBK软件对2011年年积日为94天的shao, lhaz, xian, kunm, bjfs, urum共6个IGS测站的GPS测量数据进行处理，并对处理结果进行评估。
二．处理步骤
-
安装虚拟机和Linux系统
在win7系统下安装VMware Workstation v9.0.1虚拟机，并在其上安装ubuntu 12.10。
安装完并登陆ubuntu后，在terminal中运用以下指令安装csh:
sudo apt –get install csh
安装完成后在terminal执行：
chsh
按提示输入：
/bin/csh
修改bash为csh。重启ubuntu系统，打开terminal，输入：
sudo –sH
输入密码后进入root权限，然后输入：
echo $SHELL，
返回
/bin/csh，说明SHELL修改成功。
安装库文件libX11 –dev和编译器gfortran:
sudo apt –get install libX11 –dev
sudo apt –get install gfortran
-
安装GAMIT-GLOBK软件
在/opt目录下建立gamit目录作为软件的安装目录，将下载的GAMIT-GLOBK+10.40.zip安装包文件复制到/opt/gamit目录下。打开terminal，进入gamit文件夹下，输入命令：
unzip GAMIT-GLOBK+10.40.zip
解压得到source文件夹，在terminal中进入source文件夹，运行：
./install_software
如果此时提示权限不足，则输入
chmod 777 *
为所有用户和文件设置可读写权利，然后运行./install_software命令，当到达提示
时，先进入source/libraries文件夹下用gedit打开Makefile.config，找到第283行：
# ----- for Linux from 0.0.1 to 3.0.0 -- #
OS_ID Linux 0001 3000
将其改为
# ----- for Linux from 0.0.1 to 4.0.0 -- #
OS_ID Linux 0001 4000
这是由于Linux内核版本已更新至3以上；另外在此文件的105行将
MAXSIT 60
MAXSAT 32
MAXATM 13
MAXEPC 2880
的MAXATM修改为25，MAXEPC修改为6000，关闭并保存文件后在terminal里选y继续安装，即可完成GAMIT-GLOBK软件的安装。
安装完成后，需要修改csh.cshrc文件，如下图提示所示：
进入/home/etc文件夹下，用gedit csh.cshrc打开，在endif前加入
stenv PATH /opt/gg/gamit/source/gamit/bin: /opt/gg/gamit/source/gamit/com:
/opt/gg/gamit/source/gamit/kf/bin:$PATH
stenv HELP_DIR /opt/gg/gamit/source/help/
stenv INSTITUTE GNSS CENTRER
保存并关闭文件后退出，在terminal中输入doy指令，若出现结果则表明GAMIT-GLOBK已成功安装。
-
工程文件和待处理数据准备
在任意目录下建立工程文件夹x094，文件名为四字小写，里面建立brdc，igs，rinex和年积日94四个文件夹：
需要下载的数据有rinex观测值文件，广播星历和sp3星历，分别放在rinex，brdc和igs文件夹下，可以用sh_get_rinex，sh_get_nav和sh_get_orbits自动下载，也可以通过网站http://cddis.nasa.gov/gnss_datasum.html#daily网站手动下载，本例数据从网上下载，准备好的工程文件内容如下（tree指令可通过sudo apt-get install tree安装）：
-
tables文件夹更新和准备文件制作
需要每年更新的tables文件有：
章动表nutabl.、太阳表soltal.、月亮表luntab.、跳秒表leap.sec
需要每月更新的文件有：
编码文件dcb.dat
需要每日更新的文件有：
地球自转参数ut1.、pole
在有新的接收机或者天线出现时需要更新rcvant.dat，在有新的卫星发射时需要更新svnav.dat。
这些数据可以在ftp://garner.ucsd.edu/pub/gamit/tables/下载，在所有更新完成后，在工程目录x094下运行sh_setup -yr 2011将tables表链接到工程目录中。
需要准备的文件有：
测站信息文件station.info
将tables文件夹里的station.info文件复制到rinex目录下，打开并编辑该文件，只保留以*号开始的行信息（文件头），保存并关闭该文件。使用terminal进入rinex目录，输入
sh_upd_stnfo –files *.11o
就可以生成station.info文件：
近似坐标文件lfiles.
近似坐标从观测o文件中提取，terminal进入rinex目录，输入
grep POSITION *.11o >lfiles.rnx
抓取o文件中的近似坐标值到lfiles.rnx文件中，得到测站近似坐标后，采用命令
rx2apr lfiles.rnx 2011 094
得到apr格式的文件lfiles.rnx.apr，然后通过指令
gapr_to_l lfiles.rnx.apr lfiles. 2011 094
将apr文件转换为lfiles.：
将station.info和lfiles.文件复制到工程文件夹x094下的tables目录下。
进入工程文件夹x094下的tables目录，目录控制文件process.defaults、测站列表文件sites.defaults、解算控制文件sestbl.、精度控制文件sittbl.均保持默认值。在本次处理过程中，不考虑海洋潮汐的影响，可以禁用海洋潮汐模型otl.grid文件，在sestbl.文件中设置为：Tides applied = 23; use otl.grid = N，保存并关闭sestbl.文件
-
数据处理
进入工程目录x094下，采用批处理命令
sh_gamit –d 2011 094 –expt x094 –noftp
处理过程截图如下：
-
处理结果及精度评定
最终处理结果保存在qx094a.094文件中。
运用globk进行点位平差处理，首先terminal进入工程目录x094下，运行
sh_glred –d 2011 94 –expt x094 scal –opt H G
对基线解算数据进行预处理，然后进入gsoln目录，运行
ls ../glbf/h*glx>x094.gdl
globk 6 globk_comb.prt globk_com.log x094.gdl globk_comb.cmd
处理截图如下：