GAMIT10.5大数据处理图文流程

《GAMIT10.5大数据处理图文流程》由会员分享，可在线阅读，更多相关《GAMIT10.5大数据处理图文流程（43页珍藏版）》请在装配图网上搜索。

1、word冉启顺目录一、前言2二、数据准备31.前期准备32.工程目录设置33.数据下载34.更新tables表，45.创立76.控制文件配置13准备测站列表13准备测站近似坐标文件lfile.13文件的更新167 修改配置、sittbl.与sestbl.文件19配置19配置sittbl.坐标约束文件19配置sestbl.改正模型文件20三、分布处理221 .Makexp222精细星历处理223文件检查244广播星历处理245文件检查256生成中间文件257生成bat文件268基线解算279基线解算精度评定标准30四、cosagps平差31五、建议与说明34六、相关资料35一、 前言二、数据准备

2、1.前期准备首先想好要处理什么时段、什么地方的数据本文档以处理bjfs,shao,lhaz,urum四个站点2013年2月1日的数据为例本文档使用处理环境：win10上的虚拟机2.工程目录设置在任意位置建立工程文件夹全文以test为例，并在其目录下建立如下的文件夹：3.数据下载在终端中打开test，查询2013年2月1日的年积日，和GPS周等信息代码：doy 2013 02 011 在test的目录下再建一个文件夹以032命名的文件夹，此时，test目录下有：2 根据步骤3的信息，到相关下载对应的广播星历，精细星历，观测值文件分别放在文件夹brdc,igs,rinex里面。当然还涉与到下载的数

3、据解压，.d文件到.o文件的转换等操作，其中广播星历可以不改名字即可以不用改成brdc0320.13ns表，更新准如此，下面的表格为准：还有个表文件名表含义更新频次内容更新地址 pole.或极移表每次更新pub/gamit/tables/ ut1.或地球自转参数表每次更新月亮表使用当年的表太阳表使用当年的表章动表使用当年的表跳秒表每次更新gdetic.dat,大地坐标参数文件不常更新antmod.dat,天线相位中心改正正表有新天线时更新svnav.dat,卫星列表有新卫星时更新rcvant.dat,接收机天线对照表有新天线时更新svs_exclude.info坏卫星信息文件每次更新vmf1映

4、射函数模型 使用当年的表 海潮格网数据使用当年的表i. 常见文件说明ii. 更新的准如此是：“更新频次中的红字局部每次处理都更新一下，同时结合本次处理与上次处理的年份是否一样来进展年表的选择。原本有些文件是每周更新一次，但是作为新手为了防止出错还是每次更新吧有些时候，有些文件找不到，或者里面没有，如此可以忽略。如svs_exclude.info，常常找不到我在更新时，这四个文件在服务器上没找到如下图是我更新的文件。原本还应该有2013年的各种年表，但是我之前处理过一次2013年的数据，并且我的/tables已经有了2013年的这些年表。iii. 将上述更新的数据复制到gamit安装目录下的ta

5、bles文件夹内，并替换原有文件。我的gamit安装路径是/opt/gamit，即我安装路径下的tables 的路径是：/tablesiv. 然后将更新后的tables文件夹，/tables拷贝或者设置连接，我习惯拷贝到test目录下。/tables里面放了很多数据，所以它显得特别大1.8G，往往在包含otl.grid的情况下tables也就7800M吧。此后，test的结构应是：这一步相当重要，而且常常容易出错。1） 将.o文件，.sp3文件，*.*n文件观测文件、精细星历、广播星历到032的目录下执行代码：2） 处理数据需要的各种年表等文件此时032中有如下的文件：检查并重新设置断掉的，如

6、如下图所示：如：常说的死，在032目录下，执行代码：ln -s ./tables/luntab.2013.J2000 ./luntab.最后的点号一定要加上系统报错，即“无法创建符号，这句话的意思是，032中早已经有了一个luntab.文件但是我们却看到它是个死，删除该，重新执行上面的代码。新生成，并打开查看该所指向的文件是不是2013年的luntab表重点检查6-ii的11个文件，以与soltab.,nutabl.，用类似的方法建立以下是我手动参加的局部：以下是网上的资料图：其中还要注意设置：，otl.grid，atl.grid的，我在处理数据时使用的如下：在查看原始.grid,如如下图所示

