GAMITGLOBK软件技术应用高级培训班

GAMIT/GLOBK软件技术应用主讲人：周命端博士E-mail：zmd_zry@163.comGPS测量技术工程应用高级培训班(第二期)1中国·北京2014年8月16-19日2023/5/16GAMIT/GLOBK/GMT/TRACKCosaGPSTRACKGMTGLOBKGAMIT学习22023/5/16内容提要GAMIT/GLOBK/GMT/TRACK软件安装与操作专题实操练习：应用GAMIT+CosaGPS软件进行高精度GPS工程控制网数据处理与精度评估综合实操练习：应用GAMIT/GLOBK/GMT/TRACK软件进行CORS站网数据分析高精度GPS数据处理“技术沙龙”2023/5/16第一部分:GAMIT+CosaGPS42023/5/161操作系统主要内容虚拟机(VMwareWorkstation)的介绍、下载、安装与基本使用UBUNTU操作系统的下载、安装与基本使用(cd、ls、ln、gedit、chmod命令与操作)基于VMwareWorkstation的UBUNTU操作系统的安装52023/5/161.1虚拟机(VMwareWorkstation)虚拟工作站(VMwareWorkstation)VMwareWorkstation是一款功能强大的桌面虚拟计算机软件下载地址：http:///cn最新版本为10.02(2014年4月17日发行)2023/5/1661.2UBUNTU操作系统乌班图(Ubuntu)Ubuntu是一个以桌面应用为主的Linux操作系统下载地址：http://最新版本是14.04LTS(2014年4月18日发行)ls(列出目录内容)cd(切换当前目录)[cd../][cd./]gedit(用文本编辑器打开文件)ln(创建文件链接)[ln-s]chmod(添加执行属性)[chmod+x]2023/5/1671.3操作系统安装过程基于VMwareWorkstation的UBUNTU操作系统的安装演示82023/5/162

GAMIT/GLOBK软件介绍与安装主要内容GAMIT/GLOBK软件介绍、下载方式UBUNTU操作系统下安装GAMIT/GLOBK软件的系统配置(csh、gfortran、libX11-dev)“手把手”教会GAMIT/GLOBK软件的安装GAMIT/GLOBK软件的路径配置(.bashrc)与测试(doy)92023/5/162.1GAMIT/GLOBK软件介绍与下载GAMIT/GLOBK安装包GAMIT/GLOBK是一套高精度数据处理软件。由美国麻省理工学院(MIT)和斯克里普斯海洋研究所(SIO)联合开发。下载地址：/

最新版本是10.50(2014年4月18日Release)2023/5/16102.2操作系统配置GAMIT/GLOBK安装准备工作更换源：系统→系统管理→新立得软件包管理器设置→软件库打开终端，用户名目录下执行如下命令sudoapt-getinstallupdatesudoapt-getinstallcshsudoapt-getinstallgfortransudoapt-getinstalllibX11-dev2023/5/1611连网2.3安装GAMIT/GLOBK安装安装目录下，运行：修改./libraries/Makefile.config文件./install_software#Generic(willworkonanysystemiflinksinplace)X11LIBPATH/usr/libX11INCPATH/usr/include/X11MAXSIT60MAXSAT32MAXATM13MAXEPC2880MAXSIT60MAXSAT32MAXATM25MAXEPC57602023/5/16122.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.bashrc2023/5/16133分步式“玩转”GAMIT软件主要内容分步操作GAMIT软件(分步操作的基本命令使用)标准RINEX文件准备(【算例】：某工程控制网GPS数据4个观测时段、国标D级)关键文件准备(sestbl.、、sittbl.、lfile.、lfile.apr)与详细介绍高精度GPS基线处理(“手把手”教分步操作命令)基线质量评估方法(“四看”)142023/5/163.1分步操作GAMIT软件外业数据准备文件准备数据处理(操作GAMIT)质量评估152023/5/16数据准备观测文件标准化(原始文件→RINEXo-files)采用随机格式转换软件或TEQC软件可将原始文件转换为标准的rinexo-files格式对所有的o-files文件的命名进行标准化，如sitedoyn.yyo。注意文件名字字母要小写。2023/5/1616天线高改化至DHARP(天线座底部)外业天线高量测的方式根据不同的接收机天线型号有不同的量测方法GAMIT软件里设置了不同的量高模式。建议采用天线座底部(DHARP)天线高改化至L1PHC(天线相位中心)2023/5/16172023/5/1618GPS接收机天线高改化原理、量高模式与改化方法GPS接收机天线量高模式

