GAMIT学习资料整理

GAMIT10.4安装(基于 10.04)

1、安装机VMware、Ubuntu10.04;

2、进入终端输入:sudo passwd root 为root用户创建密码,并以root用户登录,或sudo su回车;

3、系统更新、汉化;

4、安装gcc:# apt-get install gcc ;

5、安装csh: # apt-get install csh 修改bash为csh,重启;

6、安装gfortran : # apt-get install gfortran ;

7、安装libx11-dev库支持 # apt-get install libx11-dev ;

8、修改shall为bash,重启,并设置路径:

回到用户根目录,打开.bashrc ,在最后加上如下代码即可

export PATH=$PATH:/opt/gamit/gamit/bin:/opt/gamit/com:/opt/gamit/kf/bin export setenv HELP_DIR=/opt/gamit/help/

9、将gamit安装包放在目录opt/gamit/下

进入目录 :# cd /opt/gamit

10、修改install_software文件内容 : # gedit install_software , 打开install_software文件,在文件的中下部修改“usr -name libX11.a”为”usr -name libX11.so” 。(动态共享库)

11、运行install_software : # ./install_software, 按提示输入两次Y后,修改makefile.config ,在/opt/gamit/libraries里,修改Makefile.config中的一组参数

(1)MAXSIT 55 、MAXSAT 32 、MAXATM 25 、MAXEPC 5760

(2)# Specific to FC5(F6,F7,F8 )

然后,在输入两次Y完成安装;

12、安装完后,打开终端输入:doy,查看程序是否已正确安装。

(注:在root用户下安装完后,回到自己用户不能使用,

原因:需要修改软件的权限(chmod 777 软件名)

、数据准备阶段

(个人理解:这一过程可以先在下准备好,到时候再拷到的工程文件夹)

1.观测文件标准化:

1)观测文件o文件的文件名要小写,观测文件里的点名一般要大写。

2)认真检查观测文件内接收机类型、天线类型并查看GAMIT的ant.dat/rcvant.dat中是否有相应的类型,否则需要更新tables中的文件。

3)天线高及天线测量方式。

2.需要从网上下载的数据文件

1) 查询观测 0-files 文件所在的GPS周天,年积日与 UTC 之间的关系。

2)导航电文文件brdcdoyn.11n、精密星历*.sp3文件以及用到的IGS跟踪站 O-files数据。(导航电文文件下载autodoyn.11n;注意:文件名都要改成小写)

3.一些常见表文件说明

dcb.dat:码相关型接收机伪距改正参数统计表;(缺少这个文件基线处理最后一部会出错)

antmod.dat:天线相位中心参数文件;

rcvant.dat:接收机及天线名称对照表;

guess_rcvant.dat: (如果rcvant.dat中没有rinex头文件中的接收机或天线型号,将从该文件中获取,用户可以自己修改)(该文件在处理过程中会自己产生,一般没用不需要准备) otl_FES2004.grid: GAMIT10.4中的潮汐改正文件,696M

注:地球固体潮改正和海洋负荷潮改正:(详见《卫星导航定位新技术及高精度数据处理方法》p7-p8)

地球固体潮改正:摄动天体(月亮、太阳)对弹性地球的引力使地球表面产生周期性的涨落,称为固体潮现象。在小于100KM的GPS相对定位中,两个测站的固体潮影响几乎是相同的,在差分过程中可抵消,因此可不考虑此项改正。

海洋负荷潮改正:由海洋潮汐的周期性涨落所引起的。与固体潮类似,但小一个数量级。 需要每周更新的:

ut1.: 地球自转参数; ln –s ***

pole.:极移参数; ln –s ***

需要每年更新的:

leap.sec:跳秒表;

soltab.:太阳星历; ln –s ***

luntab.:月亮星力; ln –s ***

nutabl.:章动参数表; ln –s ***

gdetic.dat:各种大地坐标系参数;

pmu.usno: ln –s ***

grid.otc: ln –s ***(这个好像有问题,GAMIT10.4中不是这个文件)? svs_exclude.dat:需要剔除卫星列表;(没有的话,处理过程中会提示waring,但也可以计算)(注:应该是每次计算都需要更新的,个人理解)

vg_in:测站坐标速度表;

autcln.cmd:(注:如果没有的话,数据处理的最后一步 csh btest1.bat不能正常进行) 需要重点准备的:

sestbl:数据处理参数设置表;

sittbl:测站约束表;

station.info:测站信息表;

lfile.:站初始坐标(球坐标);(自己制作)

、文件准备阶段

1.文件结构整理

test:工程目录,在其下面包括:

1)tables/存放各种表文件;

2)igs/存放精密星历*.sp3文件;

3)rinex/存放观测数据文件*doyn.11文件;

4)brdc/存放广播星历文件brdcdoy0.11n文件;

5)doy1/ doy2/ ??/存放各个时段解算的文件及链接文件。

2.需要重点准备的4个文件

1) station.info:

记录各测站的站点、站点名称、开始观测、结束观测、天线高、 天线量测方法、接收机类型、软硬件版本号、天线类型等信息。严格按照已经给出的跟踪站格式编辑。(注意:要现在当前目录下拷入station.info文件,把内容删掉留下前面的几行,否则后面的步骤没法进行)

具体操作方法有二:

1、手动操作,按照相应文件格式对齐操作;可以从样本格式开始编辑。

2、自动提取观测文件的有关信息,生成 station.info 文件。

在各时段文件目录下操作:

