GAMIT学习资料整理.docx

GAMIT学习资料整理.docx

《GAMIT学习资料整理.docx》由会员分享，可在线阅读，更多相关《GAMIT学习资料整理.docx（14页珍藏版）》请在冰点文库上搜索。

GAMIT学习资料整理

GAMIT10.4安装（基于ubuntu10.04）

1、安装虚拟机vmware、ubuntu10.04；

2、进入终端输入：

sudopasswdroot为root用户创建密码，并以root用户登录，或sudosu回车；

3、系统更新、汉化；

4、安装gcc：

#apt-getinstallgcc；

5、安装csh：

#apt-getinstallcsh修改bash为csh，重启；

6、安装gfortran：

#apt-getinstallgfortran；

7、安装libx11-dev库支持#apt-getinstalllibx11-dev；

8、修改shall为bash，重启，并设置路径：

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

exportPATH=$PATH:

/opt/gamit/gamit/bin:

/opt/gamit/com:

/opt/gamit/kf/bin

exportsetenvHELP_DIR=/opt/gamit/help/

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

进入目录：

#cd/opt/gamit

10、修改install_software文件内容：

#geditinstall_software，打开install_software文件，在文件的中下部修改“usr-namelibX11.a”为"usr-namelibX11.so"。

（动态共享库）

11、运行install_software：

#./install_software，按提示输入两次Y后，修改make，在/opt/gamit/libraries里，修改Make中的一组参数

（1）MAXSIT55、MAXSAT32、MAXATM25、MAXEPC5760

（2）#SpecifictoFC5（F6,F7,F8）

然后，在输入两次Y完成安装；

12、安装完后，打开终端输入：

doy，查看程序是否已正确安装。

（注：

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

原因:

需要修改软件的权限（chmod777软件名）

、数据准备阶段

（个人理解：

这一过程可以先在windows下准备好，到时候再拷到linux的工程文件夹）

1.观测文件标准化：

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

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

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

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

1）v/CORS/Gpscal.html查询观测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：

（注：

如果没有的话，数据处理的最后一步cshbtest1.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–filessite1doy0.yyosite2doy0.yyo–orbtigsf–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

cdtest/rinex/all_station

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

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

geditstation.info

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

2）l-file：

记录各测站的近似坐标:

Gamit要求是球面坐标，即l-file。

将所有o文件整理在一起，然后分别运行如下命令:

grepPOSITION*.11o>l

rx2aprl2011250（注：

根据所有的o文件中的x、y、z取平均值,文件名为l）

gapr_to_lllfile.““2011250

我对l-file文件的准备（武大郭博峰）：

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

2、进入该文件夹：

cdtest/rinex/all_data

3、提取POSITION

grepPOSITION*.11o>l

4、查看l文件是否正确

geditl

5、rx2aprl2011250可得到l文件（空间直角坐标）

（注：

根据所有的o文件中的x、y、z取平均值）

6、l转到lfile.（球面坐标）

Gapr_to_lllfile.“”2011250

7、将得到的l复制到lfile.文件夹下

3）sestbl：

（该表是gamit的关键，需要仔细学习）

从tables文件夹下复制到lfiles中。

Gamit基线处理时的模型选择配置文件。

需要仔细参照技术文档和程序代码，一般情况下采用默认值，通常需要修改的地方有：

ChoiceofObservable（观测值的选择类型）、ChoiceofExperiment（基线处理类型）两项.还有高度截止角设定，禁止固体潮改正，变潮汐改正系数个数31为23.

这个模型参数的设置只要是对处理方式进行选取和误差改正模型的选取。

要是设定了误差模型改正，那在tables表里需要有这些模型改正文件，同时也需要将这些文件链接到各个观测时段文件里。

否则在处理的时候会出错的。

4）sittbl

工程网和框架网需要分开生成，例如sittbl_kjw和sittbl_gcw。

也是从tables文件夹下复制到lfiles中。

各测站的精度控制指标。

一般情况下采用默认值。

对高精度的已知坐标强约束待求点坐标松弛约束。

是否需要igs跟踪站点的文件做约束吗？

3.链接各个测站观测数据（包括brdc*.11n及igs*.sp3文件）