2023/5/1619GPS接收机天线高改化原理H为标石中心天线参考点的斜高R为天线盘半径h1为标石中心天线参考点的垂高h2为天线参考点天线座底部的垂高；h3为天线参考点天线相位中心的垂高标石中心天线相位中心的垂高：

GPS测量技术工程应用培训班2023/5/1620已知外业测量的斜高为1.564m求：

从网上查询或下载数据文件2023/5/1621IGS跟踪站o-files数据、导航电文brdc文件及精密星历SP3文件：

ftp:///pub/；ftp:///pub/gps/

ftp://igs.ensg.ign.fr/pub/。更新tables文件。可以从服务器ftp:///pub/gamit

下载更新的tables表文件其中luntab.、soltab.、nutabl.等文件按年进行更新pole.、ut1.、pmu.等文件每周都有更新222023/5/16文件结构整理：对每个工程目录做处理前的准备工作，将数据文件分类存放，便于维护和更新步骤如下：①在当前工程目录下建立的文件夹有:igs、rinex、brdc以及时段文件夹igs文件夹存放精密星历文件rinex文件夹存放o-files文件brdc文件夹存放brdc文件时段文件夹，以时段号命名，所有解算都是基于该文件夹进行的文件准备232023/5/16将更新好的完好tables表文件复制到当前工程目录下可以手动复制或命令cp–r../tables./。注意：这里尽量用复制！当然，也可以用连接ln-s../../gamit/source/tables./至此，当前工程目录文件夹结构应该为：242023/5/16修改tables文件夹下的配置文件文件：记录各测站的站点、站点名称、开始观测时间、结束观测时间、天线高、天线量测方法、接收机类型、软硬件版本号、天线类型，天线相位中心模型机构代码等信息。注意：AntHt、HtCod(标志天线量测方法DHARP、SLBGP等)和AntennaType(用来做天线高改正)，ReceiverType(用来做周跳恢复)，SessionStart、SessionStop(用来控制测段跨度，处理不好后面会提示错误)252023/5/16严格按照已经给出的跟踪站格式编辑具体操作方法有2种：手动操作：按照相应文件格式对齐操作，按照标准格式.standard编辑修改。sh_upd_stnfo自动提取观测文件的有关信息，生成，并认真核对信息是否正确。事先提供空模板后，再执行命令：

sh_upd_stnfo–ref–files*.yyo262023/5/16l-file文件：记录各测站的近似坐标Gamit要求是球面坐标，可以先在测站的O文件里手动提取近似坐标XYZ，生成l-file文件，按照格式对齐，然后保存文件，将名字重新命名为lfile.apr

，然后用gapr_to_l将lfile.apr转换到lfile.。可以用sh_rx2apr命令自动生成，利用svpos和svdiff模块单点定位生成站点坐标。也可以在测站的O文件里自动提取近似坐标：

grepPOSITION*.yyo>lfilerx2aprlfile<year><doy>gapr_to_llfile.aprlfile.“”<year><doy>272023/5/16sestbl.文件gambit基线处理时的模型选择配置文件。需要仔细参照技术文档和程序代码，一般情况下采用默认值，通常需要修改的地方有：ChoiceofObservable(观测值的选择类型)ChoiceofExperiment(基线处理类型)sittbl.文件(sittbl.short)各测站的精度控制指标，一般情况下采用默认值。高精度的已知坐标强约束，待求点坐标松弛约束。282023/5/16将tables表文件中需要的文件（包括luntab、soltab、nutabl、pole、ut1、pmu文件，还有一些误差改正模型文件等）链接到各时段文件夹。该步骤的作用是：在数据处理各个时段的工作目录下，建立该工作目录与/tables目录中部分文件的链接，为后面在工作目录下逐行输入命令进行数据解算做准备。在各时段文件夹下分别建立与tables文件夹里相应文件对应的链接，输入命令：links.dayyrdoytest（可以参考links.day命令的用法）292023/5/16然后需要手动建立以下链接，命令如下

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）302023/5/16一定要检查一下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./312023/5/16数据文件准备分步操作：在终端，输入命令：makexp，回车，然后按照提示输入相应的文件数据。提示在要求输入l-file文件时候，建议输入lfile.apr文件。批操作