7、：但是我在查看，上传是2010年8月28，大小也是3.1M，故我断定，服务器上的atl。于是我在服务器上下载了一个，然后在032里面制作了一个atl.grid的。当然也可以不用下载atl综上所述：n 更新tables常更新以下11个文件各种年表除外：但常常 guess_rcvant.dat 用不到,可以不更新n 制作特别注意以下几个：n 在更新tables,建立完连接后执行ls al 检查情况，红色的是死：发现atml.grid的已断，于是重新设置即可。至于剩下的红色死，我全是忽略了，因为我也不会设置。好似那个session.info是在正式处理过程中产生的，可以忽略总的来说设置是与其容易出错

8、，要小心又小心。在test/tables下把需要处理的测站参加该表，包括IGS .制作完lfile.文件后不要忘了将lfile.拷贝到test/tables中,同时在032中设置对应lfile.文件是测站近似点坐标越准确越好,10.4以前的版本里面存放的经纬度信息，而10.4与其之后的版本直接存放x,y,z信息GAMIT10.5 的lfile 和apr 文件是同样的格式，所以我们可以认为apr 文件就是lfile。该文件中保存的是各个站的初始坐标，至于该文件中的测站坐标要达到什么样的精度，可参考如下两条：l 如果是该站作为基线解算起算点,如此该测站的坐标做保持在5cm 以内，最好不要超过10c

9、m。l 如果该测站在基线解算中不作为起算点，如此伪距单点定位的精度即可。对于IGS站，获取其在特定框架下的准确坐标方法：方法1：我用的方法可用于非IGS站：在我的电脑上试了几次，老是最后一步出错，不推荐，但也有教程上是这个方法：l 在rinex，.apr，最后生成lfile.依次执行代码：rx2apr lfile.rnx 2013 032gapr_to_l lfile.rnx.apr lfile. * 2013 032方法2：方法3：l 伪距单点定位计算非IGS站的初始坐标实例假设shao是非IGS站：执行完之后会生成两个不同格式的测站坐标文件lfile.shao 和，然后执行l 双差定位计算

10、测站的初始坐标实例：在此例中假设shao站是未知坐标的测站，而bjfs是坐标的测站或起算站综上所述：在生成life.文件时，IGS站的坐标可以直接从框架获得，或者利用双差定位求得。而非IGS站如此可以根据要求精度的上下选择： 直接提取头信息 伪距单点定位 双差定位法 方法1常常在伪距单点定位和双差定位法中选择，然后将所须的坐标信息拷贝到lfile.中即可，并将lfile.拷贝到test/tables中,同时在032中设置对应现在，我假设bjfs，lhaz两站点是精度较高的IGS站，坐标直接从框架中获得，而shao,urum站是非IGS站，近似坐标通过伪距单点定位获得。然后将lfile.拷贝到t

11、est/tables中,同时在032中设置对应好似不需要在032中设置，但是怕出错，我还是设置了，希望没有弄巧成拙 station.info文件的更新最后要在文件在gamit全英文的官方教程版本10.6,24页里面有这样一句话：可见station.Info的可以指向如下图中选中的文件，也就是说倘假如你的test/tables中没有station.info文件，可以自己新建并指向他们。并且这个MIT,SOPAC的文件好似可以在上更新，只是登不上去。练手的话凑合着能用。在test/tables目录下：先建立一个sitelist的文件，并录入测站名字每行首字母空格，不分先后顺序然后终端打开table

12、s并执行代码：显示station.info files have 3 sites，这是错误的，一般来说总共有多少测站参与计算就会有多少个sites，于是回去检查sitelist, ，发现某某站写错了，应该是lhaz所以在操作过程中一定要小心又小心，切记不能将单词写错，同时也要不害怕出现错误，只要根据其提示，顺藤摸瓜，根本上能修改错误的地方，实在不行重来还发现，Recevier name * not found in rcvant.dat,网上百度了一下说要修改，说得又很模糊，所以不会。又去问了一下师兄，他说可以忽略，所以我也就忽略了谁知道的话告诉我一声。执行上述代码完毕后最后一行显示：，。再执

