汇报-岳彩亚-GAMIT课件

GAMIT数据处理汇报人：岳彩亚

1主要内容

●Linux操作系统

●shell脚本编写

●

GAMIT基线处理

●算例分析

●总结2Linux操作系统

1Linux简介2Linux的发展历史

3Linux系统的特点3Linux和Windows（xp、7、8）以及MaxOS一样，都是一种操作系统。整个计算机系统中可以分为四块内容：硬件、内核、系统调用、应用程序，其中内核和系统调用的组合统称为操作系统。如图所示。Linux是一种交互式操作系统。除此之外，还能使得多个用户从不同的终端同时操作主机。也就是说一个Linux系统上可以有多个用户，这些用户之间各自文件，各自的目录，各自管理着自己的任务。这样的一种系统称为分时段操作系统。

Linux简介4Linux简介硬件内核系统调用应用程序操作系统51、为了研发分时操作系统，在各大学、科研机构以及商业项目都启动了自己的计划，一个项目叫做“Multics”是当时最先进的。2、在Multics开发后期，遇到了一些困难，他的一个开发团体—贝尔实验室退出了这个项目，开发自己的操作系统，叫“UNIX”。3、九十年代，UNIX在服务器市场成为了主流的操作系统，开始的时候是免费于使用者，但当用户逐渐增多时，UNIX系统开始变的商业化，并且越来越昂贵Linux发展历史64、1991年10月，LinusTorvalds（Linux之父）开始入手，并且雇佣了一批黑客和一些对操作系统开发有兴趣的计算机专家，开始研究自己的自己的操作系统，也就是Linux系统。

LinusTorvalds（芬兰赫尔辛基大学）Linux发展历史71Linux系统中大都为开源代码并且免费提供给用户，用户可以查看源代码和修改定制以及再发布。2兼具图形界面操作（需要使用带有桌面环境的发行版）和完全的命令行操作，可以只用键盘完成一切操作，并且资源占用少、运行速度快。3Linux系统可以实现多用户共同操作同一个系统，且不会相互影响。Linux的特点8主要内容

●Linux操作系统

●

shell简单介绍

●GAMIT基线处理

●算例分析

●总结9shell介绍shell语言的分类shell在GAMIT中的用途10

shell可以分为以下几种

/bin/sh/bin/dash

/bin/bash/bin/rbash/bin/zsh

/bin/csh