sh_makexp–expttest–orbt

igsf–yr2009–doy088–sess99–srin–navbrdc0880.09n–apr

lfile.–sinfo1500005760需要准备的文件：O文件、N文件；文件；rcvant.dat文件，接收机天线类型。输出文件：数据处理(分步处理)322023/5/16精密星历处理sh_sp3fit–figs15250.sp3–oigsf–d2009088–rBERNE–t–u需要准备的文件

sp3文件；ut1.，时长变化表；

pole.，极移表；

nutabl.，章动表；

soltbl.，太阳星历表；luntab.，月亮星历表；leap.sec，跳秒；svnav.dat，卫星天线类型332023/5/16输出文件

gigsf9.088；说明：g-files文件所有卫星初始轨道坐标文件tigsf9.088；文件检查sh_check_sess–sess088–typegfile–filegigsf9.088需要准备的文件：gtest9.088；输出文件：广播星历处理makejbrdc0880.09njbrdc9.088需要准备的文件：brdc0880.09n；输出的文件：jbrdc9.088342023/5/16文件检查sh_check_sess–sess088–typejfile–filejbrdc9.088需要准备的文件jbrdc9.088输出文件形成中间文件makex

test.makex.batch需要准备的文件

test.makex.batch，以及该文件所列出的所有文件输出文件dtest9.088352023/5/16

生成bat文件输入命令

fixdrvdtest9.088需要准备的文件：dtest9.088，以及该文件所列出的所有文件输出的文件：btest9.088.bat，以及btest9.num基线解算输入命令

cshbtest9.bat需要准备文件：btest9.bat，以及btest9.num输出文件：q-、o-、h-等所有解算的结果文件。362023/5/16GAMIT基线处理结果质量评价指标检查所有测站点是否都参与计算(LC_HELP)标准均方根nrms解算结果为0.20左右若是nrms大于0.5，表明处理结果是有问题的（例如：周跳没有修复；测站的起算点坐标有问题等）。若是nrms小于0.5，则认为成功（合理）解算，否则需要检查原因，重新处理；检查各基线文件处理精度是否符合gamit的实际精度质量评估372023/5/164定制式“玩转”GAMIT软件主要内容bash编程语言介绍基于bash编程语言自主编写“定制化”批处理脚本(*.sh)【算例】：某工程控制网GPS数据4个观测时段、国标D级(*.sh文件)382023/5/162023/5/16使用GAMIT软件进行基线解算有一定程度的了解之后，可以不用每次只执行一条命令，可利用UNIX/LINX系统自带的shell编程语言，定制适合自己的批处理脚本Shell命令。值得注意的是：Shell脚本命令只对数据处理进行批处理，数据准备及文件准备建议手动操作。以某城市CORS2012年doy=201第0时段基线解算为例，定制的批处理脚本文件2010_cs.sh内容如下所示。392023/5/16#!/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'402023/5/16##step2：linkthetablesintocurrentdirectorylinks.day$yr_4$doy$expt##step3：removebrokenfilesandrestart-linkfilesintocurrentdirectoryrm

luntab.rm

soltab.rm

nutabl.rmpole.rmut1.rm

pmu.412023/5/16ln-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.422023/5/16##step4：removerinexfilesandrestrart-linkrinexfilesintocurrentdirectoryrm*.12orm*.sp3rm*.12nln../rinex/${doy}${session}/*${doy}${session}.12o./ln../brdc/$brdcfile./ln../igs/$sp3file./432023/5/16##step5：preparationfortheexperimentrm*.fatalrm

sh_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}.bat442023/5/16##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}45GAMIT软件基线解算流程图462023/5/16技术讨论：经验交流会互动环节与答疑2023/5/16475批处理“玩转”GAMIT软件

主要内容【算例】：中国大陆8个IGS站(CHAN、BJFS、URUM、GUAO、LHAZ、WUHN、SHAO、KUNM)2013年doy=001~005高精度基线处理(sh_gamit)批处理(sh_gamit)介绍与工程应用使用前提sh_gamit操作GAMIT软件的基本流程482023/5/16【算例】中国大陆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.pssh_gamitStep1：在./2013目录下分别建立igs、brdc、rinex文件夹，并存放对应数据,并在./2013目录下，终端执行：sh_setup-yr2013编辑./2013/tables/process.defaults编辑./2013/tables/sites.defaults准备准备lfile.apr和lfile.编辑sestbl.和sittbl.在./2013目录下，终端执行：-d:所处理数据所在的年、年积日-s：2013001003-copt：待压缩的文件类型-dopt：待删除的文件类型-yrext：指定待处理数据年份及年积日-pres：NO、YES和ELEVsh_gamit-exptcors-d2013001002003-presELEV-orbitIGSF-coptxkp-doptscao专题实操练习