13、行一次更新代码：在生成的station.info文件中会有如下图中红框的字符串，表示更新成功：最后在032中检查station.info的指向，如果没有的话，自己手动用代码建立一个。7 、sittbl.与sestbl.文件它里面包含目录信息，改动较小。需要在032中设置配置 对高精度的坐标强约束，待求点坐标松约束非起算点约束量可以更大,如100m。并且如果IGS站列表中没有你所处理数据中的IGS站，需要自己添加，格式一致。配置sestbl.改正模型文件配置准如此如下面的截图：下面是我的局部设置：我没有met.list,met.grid,也没有.grid,所以我设置的是N不使用我的atml.li

14、st和atl.list是死，我又不知道如何制作他们好似在麻省理工的FTP上可以下载，但是这个FTP好久都登不上去了，于是全选N，而我有atl ，于是选Y修改完了，别忘了保存，同时为了保险起见到032里面再检查一次几个配置文件是否修改、好了至此，数据准备阶段完成。下面我们开始计算。三、分布处理1 .Makexp在032目录下执行代码：makexp2精细星历处理如如下图执行代码：3文件检查执行如下图的代码：对于代码中的gigsf3.032,其中g表示文件类型，中间igsf表示轨道，第一步输入的：，3表示年的最后一位，032是年积日。4广播星历处理其中auto0320.13n没有改名字作为输入文件，

15、jauto3.032是输出文件当然也可以在将auto0320.13n改名字后，执行 make5文件检查6生成中间文件7生成bat文件8基线解算最后一步还是出错了，但是我百度了一下“RC not executed STOP FATAL Error: Stop from report_stat在这篇文章的最后一页找到解决方法，也不知道对不对，反正最后是成功解算了。在032中找到上图的两个文件并删除，然后重新执行基线解算。中间有很多省略结算结果精度还算过得去qexpta.doy:解算记录为的简略版：协方差矩阵、参数平差值9基线解算精度评定标准网上各类教程上的标准某某小异标准1：一般NRMS越小越好，

16、一般小于0.25视为解算效果较好，本次解算的NRMS为0.26，符合要求基线的X、Y方向误差X围应控制在1cm以内，Z方向的误差控制在1.5cm以内，本次解算所有基线均满足要求本次解算x,y,z方向的误差均有超过1cm,结果不是很理想，还过得去最后还有一个基线重复性检验：但是我的检验失败了，是不是要进展两天以上的数据处理才能进展重复性试验？至此可以理解为解算完毕，然后进展基线平差，也可以使用GLBOK平差，但是我选择了用COSAGPS进展的平差。如果想直接用GLBOK进展平差，参考安向东的处理流程。四、cosagps平差1. 在032目录下打开otesta.032,并在第二次出现基线信息的地方

17、参加下面红框中的一行字符串：“cosagps for gamit o-file最后保存并复制到windows桌面下。2. Cosa平差打开cosagps,新建工程，设置如下，其它参数根本默认3. 导入基线数据基线的误差有点大但为了走完流程，有精度强迫症的就凑合着看吧_4. 导入三维坐标这里就是你解算时坐标起算数据，但是要注意框架和坐标系的统一。先在点名方框内随便输入数据，然后确定，如如下图在工程目录下找到刚新建的关于起算点坐标的文件，用记事本打开，并将起算数据复制粘贴到know文件中，如如下图：5. 三维平差最后将解算出的shao,urum站的坐标与原本查询到的坐标进展比照即可。以如下图片里面

18、的坐标就是我在ubuntu里面采用：获得的，比照后发现相差不是太大，但是还是不太理想不管怎么样将软件走通了五、建议与说明六、相关资料天线查询TEQC广播星历下载(以auto开头的文件，下载后改为brdc开头)ftp:/cddisa.gsfc.nasa.gov/pub/gps/data/daily/YYYY/brdc/以brdc开头ftp:/igs.ensg.ign.fr/pub/igs/data/YYYY/DOY/以brdc开头精细星历下载用户名：anonymous，密码：email地址IGS连续站O文件下载d文件，下载后用crx2rnx 解压d文件，下载后用crx2rnx 解压crx2rnx软件压缩D文件转为O文件下载WINDOWS版LINUX版GAMIT的tables文件下载O、N文件下载MATLB TOOLS-gpsg.mit.edu/tah/GGMatlab/43 / 43