1.1   GAMIT/GLOBK软件可从网络上申请下载。该软件功能强大,用途广泛,一般包括精确定位,大气层可降水汽估计和空间电离层变化分析等。后两种用途只需要用到GAMIT模块,精确定位则还需要GLOBK模块。

1.2    注意:如果你的目标是得到ITRF下的精确测站坐标,为了达到在GLOBK模块中进行坐标转换的目标,需要引入已知ITRF坐标的站点联合解算。即在gamit模块中将周围稳定的,均匀分布的3个或者以上IGS站点引入。如果不注意这点你将前功尽弃!!!!

1.3    在求解坐标的时候,所有GPS处理软件的思想基本上都是一样的,通过基线解算得到测站间的精确相对位置,然后减少其它参数(平差)得到WGS84下的测站坐标(和卫星轨道),最后将坐标转换到某个指定的框架下(后面两步按照大学时老师教TGO的说法就是无约束平差和约束平差)。 


1.4需要掌握最基本的linux命令和语法。例如ln,ls和>>等。 ln即link,类似于windows的快捷方式。lsb-release –a查看linxu版本号,gcc –v查看gcc版本号。具体语法请看终端的命令帮助。

1.5 GAMIT/GLOBK结构。GAMIT/GLOBK安装完成后有几个重要的文件夹:gg/gamit 和gg/kf内小模块是用fortran编写的,分别为基线解算和平差(卡尔曼滤波)两个大功能。gg/com则是用cshell编写的脚本,分别将gamit或kf文件夹内的众多小模块有机的组织起来。gg/tables内包括测站接收机、天线和卫星信息等各种表文件,需要不定期更新。GAMIT/GLOBK是开源软件,允许你自己对软件进行更改,重新编译。

1.6 需要注意的更新。 LUNTAB、SOLTAB、NUTABL、LEAP.SEC的内容每年更新一次。UT1、POLE表每周更新一次。至于软件是否需要更新,请关注如下网址上的说明,包括软件改进和bug修正。如果某个模块中的bug已经严重影响到你的计算结果,请下载并重新编译,如果bug存在于lib文件中,则要求下载升级包后重新编译整个软件。http://www-gpsg.mit.edu/~simon/gtgk/down.htm

1.7运行过程需要的目录。
如果使用sh_gamit批处理,至少需要在工程目录下新建rinex,brdc和gfiles三个目录,分别存放你的观测数据o文件,广播星历n文件, 卫星轨道g文件。这样做的目的是把文件归类,最后这些文件都将被link到单天目录中,基线解算最终是在单天目录里运行的。如果你学习分步处理,目录则比 较自由,最后只要保证一点,所有必要的表文件和观测文件等都要link到单天中。