sh_upd_stnfo –files site1doy0.yyo site2doy0.yyo –orbt igsf –u –ref ../tables/station.info

运行上述命令后,就自动从观测文件 0-files中提取重要信息,生成 station.info 文件。

我对station.info文件的准备(武大郭博峰):

1、rinex文件夹下新建文件夹all_station(存放所有点)

2、工程观测文件(o-files)复制到all_station,框架网(o-files)复制到all_station

3、进入Terminal

cd test/rinex/all_station

sh_upd_stnfo –files *.11o -ref station.info(此文件为软链接文件,从~/不常更新表/station.info 链接至指定文件夹下,且只有一个表头)

4、查看station.info文件是否正确

gedit station.info

5、将得到的station.info复制到lfiles文件夹下

2) l-file:

记录各测站的近似坐标: Gamit要求是球面坐标,即l-file。将所有 o 文件整理在一起,然后分别运行如下命令:

grep POSITION *.11o > lfile.rnx

rx2apr lfile.rnx 2011 250 (注:根据所有的o文件中的x、y、z取平均值,文件名为lfile.rnx.apr)

gapr_to_l lfile.rnx.apr lfile. “ “ 2011 250

我对l-file文件的准备(武大郭博峰):

1、在rinex文件夹下新建文件夹all_data(存放所有点)

2、进入该文件夹:cd test/rinex/all_data

3、提取POSITION

grep POSITION *.11o >lfile.xyz

4、查看lfile.xyz文件是否正确

–s ../tables/soltab.2011.J2000 ./soltab.

ln –s ../tables/nutabl.2011 . /nutabl.

ln –s ../tables/pole.usno . /pole.

ln –s ../tables/ut1.usno . /ut1.

ln –s ../tables/pmu.usno . /pmu.

另外,潮汐改正文件也是很容易出错的文件。以10.4版本为例, 默认的潮汐 改正文件名为otl.grid(此文件为损坏文件)(链接前是otl._FES2004.grid )。一定要检查一下 otl.grid 是否是真实链接,出现死链接的情况需要手动链接。命令为:

ln -s ../tables/otl._FES2004.grid ./otl.grid (我在算的时候没有用这个文件,下次试试用)

(注意:链接时,如果该文件夹中已有该文件,无论链接成功与否,要重新链接时都要先删除原文件在进行链接。)

、数据处理阶段

一、批处理

1) 运行makexp 程序,生成输入文件

sh_makexp –expt test –orbt igsf –yr 2011 –doy 250 –sess 99 –srin –nav brdc2500.11n –apr lfile. –sinfo 15 00 00 5760

或者:直接运行 makexp,然后按提示输入

2) 运行sh_sp3fit 脚本,生成轨道初始根数;

(曾由于ut1.链接的不是最新的,导致错误。)

sh_sp3fit –f igs15623.sp3 –o igsf –d 2011 250 –r BERNE –t –u

3) 运行sh_check_sess 脚本,检查卫星一致性;

sh_check_sess –sess 250 –type gfile –file gigsf1.250

4) 运行makej 程序,生成卫星钟差文件;

makej brdc2500.11n jbrdc1.250

5) 运行sh_check_sess,检查卫星一致性;

sh_check_sess –sess 250 –type jfile –file jbrdc1.250

6) 运行makex 程序,生成X 文件;

makex test.makex.batch

7) 运行fixdrv 程序,生成批处理文件;

fixdrv dtest1.250 (进行这步可能发生错误,提示

FATAL :121202:2337:41.0 FIXDRV/armake: Invalid T-file interval = 10 27 2012 10 27 2012 –need either T-file or one X-file to run FIXDR 修改fixdrv/armake.f 文件中的2010为2099或9999)并在此目录下运行 make -f Makefile 进行重新编译。

8) 运行fixdrv 生成的批处理文件

csh btest1.bat 运行此命令前先把GAMIT.FATAL 和 GAMTI.WARNING 删除

二、 gamit 基线处理结果质量评价指标: (检验q文件)

A)标准均方根 nrms解算结果为 0.25左右,若是 nrms大于 0.5 ,表明处理结果是有问题的(例如:周跳没有修复;测站的起算点坐标有问题等)。若是nrms小于 0.5 ,则认为成功解算, 否则需要检查原因,重新处理;

B) 查看gamit基线解算精度,一般要求相对精度达到10-8~10-9

C) 检查所有测站点是否都参与计算

说明:如果有问题,需要查看

1.查看station.info文件: REC#、 ANT#、 HtCod

IGS跟踪站的天线高都是到DHARP的。

外业观测的需要问清楚当时量取天线的位置。

2.查看sestbl文件

1) 高度角 Elevation Cutoff(不能是0°,按规范来为15°);

2) Choice of Observable 长基线选择LC_help、短基线选择L1\L2 Independent

3) 其他小错误该正项

GAMIT相关阅读

GAMIT10.4及最新版本的GAMIT下载在公社5号FTP服务器,具体下载见

转载自:https://www.linuxidc.com/Linux/2016-02/128250.htm

声明: 除非转自他站(如有侵权,请联系处理)外,本文采用 BY-NC-SA 协议进行授权 | 智乐兔
转载请注明:转自《GAMIT学习资料整理
本文地址:https://www.zhiletu.com/archives-6002.html
关注公众号:智乐兔

赞赏

wechat pay微信赞赏alipay pay支付宝赞赏

上一篇
下一篇

相关文章

在线留言

你必须 登录后 才能留言!