GAMITGLOBK軟件技術(shù)應(yīng)用高級(jí)培訓(xùn)班

GAMIT/GLOBK軟件技術(shù)應(yīng)用主講人：周命端博士E-mail：zmd_zry@163.comGPS測(cè)量技術(shù)工程應(yīng)用高級(jí)培訓(xùn)班(第二期)1中國(guó)·北京2014年8月16-19日2023/5/16GAMIT/GLOBK/GMT/TRACKCosaGPSTRACKGMTGLOBKGAMIT學(xué)習(xí)22023/5/16內(nèi)容提要GAMIT/GLOBK/GMT/TRACK軟件安裝與操作專(zhuān)題實(shí)操練習(xí)：應(yīng)用GAMIT+CosaGPS軟件進(jìn)行高精度GPS工程控制網(wǎng)數(shù)據(jù)處理與精度評(píng)估綜合實(shí)操練習(xí)：應(yīng)用GAMIT/GLOBK/GMT/TRACK軟件進(jìn)行CORS站網(wǎng)數(shù)據(jù)分析高精度GPS數(shù)據(jù)處理“技術(shù)沙龍”2023/5/16第一部分:GAMIT+CosaGPS42023/5/161操作系統(tǒng)主要內(nèi)容虛擬機(jī)(VMwareWorkstation)的介紹、下載、安裝與基本使用UBUNTU操作系統(tǒng)的下載、安裝與基本使用(cd、ls、ln、gedit、chmod命令與操作)基于VMwareWorkstation的UBUNTU操作系統(tǒng)的安裝52023/5/161.1虛擬機(jī)(VMwareWorkstation)虛擬工作站(VMwareWorkstation)VMwareWorkstation是一款功能強(qiáng)大的桌面虛擬計(jì)算機(jī)軟件下載地址：http:///cn最新版本為10.02(2014年4月17日發(fā)行)2023/5/1661.2UBUNTU操作系統(tǒng)烏班圖(Ubuntu)Ubuntu是一個(gè)以桌面應(yīng)用為主的Linux操作系統(tǒng)下載地址：http://最新版本是14.04LTS(2014年4月18日發(fā)行)ls(列出目錄內(nèi)容)cd(切換當(dāng)前目錄)[cd../][cd./]gedit(用文本編輯器打開(kāi)文件)ln(創(chuàng)建文件鏈接)[ln-s]chmod(添加執(zhí)行屬性)[chmod+x]2023/5/1671.3操作系統(tǒng)安裝過(guò)程基于VMwareWorkstation的UBUNTU操作系統(tǒng)的安裝演示82023/5/162

GAMIT/GLOBK軟件介紹與安裝主要內(nèi)容GAMIT/GLOBK軟件介紹、下載方式UBUNTU操作系統(tǒng)下安裝GAMIT/GLOBK軟件的系統(tǒng)配置(csh、gfortran、libX11-dev)“手把手”教會(huì)GAMIT/GLOBK軟件的安裝GAMIT/GLOBK軟件的路徑配置(.bashrc)與測(cè)試(doy)92023/5/162.1GAMIT/GLOBK軟件介紹與下載GAMIT/GLOBK安裝包GAMIT/GLOBK是一套高精度數(shù)據(jù)處理軟件。由美國(guó)麻省理工學(xué)院(MIT)和斯克里普斯海洋研究所(SIO)聯(lián)合開(kāi)發(fā)。下載地址：/

最新版本是10.50(2014年4月18日Release)2023/5/16102.2操作系統(tǒng)配置GAMIT/GLOBK安裝準(zhǔn)備工作更換源：系統(tǒng)→系統(tǒng)管理→新立得軟件包管理器設(shè)置→軟件庫(kù)打開(kāi)終端，用戶(hù)名目錄下執(zhí)行如下命令sudoapt-getinstallupdatesudoapt-getinstallcshsudoapt-getinstallgfortransudoapt-getinstalllibX11-dev2023/5/1611連網(wǎng)2.3安裝GAMIT/GLOBK安裝安裝目錄下，運(yùn)行：修改./libraries/Makefile.config文件./install_software#Generic(willworkonanysystemiflinksinplace)X11LIBPATH/usr/libX11INCPATH/usr/include/X11MAXSIT60MAXSAT32MAXATM13MAXEPC2880MAXSIT60MAXSAT32MAXATM25MAXEPC57602023/5/16122.4路徑配置GAMIT/GLOBK安裝的路徑配置修改.bashrc文件用戶(hù)名目錄下：測(cè)試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分步式“玩轉(zhuǎn)”GAMIT軟件主要內(nèi)容分步操作GAMIT軟件(分步操作的基本命令使用)標(biāo)準(zhǔn)RINEX文件準(zhǔn)備(【算例】：某工程控制網(wǎng)GPS數(shù)據(jù)4個(gè)觀測(cè)時(shí)段、國(guó)標(biāo)D級(jí))關(guān)鍵文件準(zhǔn)備(sestbl.、、sittbl.、lfile.、lfile.apr)與詳細(xì)介紹高精度GPS基線處理(“手把手”教分步操作命令)基線質(zhì)量評(píng)估方法(“四看”)142023/5/163.1分步操作GAMIT軟件外業(yè)數(shù)據(jù)準(zhǔn)備文件準(zhǔn)備數(shù)據(jù)處理(操作GAMIT)質(zhì)量評(píng)估152023/5/16數(shù)據(jù)準(zhǔn)備觀測(cè)文件標(biāo)準(zhǔn)化(原始文件→RINEXo-files)采用隨機(jī)格式轉(zhuǎn)換軟件或TEQC軟件可將原始文件轉(zhuǎn)換為標(biāo)準(zhǔn)的rinexo-files格式對(duì)所有的o-files文件的命名進(jìn)行標(biāo)準(zhǔn)化，如sitedoyn.yyo。注意文件名字字母要小寫(xiě)。2023/5/1616天線高改化至DHARP(天線座底部)外業(yè)天線高量測(cè)的方式根據(jù)不同的接收機(jī)天線型號(hào)有不同的量測(cè)方法GAMIT軟件里設(shè)置了不同的量高模式。建議采用天線座底部(DHARP)天線高改化至L1PHC(天線相位中心)2023/5/16172023/5/1618GPS接收機(jī)天線高改化原理、量高模式與改化方法GPS接收機(jī)天線量高模式

