GAMIT学习资料整理.docx

GAMIT10.4安装（基于ubuntu10.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 软件名)、 数据准备阶段（个人理解：这一过程可以先在windows下准备好，到时候再拷到linux的工程文件夹）1.观测文件标准化：1）观测文件o文件的文件名要小写，观测文件里的点名一般要大写。 2）认真检查观测文件内接收机类型、天线类型并查看GAMIT的ant.dat/rcvant.dat中是否有相应的类型，否则需要更新tables中的文件。 3）天线高及天线测量方式。2.需要从网上下载的数据文件1）/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：（注：如果没有的话，数据处理的最后一步 csh btest1.bat不能正常进行）需要重点准备的：sestbl：数据处理参数设置表；sittbl：测站约束表；：测站信息表；lfile.：站初始坐标（球坐标）；（自己制作）、 文件准备阶段1.文件结构整理test：工程目录，在其下面包括：1）tables/存放各种表文件； 2）igs/存放精密星历*.sp3文件；3）rinex/存放观测数据文件*doyn.11文件； 4）brdc/存放广播星历文件brdcdoy0.11n文件；5）doy1/ doy2/ /存放各个时段解算的文件及链接文件。2.需要重点准备的4个文件1) : 记录各测站的站点、站点名称、开始观测时间、结束观测时间、天线高、 天线量测方法、接收机类型、软硬件版本号、天线类型等信息。严格按照已经给出的跟踪站格式编辑。(注意：要现在当前目录下拷入文件，把内容删掉留下前面的几行，否则后面的步骤没法进行)具体操作方法有二： 1、手动操作，按照相应文件格式对齐操作；可以从样本格式开始编辑。 2、自动提取观测文件的有关信息，生成 文件。 在各时段文件目录下操作： sh_upd_stnfo files site1doy0.yyo site2doy0.yyo orbt igsf u ref ./tables/ 运行上述命令后，就自动从观测文件 0-files中提取重要信息，生成 文件。 我对文件的准备（武大郭博峰）： 1、rinex文件夹下新建文件夹all_station(存放所有点) 2、工程观测文件(o-files)复制到all_station，框架网（o-files）复制到all_station3、进入Terminal cd test/rinex/all_station sh_upd_stnfo files *.11o -ref (此文件为软链接文件，从/不常更新表/ 链接至指定文件夹下，且只有一个表头)4、查看文件是否正确 gedit 5、将得到的复制到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文件是否正确 gedit lfile.xyz 5、rx2apr lfile.xyz 2011 250 可得到lfile.xyz.apr文件（空间直角坐标） (注：根据所有的o文件中的x、y、z取平均值)6、lfile.xyz.apr转到 lfile. (球面坐标) Gapr_to_l lfile.xyz.apr lfile. “ ” 2011 250 7、将得到的lfile.xyz.apr复制到lfile.文件夹下3）sestbl：（该表是gamit的关键，需要仔细学习）从tables文件夹下复制到lfiles中。 Gamit基线处理时的模型选择配置文件。需要仔细参照技术文档和程序代码， 一般情况下采用默认值，通常需要修改的地方有：Choice of Observable（观测值的选择类型）、Choice of Experiment（基线处理类型）两项.还有高度截止角设定，禁止固体潮改正，变潮汐改正系数个数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 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 u3) 运行sh_check_sess 脚本，检查卫星一致性；sh_check_sess sess 250 type gfile file gigsf1.250 4) 运行makej 程序，生成卫星钟差文件；makej brdc2500.11n jbrdc1.2505) 运行sh_check_sess，检查卫星一致性；sh_check_sess sess 250 type jfile file jbrdc1.2506) 运行makex 程序，生成X 文件；makex test.makex.batch7) 运行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-810-9 C) 检查所有测站点是否都参与计算说明：如果有问题，需要查看 1.查看文件： REC#、 ANT#、 HtCod IGS跟踪站的天线高都是到DHARP的。外业观测的需要问清楚当时量取天线的位置。2.查看sestbl文件 1) 高度角 Elevation Cutoff(不能是0，按规范来为15)； 2) Choice of Observable 长基线选择LC_help、短基线选择L1L2 Independent 3) 其他小错误该正项附录1：下载IGS观测数据的网站 广播星历：/pub/gps/data/daily/2010/brdc/精密星历：/freebies/gpscalendar.htmltables表：/pub/gamitIGS站点的ITRF坐标：/cgi-bin/sector.cgiIGS站点的O文件：/pub/gps/data/daily/例如/pub/gps/data/daily/2010/326/10o/广播星历和IGS站点O文件的总目录：/pub/gps/data/daily/全球IGS站点分布图：/network/netindex.html数据准备中可能用到的资源地址：各种星历表文件更新在 /pub/gamit/tables/ 全球igs分布图 /network/complete.html igs站点rinex文件下载 /pub/rinex/ igs站点广播星历下载 /pub/nav/ sp3文件下载 /pub/products/ 一个很系统的教程 /battag/GAMITwrkshp/wrkshp.html 附录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.中的几个选项说明一、Choice of Observable=1、LC AUTCLN 意味着观测值是无电离层的线性组合，当前的默认值，因为现有的接收机可同时接收L1和L2波段上的P码，所以新方案的效果总是较为理想尤其是用于在长基线时；2、由于接收机信号的交叉相关性（96年之前的接收机普遍是这样的），卫星波段在C1和P2之间波动导致了宽相模糊度较差或者错误的解算结果，所以这方面要有所注意，必要时要采取保守的方法使用电离层约束(LC HELP)；3、Choice of Observable=L1,L2 INDEPEND或者L1_ONLY用于很短的基线（For networks extending over less than a few kilometers）；4、如果只有单频接收机的数据，则选L1_RECEIVERS，L1_ONLY也可以；若测段中有单频接收机，必须在autcln命令文件中设置L1_only；5、使用大范围的基站来确定轨道，设置Choice of Observable=LC(不存在模糊度)可以节约计算时间。二、Choice of Experiment = BASELINE ； BASELINE/RELAX./ORBIT 1、BASELINE只解基线；2、RELAX.同时解基线和轨道；3、ORBIT 只确定轨道。三、Type of analysis = 1-ITER ；控制着通过model，autcln和solve等一系列过程获得最终结果的数量。较为理想的方法是Type of analysis = 1-ITER，这样可以得到两个通过model,autcln和solve的量，第一个量是为了获得一个初值用于得到接近最终结果的坐标和轨道参数值，进而整平残差编辑和减小来自于大平差的误差。第一次通过时autcln作用于前附和残差值，第二次通过时使用由solve写入M文件的平差值进行后附和编辑。四、若没有加潮汐改正则设置为：Tides applied = 23；Use otl.grid = NGAMIT10.4中的潮汐改正文件为：otl_FES2004.grid（大小696M）。(相对定位基线超过100KM才加上此项改正,个人观点)部分模型的说明(1).otl潮汐改正比较准确和成熟，一般都会使用。 (2).有研究认为vmf1映射函数目前是最好的，可比其它映射函数提高7%的重复性(但实时性较差