应用GAMIT/CosaGPS软件进行高精度GPS工程控制网数据处理与精度评估

532023/5/16主要内容基于IGS站的框架网构建与优化设计IGS站点数据(GPS观测资料、三维框架坐标)获取途径大地测量参考框架(ITRF系列框架与参考历元的基准统一及其转换)全球定位系统（GPS）测量规范(GB-T18314-2009)解读【算例】：某A级GPS变形监测网高精度数据处理中框架网构建542023/5/16利用“定制式”脚本操作GAMIT软件进行高精度GPS基线处理（可获得高精度GPS基线向量、天顶对流层延迟参数）高精度GPS基线向量质量评估【算例】：AB级GPS基线分量、长度重复性计算应用GLOBK/CosaGPS软件进行GPS网平差与精度评估数据计算报告与技术总结报告撰写552023/5/161基于IGS站的框架网优化设计与构建IGS站网2023/5/1656SOPAC网站：/

2023/5/1657IERS网站：

2023/5/16582ITRF系列框架与参考历元的基准统一2023/5/16592023/5/16603全球定位系统(GPS)测量规范解读术语和定义基本规定级别划分和测量精度布设的原则选点埋石仪器观测外业成果记录数据处理成果验收与上交资料612023/5/162023/5/16623.1基本规定GPS测量采用2000国家大地坐标系统采用GPS时间系统手簿记录宜采用世界协调时(UTC)接收机应经检定合格，在有效期内使用采用中误差为精度指标，以2倍为极限误差当需要提供80或54或其他坐标成果时，应按坐标转换方法求得。2023/5/163.2级别划分和测量精度级别划分GPS测量按照精度和用途分为A~E级5等级测量精度A级GPS网由CORS站构成B~E级的精度应不低于表2要求632023/5/1664（1）按工作原理分：码相关型接收机平方型接收机混合型接收机（2）按信号通道的类型分：多通道接收机序贯通道接收机多路复用通道接收机（3）根据接收信号的频率分：单频接收机（L1）双频接收机（L1+L2）（4）根据接收机的用途分：导航型测量型授时型3.3仪器GPS接收机的分类2023/5/16接收机选用A级网测量按有关规定执行B~E级GPS网按表4规定执行GPS接收机的检验新购置GPS接收机的检验一般性检视、通电检验、试测检验旧GPS接收机的检验一般性检视、通电检验、试测检验的准备部分652023/5/16663.4观测A级网B~E级GPS网观测的基本技术规定应符合表5的要求2023/5/163.5外业成果记录记录类型观测数据、测量手簿、其他记录记录内容原始观测数据、对应观测值的GPS时间、测站与接收机初始信息、天线高记录要求及时、认真、准确672023/5/163.6数据处理原始格式—>RINEX格式方法GPS数据文件的组织外业数据质量检验数据删除率：不宜大于10%复测基线的长度差682023/5/1669各级GPS测量的精度指标相邻点间基线长度标准差例如：a=5mmb=1ppmD=5km2023/5/1670外业数据质量检验同步环闭合差2023/5/16独立环闭合差及附和路线坐标闭合差7172GPS测量数据处理的流程四个阶段：1.数据传输2.格式转换（可选）3.基线处理4.网平差2023/5/162023/5/163.7成果验收与上交资料成果验收73上交资料2023/5/16744框架网/工程网高精度数据处理利用GAMIT/CosaGPS软件进行高精度GPS数据处理算例分析GAMIT基线处理框架网基线处理技巧与方法、质量评估方法参考框架与参考历元的统一及转换方法工程网基线处理技巧与方法、质量评估方法CosaGPS网平差框架网平差处理与精度评估方法工程网平差处理与精度评估方法2023/5/1675AB级基线分量重复性计算2023/5/16762023/5/16GPS工程控制网数据处理算例数据来源某测绘公司采用Leica接收机于2009年3月28~29日在某工程项目中布测了一个高精度GPS控制网，我们取该网中6个GPS控制点，共4个观测时段数的数据为例进行练习预期目标熟练掌握GAMIT+CosaGPS高精度GPS数据处理技术与方法GAMIT基线处理CosaGPS平差计算77工程案例“实战”演练【工程算例】：某县城GPSE级网背景资料介绍某县城GPSE级网数据处理全过程学员实操与答疑782023/5/162023/5/1679CosaGPS软件技术讨论：经验交流会互动环节与答疑2023/5/1680第二部分:TRACK/TTC812023/5/16TRACK主要内容TRACK模块介绍Track.cmd控制策略文件介绍TRACK模块操作过程(【算例】：track_example详细讲解)【工程算例】：1基准站+1流动站GPS数据,采样间隔1sTRACK模块工程应用学员实操与答疑822023/5/16TRACK：GPS差分相位动态定位程序模块，可以处理GPS动态单历元数据深入文档可以参考网站/~tah/track_example目前的版本是TRACKVersion1.21GPS动态定位程序(2009/01/11)