2023/5/1619GPS接收機(jī)天線高改化原理H為標(biāo)石中心天線參考點(diǎn)的斜高R為天線盤(pán)半徑h1為標(biāo)石中心天線參考點(diǎn)的垂高h(yuǎn)2為天線參考點(diǎn)天線座底部的垂高；h3為天線參考點(diǎn)天線相位中心的垂高標(biāo)石中心天線相位中心的垂高：

GPS測(cè)量技術(shù)工程應(yīng)用培訓(xùn)班2023/5/1620已知外業(yè)測(cè)量的斜高為1.564m求：

從網(wǎng)上查詢(xún)或下載數(shù)據(jù)文件2023/5/1621IGS跟蹤站o-files數(shù)據(jù)、導(dǎo)航電文brdc文件及精密星歷SP3文件：

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

ftp://igs.ensg.ign.fr/pub/。更新tables文件?？梢詮姆?wù)器ftp:///pub/gamit

下載更新的tables表文件其中l(wèi)untab.、soltab.、nutabl.等文件按年進(jìn)行更新pole.、ut1.、pmu.等文件每周都有更新222023/5/16文件結(jié)構(gòu)整理：對(duì)每個(gè)工程目錄做處理前的準(zhǔn)備工作，將數(shù)據(jù)文件分類(lèi)存放，便于維護(hù)和更新步驟如下：①在當(dāng)前工程目錄下建立的文件夾有:igs、rinex、brdc以及時(shí)段文件夾igs文件夾存放精密星歷文件rinex文件夾存放o-files文件brdc文件夾存放brdc文件時(shí)段文件夾，以時(shí)段號(hào)命名，所有解算都是基于該文件夾進(jìn)行的文件準(zhǔn)備232023/5/16將更新好的完好tables表文件復(fù)制到當(dāng)前工程目錄下可以手動(dòng)復(fù)制或命令cp–r../tables./。注意：這里盡量用復(fù)制！當(dāng)然，也可以用連接ln-s../../gamit/source/tables./至此，當(dāng)前工程目錄文件夾結(jié)構(gòu)應(yīng)該為：242023/5/16修改tables文件夾下的配置文件文件：記錄各測(cè)站的站點(diǎn)、站點(diǎn)名稱(chēng)、開(kāi)始觀測(cè)時(shí)間、結(jié)束觀測(cè)時(shí)間、天線高、天線量測(cè)方法、接收機(jī)類(lèi)型、軟硬件版本號(hào)、天線類(lèi)型，天線相位中心模型機(jī)構(gòu)代碼等信息。注意：AntHt、HtCod(標(biāo)志天線量測(cè)方法DHARP、SLBGP等)和AntennaType(用來(lái)做天線高改正)，ReceiverType(用來(lái)做周跳恢復(fù))，SessionStart、SessionStop(用來(lái)控制測(cè)段跨度，處理不好后面會(huì)提示錯(cuò)誤)252023/5/16嚴(yán)格按照已經(jīng)給出的跟蹤站格式編輯具體操作方法有2種：手動(dòng)操作：按照相應(yīng)文件格式對(duì)齊操作，按照標(biāo)準(zhǔn)格式.standard編輯修改。sh_upd_stnfo自動(dòng)提取觀測(cè)文件的有關(guān)信息，生成，并認(rèn)真核對(duì)信息是否正確。事先提供空模板后，再執(zhí)行命令：

sh_upd_stnfo–ref–files*.yyo262023/5/16l-file文件：記錄各測(cè)站的近似坐標(biāo)Gamit要求是球面坐標(biāo)，可以先在測(cè)站的O文件里手動(dòng)提取近似坐標(biāo)XYZ，生成l-file文件，按照格式對(duì)齊，然后保存文件，將名字重新命名為lfile.apr

，然后用gapr_to_l將lfile.apr轉(zhuǎn)換到lfile.。可以用sh_rx2apr命令自動(dòng)生成，利用svpos和svdiff模塊單點(diǎn)定位生成站點(diǎn)坐標(biāo)。也可以在測(cè)站的O文件里自動(dòng)提取近似坐標(biāo)：

