gamit105数据处理图文流程

GAMIT10.5数据处理图文流程冉启顺2016.1.20目录一、前言2二、数据准备31.前期准备32.工程目录设置33.数据下载34.更新tables表，45.创立链接76.控制文件配置136.1准备测站列表sites.defaults136.2准备测站近似坐标文件lfile.136.3 文件的更新167 修改配置process.defaults、sittbl.与sestbl.文件197.1配置process.defaults197.2配置sittbl.坐标约束文件197.3配置sestbl.改正模型文件20三、分布处理221 .Makexp222精密星历处理223文件检查244广播星历处理245文件检查256生成中间文件257生成bat文件268基线解算279基线解算精度评定标准30四、cosagps平差31五、建议与说明34六、相关资料链接35一、 前言二、数据准备1.前期准备首先想好要处理什么时段、什么地方的数据本文档以处理bjfs,shao,lhaz,urum四个站点2013年2月1日的数据为例本文档使用软件版本GAMIT10.5处理环境：win10上的虚拟机ubuntu15.102.工程目录设置在任意位置建立工程文件夹（全文以test为例），并在其目录下建立如下的文件夹：3.数据下载在终端中打开test，查询2013年2月1日的年积日，和GPS周等信息代码：doy 2013 02 011 在test的目录下再建一个文件夹以032命名的文件夹，此时，test目录下有：2 根据步骤3的信息，到相关网站下载对应的广播星历，精密星历，观测值文件分别放在文件夹brdc,igs,rinex里面。当然还涉及到下载的数据解压，.d文件到.o文件的转换等操作，其中广播星历可以不改名字（即可以不用改成brdc0320.13n）4.更新tables表，更新准则，下面的表格为准：还有个dcb.dat需要每次处理都要更新表文件名表含义更新频次内容更新地址 pole.或pole.usno极移表每次更新/pub/gamit/tables/ ut1.或ut1.usno地球自转参数表每次更新 luntab.2012.J2000月亮表使用当年的表 soltab.2012.J2000太阳表使用当年的表 nutabl.2012章动表使用当年的表 leap.sec跳秒表每次更新gdetic.dat,大地坐标参数文件不常更新antmod.dat,天线相位中心改正正表有新天线时更新svnav.dat,卫星列表有新卫星时更新rcvant.dat,接收机天线对照表有新天线时更新svs_坏卫星信息文件每次更新vmf1映射函数模型 使用当年的表/pub/GRIDS/ otl.grd海潮格网数据使用当年的表i. 常见文件说明ii. 更新的准则是：“更新频次”中的红字部分每次处理都更新一下，同时结合本次处理与上次处理的年份是否相同来进行年表的选择。原本有些文件是每周更新一次，但是作为新手为了避免出错还是每次更新吧有些时候，有些文件找不到，或者/pub/gamit/tables/里面没有，则可以忽略。如svs_，gdetic.dat常常找不到我在更新时，这四个文件在服务器上没找到下图是我更新的文件。原本还应该有2013年的各种年表，但是我之前处理过一次2013年的数据，并且我的gamit安装路径/opt/gamit/GAMIT10.5/tables已经有了2013年的这些年表。iii. 将上述更新的数据复制到gamit安装目录下的tables文件夹内，并替换原有文件。我的gamit安装路径是/opt/gamit，即我安装路径下的tables 的路径是：/opt/gamit/GAMIT10.5/tablesiv. 然后将更新后的tables文件夹，即/opt/gamit/GAMIT10.5/tables拷贝（或者设置连接，我习惯拷贝）到test目录下。因为我之前在/opt/gamit/GAMIT10.5/tables里面放了很多数据，所以它显得特别大1.8G，往往在包含otl.grid的情况下tables也就7800M吧。此后，test的结构应是：5.创立链接这一步相当重要，而且常常容易出错。1） 将.o文件，.sp3文件，*.*n文件（观测文件、精密星历、广播星历）链接到032的目录下执行代码：2） 链接处理数据需要的各种年表等文件此时032中有如下的文件：检查并重新设置断掉的链接，如下图所示：如：luntab.链接已断（常说的死链接），在032目录下，执行代码：ln -s ./tables/luntab.2013.J2000 ./luntab.（最后的点号一定要加上）系统报错，即“无法创建符号链接，luntab.早已存在”，这句话的意思是，032中早已经有了一个luntab.文件（但是我们却看到它是个死链接），删除该链接，重新执行上面的代码。新链接生成，并打开查看该链接所指向的文件是不是2013年的luntab表重点检查6-ii的11个文件链接，以及soltab.,nutabl.，用类似的方法建立链接以下是我手动加入的部分链接：以下是网上的资料图：其中还要注意设置：map.grid，otl.grid，atl.grid的链接，我在处理数据时使用的如下链接：在查看原始GAMIT10.5的解压文件中的tables文件夹的链接atl.grid,如下图所示：但是我在查看/pub/GRIDS/发现它原本有一个atl.grid，上传是2010年8月28，大小也是3.1M，故我断定，服务器上的atl.grid就是tables里面的ANU_grid.atl。于是我在服务器上下载了一个3.1M的atl.grid拷贝到test/tables中，然后在032里面制作了一个atl.grid的链接。（当然也可以不用下载atl.grid,直接制作指向ANU_grid.atl的链接）综上所述：n 更新tables常更新以下11个文件（各种年表除外）：但常常 guess_rcvant.dat 用不到,可以不更新n 制作链接特别注意以下几个链接：n 在更新tables,建立完连接后执行ls al 检查链接情况，红色的是死链接：发现atml.grid的链接已断，于是重新设置即可。至于剩下的红色死链接，我全是忽略了，因为我也不会设置。（好像那个是在正式处理过程中产生的，可以忽略）总的来说设置链接是及其容易出错，要小心又小心。6.控制文件配置6.1准备测站列表sites.defaults在test/tables下把需要处理的测站加入该表，包括IGS测站和本地测站，并确定工程名。然后将site.defaults在032下设置一个链接 6.2准备测站近似坐标文件lfile.制作完lfile.文件后不要忘了将lfile.拷贝到test/tables中,同时在032中设置对应链接lfile.文件是测站近似点坐标（越精确越好,10.4以前的版本里面存放的经纬度信息，而10.4及其之后的版本直接存放x,y,z信息）GAMIT10.5 的lfile 和apr 文件是同样的格式，所以我们可以认为apr 文件就是lfile。该文件中保存的是各个站的初始坐标，至于该文件中的测站坐标要达到什么样的精度，可参考如下两条：l 如果是该站作为基线解算起算点,则该测站的坐标做保持在5cm 以内，最好不要超过10cm。l 如果该测站在基线解算中不作为起算点，则伪距单点定位的精度即可。对于IGS站，获取其在特定框架下的精确坐标方法：首先要确保test/rinex下有itrf08.apr方法1：我用的方法（可用于非IGS站）：在我的电脑上试了几次，老是最后一步出错，不推荐，但也有教程上是这个方法：l 在rinex文件夹下先生成lfile.rnx文件，然后生成lfile.rnx.apr，最后生成lfile.依次执行代码：grep POSITION *.13o lfile.rnxrx2apr lfile.rnx 2013 032gapr_to_l lfile.rnx.apr lfile. * 2013 032方法2：方法3：l 伪距单点定位计算非IGS站的初始坐标实例（假设shao是非IGS站）：执行完之后会生成两个不同格式的测站坐标文件lfile.shao 和shao.apr，然后执行l 双差定位计算测站的初始坐标实例：在此例中假设shao站是未知坐标的测站，而bjfs是已知坐标的测站（或起算站）综上所述：在生成life.文件时，IGS站的坐标可以直接从框架获得，或者利用双差定位求得。而非IGS站则可以根据要求精度的高低选择： 直接提取头信息 伪距单点定位 双差定位法 方法1常常在伪距单点定位和双差定位法中选择，然后将所须的坐标信息拷贝到lfile.中即可，并将lfile.拷贝到test/tables中,同时在032中设置对应链接现在，我假设bjfs，lhaz两站点是精度较高的IGS站，坐标直接从框架中获得，而shao,urum站是非IGS站，近似坐标通过伪距单点定位获得。最后按照特定的格式均复制到lfile.文件中然后将lfile.拷贝到test/tables中,同时在032中设置对应链接好像不需要在032中设置链接，但是怕出错，我还是设置了，希望没有弄巧成拙6.3 文件的更新最后要在032中检查链接是否指向了更新后的文件在gamit全英文的官方教程（版本10.6,24页）里面有这样一句话：可见station.Info的链接可以指向下图中选中的文件，也就是说倘若你的test/tables中没有文件，可以自己新建并指向他们。并且这个MIT,SOPAC的文件好像可以在上更新，只是登不上去。练手的话凑合着能用。在test/tables目录下：先建立一个sitelist的文件，并录入测站名字（每行首字母空格），不分先后顺序然后终端打开tables并执行代码：显示 files have 3 sites，这是错误的，一般来说总共有多少测站参与计算就会有多少个sites，于是回去检查sitelist, ，发现拉萨站写错了，应该是lhaz（所以在操作过程中一定要小心又小心，切记不能将单词写错，同时也要不害怕出现错误，只要根据其提示，顺藤摸瓜，基本上能修改错误的地方，实在不行重来）还发现，Recevier name * not found in rcvant.dat,网上百度了一下说要修改，说得又很模糊，所以不会。又去问了一下师兄，他说可以忽略，所以我也就忽略了（谁知道的话告诉我一声）。执行上述代码结束后最后一行显示：生成了一个.new，此时需要将其改名为，同时删除原有的的链接。再执行一次更新代码：在生成的文件中会有下图中红框的字符串，表示更新成功：最后在032中检查链接的指向，如果没有的话，自己手动用代码建立一个。7 修改配置process.defaults、sittbl.与sestbl.文件7.1配置process.defaults它里面包含目录信息，改动较小。需要在032中设置链接7.2配置sittbl.坐标约束文件 对高精度的已知坐标强约束，待求点坐标松约束非起算点约束量可以更大,如100m。并且如果IGS站列表中没有你所处理数据中的IGS站，需要自己添加，格式一致。7.3配置sestbl.改正模型文件配置准则如下面的截图：下面是我的部分设置：我没有met.list,met.grid,也没有map.list和map.grid,所以我设置的是N（不使用）我的atml.list和atl.list是死链接，我又不知道如何制作他们（好像在麻省理工的FTP上可以下载，但是这个FTP好久都登不上去了），于是全选N，而我有atl.grid和atml.grid ，于是选Y修改完了，别忘了保存，同时为了保险起见到032里面再检查一次几个配置文件是否修改、链接好了至此，数据准备阶段完成。下面我们开始计算。三、分布处理1 .Makexp在032目录下执行代码：makexp2精密星历处理如下图执行代码：3文件检查执行下图的代码：对于代码中的gigsf3.032,其中g表示文件类型，中间igsf表示轨道，第一步输入的：，3表示年的最后一位，032是年积日。4广播星历处理其中auto0320.13n没有改名字作为输入文件，jauto3.032是输出文件当然也可以在将auto0320.13n改名字后，执行 makej brdc0320.13n jbrdc3.0325文件检查6生成中间文件7生成bat文件8基线解算最后一步还是出错了，但是我百度了一下“RC not executed STOP FATAL Error: Stop from report_stat”在这篇文章的最后一页找到解决方法，也不知道对不对，反正最后是成功解算了。/p-2.html在032中找到上图的两个文件并删除，然后重新执行基线解算。中间有很多省略结算结果精度还算过得去但是我往往以otesta.032为准qexpta.doy:解算记录oexpta.doy为qexpta.doy的简略版hexpta.yrdoy：协方差矩阵、参数平差值9基线解算精度评定标准（网上各类教程上的标准大同小异）标准1：一般NRMS越小越好，一般小于0.25视为解算效果较好，本次解算的NRMS为0.26，符合要求基线的X、Y方向误差范围应控制在1cm以内，Z方向的误差控制在1.5cm以内，本次解算所有基线均满足要求本次解算x,y,z方向的误差均有超过1cm,结果不是很理想，还过得去最后还有一个基线重复性检验：但是我的检验失败了，是不是要进行两天以上的数据处理才能进行重复性试验？至此可以理解为解算完毕，然后进行基线平差，也可以使用GLBOK平差，但是我选择了用COSAGPS进行的平差。如果想直接用GLBOK进行平差，参考安向东的处理流程。四、cosagps平差1. 在032目录下打开otesta.032,并在第二次出现基线信息的地方加入下面红框中的一行字符串：“cosagps for gamit o-file”最后保存并复制到windows桌面下。2. Cosa平差打开cosagps,新建工程，设置如下，其它参数基本默认3. 导入基线数据导入从ubuntu里面复制出来的otesta.032基线的误差有点大（但为了走完流程，有精度强迫症的就凑合着看吧_）4. 导入三维已知坐标这里就是你解算时坐标起算数据，但是要注意框架和坐标系的统一。先在点名方框内随便输入数据，然后确定，如下图在工程目录下找到刚新建的关于起算点坐标的文件，用记事本打开，并将已知起算数据复制粘贴到know文件中，如下图：5. 三维平差最后将解算出的shao,urum站的坐标与原本查询到的坐标进行对比即可。以下图片里面的坐标就是我在ubuntu里面采用：获得的，对比后发现相差不是太大，但是还是不太理想（不管怎么样将软件走通了）五、建议与说明六、相关资料链接天线查询/ANTCAL/TEQC/software/teqc/广播星历下载/pub/rinex/YYYY/DOY/(以auto开头的文件，下载后改为brdc开头)/pub/gps/data/daily/YYYY/brdc/（以brdc开头）ftp:/igs.ensg.ign.fr/pub/igs/data/YYYY/DOY/（以brdc开头）精密星历下载/p