在各时段文件夹下建立观测数据连接，在各时段文件夹下输入命令：

（由于此法将GAMIT装在/OPT下，而test文件夹在/home下，为两个不同的文件系统，因此要用符号链接。

而我自己的方法将GAMIT装在/home/jiang下为相同的文件系统。

）

ln–s../rinex/*doy？

.yyo./

ln–s../igs/*igsweek.sp3./

ln–s../brdc/*doy？

.yyn./

4.将tables表文件中需要的文件（包括luntab、soltab、nutabl、pole、ut1、pmu文件，还有一些误差改正模型文件等）链接到各时段文件夹。

输入以下命令建立链接：

（注：

ln–s是符号连接的意思）

ln–s../tables/luntab.2011.J2000./luntab.

ln–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–expttest–orbtigsf–yr2011–doy250–sess99–srin–navbrdc2500.11n–aprlfile.–sinfo1500005760

或者：

直接运行makexp，然后按提示输入

2）运行sh_sp3fit脚本，生成轨道初始根数；

（曾由于ut1.链接的不是最新的，导致错误。

）

sh_sp3fit–figs15623.sp3–oigsf–d2011250–rBERNE–t–u

3）运行sh_check_sess脚本，检查卫星一致性；

sh_check_sess–sess250–typegfile–50

4）运行makej程序，生成卫星钟差文件；

makejbrdc2500.11njbrdc1.250

5）运行sh_check_sess，检查卫星一致性；

sh_check_sess–sess250–typejfile–

6）运行makex程序，生成X文件；

makextest.makex.batch

7）运行fixdrv程序，生成批处理文件；

fixdrvdtest1.250（进行这步可能发生错误，提示FATAL:

121202:

2337:

41.0FIXDRV/armake:

InvalidT-=1027201210272012--needeitherT-oneX-runFIXDR修改fixdrv/armake.f文件中的2010为2099或9999）并在此目录下运行make-fMakefile进行重新编译。

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

cshbtest1.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）高度角ElevationCutoff（不能是0°，按规范来为15°）；

2）ChoiceofObservable长基线选择LC_help、短基线选择L1\L2Independent

3）其他小错误该正项

附录1：

下载IGS观测数据的网站

garner.ucsd.edu

everest.mit.edu

cddis.gsfc.nasa.gov

广播星历：

精密星历：

tables表：

IGS站点的ITRF坐标：

IGS站点的O文件：

例如

广播星历和IGS站点O文件的总目录：

全球IGS站点分布图：

数据准备中可能用到的资源地址：

各种星历表文件更新在

全球igs分布图

igs站点rinex文件下载

igs站点广播星历下载

sp3文件下载

一个很系统的教程

附录2:

GAMIT文件说明

A-文件：

ASCII码的T文件。

B-文件：

控制批处理的数据文件。

C-文件：

残差（O-C）和偏导数文件。

D-文件：

时段和接收机开关文件。

E-文件：

广播星历文件或RINEX的导航文件。

G-文件：

初始轨道状态向量和光压参数值文件。

H-文件：

GAMIT解输出的方差协方差文件，是GLOBK的输入文件。

I-文件：

接收机多项式钟差文件。

J-文件：

卫星多项式钟差系数文件。

K-文件：

从伪距得到的接收机钟差文件。

L-文件：

测站坐标文件。

M-文件：

控制C-文件和SOLVE模块和编辑模块的文件。

N-文件：

autcln.sum.posfit生成的数据加权文件。

O-文件：

用于事后分析的解文件（强约束）。

P-文件：

MODEL运行过程记录文件。

Q-文件：

分析结果文件。

S-文件：

测站坐标和天线文件（已不使用）。

T-文件：

表格化的星历文件。

U-文件：

测站海潮文件。

V-文件：

SINCLN，DBLCLN和SCANRMS的编辑输出文件。

W-文件：

气象数据文件。

X-文件：

预处理后的观测数据文件。

Y-文件：

卫星YAW偏航文件。

Z-文件：

水汽辐射计数据文件。

Sestbl.中的几个选项说明

一、ChoiceofObservable=

1、LCAUTCLN意味着观测值是无电离层的线性组合，当前的默认值，因为现有的接收机可同时接收L1和L2波段上的P码，所以新方案的效果总是较为理想尤其是用于在长基线时；2、由于接收机信号的交叉相关性（96年之前的接收机普遍是这样的），卫星波段在C1和P2之间波动导致了宽相模糊度较差或者错误的解算结果，所以这方面要有所注意，必要时要采取保守的方法使用电离层约束（LCHELP）；

3、ChoiceofObservable=L1,L2INDEPEND或者L1_ONLY用于很短的基线（Fornetworksextendingoverlessthanafewkilometers）；

4、如果只有单频接收机的数据，则选L1_RECEIVERS，L1_ONLY也可以；若测段中有单频接收机，必须在autcln命令文件中设置L1_only；

5、使用大范围的基站来确定轨道，设置ChoiceofObservable=LC（不存在模糊度）可以节约计算时间。

二、ChoiceofExperiment=BASELINE；BASELINE/RELAX./ORBIT

1、BASELINE只解基线；

2、RELAX.同时解基线和轨道；

3、ORBIT只确定轨道。

三、Typeofanalysis=1-ITER；

控制着通过model，autcln和solve等一系列过程获得最终结果的数量。

较为理想的方法是Typeofanalysis=1-ITER，这样可以得到两个通过model,autcln和solve的量，第一个量是为了获得一个初值用于得到接近最终结果的坐标和轨道参数值，进而整平残差编辑和减小来自于大平差的误差。

第一次通过时autcln作用于前附和残差值，第二次通过时使用由solve写入M文件的平差值进行后附和编辑。

四、

若没有加潮汐改正则设置为：

Tidesapplied=23；Useotl.grid=N

GAMIT10.4中的潮汐改正文件为：

otl_FES2004.grid（大小696M）。

（相对定位基线超过100KM才加上此项改正,个人观点）

部分模型的说明

（1）.otl潮汐改正比较准确和成熟，一般都会使用。

（2）.有研究认为vmf1映射函数目前是最好的，可比其它映射函数提高7%的重复性（但实时性较差，34h的延迟）（改进的GMF和它近似，但无时延问题，两者可选用）。

对流层相关知识详见《GPS测量与数据处理》（第二版）P113-116

（3）.atml大气（非潮汐）荷载模型只对高程方向影响较大，可消除大部分周期性跳动，其可靠性还需要进一步证实，因此应用不广。

（4）.atl大气潮汐荷载模型和met气象模型是软件中保留的模型，目前（截至2007年7月）尚未开发。

（5）.现在的GAMIT/GLOBK软件把改正模型中的栅格文件都分离出来，需要到[url=][/url]（不需要权限）下载。

并做好tables中的各种link。

*.Sp3文件说明

下面对sp3文件做个简单格式说明。

文件部分截图如下：

第一行的第二列为版本，这里为a，后续高版本将依次按照字母顺序命名；

           第三列P表示坐标，注意有的版本包含速率的时候该字符为V，接着为时间（4-7列（年）、9-10列（月）、12-13列（日）、15-16列（时）、18-19列（分）、21-31列（秒））；

          96代表历元数（33-39列）；

          IGS00代表坐标系统（47-51列）。

    第二行1148表示GPS周（4-7列）；

           0.00000000周内秒（9-23列）；

            900.00000000历元间隔（25-38列）；

           52280儒略日（40-44列）；

            0.00儒略日小数秒（46-60）。

    第三行到第七行为卫星编码，其中第三行的28为卫星数（5-6列）。

    第八行至第十二行为精度。

    十三至十八行不说了。

    19到22行为注释行。

    23行观测历元：

4-7列（年）、9-10列（月）、12-13列（日）、15-16列（时）、18-19列（分）、21-31列（秒）。

    24行第一列为P（位置，单位km），2-4列卫星号，

          5-18列 x坐标，

          19-32列 y坐标 

          33-46列  z坐标（注意当第一列标识为V表示速率时候，其单位为分米/秒（dm/sec））

           47-60列钟差（单位微妙，当第一列标识为V时表示时间变化率，其单位为毫秒/秒（microsec/sec））。

注意如果数据行第75列有’E‘字符，表示卫星坐标或者钟记录存在事件。