汇报—岳彩亚---GAMIT课件

GAMIT数据处理汇报人：岳彩亚 1 GAMIT数据处理1主要内容 Linux操作系统 shell脚本编写 GAMIT基线处理 算例分析 总结2主要内容 2Linux操作系统 1 Linux简介 2 Linux的发展历史 3 Linux系统的特点3Linux操作系统 Linux和Windows（xp、7、8）以及MaxOS一样，都是一种操作系统。整个计算机系统中可以分为四块内容：硬件、内核、系统调用、应用程序，其中内核和系统调用的组合统称为操作系统。如图所示。Linux是一种交互式操作系统。除此之外，还能使得多个用户从不同的终端同时操作主机。也就是说一个Linux系统上可以有多个用户，这些用户之间各自文件，各自的目录，各自管理着自己的任务。这样的一种系统称为分时段操作系统。Linux简介4 Linux和Windows（xp、7、8）以及MaxLinuxLinux简介简介硬件内核系统调用应用程序操作系统5Linux简介硬件内核系统调用应用程序操作系统51、为了研发分时操作系统，在各大学、科研机构以及商业项目都启动了自己的计划，一个项目叫做“Multics”是当时最先进的。2、在Multics开发后期，遇到了一些困难，他的一个开发团体贝尔实验室退出了这个项目，开发自己的操作系统，叫“UNIX”。3、九十年代，UNIX在服务器市场成为了主流的操作系统，开始的时候是免费于使用者，但当用户逐渐增多时，UNIX系统开始变的商业化，并且越来越昂贵Linux发展历史6 1、为了研发分时操作系统，在各大学、科研机构以及商4、1991年10月，LinusTorvalds（Linux之父）开始入手，并且雇佣了一批黑客和一些对操作系统开发有兴趣的计算机专家，开始研究自己的自己的操作系统，也就是Linux系统。LinusTorvalds（芬兰赫尔辛基大学）Linux发展历史7 4、1991年10月，Linus Torvalds（1Linux系统中大都为开源代码并且免费提供给用户，用户可以查看源代码和修改定制以及再发布。2兼具图形界面操作（需要使用带有桌面环境的发行版）和完全的命令行操作，可以只用键盘完成一切操作，并且资源占用少、运行速度快。3Linux系统可以实现多用户共同操作同一个系统，且不会相互影响。Linux的特点8 1 Linux系统中大都为开源代码并且免费提供给主要内容 Linux操作系统 shell简单介绍 GAMIT基线处理 算例分析 总结9主要内容 9shell介绍shell语言的分类shell在GAMIT中的用途10shell介绍10shell可以分为以下几种/bin/sh/bin/dash/bin/bash/bin/rbash/bin/zsh/bin/cshshell介绍11 查看当前系统默认的shell编译语言命令：%echo$SHELL查看当前系统支持的shell编译语言命令：%more/etc/shells改变默认编译语言的命令：%chsh-s/bin/*sh(bash)，重启后即改变成功。shell语言分类12查看当前系统默认的shell编译语言命令：shell语言分类GAMIT是由一些程序模块组成，并且这些程序模块可以实现调用，他们在GAMIT数据处理过程中起到的作用互不相同。使用shell脚本的编写可以将不同的程序模块组合起来实现不同的作用。GAMIT自身带了一些批处理命令，这些命令只是达到固定的结果，通过编写shell脚本语言可以实现新的批处理过程，节约了时间，提高效率。shell在GAMIT用途13 GAMIT是由一些程序模块组成，并且这些程序模块可以主要内容 Linux操作系统 shell简单介绍 GAMIT基线处理 算例分析 总结14主要内容 14GAMIT基线处理：参数文件更新数据处理软件分类分步数据处理流程 自动批处理 结果评定15GAMIT基线处理：针对大区域GPS高精度定位解算，通常采用专业机构开发的高精度科研软件进行处理。目前被广泛使用且较为成熟的GPS数据处理软件可以分为两类：一、基于非差模式的定规定位软件美国航空航天局喷气动力实验室开发的GIPSY/OASIS德国地学研究中心开发的EPOS中国武汉大学GPS研究中心开发的PANDA数据处理软件分类：16 针对大区域GPS高精度定位解算，通常采用专业机二、基于双差模式的定轨定位软件瑞士伯尔尼大学天文学院开发的BERNESE美 国 麻 省 理 工 学 院 和 美 国 斯 克 利 普 斯 海 洋 研 究 所 开 发 的GAMIT/GLOBK数据处理软件分类：17二、基于双差模式的定轨定位软件数据处理软件分类：17每年更新一次：nutabl.soltal.luntab.leap.sec每月更新一次：dcb.dat(编码文件)每日更新一次：ut1（地球自转参数）pole.当出现新的接收机：rcvant.dat新的卫星发射时：svnav.datftp:/garger.usd.edu/pub/gamit/tables（使用迅雷下载）参数文件更新：18每年更新一次：nutabl.soltal.一、创建相关文件夹：下面是创建文件夹的名字，是与控制文件process.defaults中默认的路径是相对应的，如果想修改，必须也把控制文件中默认的名字做相应的修改。处理基线需用到的文件夹：文件夹igs：用来存放SP3精密星历。文件夹brdc：用来存放广播星历，也就是n文件。文件夹rinex：用来存放观测文件，o文件。文件夹met：用来存放气象观测数据文件，一般研究对象是气象方面的时能用的到。文件夹“年纪日”：如果处理的是2012年，年积日001天的数据，就将文件改为001文件夹tables：用来存放控制文件，处理用到的参数文件，即所谓的全局行文件。sh_links.tables-yr|links.tables数据处理流程19一、创建相关文件夹：数据处理流程19与网平差相关的文件夹hfile文件夹：存放文本形式的h文件glbf文件夹：存放二进制形式的h文件gslon文件夹：进行平常操作tables：控制文件，可以直接使用基线解算的tables数据处理流程20 与网平差相关的文件夹数据处理流程20二、IGS核心站数据的获取广播星历：ftp:/garner.ucsd.edu/pub/rinex/YYYY/DOY/(以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开头）sh_get_brdc-archivesopac-yr$year-doy$njr-ndays$ts数据处理流程21二、IGS核心站数据的获取数据处理流程21观测文件（o文件）获取：ftp:/garner.ucsd.edu/pub/rinex/YYYY/DOY/（d文件，下载后用crx2rnx解压）ftp:/cddisa.gsfc.nasa.gov/pub/gps/data/daily/YYYY/DOY/YYo/ftp:/igs.ensg.ign.fr/pub/igs/data/YYYY/DOY/（d文 件，下 载 后 用crx2rnx解压）sh_get_rinex-archive sopac-yr$year-doy$njr-ndays$ts-sites$zm(bjfswuhn)数据处理流程22 观测文件（o文件）获取：数据处理流程22精密星历文件:ftp:/garner.ucsd.edu/pub/products/wwww/ftp:/cddisa.gsfc.nasa.gov/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数据处理流程23三、文件配置（1）制作测站信息文件station.info：拷贝GAMIT/GLOBK安装目录tables/文件夹下的station.info文件至解算目录的tables/文件夹，通过cd命令进入该文件夹（后者），执行：(执行之前一定要先设置station.info的文件的头)%sh_upd_stnfo-files./rinex/*o数据处理流程24三、文件配置数据处理流程24（2）制作台站先验坐标 lfile.执行：%grep POSITION./rinex/*o lfile.rnx（提取o文件中的测站坐标）%rx2apr lfile.rnx$year$doy（转化成GAMIT识别的文件）%gapr_to_l lfile.rnx.apr lfile.$year$doy（将空间坐标转换成球面坐标）这样，就得到了测站信息文件 lfile.（rinex/文件夹下），将其拷贝到tables/文件夹即可。数据处理流程25（2）制作台站先验坐标 lfile.数据处理流程25（3）文件配置数据处理流程sestbl.（策略参数文件）sistbl.（测站约束控制文件）site.default（测站控制文件）process.default（指定文件）26（3）文件配置数据处理流程sestbl.（策略参数文件）s （4）数据连接：cd命令进入$doy（001）文件夹，执行命令：%ln-s./rinex/*o./#链接观测数据%ln-s./igs/*sp3./#链接精密星历%ln-s./brdc/*n./#链接广播星历%ln-s./met/*m./#链接气象数据%sh_links.tables J2000$year#链接表文件数据处理流程33 （4）数据连接：数据处理流程33（5）解算：%makexp#数据预处理%sh_sp3fit-f*sp3-o igsf-d$year$doy#使用精密星历生成轨道参数%makej brdc*n jbrdc$y.$doy#生成卫星钟文件%makex$expt.makex.batch#创建观测数据%fixdrv d$expt$y.$doy#启动数据驱动程序%csh b*bat#执行数据解算数据处理流程34数据处理流程34（1）创建文件：RINEX文件夹（rinex igs brdc）MET文件夹 STABLES文件夹（2）解算文件配置：sestbl.解算策略配置文件 sittbl.测站约束控制文件 process.defaults 参数指定 sites.defaults站点控制文件 *.apr固定坐标（可以用L文件）station.info测站信息文件 自动批处理35（1）创建文件：自动批处理35（3）开始解算：%sh_gamit-expt$expt-d$year days sh_gamit自动批处理程序会自动创建其他所需文件夹，并进行数据解算。解算结果将保存到年积日文件夹。执行这个命令后，会自动创建一些文件夹，创建的这些文件夹被使用过后有的留了下来，有的被删掉。自动批处理36自动批处理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 在输出文件中，在q$expta.$doy($expt主要内容 Linux操作系统 shell简单介绍 GAMIT基线处理 GLOBK网平差 算例分析 总结38主要内容 38文件准备参数配置平差过程精度评定GLOBK网平差步骤：39GLOBK网平差步骤：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 文件准备：40 globk_comb.cmd 文 件 对 平 差 过 程 进 行 控 制，编 辑 globk_comb.cmd 文件以满足任务要求，主要有以下项：参数配置：41 globk_comb.cmd 文件对平差过程进行控 glorg_comb.cmd 文件对固定点进行控制，以将平差结果纳入到一定的坐标框架。编辑 glorg_comb.cmd 文件以满足任务要求，主要有以下项参数配置：42 glorg_comb.cmd 文件对固定点进行控制，分步处理：（1）将GAMIT的H文件（或其他解算数据）转为GLOBK认可的二进制的H文件：%cd hfile/%htoglb./glbf/./tables/svnav.dat h?a*将文本形式的H文件转化为二进制形式（扩展名为glr或glx），并存放到glbf/文件夹。（2）创建H文件列表：%cd./gsoln/%ls./glbf/h*glx$expt_glx.gdl将所有需解算的H文件的路径和文件名列入gdl文件。（3）创建待平差测站的先验坐标apr文件：文件应包括测站的三维坐标+站速度+历元也可使用现有坐标框架的值，如itrf08_comb.apr（以该坐标框架为例）平差过程：43分步处理：平差过程：43（4）在gsoln/文件夹下执行：%glred 6 glred.opt glred.log$expt_glx.gdl globk_comb.cmd生成测站坐标的时间序列，并进行测站重复性检验。（5）执行命令，进行自由网平差：%globk 6 globk.prt globk.log$expt_glx.gdl globk_comb.cmd（6）执行命令：%glorg 6 glorg.org glorg.log$expt_glx.gdl glorg_comb.cmd进行约束平差，将测站纳入到所选的坐标框架。平差过程：44（4）在gsoln/文件夹下执行：平差过程：44自动批处理：（1）配置文件：主要是配置sites.defaults globk_comb.cmd（如上所述）glorg_comb.cmd（如上所述）（2）执行批处理命令：%sh_glred-d$year$doy-expt$expt -opt H G E平差过程：45自动批处理：平差过程：45结果输出文件：log文件，解算过程记录 prt文件，自由网平差输出结果 org文件，约束网平差输出结果 SUM文件，解算结果总结网平差阶段的精度评定：RMS值，均方根误差 测站坐标重复性图 已知点坐标符合度（内符合和外符合）精度评定：46结果输出文件：精度评定：46主要内容 Linux操作系统 shell简单介绍 GAMIT基线处理 GLOBK网平差 算例分析 总结47主要内容 47 取中国境内的六个IGS核心站，分别为BJFS、WUHN、CHAN、SHAO、KUNM、URUM。这六个站的分布如图所示。（GUAO、XIAN）（TAIW、TNML）算例分析48 取中国境内的六个IGS核心站，分别为BJFS、WUH 其中以CHAN、SHAO、KUNM、URUM为约束点，对BJFS、WUHN进行约束，实现这种处理方式，主要是在sisbl.中设置，如图二。算例分析49 其中以CHAN、SHAO、KUNM、URUM为约束点，精度评定：1、q文件中的nram值。算例分析50精度评定：算例分析50 2、基线误差不能大于约束值的2倍：算例分析51 2、基线误差不能大于约束值的2倍：算例分析51 3、控制站点的情况：算例分析52 3、控制站点的情况：算例分析52 1、基线解算结果：算例分析53 1、基线解算结果：算例分析53 1、平差解算结果：(globk.org)算例分析54 1、平差解算结果：(globk.org)算例分析541、GAMIT数据处理软件：（1）版本 （2）掌握重点2、LINUX系统：3、shell语言：shell语言是一种Linux系统中的脚步编写语言，学习shell语言时主要的是在学习里面的命令。在使用shell编写脚本时就是一些简单的命令组合，如果只是为了使用GAMIT而学习，可以简单的掌握一些常用命令即可。总结：注册一个账号Linux入门基础（Java、C、C+）实验楼第一家以实验为核心的IT在线教育平台实验楼551、GAMIT数据处理软件：总结：注册一个账号Linux入门谢谢观看56谢谢观看56