《对GAMIT GLOBK的基本认识》由会员分享，可在线阅读，更多相关《对GAMIT GLOBK的基本认识（3页珍藏版）》请在金锄头文库上搜索。
1、1.1 GAMIT/GLOBK 软件可从网络上申请下载。该软件功能强大，用途广泛，一般包括精确定位，大气层可降水汽估计和空间电离层变 化分析等。后两种用途只需要用到GAMIT模块，精确定位则还需要GLOBK模块。1.2 注意：如果你的目标是得到 ITRF 下的精确测站坐标，为了达到在 GLOBK 模块中进行坐标转换的目标,需要引入已知 ITRF 坐标的 站点联合解算。即在gamit模块中将周围稳定的，均匀分布的3个或者以上IGS站点引入。如果不注意这点你将前功尽弃！ ！1.3在求解坐标的时候，所有GPS处理软件的思想基本上都是一样的，通过基线解算得到测站间的精确相对位置，然后减少其它参数（平
2、差）得到WGS84下的测站坐标（和卫星轨道），最后将坐标转换到某个指定的框架下（后面两步按照大学时老师教TGO的说法就是无 约束平差和约束平差）。1.4需要掌握最基本的linux命令和语法。例如In, Is和等。In即link，类似于windows的快捷方式。lsb-release -a查看linxu 版本号，gcc -v查看gcc版本号。具体语法请看终端的命令帮助。1.5 GAMIT/GLOBK结构。GAMIT/GLOBK安装完成后有几个重要的文件夹：gg/gamit和gg/kf内小模块是用fortran编写的，分 别为基线解算和平差（卡尔曼滤波）两个大功能。gg/com则是用cshell编
3、写的脚本，分别将gamit或kf文件夹内的众多小模块有机 的组织起来。gg/tables内包括测站接收机、天线和卫星信息等各种表文件，需要不定期更新。GAMIT/GLOBK是开源软件，允许你自 己对软件进行更改，重新编译。1.6需要注意的更新。LUNTAB、SOLTAB、NUTABL、LEAP.SEC的内容每年更新一次。UT1、POLE表每周更新一次。至于软件是 否需要更新，请关注如下网址上的说明，包括软件改进和bug修正。如果某个模块中的bug已经严重影响到你的计算结果，请下载并重 新编译 如果bug存在于lib文件中，则要求下载升级包后重新编译整个软件。http:/www-gpsg.mit
4、.edu/simon/gtgk/down.htm1.7 运行过程需要的目录。如果使用sh_gamit批处理，至少需要在工程目录下新建rinex，brdc和gfiles三个目录，分别存放你的观测数据o文件，广播星历n 文件，卫星轨道g文件。这样做的目的是把文件归类，最后这些文件都将被link到单天目录中，基线解算最终是在单天目录里运行的。 如果你学习分步处理，目录则比较自由，最后只要保证一点，所有必要的表文件和观测文件等都要link到单天中。1.8 部分模型的说明。（1）.otl 潮汐改正比较准确和成熟，一般都会使用。（2）.有研究认为vmf1映射函数目前是最好的，可比其它映射函数提高7%的重复
5、性（改进的GMF和它近似，两者可选用）。（3）.atml 大气（非潮汐）荷载模型只对高程方向影响较大，可消除大部分周期性跳动，其可靠性还需要进一步证实，因此应用不广。（4）.atl大气潮汐荷载模型和met气象模型是软件中保留的模型，目前（截至2007年7月）尚未开发。（5）. 现 在 的 GAMIT/GLOBK 软 件 把 改 正 模 型 中 的 栅 格 文 件 都 分 离 出 来 ， 需 要 到 url=ftp:/erverest.mit.eduftp:/erverest.mit.edu/url（不需要权限）下载。并做好 tables 中的各种 link。1.9与星历有关的四个文件关于brd
6、c文件、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:/igscbjpl.nasa.gov/network/complete.h
7、tmligs 站点 rinex 文件下载 http:/garner.ucsd.edu/pub/rinex/igs 站点广播星历下载 http:/garner.ucsd.edu/pub/nav/sp3 文 件 下 载http:/garner.ucsd.edu/pub/p roducts/一 个 很 系 统 的 教 程http:/seismo.be rkeley.edu/battag/GAMITw rkshp/w rkshp.html确保成功安装GAMIT/BLOCK;（需要交流具体细节者请按上面的方式联系我） 数据准备中可能用到的资源地址：各种星历表文件更新在 http:/garner.ucsd
8、.edu/pub/gamit/tables/全球 igs 分布图 http:/igscbjpl.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/tab
9、les 和单天目录 mywk/19983532. 把rinex格式的观测数据放入单天目录，例如mhcb3530.98n.Z解压 uncompress mhcb3530.98n.Zdoy命令换算日期和gps周和此周第几天， 下载精密星历放入单天目录，例如 igs09886.sp3.gz解压 gunzip igs09886.sp3.gz3. 在 tables 目录下运行 links.tables J2000 1998，单天 353 里运行 links.day 1998 353 lake参照.log文件和cvan.dat以及o观测文件编辑单天里文件station.info的内容，务必准确，例如*SI
10、TE 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
11、 0 0.0687 DHBCR 0.0000 0.0000 ASHZ12 8.24 ATDMCC4. 用 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输入工程名、轨道名，都是lake6. sh_sp3fit -f igs
12、12481.sp3 -o lake如果station.info没有指定轨道名，此处需要o参数sh_check_sess -sess 342 -type gfile -file glake3.342makej farb3420.03n jfarb3.342 生成卫星钟文件sh_check_sess -sess 342 -type jfile -file jfarb3.342makex lake.makex.batchfixdrv dlake3.342 生成所有准备文件7. 运行批处理 csh blake3.bat得到Ofile、Q-file文件以及用于后续GLOBK软件（运用卡尔曼滤波进行全网整
13、体平差的网平差软件）所需的Hfile文件注意：为了达到在下一步GLOBK模块中进行坐标转换的目标，需要引入已知点联合解算。例如全部固定到ITRF坐标，需要在gamit模 块中将周围稳定的，均匀分布的3个或者以上IGS站点引入联合解算。用 GLOBK 模块进行平差和坐标转换，对于区域网分析，不引入全球解的情况&建立目录glbf和soln（均与gamit中的单天目录平行），用htoglb将h解转换成二进制格式，用ls生成gdl文件，准备好需要的 apr 文件， globk 命令文件和 glorg 命令文件，并做好命令文件的相关设置。注意：globk命令进行平差，最后得到每天的坐标.prt文件，上述
14、坐标结果是在WGS84下的坐标，其内部几何形很好，但是每天坐标 框架的浮动都很大（cm级的浮动范围下即可求出整周模糊度），需要把每天的解都统一到同一框架中，这就是引入IGS站点的作用。IGS 站点在ITRF中的坐标是通过各种技术精确计算出来的，由某个历元的坐标和运动速度组成。例如在ITRFOO框架中，某个IGS站的坐 标=给定初始坐标+速度x时间差glorg命令可以通过IGS站计算出单天解对应的7个（或者3个）转换参数，并且利用转换参数将其 他待求点的 ITRF 坐标求出来。9. 运行glred命令，它会自动激发globk和glorg。前者命令进行平差，得到每天的坐标.prt文件。后者得到IT
15、RF坐标.org文件。10. 现在，用plotcrd把时间序列输出到屏幕吧，发现有粗差后剔除它，或者返回到gamit重新解算这天的基线。11如果你要的是周解或其他平均解，请在剔除粗差后运行globk命令将单天解联合解算，后续步骤类似上面介绍。sh_gamit 批处理1.如不是IGS站，将处理的站信息加入到软件的station.info中,将作为参考转换站的ITRF2000坐标和速度加入到itrfOO.apr中2建立工程名，其下建立子目录rinex ,brdc ,gfiles.分别放入（或者链接，下同）rinex文件（文件名中站名小写字母），广播星历文件（前 四个字符必须是brdc，而不是aut
16、o等其它）,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 100process.de