grepPOSITION*.yyo>lfilerx2aprlfile<year><doy>gapr_to_llfile.aprlfile.“”<year><doy>272023/5/16sestbl.文件gambit基線處理時(shí)的模型選擇配置文件。需要仔細(xì)參照技術(shù)文檔和程序代碼，一般情況下采用默認(rèn)值，通常需要修改的地方有：ChoiceofObservable(觀測(cè)值的選擇類(lèi)型)ChoiceofExperiment(基線處理類(lèi)型)sittbl.文件(sittbl.short)各測(cè)站的精度控制指標(biāo)，一般情況下采用默認(rèn)值。高精度的已知坐標(biāo)強(qiáng)約束，待求點(diǎn)坐標(biāo)松弛約束。282023/5/16將tables表文件中需要的文件（包括luntab、soltab、nutabl、pole、ut1、pmu文件，還有一些誤差改正模型文件等）鏈接到各時(shí)段文件夾。該步驟的作用是：在數(shù)據(jù)處理各個(gè)時(shí)段的工作目錄下，建立該工作目錄與/tables目錄中部分文件的鏈接，為后面在工作目錄下逐行輸入命令進(jìn)行數(shù)據(jù)解算做準(zhǔn)備。在各時(shí)段文件夾下分別建立與tables文件夾里相應(yīng)文件對(duì)應(yīng)的鏈接，輸入命令：links.dayyrdoytest（可以參考links.day命令的用法）292023/5/16然后需要手動(dòng)建立以下鏈接，命令如下

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.另外，潮汐改正文件也是很容易出錯(cuò)的文件。以gamit10.50版本為例，Gamit默認(rèn)的潮汐改正文件名為otl.grid（鏈接前是otl_FES2004.grid）302023/5/16一定要檢查一下otl.grid是否是真實(shí)鏈接，出現(xiàn)死鏈接的情況需要手動(dòng)鏈接。命令為(在tables文件夾目錄下)：