使用方法track-f<commandfile>-a<ambiguityfile>-d<day>-w<week>-s<S01><S02>..<S10>2023/5/1683<commandfile>是一个必须的含有命令列表来控制程序运行的文件<ambiguityfile>是一个可选的含有一个改进设置整周偏差参数和配置的文件<day>是字符串，用于替代命令文件行的<day>，例如：如果给定选项-d222，则bas1<day>0.03o将变为bas12220.03o<week>是字符串，用于替代命令文件的任何<week>，为nav_file名字提供sp3文件的周数<S01>,<S02>..<S10>是10个字符串，在命令文件中用于替代，例如<S01>字符串，在命令文件中将替代第一个字符串，<S02>字符串，在命令文件中将替代第二个字符串，依次类推2023/5/1684track.cmd控制策略文件2023/5/16852023/5/16862023/5/16872023/5/16882023/5/1689TRACK模块使用案例分析算例数据流动站QXJA观测时间：2012年7月21日0~12日基准站：选用IGS站LHAZ精密星历：igr16976.sp3执行命令：track-ftrack.cmd-d203_12-w16976>!TRAK_203_12.out2023/5/16902023/5/1691TTC软件第三部分:GMT922023/5/16主要内容GMT软件介绍、下载方式“手把手”教会UBUNTU操作系统下GMT软件的安装方法GMT软件的基本命令与使用方法应用GMT软件实现大地测量成果图形表达932023/5/16GMT工程应用常用命令常用参数安装GMT942023/5/161GMT的安装通用制图工具(TheGenericMappingTools,GMT)是一个被地理学界广泛使用的绘图工具，可以完成海岸线、国界、河流等的绘制GMT绘图工具的获取2023/5/1695GMT安装准备工作在安装GMT之前，需要先安装如下软件NetCDF库，版本3.4甚至更高(可以从网站下载)C编译器(一般Linux系统都已安装，如果没有可以从网站下载)此外，还需要安装ps文档阅读器，如ghostview等2023/5/1696GMT安装将整个GMT安装文件夹(gmt)拷贝到：/home/user/下或任意目录下打开终端，进入到gmt文件目录下(cd/home/user/gmt)，然后输入：sudoshinstall_gmtGMTparam.txt2023/5/1697安装完成后配置环境以管理员身份(sudopasswdroot)打开/etc下的environment文件，在PATH行后面追加：/usr/local/gmt/bin:$PATH再最后行添加：ETCDFHOME=“/usr/local/netcdf”注意：第一个不用引号，因为已经存有引号，第二个需要加上引号。然后，保存退出，重启系统。在终端测试是否安装成功：输入psxy等进行测试2023/5/1698GMT软件安装演示过程992023/5/162Linux相关知识点重定向：使用bash可以方便地用<和>实现输出输入的重定向cmd>filecmd>>filecmd<filecmd<<textcmd<<<wordcmd<>filecmd>|file1002023/5/16管道符(|)：利用Linux所提供的管道符“|”将两个命令隔开，管道符左边命令的输出就会作为管道符右边命令的输入利用一个管道#rpm–qa|greplicq利用多个管道#cat/etc/passwd|grep/bin/bash|wc-11012023/5/163GMT常用参数参数含义-B定义底图和坐标轴的刻度、注记和标签-H输入/输出表的头记录行数-J选择所用的地图投影或坐标转换-K允许在当前图件中追加图画-O允许将当前图件追加已存在的图件中-P将图件页面设置成纵向模式(默认是横向)-R定义地图/图件的区域-U添加时间戳，默认位于页面的左下角-V备注默认，报告运行的进度-X设置图件在页面中的起始x坐标-Y设置图件在页面中的起始y坐标-b设定输入/输出文件采用二进制格式-c设定图件的备份数-f设定每列数据的具体格式-：假定地理数据的格式为(纬度、经度)，而不是(经度、纬度)1022023/5/16数据域或者地图区域：-R选项

