论文部分内容阅读
[摘 要]GAMIT是美国麻省理工学院( MIT)与斯克里普斯海洋研究所( SIO)联合研制的一款GPS处理软件,具有处理速度快、精度高、开源的优势。文章主要介绍了该软件的安装方法并结合全球IGS跟踪站的观测数据来进行基线平差,并对实验结果进行评定。
[关键词]GAMIT;跟踪站;GPS基线;均方差
中图分类号:P228 文献标识号:A 文章编号:2306-1499(2014)08-0179-02
1.引言
随着GPS 技术在大地测量、精密工程测量、地壳运动监测、地球动力学等领域的广泛应用, GPS数据精密处理技术越来越受到重视和发展。而国内或者一些商业软件并不能满足精度的需求。如今,国内外研制出了很多优秀的GPS 数据精密处理软件, 其中包括瑞士伯尔尼大学天文研究所研制的BERNESE 软件, 美国德克萨斯大学的TEXGAP软件,英国的GAS软件等, 以及本文将要介绍的GAMIT 软件。
2.GAMIT软件功能简介
GAMlI软件主要由以下几个模块组成:ARC(轨道积分)、MODEL(组成观测方程)、SINCLN(单差自动修复周跳)、DBLCLN(双差自动修复周跳)、CVIEW(人工交互式修复周跳)、SOLVE(利用双差观测按最小二乘法求解参数)、FXDRV ( 生成数据处理的程序) 、GLOBK( 运用卡尔曼滤波进行网平差模块) 等。GAMIT软件处理双差观测量.采用最小二乘算法进行参数估计,采用双差观测量的优点是可以完全消除卫星钟差和接收机钟差的影响,同时也可以明显减弱诸如轨道误差、大气折射误差等系统性误差的影响。
3.GAMIT软件的安装准备
3.1GAMIT软件所需编译器的安装
在安装GAMIT之前需要对软件使用环境进行配置,需要安装以下软件,gcc-4.4,gfortran-4.4,libx11-dev,csh及tcsh,gmt,gftp,teqc。在安装的过程中需要在root环境下,在登陆的时候选择root用户再进行安装,在终端里直接使用命令apt-get install gcc-4.4,apt-get install libx11-dev,apt-get install csh,apt-get install tcsh,apt-get install gmt,apt-get install gftp,apt-get install gfortran-4.4,以上安装不分先后顺序,系统会自动下载所需的安装包自动安装。
3.2GAMIT软件的安装
从麻省理工学院的服务器(ftp://chandler.mit.edu)下的/updates/source目录下,下载GAMIT安装所需软件包,其中包括安装批处理文件install_ software、更新批处理文件install_updates和几个以.tar或.tar.gz结尾的压缩文件。
在/opt目录下创建gamit文件夹,将以上下载的压缩包全部拷入gamit文件夹下,在终端中定位到该目录下,执行命令./ installsoftware,按照提示就可以进行自动安装,在回答2个Y之后,对安装程序的配置文件Makefile.config进行修改,将其中的3000改为你第一次报错是提示的数值,不同的安装环境会提示不同的数值。
3.4GAMIT软件路径配置
上述安装过程完成之后,GAMIT软件并不能正常运行,需对.bashrc文件进行配置。将路径/opt/gamit/com、/opt/gamit/gamit/bin、/opt/gamit/ kf/bin,加入到ubuntu的PATH中。打开.bashrc文件在其最后面添加如下:
#set gamitpath
gamitpath=/opt/gamit
export PATH=$PATH:${gamitpath}/com/:${gamitpath}/gamit/ bin:${gamitpath}/kf/bin
HELP_DIR=${gamitpath}/help/
export HELP_DIR
############GMT##############
#set GMT path
gmtpath=/user/lib/gmt
PATH=$PATH:${gamitpath}/bin
配置完后,在root用户下的安装与配置便完成了,可以再终端中输入doy命令来测试安装是否成功,如果输入没有任何反应则需要重新安装。在以上所有步骤都完成后还需要再进行修改opt/gamit/gamit/fixdrv文件夹下的armake.f,将2010改为9999再打开端口运行make -f Makefile命令,否则会出现2010年之后的数据无法处理的現象。
4.GAMIT处理数据的流程
首先在/home目录下创建工作目录40,以2013年第40年积日为例进行数据处理,
在该目录下使用sh_setup -yr 2013链接tables文件夹,但在此之前要更新tables文件夹下的文件。
4.1准备工作:
下载所需的精密星历sp3文件、导航文件n文件和观测站的o文件,解压并拷贝到tables目录下。手动修改这些文件:测站点近似坐标文件lfile.,测站文件sites.default需要将所测的站点名称按一定格式加入进去,采用4个字节,将工程名称“expt”替换成ch13,从o文件来更新测站信息文件station.info,设置测站信息控制文件sestbl.,基线解算时,处理模式选用松驰解RELAX,电离层折射影响选用LC_HELP 观测值消除, 对流层折射根据标准大气模型选用萨斯坦莫宁( Saastamoinen)模型改正,设置sittbl.(测站精度控制约束控制文件)。 4.2单天处理方式
(1)执行makexp命令时根据提示输入相关信息,包括: 项目名、轨道(接上页)名、处理数据的年份、天数 、session number、近似坐标文件、导航文件名, 采样间隔、起始时间、最大历元数。这个过程将会产生“session. info”文件和“d”文件,需要注意的是如果处理的是10年之前的数据,轨道名输入igs0,;如果处理的数据是10年之后的轨道名则输入igs1,否则会出来生成的**.bat文件不能将所有的文件串联起来。
(2)执行sh_sp3fit -f igs17266.sp3 -d 2013 40 –t会产生“g”文件和“t”文件。
(3)makej auto0400.13n jauto3.040会产生卫星钟差文件。
(4) makex ch13.makex.batch会产生接收机时钟文件和观测文件。
(5)fixdrv dch133.040会产生批处理文件以.bat结尾的文件b文件。ARC,MODEL,AUTCLN,CFMRG和SOLVE这些处理模块将会自动连接起来。
(6) csh bch133.bat该过程无需人与干涉,会自动完成解算。AUTCLN程序可对相位观测的周跳(CYCLESUP)和粗差(OUTUER)进行自动剔除,SOLVE程序完成最小二乘法的解析,并将结果输出到生成的q文件中。
4.3连续多天的批量处理方式
(1)需要在工程目录下建立igs、brdc、rinex三个文件夹,分别存放精密星历文件、导航文件(要把导航文件名称的auto改为brdc)、观测文件(在解压的时候使用sh_crx2rnx -c n -d y -f *.Z),链接tables文件夹。
(2)在igs文件夾下需要链接以下文件,命令如下:
ln -s ../tables/nutabl. ln -s ../tables/ut1. ln -s ../tables/pole.
ln -s ../tables/leap.sec ln -s ../tables/svnav.dat ln -s ../tables/luntab. ln -s ../ tables/soltab.
(3)最后使用sh_gamit -s 2013 10 12 -expt pgga -orbit IGSF -noftp
4.4实例应用
通过对全球跟踪站的BJFS(北京房山)、CHAN(长春)、WUHN(武汉)、SHAO(上海)的2013年第40个年积日的观测数据为例进行解算,结果如下:
GAMIT解算的评估因子是均方差nrms,本例中的nrms=0.20,解算质量可靠。
5.结论
GAMIT对文件的位置和命名的要求很严格,在实际操作中要注意这些,数据处理结果表明, GAMIT 软件所获得的基线解是高精度的, 对于超长基线的相对精度能达到10-9。可以满足精密工程的应用以及大地板块监测等方面。
参考文献
[1]T.A.Herring, R.W.King, S.C.McClusky. GAMIT Reference Manual release 10(5).
[2]李征航,黄劲松.GPS测量与数据处理[M].武汉:武汉大学出版社,2004.
[3]T.A.Herring;R.W.King;S.C.McClusky GAMIT Reference Manual release 10(4).
[4]梁光民,王隆杰.LINUX操作系统实用教程,2004.
[关键词]GAMIT;跟踪站;GPS基线;均方差
中图分类号:P228 文献标识号:A 文章编号:2306-1499(2014)08-0179-02
1.引言
随着GPS 技术在大地测量、精密工程测量、地壳运动监测、地球动力学等领域的广泛应用, GPS数据精密处理技术越来越受到重视和发展。而国内或者一些商业软件并不能满足精度的需求。如今,国内外研制出了很多优秀的GPS 数据精密处理软件, 其中包括瑞士伯尔尼大学天文研究所研制的BERNESE 软件, 美国德克萨斯大学的TEXGAP软件,英国的GAS软件等, 以及本文将要介绍的GAMIT 软件。
2.GAMIT软件功能简介
GAMlI软件主要由以下几个模块组成:ARC(轨道积分)、MODEL(组成观测方程)、SINCLN(单差自动修复周跳)、DBLCLN(双差自动修复周跳)、CVIEW(人工交互式修复周跳)、SOLVE(利用双差观测按最小二乘法求解参数)、FXDRV ( 生成数据处理的程序) 、GLOBK( 运用卡尔曼滤波进行网平差模块) 等。GAMIT软件处理双差观测量.采用最小二乘算法进行参数估计,采用双差观测量的优点是可以完全消除卫星钟差和接收机钟差的影响,同时也可以明显减弱诸如轨道误差、大气折射误差等系统性误差的影响。
3.GAMIT软件的安装准备
3.1GAMIT软件所需编译器的安装
在安装GAMIT之前需要对软件使用环境进行配置,需要安装以下软件,gcc-4.4,gfortran-4.4,libx11-dev,csh及tcsh,gmt,gftp,teqc。在安装的过程中需要在root环境下,在登陆的时候选择root用户再进行安装,在终端里直接使用命令apt-get install gcc-4.4,apt-get install libx11-dev,apt-get install csh,apt-get install tcsh,apt-get install gmt,apt-get install gftp,apt-get install gfortran-4.4,以上安装不分先后顺序,系统会自动下载所需的安装包自动安装。
3.2GAMIT软件的安装
从麻省理工学院的服务器(ftp://chandler.mit.edu)下的/updates/source目录下,下载GAMIT安装所需软件包,其中包括安装批处理文件install_ software、更新批处理文件install_updates和几个以.tar或.tar.gz结尾的压缩文件。
在/opt目录下创建gamit文件夹,将以上下载的压缩包全部拷入gamit文件夹下,在终端中定位到该目录下,执行命令./ installsoftware,按照提示就可以进行自动安装,在回答2个Y之后,对安装程序的配置文件Makefile.config进行修改,将其中的3000改为你第一次报错是提示的数值,不同的安装环境会提示不同的数值。
3.4GAMIT软件路径配置
上述安装过程完成之后,GAMIT软件并不能正常运行,需对.bashrc文件进行配置。将路径/opt/gamit/com、/opt/gamit/gamit/bin、/opt/gamit/ kf/bin,加入到ubuntu的PATH中。打开.bashrc文件在其最后面添加如下:
#set gamitpath
gamitpath=/opt/gamit
export PATH=$PATH:${gamitpath}/com/:${gamitpath}/gamit/ bin:${gamitpath}/kf/bin
HELP_DIR=${gamitpath}/help/
export HELP_DIR
############GMT##############
#set GMT path
gmtpath=/user/lib/gmt
PATH=$PATH:${gamitpath}/bin
配置完后,在root用户下的安装与配置便完成了,可以再终端中输入doy命令来测试安装是否成功,如果输入没有任何反应则需要重新安装。在以上所有步骤都完成后还需要再进行修改opt/gamit/gamit/fixdrv文件夹下的armake.f,将2010改为9999再打开端口运行make -f Makefile命令,否则会出现2010年之后的数据无法处理的現象。
4.GAMIT处理数据的流程
首先在/home目录下创建工作目录40,以2013年第40年积日为例进行数据处理,
在该目录下使用sh_setup -yr 2013链接tables文件夹,但在此之前要更新tables文件夹下的文件。
4.1准备工作:
下载所需的精密星历sp3文件、导航文件n文件和观测站的o文件,解压并拷贝到tables目录下。手动修改这些文件:测站点近似坐标文件lfile.,测站文件sites.default需要将所测的站点名称按一定格式加入进去,采用4个字节,将工程名称“expt”替换成ch13,从o文件来更新测站信息文件station.info,设置测站信息控制文件sestbl.,基线解算时,处理模式选用松驰解RELAX,电离层折射影响选用LC_HELP 观测值消除, 对流层折射根据标准大气模型选用萨斯坦莫宁( Saastamoinen)模型改正,设置sittbl.(测站精度控制约束控制文件)。 4.2单天处理方式
(1)执行makexp命令时根据提示输入相关信息,包括: 项目名、轨道(接上页)名、处理数据的年份、天数 、session number、近似坐标文件、导航文件名, 采样间隔、起始时间、最大历元数。这个过程将会产生“session. info”文件和“d”文件,需要注意的是如果处理的是10年之前的数据,轨道名输入igs0,;如果处理的数据是10年之后的轨道名则输入igs1,否则会出来生成的**.bat文件不能将所有的文件串联起来。
(2)执行sh_sp3fit -f igs17266.sp3 -d 2013 40 –t会产生“g”文件和“t”文件。
(3)makej auto0400.13n jauto3.040会产生卫星钟差文件。
(4) makex ch13.makex.batch会产生接收机时钟文件和观测文件。
(5)fixdrv dch133.040会产生批处理文件以.bat结尾的文件b文件。ARC,MODEL,AUTCLN,CFMRG和SOLVE这些处理模块将会自动连接起来。
(6) csh bch133.bat该过程无需人与干涉,会自动完成解算。AUTCLN程序可对相位观测的周跳(CYCLESUP)和粗差(OUTUER)进行自动剔除,SOLVE程序完成最小二乘法的解析,并将结果输出到生成的q文件中。
4.3连续多天的批量处理方式
(1)需要在工程目录下建立igs、brdc、rinex三个文件夹,分别存放精密星历文件、导航文件(要把导航文件名称的auto改为brdc)、观测文件(在解压的时候使用sh_crx2rnx -c n -d y -f *.Z),链接tables文件夹。
(2)在igs文件夾下需要链接以下文件,命令如下:
ln -s ../tables/nutabl. ln -s ../tables/ut1. ln -s ../tables/pole.
ln -s ../tables/leap.sec ln -s ../tables/svnav.dat ln -s ../tables/luntab. ln -s ../ tables/soltab.
(3)最后使用sh_gamit -s 2013 10 12 -expt pgga -orbit IGSF -noftp
4.4实例应用
通过对全球跟踪站的BJFS(北京房山)、CHAN(长春)、WUHN(武汉)、SHAO(上海)的2013年第40个年积日的观测数据为例进行解算,结果如下:
GAMIT解算的评估因子是均方差nrms,本例中的nrms=0.20,解算质量可靠。
5.结论
GAMIT对文件的位置和命名的要求很严格,在实际操作中要注意这些,数据处理结果表明, GAMIT 软件所获得的基线解是高精度的, 对于超长基线的相对精度能达到10-9。可以满足精密工程的应用以及大地板块监测等方面。
参考文献
[1]T.A.Herring, R.W.King, S.C.McClusky. GAMIT Reference Manual release 10(5).
[2]李征航,黄劲松.GPS测量与数据处理[M].武汉:武汉大学出版社,2004.
[3]T.A.Herring;R.W.King;S.C.McClusky GAMIT Reference Manual release 10(4).
[4]梁光民,王隆杰.LINUX操作系统实用教程,2004.