GAMIT GLOBK软件技术应用高级培训班ppt课件

2020 4 7 1 GAMIT GLOBK软件技术应用主讲人 周命端博士E mail zmd zry GPS测量技术工程应用高级培训班 第二期 1 中国 北京2014年8月16 19日 GAMIT GLOBK GMT TRACK 2 2020 4 7 2020 4 7 3 内容提要 GAMIT GLOBK GMT TRACK软件安装与操作专题实操练习 应用GAMIT CosaGPS软件进行高精度GPS工程控制网数据处理与精度评估综合实操练习 应用GAMIT GLOBK GMT TRACK软件进行CORS站网数据分析高精度GPS数据处理 技术沙龙 第一部分 GAMIT CosaGPS 4 2020 4 7 2020 4 7 1操作系统 主要内容虚拟机 VMwareWorkstation 的介绍 下载 安装与基本使用UBUNTU操作系统的下载 安装与基本使用 cd ls ln gedit chmod命令与操作 基于VMwareWorkstation的UBUNTU操作系统的安装 5 1 1虚拟机 VMwareWorkstation 虚拟工作站 VMwareWorkstation VMwareWorkstation是一款功能强大的桌面虚拟计算机软件下载地址 最新版本为10 02 2014年4月17日发行 2020 4 7 6 1 2UBUNTU操作系统 乌班图 Ubuntu Ubuntu是一个以桌面应用为主的Linux操作系统下载地址 最新版本是14 04LTS 2014年4月18日发行 ls 列出目录内容 cd 切换当前目录 cd cd gedit 用文本编辑器打开文件 ln 创建文件链接 ln s chmod 添加执行属性 chmod x 2020 4 7 7 1 3操作系统安装过程 基于VMwareWorkstation的UBUNTU操作系统的安装演示 8 2020 4 7 2GAMIT GLOBK软件介绍与安装 主要内容GAMIT GLOBK软件介绍 下载方式UBUNTU操作系统下安装GAMIT GLOBK软件的系统配置 csh gfortran libX11 dev 手把手 教会GAMIT GLOBK软件的安装GAMIT GLOBK软件的路径配置 bashrc 与测试 doy 9 2020 4 7 2 1GAMIT GLOBK软件介绍与下载 GAMIT GLOBK安装包GAMIT GLOBK是一套高精度数据处理软件 由美国麻省理工学院 MIT 和斯克里普斯海洋研究所 SIO 联合开发 下载地址 ftp chandler mit edu 最新版本是10 50 2014年4月18日Release 2020 4 7 10 2 2操作系统配置 GAMIT GLOBK安装准备工作更换源 系统 系统管理 新立得软件包管理器 设置 软件库打开终端 用户名目录下执行如下命令 sudoapt getinstallupdatesudoapt getinstallcshsudoapt getinstallgfortransudoapt getinstalllibX11 dev 2020 4 7 11 连网 2 3安装 GAMIT GLOBK安装安装目录下 运行 修改 libraries Makefile config文件 install software Generic willworkonanysystemiflinksinplace X11LIBPATH usr libX11INCPATH usr include X11 MAXSIT60MAXSAT32MAXATM13MAXEPC2880 MAXSIT60MAXSAT32MAXATM25MAXEPC5760 2020 4 7 12 2 4路径配置 GAMIT GLOBK安装的路径配置修改 bashrc文件用户名目录下 测试doy 检查是否安装成功 exportPATH PATH home mdzhou gamit source20120815 gamit bin home mdzhou gamit source20120815 com home mdzhou gamit source20120815 kf binexportsetenvHELP DIR home mdzhou gamit source20120815 help gedit bashrc 2020 4 7 13 3分步式 玩转 GAMIT软件 主要内容分步操作GAMIT软件 分步操作的基本命令使用 标准RINEX文件准备 算例 某工程控制网GPS数据4个观测时段 国标D级 关键文件准备 sestbl station info sittbl lfile lfile apr 与详细介绍高精度GPS基线处理 手把手 教分步操作命令 基线质量评估方法 四看 14 2020 4 7 3 1分步操作GAMIT软件 15 2020 4 7 数据准备 观测文件标准化 原始文件 RINEXo files 采用随机格式转换软件或TEQC软件可将原始文件转换为标准的rinexo files格式对所有的o files文件的命名进行标准化 如sitedoyn yyo 注意文件名字字母要小写 2020 4 7 16 天线高改化至DHARP 天线座底部 外业天线高量测的方式根据不同的接收机天线型号有不同的量测方法GAMIT软件里设置了不同的量高模式 建议采用天线座底部 DHARP 天线高改化至L1PHC 天线相位中心 2020 4 7 17 2020 4 7 18 GPS接收机天线高改化原理 量高模式与改化方法 GPS接收机天线量高模式 2020 4 7 19 GPS接收机天线高改化原理H为标石中心 天线参考点的斜高R为天线盘半径h1为标石中心 天线参考点的垂高h2为天线参考点 天线座底部的垂高 h3为天线参考点 天线相位中心的垂高标石中心 天线相位中心的垂高 GPS测量技术工程应用培训班 2020 4 7 20 已知外业测量的斜高为1 564m求 从网上查询或下载数据文件 2020 4 7 21 2020 4 7 IGS跟踪站o files数据 导航电文brdc文件及精密星历SP3文件 ftp garner ucsd edu pub ftp cddis gsfc nasa gov pub gps ftp igs ensg ign fr pub 更新tables文件 可以从服务器ftp garner ucsd edu pub gamit下载更新的tables表文件其中luntab soltab nutabl 等文件按年进行更新pole ut1 pmu 等文件每周都有更新 22 2020 4 7 文件结构整理 对每个工程目录做处理前的准备工作 将数据文件分类存放 便于维护和更新步骤如下 在当前工程目录下建立的文件夹有 igs rinex brdc以及时段文件夹igs文件夹存放精密星历文件rinex文件夹存放o files文件brdc文件夹存放brdc文件时段文件夹 以时段号命名 所有解算都是基于该文件夹进行的 文件准备 23 2020 4 7 将更新好的完好tables表文件复制到当前工程目录下可以手动复制或命令cp r tables 注意 这里尽量用复制 当然 也可以用连接ln s gamit source tables 至此 当前工程目录文件夹结构应该为 24 2020 4 7 修改tables文件夹下的配置文件station info文件 记录各测站的站点 站点名称 开始观测时间 结束观测时间 天线高 天线量测方法 接收机类型 软硬件版本号 天线类型 天线相位中心模型机构代码等信息 注意 AntHt HtCod 标志天线量测方法DHARP SLBGP等 和AntennaType 用来做天线高改正 ReceiverType 用来做周跳恢复 SessionStart SessionStop 用来控制测段跨度 处理不好后面会提示错误 25 2020 4 7 严格按照已经给出的跟踪站格式编辑具体操作方法有2种 手动操作 按照相应文件格式对齐操作 按照标准格式station info standard编辑修改 sh upd stnfo自动提取观测文件的有关信息 生成station info 并认真核对信息是否正确 事先提供station info空模板后 再执行命令 sh upd stnfo refstation info files yyo 26 2020 4 7 l file文件 记录各测站的近似坐标Gamit要求是球面坐标 可以先在测站的O文件里手动提取近似坐标XYZ 生成l file文件 按照格式对齐 然后保存文件 将名字重新命名为lfile apr 然后用gapr to l将lfile apr转换到lfile 可以用sh rx2apr命令自动生成 利用svpos和svdiff模块单点定位生成站点坐标 也可以在测站的O文件里自动提取近似坐标 grepPOSITION yyo lfilerx2aprlfile gapr to llfile aprlfile 27 2020 4 7 sestbl 文件gambit基线处理时的模型选择配置文件 需要仔细参照技术文档和程序代码 一般情况下采用默认值 通常需要修改的地方有 ChoiceofObservable 观测值的选择类型 ChoiceofExperiment 基线处理类型 sittbl 文件 sittbl short 各测站的精度控制指标 一般情况下采用默认值 高精度的已知坐标强约束 待求点坐标松弛约束 28 2020 4 7 将tables表文件中需要的文件 包括luntab soltab nutabl pole ut1 pmu文件 还有一些误差改正模型文件等 链接到各时段文件夹 该步骤的作用是 在数据处理各个时段的工作目录下 建立该工作目录与 tables目录中部分文件的链接 为后面在工作目录下逐行输入命令进行数据解算做准备 在各时段文件夹下分别建立与tables文件夹里相应文件对应的链接 输入命令 links dayyrdoytest 可以参考links day命令的用法 29 2020 4 7 然后需要手动建立以下链接 命令如下ln s tables luntab 2009 J2000 luntab ln s tables soltab 2009 J2000 soltab ln s tables nutabl 2009 nutabl ln s tables pole usno pole ln s tables ut1 usno ut1 ln s tables pmu usno pmu 另外 潮汐改正文件也是很容易出错的文件 以gamit10 50版本为例 Gamit默认的潮汐改正文件名为otl grid 链接前是otl FES2004 grid 30 2020 4 7 一定要检查一下otl grid是否是真实链接 出现死链接的情况需要手动链接 命令为 在tables文件夹目录下 ln s otl FES2004 grid otl grid链接各个测站观测数据在各时段文件夹下建立观测数据连接 在各时段文件夹下输入命令 ln s rinex doy yyo ln s igs igsweek sp3 ln s brdc doy yyn 31 2020 4 7 数据文件准备分步操作 在终端 输入命令 makexp 回车 然后按照提示输入相应的文件数据 提示在要求输入l file文件时候 建议输入lfile apr文件 批操作sh makexp expttest orbtigsf yr2009 doy088 sess99 srin navbrdc0880 09n aprlfile sinfo1500005760需要准备的文件 O文件 N文件 station info文件 rcvant dat文件 接收机天线类型 输出文件 session info 数据处理 分步处理 32 2020 4 7 精密星历处理sh sp3fit figs15250 sp3 oigsf d2009088 rBERNE t u需要准备的文件sp3文件 ut1 时长变化表 pole 极移表 nutabl 章动表 soltbl 太阳星历表 luntab 月亮星历表 leap sec 跳秒 svnav dat 卫星天线类型 33 2020 4 7 输出文件gigsf9 088 说明 g files文件所有卫星初始轨道坐标文件tigsf9 088 文件检查sh check sess sess088 typegfile filegigsf9 088需要准备的文件 gtest9 088 输出文件 session info广播星历处理makejbrdc0880 09njbrdc9 088需要准备的文件 brdc0880 09n 输出的文件 jbrdc9 088 34 2020 4 7 文件检查sh check sess sess088 typejfile filejbrdc9 088需要准备的文件jbrdc9 088输出文件session info形成中间文件makextest makex batch需要准备的文件test makex batch 以及该文件所列出的所有文件输出文件dtest9 088 35 2020 4 7 生成bat文件输入命令fixdrvdtest9 088需要准备的文件 dtest9 088 以及该文件所列出的所有文件输出的文件 btest9 088 bat 以及btest9 num基线解算输入命令cshbtest9 bat需要准备文件 btest9 bat 以及btest9 num输出文件 q o h 等所有解算的结果文件 36 2020 4 7 GAMIT基线处理结果质量评价指标检查所有测站点是否都参与计算 LC HELP 标准均方根nrms解算结果为0 20左右若是nrms大于0 5 表明处理结果是有问题的 例如 周跳没有修复 测站的起算点坐标有问题等 若是nrms小于0 5 则认为成功 合理 解算 否则需要检查原因 重新处理 检查各基线文件处理精度是否符合gamit的实际精度 质量评估 37 4定制式 玩转 GAMIT软件 主要内容bash编程语言介绍基于bash编程语言自主编写 定制化 批处理脚本 sh 算例 某工程控制网GPS数据4个观测时段 国标D级 sh文件 38 2020 4 7 2020 4 7 使用GAMIT软件进行基线解算有一定程度的了解之后 可以不用每次只执行一条命令 可利用UNIX LINX系统自带的shell编程语言 定制适合自己的批处理脚本Shell命令 值得注意的是 Shell脚本命令只对数据处理进行批处理 数据准备及文件准备建议手动操作 以某城市CORS2012年doy 201第0时段基线解算为例 定制的批处理脚本文件2010 cs sh内容如下所示 39 2020 4 7 bin csh f Step1 customizedyourexperimentsetexpt njcs setyr 4 2012 setyr 2 12 setyr 1 2 setdoy 201 setsp3type igsr setsp3file igr16974 sp3 setbrdcfile brdc2010 12n setsession 0 40 2020 4 7 step2 linkthetablesintocurrentdirectorylinks day yr 4 doy expt step3 removebrokenfilesandrestart linkfilesintocurrentdirectoryrmluntab rmsoltab rmnutabl rmpole rmut1 rmpmu 41 2020 4 7 ln s tables luntab 2012 J2000 luntab ln s tables soltab 2012 J2000 soltab ln s tables nutabl 2012 nutabl ln s tables pole usno pole ln s tables ut1 usno ut1 ln s tables pmu usno pmu 42 2020 4 7 step4 removerinexfilesandrestrart linkrinexfilesintocurrentdirectoryrm 12orm sp3rm 12nln rinex doy session doy session 12o ln brdc brdcfile ln igs sp3file 43 2020 4 7 step5 preparationfortheexperimentrm fatalrmsession infosh makexp expt expt orbt sp3type yr yr 4 doy doy sess99 srin nav brdcfile aprl expt yr 1 doy sinfo3000002880sh sp3fit f sp3file o sp3type d yr 4 doy rBERNE t ush check sess sess doy typegfile fileg sp3type yr 1 doymakej brdcfilejbrdc yr 1 doysh check sess sess doy typejfile filejbrdc yr 1 doymakex expt makex batchfixdrvd expt yr 1 doycshb expt yr 1 bat 44 2020 4 7 step6 resultsofbaselinesolutionsif e expt kjw mkdir expt kjwcp r q expt a doy expt kjwmv expt kjw q expt a doy expt kjw q doy session a doy cp r o expt a doy expt kjwmv expt kjw o expt a doy expt kjw o doy session a doy cp r h expt a yr 2 doy expt kjwmv expt kjw h expt a yr 2 doy expt kjw h doy session a doy 45 2020 4 7 GAMIT软件基线解算流程图 46 技术讨论 经验交流会 互动环节与答疑 2020 4 7 47 5批处理 玩转 GAMIT软件 主要内容 算例 中国大陆8个IGS站 CHAN BJFS URUM GUAO LHAZ WUHN SHAO KUNM 2013年doy 001 005高精度基线处理 sh gamit 批处理 sh gamit 介绍与工程应用使用前提sh gamit操作GAMIT软件的基本流程 48 2020 4 7 算例 中国大陆8个IGS站 GMT绘图命令 bin cshpsxyChina Province dat Jl103 35 28 40 1 28000000 R64 138 18 54 m K China IGS pscatchina igs txt awk print 5 4 psxy R J O K Sc0 2 Gred China IGS pscatchina igs txt awk print 5 4 12 0 7 BR 1 pstext R J O Gblue V China IGS ps sh gamit Step1 在 2013目录下分别建立igs brdc rinex文件夹 并存放对应数据 并在 2013目录下 终端执行 sh setup yr2013编辑 2013 tables process defaults编辑 2013 tables sites defaults准备station info准备lfile apr和lfile 编辑sestbl 和sittbl 在 2013目录下 终端执行 d 所处理数据所在的年 年积日 s 2013001003 copt 待压缩的文件类型 dopt 待删除的文件类型 yrext 指定待处理数据年份及年积日 pres NO YES和ELEV sh gamit exptcors d2013001002003 presELEV orbitIGSF coptxkp doptscao 专题实操练习 应用GAMIT CosaGPS软件进行高精度GPS工程控制网数据处理与精度评估 53 2020 4 7 主要内容 基于IGS站的框架网构建与优化设计IGS站点数据 GPS观测资料 三维框架坐标 获取途径大地测量参考框架 ITRF系列框架与参考历元的基准统一及其转换 全球定位系统 GPS 测量规范 GB T18314 2009 解读 算例 某A级GPS变形监测网高精度数据处理中框架网构建 54 2020 4 7 利用 定制式 脚本操作GAMIT软件进行高精度GPS基线处理 可获得高精度GPS基线向量 天顶对流层延迟参数 高精度GPS基线向量质量评估 算例 AB级GPS基线分量 长度重复性计算应用GLOBK CosaGPS软件进行GPS网平差与精度评估数据计算报告与技术总结报告撰写 55 2020 4 7 1基于IGS站的框架网优化设计与构建 IGS站网 2020 4 7 56 SOPAC网站 http sopac ucsd edu 2020 4 7 57 IERS网站 www iers org 2020 4 7 58 2ITRF系列框架与参考历元的基准统一 2020 4 7 59 2020 4 7 60 3全球定位系统 GPS 测量规范解读 术语和定义基本规定级别划分和测量精度布设的原则选点埋石仪器观测外业成果记录数据处理成果验收与上交资料 61 2020 4 7 2020 4 7 62 3 1基本规定 GPS测量采用2000国家大地坐标系统采用GPS时间系统手簿记录宜采用世界协调时 UTC 接收机应经检定合格 在有效期内使用采用中误差为精度指标 以2倍为极限误差当需要提供80或54或其他坐标成果时 应按坐标转换方法求得 2020 4 7 3 2级别划分和测量精度 级别划分GPS测量按照精度和用途分为A E级5等级测量精度A级GPS网由CORS站构成B E级的精度应不低于表2要求 63 2020 4 7 64 1 按工作原理分 码相关型接收机平方型接收机混合型接收机 2 按信号通道的类型分 多通道接收机序贯通道接收机多路复用通道接收机 3 根据接收信号的频率分 单频接收机 L1 双频接收机 L1 L2 4 根据接收机的用途分 导航型测量型授时型 3 3仪器 GPS接收机的分类 2020 4 7 接收机选用A级网测量按有关规定执行B E级GPS网按表4规定执行GPS接收机的检验新购置GPS接收机的检验一般性检视 通电检验 试测检验旧GPS接收机的检验一般性检视 通电检验 试测检验的准备部分 65 2020 4 7 66 3 4观测 A级网B E级GPS网观测的基本技术规定应符合表5的要求 2020 4 7 3 5外业成果记录 记录类型观测数据 测量手簿 其他记录记录内容原始观测数据 对应观测值的GPS时间 测站与接收机初始信息 天线高记录要求及时 认真 准确 67 2020 4 7 3 6数据处理 原始格式 RINEX格式方法GPS数据文件的组织外业数据质量检验数据删除率 不宜大于10 复测基线的长度差 68 2020 4 7 69 各级GPS测量的精度指标 相邻点间基线长度标准差例如 a 5mmb 1ppmD 5km 2020 4 7 70 外业数据质量检验 同步环闭合差 2020 4 7 独立环闭合差及附和路线坐标闭合差 71 72 GPS测量数据处理的流程 四个阶段 1 数据传输2 格式转换 可选 3 基线处理4 网平差 2020 4 7 2020 4 7 3 7成果验收与上交资料 成果验收 73 上交资料 2020 4 7 74 4框架网 工程网高精度数据处理 利用GAMIT CosaGPS软件进行高精度GPS数据处理算例分析GAMIT基线处理框架网基线处理技巧与方法 质量评估方法参考框架与参考历元的统一及转换方法工程网基线处理技巧与方法 质量评估方法CosaGPS网平差框架网平差处理与精度评估方法工程网平差处理与精度评估方法 2020 4 7 75 AB级基线分量重复性计算 2020 4 7 76 2020 4 7 GPS工程控制网数据处理算例 数据来源某测绘公司采用Leica接收机于2009年3月28 29日在某工程项目中布测了一个高精度GPS控制网 我们取该网中6个GPS控制点 共4个观测时段数的数据为例进行练习预期目标熟练掌握GAMIT CosaGPS高精度GPS数据处理技术与方法GAMIT基线处理CosaGPS平差计算 77 工程案例 实战 演练 工程算例 某县城GPSE级网背景资料介绍某县城GPSE级网数据处理全过程学员实操与答疑 78 2020 4 7 2020 4 7 79 CosaGPS软件 技术讨论 经验交流会 互动环节与答疑 2020 4 7 80 第二部分 TRACK TTC 81 2020 4 7 TRACK 主要内容TRACK模块介绍Track cmd控制策略文件介绍TRACK模块操作过程 算例 track example详细讲解 工程算例 1基准站 1流动站GPS数据 采样间隔1sTRACK模块工程应用学员实操与答疑 82 2020 4 7 TRACK GPS差分相位动态定位程序模块 可以处理GPS动态单历元数据深入文档可以参考网站http geoweb mit edu tah track example目前的版本是TRACKVersion1 21GPS动态定位程序 2009 01 11 使用方法track f a d w s 2020 4 7 83 是一个必须的含有命令列表来控制程序运行的文件是一个可选的含有一个改进设置整周偏差参数和配置的文件是字符串 用于替代命令文件行的 例如 如果给定选项 d222 则bas10 03o将变为bas12220 03o是字符串 用于替代命令文件的任何 为nav file名字提供sp3文件的周数 是10个字符串 在命令文件中用于替代 例如字符串 在命令文件中将替代第一个字符串 字符串 在命令文件中将替代第二个字符串 依次类推 2020 4 7 84 track cmd控制策略文件 2020 4 7 85 2020 4 7 86 2020 4 7 87 2020 4 7 88 2020 4 7 89 TRACK模块使用案例分析算例数据流动站QXJA观测时间 2012年7月21日0 12日基准站 选用IGS站LHAZ精密星历 igr16976 sp3执行命令 track ftrack cmd d203 12 w16976 TRAK 203 12 out 2020 4 7 90 2020 4 7 91 TTC软件 第三部分 GMT 92 2020 4 7 主要内容GMT软件介绍 下载方式 手把手 教会UBUNTU操作系统下GMT软件的安装方法GMT软件的基本命令与使用方法应用GMT软件实现大地测量成果图形表达 93 2020 4 7 GMT 94 2020 4 7 1GMT的安装 通用制图工具 TheGenericMappingTools GMT 是一个被地理学界广泛使用的绘图工具 可以完成海岸线 国界 河流等的绘制GMT绘图工具的获取http gmt soest hawaii edu 2020 4 7 95 GMT安装准备工作在安装GMT之前 需要先安装如下软件NetCDF库 版本3 4甚至更高 可以从www unidata ucar edu网站下载 C编译器 一般Linux系统都已安装 如果没有可以从www gnu org网站下载 此外 还需要安装ps文档阅读器 如ghostview等 2020 4 7 96 GMT安装将整个GMT安装文件夹 gmt 拷贝到 home user 下或任意目录下打开终端 进入到gmt文件目录下 cd home user gmt 然后输入 sudoshinstall gmtGMTparam txt 2020 4 7 97 安装完成后配置环境以管理员身份 sudopasswdroot 打开 etc下的environment文件 在PATH行后面追加 usr local gmt bin PATH再最后行添加 ETCDFHOME usr local netcdf 注意 第一个不用引号 因为已经存有引号 第二个需要加上引号 然后 保存退出 重启系统 在终端测试是否安装成功 输入psxy等进行测试 2020 4 7 98 GMT软件安装演示过程 99 2020 4 7 2Linux相关知识点 重定向 使用bash可以方便地用实现输出输入的重定向cmd filecmd filecmdfilecmd file 100 2020 4 7 管道符 利用Linux所提供的管道符 将两个命令隔开 管道符左边命令的输出就会作为管道符右边命令的输入利用一个管道 rpm qa greplicq利用多个管道 cat etc passwd grep bin bash wc 1 101 2020 4 7 3GMT常用参数 102 2020 4 7 数据域或者地图区域 R选项 R参数用来定义作图的地图或者是感兴趣数据所在区域 通常采用以下三种方式进行定义 Rxmin xmax ymin ymax Rxlleft ylleft xuright yurightr Rgridfile 103 2020 4 7 坐标变换和地图投影 J选项 Ji parameters scale 其中i为一小写字母 用来选择具体的地图投影方式 parameters为0或者是一组以斜画线分割的其他投影参数 scale为地图比列尺 按距离 度或者1 xxxxx的格式给出 JI parameters width I是一个大写字母 用来选择具体的投影方式 parameters为0或者是一组以斜画线分割的其他投影参数 width为地图的宽度 104 2020 4 7 地图框架和坐标轴注记 B选项使用格式为 B p s xinfo yinfo zinfo 标题 W w E e S s N n Z z 可以通过大小写开关来选择需要作图的边界 或轴 xinfo yinfo和zinfo的格式如下 a tick m c ftick m c gtick m c l p axislabel unitlabel 其中a f g分别为坐标轴标尺的大 中和小的间隔 axislabel 为坐标轴的标签以及 unitlabel 为坐标轴的单位 105 2020 4 7 4GMT常用命令 创建图件底图命令 psbasemappsbasemap R95 5 108 20r Ju46 1 10000000 B3g3 UTM utm ps 106 2020 4 7 绘制陆地 水体 海岸线 边界和河流命令 pscoastpscoast R 30 10 60 65 Jm1c B5 Gp100 28 iceland ps 107 2020 4 7 绘制线 多边形和标志的命令 psxy绘制中国边界 CHINA Boundary dat psxyChina Boundary dat Jl103 35 28 40 1 28000000 R64 138 18 54 m China Boundary ps 108 2020 4 7 绘制中国省界 CHINA Prov dat psxyChina Province dat Jl103 35 28 40 1 28000000 R64 138 18 54 m China Province ps 109 2020 4 7 绘制速度矢量 十字和楔形标志命令 psvelopsvelop中 S的含义如下 Sevelscale confidence fontsize 绘制 N E 的速度椭圆 Snbarscale 绘制各向异形柱 Swwedge scale wedge mag 绘制旋转扇形 Sxcross scale 绘制应变十字 110 2020 4 7 这条命令将生成具有边界红色内部绿色的箭头和边界红色椭圆的速度矢量图 数据的置信期间为39 psvelotest ps 111 2020 4 7 Long Lat EvelNvelEsigNsigCorENSITE deg deg mm yr mm yr 0 800460 54x6 8533000 53x30046460 5 5 564640 56x450 6464 0 5 6x40 56 464 0 56x 4 112 2020 4 7 113 2020 4 7 绘制震源机制解命令 psmeca其中的 S参数的具体含义如下 Sascale fontsize offset u 绘制AkiandRichards格式的震源机制 Scscale fontsize offset u 绘制HarvardCMT格式的震源机制 Sm d zscale fontsize offset u 绘制地震矩张量 Spscale fontsize offset u 绘制分离格式的震源机制解 Sx y tscale fontsize offset u 绘制主轴 114 2020 4 7 下面这条命令将绘制一个正断层的CMT震源机制图psmeca R239 240 34 35 2 Jm4 Sc0 4 H1test pslonlatdepthstrdipslipstdipslipmantexpplonplat239 38434 5561218018 88072 905 5000239 82635 0001218018 88072 905 5000END 115 2020 4 7 5工程应用 GMT绘制速度场 分步操作 生成cpt文件makecpt Cglobe T 10000 10000 200 Z V tmp cpt选定范围与投影psbasemap R73 136 3 54 Jq1 50000000 B10g5SWen K P X3 5 Y10 5 velf ps由全球地形文件 global grd 绘制地图grdimageglobal grd R73 136 3 54 Jq O K Q Ctmp cpt velf ps 2020 4 7 116 绘制海岸线 国界 常用河流信息pscoast R73 136 3 54 Jq O K Dh A300 I1 W0 1p N1 1p I1 1p P velf ps由速度文件 vel data txt 绘制速度场psvelovel data txt H1 Gblue Se0 02 0 4 3 Jq R A0 01 0 02 0 02 P V O velf ps 2020 4 7 117 利用GMT软件绘制速度场效果图 2020 4 7 118 利用GMT软件绘制GPS速度场 脚本 bin csh 设定该脚本所调用的shell 该程序调用的是csh gmtsetBASEMAP TYPEPLAIN 设定地图地图样式为PLAIN 另一个选项是FANCY setrange 70 140 10 60 设定地图的坐标范围 setprojection q96 1 32000000 设定地图的投影格式和比例尺大小 119 2020 4 7 setticks 10f5 10f5 设定坐标注记的方式 setpsfile vel ps 设定输出文件的名称 makecpt Cglobe T 10000 10000 200 Z V tmp cpt 设定自定义颜色查找表 psbasemap R range J projection K B ticks WSen X2 Y2 psfile 绘制底图 采用前述的投影格式 范围 120 2020 4 7 grdimageglobal grd Ctmp cpt Jq R O K psfile 绘制二进制格网数据Asia 3m grd 采用系统自带的颜色查找表GMT topo cpt pscoast R Jq O K Di I1 N1 psfile 绘制基本的边界和海岸线 awk print 1 2 3 4 5 6 7 0 vel data txt psvelo R W4 0 255 0 G0 0 255 Se0 03 0 68 12 Jq A0 02 0 24 0 06 O K psfile 绘制GPS速度场 使用awk命令将原始数据重排列为psvelo命令采用的格式 W为画笔属性 G为矢量填色 A为箭头属性 121 2020 4 7 利用GMT软件绘制GPS速度场效果图 122 2020 4 7 第四部分 GLOBK CosaGPS 123 2020 4 7 分步式 批处理 玩转 GLOBK软件 主要内容分步操作GLOBK软件 分步操作的基本命令与使用方法 globk com cmd控制文件介绍glorg com cmd控制文件介绍应用GLOBK软件进行测站坐标时间序列分析应用GLOBK软件进行测站速度场 水平和垂直 分析批处理 sh globk sh glorg 介绍sh globk sh glorg操作GLOBK软件的基本流程 算例 GLOBK10 50软件自带算例详细讲解 124 2020 4 7 GLOBK 125 2020 4 7 1GLOBK软件的安装 安装GLOBK 输入y 126 2020 4 7 2GLOBK软件的模块 格式转换模块 htoglb 将处理GPS VLBI SLR数据的分析软件的解算文件转换成GLOBK软件所需要的二进制文件 h 文件 目前支持的几类文件 GAMIT软件h 文件关于GPS 或其他空间大地测量技术 SINEX格式文件FONDA软件h 文件 127 2020 4 7 JPL机构提供的Stacov文件包含站坐标和速度场的SLR GSFC文件包含站坐标和速度场的VLBI GSFC文件运算模块 glred globk glorg glred模块通过调用globk模块分析单天解 对于分析基准站网 CORS 非常适用 其生成的解文件可以用来形成时间序列globk模块是GLOBK软件的主模块 实现该软件的功能 128 2020 4 7 glorg模块可以为平差结果定义参考框架 具体通过固定 或约束 站坐标和速度由坐标转换来实现glorg模块可以单独运行 也可以被globk glred模块调用cmd file是glorg和globk的控制文件 里面包含解算的策略等GMT图形应用模块sh plotcrd 129 2020 4 7 sh globk scattermultibasesh plotvel主要功能是利用GMT软件绘制时间序列 速度场等图形 可用于分析数据质量和测站的地壳运动等情况其他辅助模块一类为GLOBK模块服务的 glist glsave一类是用于框架之间和板块运动分析的 corcom cvframe velrot 130 2020 4 7 3常用文件与格式说明 输入文件GPS SLR VLBI和SINEX文件随机特征可由apr XXX和mar XXX表述数据文件 二进制h 文件控制文件 指令cmd文件GLOBK结果文件中的NEU坐标GLOBK的输出文件一般为 prt和 org在给出ITRF坐标同时还给出了新的NEU坐标 它与站坐标定义的NEU不同 这种坐标类似于平面坐标 属于圆锥投影 131 2020 4 7 运行GLOBK软件时的注意事项GLOBK是基于线性模型的 测站坐标改正值大于10米或轨道参数的改正值大于100米时需要进行数据的再处理GLOBK不能解决因周跳未得到完全探测 数据质量差或大气层延迟模型误差所带来的问题GLOBK不能进行整周模糊度的解算 前期数据处理解算必须完成整周模糊度的解算 132 2020 4 7 4GLOBK软件的使用 GLOBK数据处理流程将ASCII格式的h file转换成可被globk读取的二进制h file 然后运行glred globk以获得测站坐标的时间序列通过时间序列分析 确定具有异常域的特定站或特定历元 在earthquakefile中 运用rename命令删除具有异常域的特定站的特定历元或直接删除对应的h file 133 运行globk将单时段解的h file合并成一个h file 其代表在所选择的时间跨度里测站的平均坐标使用合并后的h file 再次运行glred globk获得时间序列 而运行globk glorg则可获得测站速度利用GLOBK进行网平差 分布操作 建立工作目录 比如 globk test 在该工作目录下还应包括以下三个子目录 134 2020 4 7 glbf 用于存储经htoglb命令运行所生成的二进制H 文件soln 存储globk的控制文件 globk com cmd glorg com cmd ASCII文本H 文件 GAMIT基线结果文件以及经平差处理后的平差结果文件 prt 并在此目录下运行GLOBK软件tables pmu unso 测站先验坐标文件 apr 经htoglb生成的 svs参数文件 2020 4 7 135 生成二进制H 文件将svnav dat文件拷贝到soln文件夹下将GAMIT生成的文本H 文件拷贝至soln文件夹下在soln工作目录下 运行命令 生成 gdl文件列表 通常以gdl为后缀 globaldirect