-R参数用来定义作图的地图或者是感兴趣数据所在区域，通常采用以下三种方式进行定义：

-Rxmin/xmax/ymin/ymax

-Rxlleft/ylleft/xuright/yurightr

-Rgridfile1032023/5/16坐标变换和地图投影：-J选项-Ji[parameters/]scale，其中i为一小写字母，用来选择具体的地图投影方式，parameters为0或者是一组以斜画线分割的其他投影参数，scale为地图比列尺，按距离/度或者1：xxxxx的格式给出-JI[parameters/]width，I是一个大写字母，用来选择具体的投影方式，parameters为0或者是一组以斜画线分割的其他投影参数，width为地图的宽度1042023/5/16地图框架和坐标轴注记：-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”为坐标轴的单位1052023/5/164GMT常用命令创建图件底图命令：psbasemappsbasemap-R95/5/108/20r-Ju46/1:10000000-B3g3:.UTM:>utm.ps1062023/5/16绘制陆地、水体、海岸线、边界和河流命令：pscoastpscoast-R-30/-10/60/65-Jm1c-B5-Gp100/28>iceland.ps1072023/5/16绘制线、多边形和标志的命令：psxy绘制中国边界(CHINA_Boundary.dat)psxyChina_Boundary.dat-Jl103/35/28/40/1:28000000-R64/138/18/54-m>China_Boundary.ps1082023/5/16绘制中国省界(CHINA_Prov.dat)psxyChina_Province.dat-Jl103/35/28/40/1:28000000-R64/138/18/54-m>China_Province.ps1092023/5/16绘制速度矢量、十字和楔形标志命令：psvelopsvelop中-S的含义如下：

-Sevelscale/confidence/fontsize：绘制(N，E)的速度椭圆-Snbarscale：绘制各向异形柱-Swwedge_scale/wedge_mag：绘制旋转扇形-Sxcross_scale：绘制应变十字1102023/5/16这条命令将生成具有边界红色内部绿色的箭头和边界红色椭圆的速度矢量图，数据的置信期间为39%psvelo<<END-H2-R-10/10/-10/10-W0.25p,red-Ggreen-L-Se0.2/0.39/18-B:"time(s)":2g2/:"year(n)":2g2:."figure":WeSn-Jx0.4/0.4-A0.1/0.3/0.3-P-V>test.ps1112023/5/16Long.Lat.Evel

Nvel

Esig

Nsig

CorENSITE(deg)(deg)(mm/yr)(mm/yr)0-800460.54x6-8533000.53x30046460.5-5-564640.56x450-6464-0.5-6x40-56-464-0.56x-41122023/5/161132023/5/16绘制震源机制解命令：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]]]：绘制主轴1142023/5/16下面这条命令将绘制一个正断层的CMT震源机制图psmeca-R239/240/34/35.2-Jm4-Sc0.4-H1<<END>test.pslonlatdepthstrdipslipstdipslipmantexpplonplat239.38434.5561218018-88072-905.5000239.82635.0001218018-88072-905.5000END1152023/5/165工程应用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)绘制地图grdimage

global.grd-R73/136/3/54-Jq-O-K-Q-Ctmp.cpt>>velf.ps

2023/5/16116绘制海岸线、国界、常用河流信息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