shell介绍11查看当前系统默认的shell编译语言命令：%echo$SHELL查看当前系统支持的shell编译语言命令：%more/etc/shells改变默认编译语言的命令：%chsh-s/bin/*sh(bash)，重启后即改变成功。shell语言分类12GAMIT是由一些程序模块组成，并且这些程序模块可以实现调用，他们在GAMIT数据处理过程中起到的作用互不相同。使用shell脚本的编写可以将不同的程序模块组合起来实现不同的作用。GAMIT自身带了一些批处理命令，这些命令只是达到固定的结果，通过编写shell脚本语言可以实现新的批处理过程，节约了时间，提高效率。shell在GAMIT用途13主要内容

●Linux操作系统

●shell简单介绍

●

GAMIT基线处理

●算例分析

●总结14GAMIT基线处理：

参数文件更新数据处理软件分类分步数据处理流程

自动批处理结果评定15针对大区域GPS高精度定位解算，通常采用专业机构开发的高精度科研软件进行处理。目前被广泛使用且较为成熟的GPS数据处理软件可以分为两类：一、基于非差模式的定规定位软件美国航空航天局喷气动力实验室开发的GIPSY/OASIS德国地学研究中心开发的EPOS中国武汉大学GPS研究中心开发的PANDA数据处理软件分类：16二、基于双差模式的定轨定位软件瑞士伯尔尼大学天文学院开发的BERNESE美国麻省理工学院和美国斯克利普斯海洋研究所开发的GAMIT/GLOBK数据处理软件分类：17每年更新一次：nutabl.soltal.luntab.leap.sec每月更新一次：dcb.dat(编码文件)每日更新一次：ut1（地球自转参数）pole.当出现新的接收机：rcvant.dat新的卫星发射时：svnav.datftp:///pub/gamit/tables（使用迅雷下载）参数文件更新：18一、创建相关文件夹：下面是创建文件夹的名字，是与控制文件process.defaults中默认的路径是相对应的，如果想修改，必须也把控制文件中默认的名字做相应的修改。处理基线需用到的文件夹：文件夹igs：用来存放SP3精密星历。文件夹brdc：用来存放广播星历，也就是n文件。文件夹rinex：用来存放观测文件，o文件。文件夹met：用来存放气象观测数据文件，一般研究对象是气象方面的时能用的到。文件夹“年纪日”：如果处理的是2012年，年积日001天的数据，就将文件改为001文件夹tables：用来存放控制文件，处理用到的参数文件，即所谓的全局行文件。sh_links.tables-yr~|links.tables数据处理流程19与网平差相关的文件夹hfile文件夹：存放文本形式的h文件glbf文件夹：存放二进制形式的h文件gslon文件夹：进行平常操作

tables：控制文件，可以直接使用基线解算的tables数据处理流程20二、IGS核心站数据的获取广播星历：/pub/rinex/YYYY/DOY/(以auto开头的文件，下载后改为brdc开头)/pub/gps/data/daily/YYYY/brdc/（以brdc开头）ftp://igs.ensg.ign.fr/pub/igs/data/YYYY/DOY/（以brdc开头）

sh_get_brdc-archivesopac-yr$year-doy$njr-ndays$ts······数据处理流程21观测文件（o文件）获取：/pub/rinex/YYYY/DOY/（d文件，下载后用crx2rnx解压）/pub/gps/data/daily/YYYY/DOY/YYo/ftp://igs.ensg.ign.fr/pub/igs/data/YYYY/DOY/（d文件，下载后用crx2rnx解压）

sh_get_rinex-archivesopac-yr$year-doy$njr-ndays$ts-sites$zm(bjfswuhn)······数据处理流程22精密星历文件:/pub/products/wwww//pub/gps/products/WWWW/ftp://igs.ensg.ign.fr/pub/igs/products/WWWW/sh_get_orbits-archivesopac-yr$year-doy$njr-ndays$ts······将下载后的*sp3文件，*n文件，*o文件拷贝到对应的创建文件夹中。数据处理流程23三、文件配置（1）制作测站信息文件：

拷贝GAMIT/GLOBK安装目录tables/文件夹下的文件至解算目录的tables/文件夹，通过cd命令进入该文件夹（后者），执行：(执行之前一定要先设置的文件的头)

%sh_upd_stnfo-files../rinex/*o数据处理流程24（2）制作台站先验坐标lfile.执行：

%grepPOSITION../rinex/*o>lfile.rnx（提取o文件中的测站坐标）

%rx2aprlfile.rnx$year$doy（转化成GAMIT识别的文件）

%gapr_to_llfile.rnx.aprlfile.$year$doy（将空间坐标转换成球面坐标）这样，就得到了测站信息文件lfile.（rinex/文件夹下），将其拷贝到tables/文件夹即可。数据处理流程25（3）文件配置数据处理流程sestbl.（策略参数文件）sistbl.（测站约束控制文件）site.default（测站控制文件）process.default（指定文件）26sestbl.27sestbl.28总结：

参数控制文件为数据处理方案的核心控制文件，其中主要的参数包括处理模式、观测量的选取、批处理迭代方案、对流层误差模型、天顶延迟个数、光压模型、数据筛选、ARC参考系统、惯性框架、L坐标更新。sestbl.29sistbl.30siste.default删除“#”，输入：chan$exptftprnxxstinfoglrepuglrepsglts31process.default32（4）数据连接：

cd命令进入$doy（001）文件夹，执行命令：

%ln-s../rinex/*o./ #链接观测数据

%ln-s../igs/*sp3./ #链接精密星历

%ln-s../brdc/*n./ #链接广播星历

%ln-s../met/*m./ #链接气象数据

%sh_links.tablesJ2000$year

#链接表文件数据处理流程33（5）解算：

%makexp

#数据预处理

%sh_sp3fit-f*sp3-oigsf-d$year$doy#使用精密星历生成轨道参数

%makejbrdc*njbrdc$y.$doy #生成卫星钟文件

%makex$expt.makex.batch #创建观测数据

%fixdrvd$expt$y.$doy

#启动数据驱动程序

%cshb*bat

#执行数据解算数据处理流程34（1）创建文件：

RINEX文件夹（rinexigsbrdc）MET文件夹STABLES文件夹（2）解算文件配置：

sestbl.

解算策略配置文件

sittbl.

测站约束控制文件

process.defaults

参数指定

sites.defaults

站点控制文件

*.apr 固定坐标（可以用L文件）

测站信息文件

自动批处理35（3）开始解算：

%sh_gamit-expt$expt-d$yeardays

sh_gamit自动批处理程序会自动创建其他所需文件夹，并进行数据解算。解算结果将保存到年积日文件夹。执行这个命令后，会自动创建一些文件夹，创建的这些文件夹被使用过后有的留了下来，有的被删掉。

自动批处理36在输出文件中，在q$expta.$doy($expt为工程名。$doy为解算的年积日，如下面的算例用的是2012年年积日为001天，工程名为gycy因此文件名为qgycya.001。)和o$expta.$doy，o*是q*的缩略版，在autcln.post.sum中可以看控制网点的相对质量。主要指标：（1）查看q*文件中的，后验nrms值，最好在0.25附近，一般不可以超过0.3。（2）基线分量误差的大小，参数的改正量不应大于约束量的2倍。

（3）*.sum总结文件中的rms值，看控制网点的质量。（4）基线向量的重复性：

结果分析37主要内容

●Linux操作系统

●shell简单介绍

●GAMIT基线处理●

GLOBK网平差

●算例分析

●总结38文件准备参数配置平差过程精度评定GLOBK网平差步骤：39

拷贝GAMIT/GLOBK安装目录下的tables/文件夹下的pmu.usno（地球自转参数）+itrf08_comb.apr+itrf08_comb.eq（地震信息）文件到解算目录下的tables/文件夹。

拷贝GAMIT/GLOBK安装目录下的tables/文件夹下的globk_comb.cmd+glorg_comb.cmd文件到解算目录下的gsoln/文件夹。

对globk_comb.cmd文件和glorg_comb.cmd文件进行编辑以适合你的解算任务。文件准备：40

globk_comb.cmd文件对平差过程进行控制，编辑globk_comb.cmd文件以满足任务要求，主要有以下项：参数配置：参数意义说明eq_file地震信息文件目录（目录全部指到文件存放处）apr_file先验坐标文件目录in_pmu地球自转参数文件目录apr_neu测站坐标值约束41

glorg_comb.cmd文件对固定点进行控制，以将平差结果纳入到一定的坐标框架。编辑glorg_comb.cmd文件以满足任务要求，主要有以下项参数配置：参数意义说明apr_file先验坐标文件目录stab_siteclear清除之前的固定站标记(有这几个字就行)stab_site固定站点（在“上”个后添加，在将那些固定的站写在后面即可。“这两行的意思是先将itrf08z8888中的固定坐标全部删除，然后再从下一行中指定固定坐标”）42分步处理：（1）将GAMIT的H文件（或其他解算数据）转为GLOBK认可的二进制的H文件：%cdhfile/%htoglb../glbf/../tables/svnav.dath????a*将文本形式的H文件转化为二进制形式（扩展名为glr或glx），并存放到glbf/文件夹。（2）创建H文件列表：%cd

../gsoln/%ls../glbf/h*glx>$expt_glx.gdl将所有需解算的H文件的路径和文件名列入gdl文件。（3）创建待平差测站的先验坐标apr文件：文件应包括测站的三维坐标+站速度+历元也可使用现有坐标框架的值，如itrf08_comb.apr（以该坐标框架为例）平差过程：43（4）在gsoln/文件夹下执行：%glred6glred.optglred.log$expt_glx.gdlglobk_comb.cmd生成测站坐标的时间序列，并进行测站重复性检验。（5）执行命令，进行自由网平差：%globk6globk.prtglobk.log$expt_glx.gdlglobk_comb.cmd（6）执行命令：%glorg6glorg.log$expt_glx.gdlglorg_comb.cmd进行约束平差，将测站纳入到所选的坐标框架。平差过程：44自动批处理：（1）配置文件：主要是配置sites.defaults

globk_comb.cmd（如上所述）

glorg_comb.cmd（如上所述）（2）执行批处理命令：%

sh_glred-d$year$doy-expt$expt-optHGE平差过程：45结