1.8   部分模型的说明。 
(1).otl潮汐改正比较准确和成熟,一般都会使用。 
(2).有研究认为vmf1映射函数目前是最好的,可比其它映射函数提高7%的重复性(改进的GMF和它近似,两者可选用)。 
(3).atml大气(非潮汐)荷载模型只对高程方向影响较大,可消除大部分周期性跳动,其可靠性还需要进一步证实,因此应用不广。 
(4).atl大气潮汐荷载模型和met气象模型是软件中保留的模型,目前(截至2007年7月)尚未开发。 
(5).现在的GAMIT/GLOBK软件把改正模型中的栅格文件都分离出来,需要到[url=ftp://erverest.mit.edu]ftp://erverest.mit.edu[/url](不需要权限)下载。并做好tables中的各种link。

1.9     与星历有关的四个文件关于brdc文件、sp3文件、g文件和t文件的说明。brdc文件是预报的星历文件,gamit中仅用来拟合卫星钟 差。sp3是事后精密星历文件,每30秒间隔的卫星位置。g文件是根据sp3拟合(积分)的圆形轨道参数文件,t文件是gamit(前三个文件格式是统一通用的,仅t文件是gamit根据自身需要而生成的)根据观测文件和g文件求出的GPS接收机采样时刻的卫星位置。

1.10    数据准备中可能用到的资源地址: 
数据准备中可能用到的资源地址:
各种星历表文件更新在 http://garner.ucsd.edu/pub/gamit/tables/
全球igs分布图 http://igscb.jpl.nasa.gov/network/complete.html
igs站点rinex文件下载 http://garner.ucsd.edu/pub/rinex/
igs站点广播星历下载 http://garner.ucsd.edu/pub/nav/
sp3文件下载 http://garner.ucsd.edu/pub/products/
一个很系统的教程 http://seismo.berkeley.edu/~battag/GAMITwrkshp/wrkshp.html 
步骤: 
1.建立目录mywk/tables 和单天目录 mywk/1998353
2.把rinex格式的观测数据放入单天目录,例如 mhcb3530.98n.Z
解压uncompress mhcb3530.98n.Z
doy 命令换算日期和gps周和此周第几天,
下载精密星历放入单天目录,例如 igs09886.sp3.gz
解压gunzip igs09886.sp3.gz
3.在tables目录下运行links.tables J2000 1998,单天353里运行 links.day 1998 353 lake 
参照.log文件和cvan.dat以及o观测文件 编辑单天里文件station.info的内容,务必准确 ,例如
*SITE Station Name Session Start Session Stop Ant Ht HtCod Ant N Ant E RcvCod SwVer AntCod FARB FARB 1998 353 0 0 30 1998 354 0 0 0 0.0333 DHBCR 0.0000 0.0000 ASHZ12 8.25 ATDMCC HOPB HOPB 1998 353 0 0 30 1998 354 0 0 0 0.0341 DHBCR 0.0000 0.0000 ASHZ12 8.36 ATDMCC MHCB MHCB 1998 353 0 0 30 1998 354 0 0 0 0.0687 DHBCR 0.0000 0.0000 ASHZ12 8.24 ATDMCC
4.用sh_rx2apr -site farb3530.97o -nav farb3530.97n生成先验坐标文件lfile.farb farb.apr 
用sh_rx2apr -site hopb3530.97o -nav hopb3530.97n生成先验坐标文件lfile.hopb hopb.apr 
把各个l文件的坐标录入到l工程文件
5.删掉session.info,rm session.info,运行makexp输入工程名、轨道名,都是lake
6.sh_sp3fit -f igs12481.sp3 -o lake 如果station.info没有指定轨道名,此处需要o 参数 
sh_check_sess -sess 342 -type gfile -file glake3.342 
makej farb3420.03n jfarb3.342 生成卫星钟文件 
sh_check_sess -sess 342 -type jfile -file jfarb3.342 
makex lake.makex.batch 
fixdrv dlake3.342 生成所有准备文件
7.运行批处理csh blake3.bat 
得到O-file、Q-file 文件以及用于后续GLOBK 软件(运用卡尔曼滤波进行全网整体平差的网平差软件) 所需的 H-file 文件
注意:为了达到在下一步GLOBK模块中进行坐标转换的目标,需要引入已知点联合解算。例如全部固定到ITRF坐标,需要在gamit模块中将周围稳定的,均匀分布的3个或者以上IGS站点引入联合解算。
用GLOBK模块进行平差和坐标转换,对于区域网分析,不引入全球解的情况 
8.建立目录glbf和soln(均与gamit中的单天目录平行), 用htoglb将h解转换成二进制格式,用ls 生成gdl文件,准备好需要的apr文件,globk命令文件和glorg命令文件,并做好命令文件的相关设置。 
注意: globk命令进行平差,最后得到每天的坐标.prt文件,上述坐标结果是在WGS84下的坐标,其内部几何形很好,但是每天坐标框架的浮动都很大(cm 级的浮动范围下即可求出整周模糊度),需要把每天的解都统一到同一框架中,这就是引入IGS站点的作用。IGS站点在ITRF中的坐标是通过各种技术精确 计算出来的,由某个历元的坐标和运动速度组成。例如在ITRF00框架中,某个IGS站的坐标=给定初始坐标+速度×时间差 glorg命令可以通过IGS站计算出单天解对应的7个(或者3个)转换参数,并且利用转换参数将其他待求点的ITRF坐标求出来。

9.运行glred命令,它会自动激发globk和glorg。前者命令进行平差,得到每天的坐标.prt文件。后者得到ITRF坐标.org文件。 
10.现在,用plotcrd把时间序列输出到屏幕吧,发现有粗差后剔除它,或者返回到gamit重新解算这天的基线。 
11.如果你要的是周解或其他平均解,请在剔除粗差后运行globk命令将单天解联合解算,后续步骤类似上面介绍。

sh_gamit 批处理 
1.如不是IGS站,将处理的站信息加入到软件的station.info中,将作为参考转换站的ITRF2000坐标和速度加入到itrf00.apr中 
2.建立 工程名,其下建立子目录rinex ,brdc ,gfiles.分别放入(或者链接,下同)rinex文件(文件名中站名小写字母),广播星历文件(前四个字符必须是brdc,而不是auto等其它),g文件。 sh_gamit提供了完善的寻找功能,或者工程名下还有sp3文件夹,当找不到某天g文件时,再到sp3里去找,再找不到则从网络上下载,最后下载不成功则提示错误并退出。 
3.在工程名下运行sh_setup yr 2007 得到一个目录tables. 现在目录为rinex ,brdc, gfiles, tables. 
编辑tables内的sites.default,所有站均xstinfo,不再用rinex来更新station.info. 在GLOBK中用到(used)的站glrepu,为参考转换(stabilize)的站设置glreps。 
编辑sittbl.将参考站紧约束0.005 0.005 0.01,区域站松约束100 100 100 
process.default,sestbl.一般不用编辑。 
4.运行例如sh_gamit –s 2007 88 102 -orbt lake -expt expt 即可得到2007年88-102天的h单天解。
sh_glred 批处理检验重复性和求平均解
1.如果是没有联合全球h文件,则需要自己建立目录gsoln(与单天平行),构建globk_comb.cmd  和glorg_comb.cmd
2.在工程名下运行命令即可,例如检验重复性sh_glred -s 2007 110 2007 130 -expt haha -opt H G E >!rep.log会得到一个org结果文件和时间序列图。
3.求平均解 sh_glred -s 2007 110 2007 130 -comb 21 -expt haha -opt H G >!rep.log注意这条命令增加了comb联合多少天的解,没有了E,因为平均解只有一个,不存在画时间序列的说法

05-11 11:09