2023/5/16117利用GMT软件绘制速度场效果图2023/5/16118利用GMT软件绘制GPS速度场(脚本)#!/bin/csh#设定该脚本所调用的shell，该程序调用的是csh。gmtsetBASEMAP_TYPEPLAIN#设定地图地图样式为PLAIN，另一个选项是FANCY。setrange=70/140/10/60#设定地图的坐标范围。setprojection=q96/1:32000000#设定地图的投影格式和比例尺大小。1192023/5/16setticks=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}#绘制底图，采用前述的投影格式、范围。1202023/5/16grdimageglobal.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为箭头属性。1212023/5/16利用GMT软件绘制GPS速度场效果图1222023/5/16第四部分:GLOBK/CosaGPS1232023/5/16分步式/批处理“玩转”GLOBK软件主要内容分步操作GLOBK软件(分步操作的基本命令与使用方法)globk_com.cmd控制文件介绍glorg_com.cmd控制文件介绍应用GLOBK软件进行测站坐标时间序列分析应用GLOBK软件进行测站速度场(水平和垂直)分析批处理(sh_globk、sh_glorg)介绍sh_globk/sh_glorg操作GLOBK软件的基本流程【算例】：GLOBK10.50软件自带算例详细讲解1242023/5/16GLOBK安装工程应用软件使用常用文件与格式说明软件模块GLOBK1252023/5/161GLOBK软件的安装安装GLOBK,输入y1262023/5/162GLOBK软件的模块格式转换模块(htoglb)将处理GPS,VLBI,SLR数据的分析软件的解算文件转换成GLOBK软件所需要的二进制文件(h-文件)。目前支持的几类文件：GAMIT软件h-文件关于GPS(或其他空间大地测量技术)SINEX格式文件FONDA软件h-文件1272023/5/16JPL机构提供的Stacov文件包含站坐标和速度场的SLR/GSFC文件包含站坐标和速度场的VLBI/GSFC文件运算模块(glred,globk,glorg)glred模块通过调用globk模块分析单天解，对于分析基准站网(CORS)非常适用，其生成的解文件可以用来形成时间序列globk模块是GLOBK软件的主模块，实现该软件的功能1282023/5/16glorg模块可以为平差结果定义参考框架，具体通过固定(或约束)站坐标和速度由坐标转换来实现glorg模块可以单独运行，也可以被globk、glred模块调用cmd_file是glorg和globk的控制文件，里面包含解算的策略等GMT图形应用模块sh_plotcrd1292023/5/16sh_globk_scattermultibasesh_plotvel主要功能是利用GMT软件绘制时间序列、速度场等图形，可用于分析数据质量和测站的地壳运动等情况其他辅助模块一类为GLOBK模块服务的：glist、glsave一类是用于框架之间和板块运动分析的：corcom，cvframe，velrot1302023/5/163常用文件与格式说明输入文件GPS,SLR,VLBI和SINEX文件随机特征可由apr_XXX和mar_XXX表述数据文件：二进制h-文件控制文件：指令cmd文件GLOBK结果文件中的NEU坐标GLOBK的输出文件一般为*.prt和*.org在给出ITRF坐标同时还给出了新的NEU坐标，它与站坐标定义的NEU不同，这种坐标类似于平面坐标，属于圆锥投影1312023/5/16运行GLOBK软件时的注意事项GLOBK是基于线性模型的：测站坐标改正值大于10米或轨道参数的改正值大于100米时需要进行数据的再处理GLOBK不能解决因周跳未得到完全探测、数据质量差或大气层延迟模型误差所带来的问题GLOBK不能进行整周模糊度的解算：前期数据处理解算必须完成整周模糊度的解算1322023/5/164GLOBK软件的使用GLOBK数据处理流程将ASCII格式的h-file转换成可被globk读取的二进制h-file，然后运行glred/globk以获得测站坐标的时间序列通过时间序列分析，确定具有异常域的特定站或特定历元，在earthquakefile中，运用rename命令删除具有异常域的特定站的特定历元或直接删除对应的h-file133运行globk将单时段解的h-file合并成一个h-file，其代表在所选择的时间跨度里测站的平均坐标使用合并后的h-file，再次运行glred/globk获得时间序列，而运行globk/glorg则可获得测站速度利用GLOBK进行网平差(分布操作)建立工作目录，比如：globk_test，在该工作目录下还应包括以下三个子目录：1342023/5/16glbf：用于存储经htoglb命令运行所生成的二进制H-文件soln：存储globk的控制文件(globk_com.cmd、glorg_com.cmd)、ASCII文本H-文件、GAMIT基线结果文件以及经平差处理后的平差结果文件(*.prt)，并在此目录下运行GLOBK软件tables：pmu.unso、测站先验坐标文件(*.apr)、经htoglb生成的*.svs参数文件2023/5/16135生成二进制H-文件将svna