ln-s./otl_FES2004.grid./otl.grid鏈接各個(gè)測(cè)站觀測(cè)數(shù)據(jù)在各時(shí)段文件夾下建立觀測(cè)數(shù)據(jù)連接，在各時(shí)段文件夾下輸入命令：ln–s../rinex/*doy？.yyo./ln–s../igs/*igsweek.sp3./ln–s../brdc/*doy？.yyn./312023/5/16數(shù)據(jù)文件準(zhǔn)備分步操作：在終端，輸入命令：makexp，回車(chē)，然后按照提示輸入相應(yīng)的文件數(shù)據(jù)。提示在要求輸入l-file文件時(shí)候，建議輸入lfile.apr文件。批操作

sh_makexp–expttest–orbt

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

lfile.–sinfo1500005760需要準(zhǔn)備的文件：O文件、N文件；文件；rcvant.dat文件，接收機(jī)天線類(lèi)型。輸出文件：數(shù)據(jù)處理(分步處理)322023/5/16精密星歷處理sh_sp3fit–figs15250.sp3–oigsf–d2009088–rBERNE–t–u需要準(zhǔn)備的文件

sp3文件；ut1.，時(shí)長(zhǎng)變化表；

pole.，極移表；

nutabl.，章動(dòng)表；

soltbl.，太陽(yáng)星歷表；luntab.，月亮星歷表；leap.sec，跳秒；svnav.dat，衛(wèi)星天線類(lèi)型332023/5/16輸出文件

gigsf9.088；說(shuō)明：g-files文件所有衛(wèi)星初始軌道坐標(biāo)文件tigsf9.088；文件檢查sh_check_sess–sess088–typegfile–filegigsf9.088需要準(zhǔn)備的文件：gtest9.088；輸出文件：廣播星歷處理makejbrdc0880.09njbrdc9.088需要準(zhǔn)備的文件：brdc0880.09n；輸出的文件：jbrdc9.088342023/5/16文件檢查sh_check_sess–sess088–typejfile–filejbrdc9.088需要準(zhǔn)備的文件jbrdc9.088輸出文件形成中間文件makex

test.makex.batch需要準(zhǔn)備的文件

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

生成bat文件輸入命令

fixdrvdtest9.088需要準(zhǔn)備的文件：dtest9.088，以及該文件所列出的所有文件輸出的文件：btest9.088.bat，以及btest9.num基線解算輸入命令

cshbtest9.bat需要準(zhǔn)備文件：btest9.bat，以及btest9.num輸出文件：q-、o-、h-等所有解算的結(jié)果文件。362023/5/16GAMIT基線處理結(jié)果質(zhì)量評(píng)價(jià)指標(biāo)檢查所有測(cè)站點(diǎn)是否都參與計(jì)算(LC_HELP)標(biāo)準(zhǔn)均方根nrms解算結(jié)果為0.20左右若是nrms大于0.5，表明處理結(jié)果是有問(wèn)題的（例如：周跳沒(méi)有修復(fù)；測(cè)站的起算點(diǎn)坐標(biāo)有問(wèn)題等）。若是nrms小于0.5，則認(rèn)為成功（合理）解算，否則需要檢查原因，重新處理；檢查各基線文件處理精度是否符合gamit的實(shí)際精度質(zhì)量評(píng)估372023/5/164定制式“玩轉(zhuǎn)”GAMIT軟件主要內(nèi)容bash編程語(yǔ)言介紹基于bash編程語(yǔ)言自主編寫(xiě)“定制化”批處理腳本(*.sh)【算例】：某工程控制網(wǎng)GPS數(shù)據(jù)4個(gè)觀測(cè)時(shí)段、國(guó)標(biāo)D級(jí)(*.sh文件)382023/5/162023/5/16使用GAMIT軟件進(jìn)行基線解算有一定程度的了解之后，可以不用每次只執(zhí)行一條命令，可利用UNIX/LINX系統(tǒng)自帶的shell編程語(yǔ)言，定制適合自己的批處理腳本Shell命令。值得注意的是：Shell腳本命令只對(duì)數(shù)據(jù)處理進(jìn)行批處理，數(shù)據(jù)準(zhǔn)備及文件準(zhǔn)備建議手動(dòng)操作。以某城市CORS2012年doy=201第0時(shí)段基線解算為例，定制的批處理腳本文件2010_cs.sh內(nèi)容如下所示。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技術(shù)討論：經(jīng)驗(yàn)交流會(huì)互動(dòng)環(huán)節(jié)與答疑2023/5/16475批處理“玩轉(zhuǎn)”GAMIT軟件

主要內(nèi)容【算例】：中國(guó)大陸8個(gè)IGS站(CHAN、BJFS、URUM、GUAO、LHAZ、WUHN、SHAO、KUNM)2013年doy=001~005高精度基線處理(sh_gamit)批處理(sh_gamit)介紹與工程應(yīng)用使用前提sh_gamit操作GAMIT軟件的基本流程482023/5/16【算例】中國(guó)大陸8個(gè)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文件夾，并存放對(duì)應(yīng)數(shù)據(jù),并在./2013目錄下，終端執(zhí)行：sh_setup-yr2013編輯./2013/tables/process.defaults編輯./2013/tables/sites.defaults準(zhǔn)備準(zhǔn)備lfile.apr和lfile.編輯sestbl.和sittbl.在./2013目錄下，終端執(zhí)行：-d:所處理數(shù)據(jù)所在的年、年積日-s：2013001003-copt：待壓縮的文件類(lèi)型-dopt：待刪除的文件類(lèi)型-yrext：指定待處理數(shù)據(jù)年份及年積日-pres：NO、YES和ELEVsh_gamit-exptcors-d2013001002003-presELEV-orbitIGSF-coptxkp-doptscao專(zhuān)題實(shí)操練習(xí)

應(yīng)用GAMIT/CosaGPS軟件進(jìn)行高精度GPS工程控制網(wǎng)數(shù)據(jù)處理與精度評(píng)估

532023/5/16主要內(nèi)容基于IGS站的框架網(wǎng)構(gòu)建與優(yōu)化設(shè)計(jì)IGS站點(diǎn)數(shù)據(jù)(GPS觀測(cè)資料、三維框架坐標(biāo))獲取途徑大地測(cè)量參考框架(ITRF系列框架與參考?xì)v元的基準(zhǔn)統(tǒng)一及其轉(zhuǎn)換)全球定位系統(tǒng)（GPS）測(cè)量規(guī)范(GB-T18314-2009)解讀【算例】：某A級(jí)GPS變形監(jiān)測(cè)網(wǎng)高精度數(shù)據(jù)處理中框架網(wǎng)構(gòu)建542023/5/16利用“定制式”腳本操作GAMIT軟件進(jìn)行高精度GPS基線處理（可獲得高精度GPS基線向量、天頂對(duì)流層延遲參數(shù)）高精度GPS基線向量質(zhì)量評(píng)估【算例】：AB級(jí)GPS基線分量、長(zhǎng)度重復(fù)性計(jì)算應(yīng)用GLOBK/CosaGPS軟件進(jìn)行GPS網(wǎng)平差與精度評(píng)估數(shù)據(jù)計(jì)算報(bào)告與技術(shù)總結(jié)報(bào)告撰寫(xiě)552023/5/161基于IGS站的框架網(wǎng)優(yōu)化設(shè)計(jì)與構(gòu)建IGS站網(wǎng)2023/5/1656SOPAC網(wǎng)站：/

2023/5/1657IERS網(wǎng)站：

2023/5/16582ITRF系列框架與參考?xì)v元的基準(zhǔn)統(tǒng)一2023/5/16592023/5/16603全球定位系統(tǒng)(GPS)測(cè)量規(guī)范解讀術(shù)語(yǔ)和定義基本規(guī)定級(jí)別劃分和測(cè)量精度布設(shè)的原則選點(diǎn)埋石儀器觀測(cè)外業(yè)成果記錄數(shù)據(jù)處理成果驗(yàn)收與上交資料612023/5/162023/5/16623.1基本規(guī)定GPS測(cè)量采用2000國(guó)家大地坐標(biāo)系統(tǒng)采用GPS時(shí)間系統(tǒng)手簿記錄宜采用世界協(xié)調(diào)時(shí)(UTC)接收機(jī)應(yīng)經(jīng)檢定合格，在有效期內(nèi)使用采用中誤差為精度指標(biāo)，以2倍為極限誤差當(dāng)需要提供80或54或其他坐標(biāo)成果時(shí)，應(yīng)按坐標(biāo)轉(zhuǎn)換方法求得。2023/5/163.2級(jí)別劃分和測(cè)量精度級(jí)別劃分GPS測(cè)量按照精度和用途分為A~E級(jí)5等級(jí)測(cè)量精度A級(jí)GPS網(wǎng)由CORS站構(gòu)成B~E級(jí)的精度應(yīng)不低于表2要求632023/5/1664（1）按工作原理分：碼相關(guān)型接收機(jī)平方型接收機(jī)混合型接收機(jī)（2）按信號(hào)通道的類(lèi)型分：多通道接收機(jī)序貫通道接收機(jī)多路復(fù)用通道接收機(jī)（3）根據(jù)接收信號(hào)的頻率分：?jiǎn)晤l接收機(jī)（L1）雙頻接收機(jī)（L1+L2）（4）根據(jù)接收機(jī)的用途分：導(dǎo)航型測(cè)量型授時(shí)型3.3儀器GPS接收機(jī)的分類(lèi)2023/5/16接收機(jī)選用A級(jí)網(wǎng)測(cè)量按有關(guān)規(guī)定執(zhí)行B~E級(jí)GPS網(wǎng)按表4規(guī)定執(zhí)行GPS接收機(jī)的檢驗(yàn)新購(gòu)置GPS接收機(jī)的檢驗(yàn)一般性檢視、通電檢驗(yàn)、試測(cè)檢驗(yàn)舊GPS接收機(jī)的檢驗(yàn)一般性檢視、通電檢驗(yàn)、試測(cè)檢驗(yàn)的準(zhǔn)備部分652023/5/16663.4觀測(cè)A級(jí)網(wǎng)B~E級(jí)GPS網(wǎng)觀測(cè)的基本技術(shù)規(guī)定應(yīng)符合表5的要求2023/5/163.5外業(yè)成果記錄記錄類(lèi)型觀測(cè)數(shù)據(jù)、測(cè)量手簿、其他記錄記錄內(nèi)容原始觀測(cè)數(shù)據(jù)、對(duì)應(yīng)觀測(cè)值的GPS時(shí)間、測(cè)站與接收機(jī)初始信息、天線高記錄要求及時(shí)、認(rèn)真、準(zhǔn)確672023/5/163.6數(shù)據(jù)處理原始格式—>RINEX格式方法GPS數(shù)據(jù)文件的組織外業(yè)數(shù)據(jù)質(zhì)量檢驗(yàn)數(shù)據(jù)刪除率：不宜大于10%復(fù)測(cè)基線的長(zhǎng)度差682023/5/1669各級(jí)GPS測(cè)量的精度指標(biāo)相鄰點(diǎn)間基線長(zhǎng)度標(biāo)準(zhǔn)差例如：a=5mmb=1ppmD=5km2023/5/1670外業(yè)數(shù)據(jù)質(zhì)量檢驗(yàn)同步環(huán)閉合差2023/5/16獨(dú)立環(huán)閉合差及附和路線坐標(biāo)閉合差7172GPS測(cè)量數(shù)據(jù)處理的流程四個(gè)階段：1.數(shù)據(jù)傳輸2.格式轉(zhuǎn)換（可選）3.基線處理4.網(wǎng)平差2023/5/162023/5/163.7成果驗(yàn)收與上交資料成果驗(yàn)收73上交資料2023/5/16744框架網(wǎng)/工程網(wǎng)高精度數(shù)據(jù)處理利用GAMIT/CosaGPS軟件進(jìn)行高精度GPS數(shù)據(jù)處理算例分析GAMIT基線處理框架網(wǎng)基線處理技巧與方法、質(zhì)量評(píng)估方法參考框架與參考?xì)v元的統(tǒng)一及轉(zhuǎn)換方法工程網(wǎng)基線處理技巧與方法、質(zhì)量評(píng)估方法CosaGPS網(wǎng)平差框架網(wǎng)平差處理與精度評(píng)估方法工程網(wǎng)平差處理與精度評(píng)估方法2023/5/1675AB級(jí)基線分量重復(fù)性計(jì)算2023/5/16762023/5/16GPS工程控制網(wǎng)數(shù)據(jù)處理算例數(shù)據(jù)來(lái)源某測(cè)繪公司采用Leica接收機(jī)于2009年3月28~29日在某工程項(xiàng)目中布測(cè)了一個(gè)高精度GPS控制網(wǎng)，我們?nèi)≡摼W(wǎng)中6個(gè)GPS控制點(diǎn)，共4個(gè)觀測(cè)時(shí)段數(shù)的數(shù)據(jù)為例進(jìn)行練習(xí)預(yù)期目標(biāo)熟練掌握GAMIT+CosaGPS高精度GPS數(shù)據(jù)處理技術(shù)與方法GAMIT基線處理CosaGPS平差計(jì)算77工程案例“實(shí)戰(zhàn)”演練【工程算例】：某縣城GPSE級(jí)網(wǎng)背景資料介紹某縣城GPSE級(jí)網(wǎng)數(shù)據(jù)處理全過(guò)程學(xué)員實(shí)操與答疑782023/5/162023/5/1679CosaGPS軟件技術(shù)討論：經(jīng)驗(yàn)交流會(huì)互動(dòng)環(huán)節(jié)與答疑2023/5/1680第二部分:TRACK/TTC812023/5/16TRACK主要內(nèi)容TRACK模塊介紹Track.cmd控制策略文件介紹TRACK模塊操作過(guò)程(【算例】：track_example詳細(xì)講解)【工程算例】：1基準(zhǔn)站+1流動(dòng)站GPS數(shù)據(jù),采樣間隔1sTRACK模塊工程應(yīng)用學(xué)員實(shí)操與答疑822023/5/16TRACK：GPS差分相位動(dòng)態(tài)定位程序模塊，可以處理GPS動(dòng)態(tài)單歷元數(shù)據(jù)深入文檔可以參考網(wǎng)站/~tah/track_example目前的版本是TRACKVersion1.21GPS動(dòng)態(tài)定位程序(2009/01/11)

使用方法track-f<commandfile>-a<ambiguityfile>-d<day>-w<week>-s<S01><S02>..<S10>2023/5/1683<commandfile>是一個(gè)必須的含有命令列表來(lái)控制程序運(yùn)行的文件<ambiguityfile>是一個(gè)可選的含有一個(gè)改進(jìn)設(shè)置整周偏差參數(shù)和配置的文件<day>是字符串，用于替代命令文件行的<day>，例如：如果給定選項(xiàng)-d222，則bas1<day>0.03o將變?yōu)閎as12220.03o<week>是字符串，用于替代命令文件的任何<week>，為nav_file名字提供sp3文件的周數(shù)<S01>,<S02>..<S10>是10個(gè)字符串，在命令文件中用于替代，例如<S01>字符串，在命令文件中將替代第一個(gè)字符串，<S02>字符串，在命令文件中將替代第二個(gè)字符串，依次類(lèi)推2023/5/1684track.cmd控制策略文件2023/5/16852023/5/16862023/5/16872023/5/16882023/5/1689TRACK模塊使用案例分析算例數(shù)據(jù)流動(dòng)站QXJA觀測(cè)時(shí)間：2012年7月21日0~12日基準(zhǔn)站：選用IGS站LHAZ精密星歷：igr16976.sp3執(zhí)行命令：track-ftrack.cmd-d203_12-w16976>!TRAK_203_12.out2023/5/16902023/5/1691TTC軟件第三部分:GMT922023/5/16主要內(nèi)容GMT軟件介紹、下載方式“手把手”教會(huì)UBUNTU操作系統(tǒng)下GMT軟件的安裝方法GMT軟件的基本命令與使用方法應(yīng)用GMT軟件實(shí)現(xiàn)大地測(cè)量成果圖形表達(dá)932023/5/16GMT工程應(yīng)用常用命令常用參數(shù)安裝GMT942023/5/161GMT的安裝通用制圖工具(TheGenericMappingTools,GMT)是一個(gè)被地理學(xué)界廣泛使用的繪圖工具，可以完成海岸線、國(guó)界、河流等的繪制GMT繪圖工具的獲取2023/5/1695GMT安裝準(zhǔn)備工作在安裝GMT之前，需要先安裝如下軟件NetCDF庫(kù)，版本3.4甚至更高(可以從網(wǎng)站下載)C編譯器(一般Linux系統(tǒng)都已安裝，如果沒(méi)有可以從網(wǎng)站下載)此外，還需要安裝ps文檔閱讀器，如ghostview等2023/5/1696GMT安裝將整個(gè)GMT安裝文件夾(gmt)拷貝到：/home/user/下或任意目錄下打開(kāi)終端，進(jìn)入到gmt文件目錄下(cd/home/user/gmt)，然后輸入：sudoshinstall_gmtGMTparam.txt2023/5/1697安裝完成后配置環(huán)境以管理員身份(sudopasswdroot)打開(kāi)/etc下的environment文件，在PATH行后面追加：/usr/local/gmt/bin:$PATH再最后行添加：ETCDFHOME=“/usr/local/netcdf”注意：第一個(gè)不用引號(hào)，因?yàn)橐呀?jīng)存有引號(hào)，第二個(gè)需要加上引號(hào)。然后，保存退出，重啟系統(tǒng)。在終端測(cè)試是否安裝成功：輸入psxy等進(jìn)行測(cè)試2023/5/1698GMT軟件安裝演示過(guò)程992023/5/162Linux相關(guān)知識(shí)點(diǎn)重定向：使用bash可以方便地用<和>實(shí)現(xiàn)輸出輸入的重定向cmd>filecmd>>filecmd<filecmd<<textcmd<<<wordcmd<>filecmd>|file1002023/5/16管道符(|)：利用Linux所提供的管道符“|”將兩個(gè)命令隔開(kāi)，管道符左邊命令的輸出就會(huì)作為管道符右邊命令的輸入利用一個(gè)管道#rpm–qa|greplicq利用多個(gè)管道#cat/etc/passwd|grep/bin/bash|wc-11012023/5/163GMT常用參數(shù)參數(shù)含義-B定義底圖和坐標(biāo)軸的刻度、注記和標(biāo)簽-H輸入/輸出表的頭記錄行數(shù)-J選擇所用的地圖投影或坐標(biāo)轉(zhuǎn)換-K允許在當(dāng)前圖件中追加圖畫(huà)-O允許將當(dāng)前圖件追加已存在的圖件中-P將圖件頁(yè)面設(shè)置成縱向模式(默認(rèn)是橫向)-R定義地圖/圖件的區(qū)域-U添加時(shí)間戳，默認(rèn)位于頁(yè)面的左下角-V備注默認(rèn)，報(bào)告運(yùn)行的進(jìn)度-X設(shè)置圖件在頁(yè)面中的起始x坐標(biāo)-Y設(shè)置圖件在頁(yè)面中的起始y坐標(biāo)-b設(shè)定輸入/輸出文件采用二進(jìn)制格式-c設(shè)定圖件的備份數(shù)-f設(shè)定每列數(shù)據(jù)的具體格式-：假定地理數(shù)據(jù)的格式為(緯度、經(jīng)度)，而不是(經(jīng)度、緯度)1022023/5/16數(shù)據(jù)域或者地圖區(qū)域：-R選項(xiàng)

-R參數(shù)用來(lái)定義作圖的地圖或者是感興趣數(shù)據(jù)所在區(qū)域，通常采用以下三種方式進(jìn)行定義：

-Rxmin/xmax/ymin/ymax

-Rxlleft/ylleft/xuright/yurightr

-Rgridfile1032023/5/16坐標(biāo)變換和地圖投影：-J選項(xiàng)-Ji[parameters/]scale，其中i為一小寫(xiě)字母，用來(lái)選擇具體的地圖投影方式，parameters為0或者是一組以斜畫(huà)線分割的其他投影參數(shù)，scale為地圖比列尺，按距離/度或者1：xxxxx的格式給出-JI[parameters/]width，I是一個(gè)大寫(xiě)字母，用來(lái)選擇具體的投影方式，parameters為0或者是一組以斜畫(huà)線分割的其他投影參數(shù)，width為地圖的寬度1042023/5/16地圖框架和坐標(biāo)軸注記：-B選項(xiàng)使用格式為-B[p|s]xinfo[/yinfo[/zinfo]][:.“標(biāo)題”:][W|w][E|e][S|s][N|n][Z|z[+]]，可以通過(guò)大小寫(xiě)開(kāi)關(guān)來(lái)選擇需要作圖的邊界(或軸)xinfo，yinfo和zinfo的格式如下：[a]tick[m|c][ftick[m|c][gtick[m|c]][l|p][:“axislabel”:][:,”unitlabel”:]，其中a、f、g分別為坐標(biāo)軸標(biāo)尺的大、中和小的間隔，“axislabel”為坐標(biāo)軸的標(biāo)簽以及“unitlabel”為坐標(biāo)軸的單位1052023/5/164GMT常用命令創(chuàng)建圖件底圖命令：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繪制線、多邊形和標(biāo)志的命令：psxy繪制中國(guó)邊界(CHINA_Boundary.dat)psxyChina_Boundary.dat-Jl103/35/28/40/1:28000000-R64/138/18/54-m>China_Boundary.ps1082023/5/16繪制中國(guó)省界(CHINA_Prov.dat)psxyChina_Province.dat-Jl103/35/28/40/1:28000000-R64/138/18/54-m>China_Province.ps1092023/5/16繪制速度矢量、十字和楔形標(biāo)志命令：psvelopsvelop中-S的含義如下：

-Sevelscale/confidence/fontsize：繪制(N，E)的速度橢圓-Snbarscale：繪制各向異形柱-Swwedge_scale/wedge_mag：繪制旋轉(zhuǎn)扇形-Sxcross_scale：繪制應(yīng)變十字1102023/5/16這條命令將生成具有邊界紅色內(nèi)部綠色的箭頭和邊界紅色橢圓的速度矢量圖，數(shù)據(jù)的置信期間為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繪制震源機(jī)制解命令：psmeca其中的-S參數(shù)的具體含義如下：

-Sascale[/fontsize[/offset[u]]]：繪制AkiandRichards格式的震源機(jī)制-Scscale[/fontsize[/offset[u]]]：繪制HarvardCMT格式的震源機(jī)制-Sm|d|zscale[/fontsize[/offset[u]]]：繪制地震矩張量-Spscale[/fontsize[/offset[u]]]：繪制分離格式的震源機(jī)制解-Sx|y|tscale[/fontsize[/offset[u]]]：繪制主軸1142023/5/16下面這條命令將繪制一個(gè)正斷層的CMT震源機(jī)制圖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工程應(yīng)用GMT繪制速度場(chǎng)(分步操作)生成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繪制海岸線、國(guó)界、常用河流信息pscoast-R73/136/3/54-Jq-O-K-Dh-A300-I1-W0.1p-N1/1p-I1/1p-P>>velf.ps

由速度文件(vel_data.txt)繪制速度場(chǎng)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軟件繪制速度場(chǎng)效果圖2023/5/16118利用GMT軟件繪制GPS速度場(chǎng)(腳本)#!/bin/csh#設(shè)定該腳本所調(diào)用的shell，該程序調(diào)用的是csh。gmtsetBASEMAP_TYPEPLAIN#設(shè)定地圖地圖樣式為PLAIN，另一個(gè)選項(xiàng)是FANCY。setrange=70/140/10/60#設(shè)定地圖的坐標(biāo)范圍。setprojection=q96/1:32000000#設(shè)定地圖的投影格式和比例尺大小。1192023/5/16setticks=10f5/10f5#設(shè)定坐標(biāo)注記的方式。setpsfile=vel.ps#設(shè)定輸出文件的名稱(chēng)。makecpt-Cglobe-T-10000/10000/200-Z-V>tmp.cpt#設(shè)定自定義顏色查找表。psbasemap-R${range}-J${projection}-K-B${ticks}WSen-X2-Y2>${psfile}#繪制底圖，采用前述的投影格式、范圍。1202023/5/16grdimageglobal.grd-Ctmp.cpt-Jq-R-O-K>>${psfile}#繪制二進(jìn)制格網(wǎng)數(shù)據(jù)Asia_3m.grd，采用系統(tǒng)自帶的顏色查找表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速度場(chǎng)，使用awk命令將原始數(shù)據(jù)重排列為psvelo命令采用的格式，-W為畫(huà)筆屬性，-G為矢量填色，-A為箭頭屬性。1212023/5/16利用GMT軟件繪制GPS速度場(chǎng)效果圖1222023/5/16第四部分:GLOBK/CosaGPS1232023/5/16分步式/批處理“玩轉(zhuǎn)”GLOBK軟件主要內(nèi)容分步操作GLOBK軟件(分步操作的基本命令與使用方法)globk_com.cmd控制文件介紹glorg_com.cmd控制文件介紹應(yīng)用GLOBK軟件進(jìn)行測(cè)站坐標(biāo)時(shí)間序列分析應(yīng)用GLOBK軟件進(jìn)行測(cè)站速度場(chǎng)(水平和垂直)分析批處理(sh_globk、sh_glorg)介紹sh_globk/sh_glorg操作GLOBK軟件的基本流程【算例】：GLOBK10.50軟件自帶算例詳細(xì)講解1242023/5/16GLOBK安裝工程應(yīng)用軟件使用常用文件與格式說(shuō)明軟件模塊GLOBK1252023/5/161GLOBK軟件的安裝安裝GLOBK,輸入y1262023/5/162GLOBK軟件的模塊格式轉(zhuǎn)換模塊(htoglb)將處理GPS,VLBI,SLR數(shù)據(jù)的分析軟件的解算文件轉(zhuǎn)換成GLOBK軟件所需要的二進(jìn)制文件(h-文件)。目前支持的幾類(lèi)文件：GAMIT軟件h-文件關(guān)于GPS(或其他空間大地測(cè)量技術(shù))SINEX格式文件FONDA軟件h-文件1272023/5/16JPL機(jī)構(gòu)提供的Stacov文件包含站坐標(biāo)和速度場(chǎng)的SLR/GSFC文件包含站坐標(biāo)和速度場(chǎng)的VLBI/GSFC文件運(yùn)算模塊(glred,globk,glorg)glred模塊通過(guò)調(diào)用globk模塊分析單天解，對(duì)于分析基準(zhǔn)站網(wǎng)(CORS)非常適用，其生成的解文件可以用來(lái)形成時(shí)間序列g(shù)lobk模塊是GLOBK軟件的主模塊，實(shí)現(xiàn)該軟件的功能1282023/5/16glorg模塊可以為平差結(jié)果定義參考框架，具體通過(guò)固定(或約束)站坐標(biāo)和速度由坐標(biāo)轉(zhuǎn)換來(lái)實(shí)現(xiàn)glorg模塊可以單獨(dú)運(yùn)行，也可以被globk、glred模塊調(diào)用cmd_file是glorg和globk的控制文件，里面包含解算的策略等GMT圖形應(yīng)用模塊sh_plotcrd1292023/5/16sh_globk_scattermultibasesh_plotvel主要功能是利用GMT軟件繪制時(shí)間序列、速度場(chǎng)等圖形，可用于分析數(shù)據(jù)質(zhì)量和測(cè)站的地殼運(yùn)動(dòng)等情況其他輔助模塊一類(lèi)為GLOBK模塊服務(wù)的：glist、glsave一類(lèi)是用于框架之間和板塊運(yùn)動(dòng)分析的：corcom，cvframe，velrot1302023/5/163常用文件與格式說(shuō)明輸入文件GPS,SLR,VLBI和SINEX文件隨機(jī)特征可由apr_XXX和mar_XXX表述數(shù)據(jù)文件：二進(jìn)制h-文件控制文件：指令cmd文件GLOBK結(jié)果文件中的NEU坐標(biāo)GLOBK的輸出文件一般為*.prt和*.org在給出ITRF坐標(biāo)同時(shí)還給出了新的NEU坐標(biāo)，它與站坐標(biāo)定義的NEU不同，這種坐標(biāo)類(lèi)似于平面坐標(biāo)，屬于圓錐投影1312023/5/16運(yùn)行GLOBK軟件時(shí)的注意事項(xiàng)GLOBK是基于線性模型的：測(cè)站坐標(biāo)改正值大于10米或軌道參數(shù)的改正值大于100米時(shí)需要進(jìn)行數(shù)據(jù)的再處理GLOBK不能解決因周跳未得到完全探測(cè)、數(shù)據(jù)質(zhì)量差或大氣層延遲模型誤差所帶來(lái)的問(wèn)題GLOBK不能進(jìn)行整周模糊度的解算：前期數(shù)據(jù)處理解算必須完成整周模糊度的解算1322023/5/164GLOBK軟件的使用GLOBK數(shù)據(jù)處理流程將ASCII格式的h-file轉(zhuǎn)換成可被globk讀取的二進(jìn)制h-file，然后運(yùn)行g(shù)lred/globk以獲得測(cè)站坐標(biāo)的時(shí)間序列通過(guò)時(shí)間序列分析，確定具有異常域的特定站或特定歷元，在earthquakefile中，運(yùn)用rename命令刪除具有異常域的特定站的特定歷元或直接刪除對(duì)應(yīng)的h-file133運(yùn)行g(shù)lobk將單時(shí)段解的h-file合并成一個(gè)h-file，其代表在所選擇的時(shí)間跨度里測(cè)站的平均坐標(biāo)使用合并后的h-file，再次運(yùn)行g(shù)lred/globk獲得時(shí)間序列，而運(yùn)行g(shù)lobk/glorg則可獲得測(cè)站速度利用GLOBK進(jìn)行網(wǎng)平差(分布操作)建立工作目錄，比如：globk_test，在該工作目錄下還應(yīng)包括以下三個(gè)子目錄：1342023/5/16glbf：用于存儲(chǔ)經(jīng)htoglb命令運(yùn)行所生成的二進(jìn)制H-文件soln：存儲(chǔ)globk的控制文件(globk_com.cmd、glorg_com.cmd)、ASCII文本H-文件、GAMIT基線結(jié)果文件以及經(jīng)平差處理后的平差結(jié)果文件(*.prt)，并在此目錄下運(yùn)行GLOBK軟件tables：pmu.unso、測(cè)站先驗(yàn)坐標(biāo)文件(*.apr)、經(jīng)htoglb生成的*.svs參數(shù)文件2023/5/16135生成二進(jìn)制H-文件將svna