匯報(bào)-岳彩亞-GAMIT課件

GAMIT數(shù)據(jù)處理匯報(bào)人：岳彩亞

1主要內(nèi)容

●Linux操作系統(tǒng)

●shell腳本編寫

●

GAMIT基線處理

●算例分析

●總結(jié)2Linux操作系統(tǒng)

1Linux簡(jiǎn)介2Linux的發(fā)展歷史

3Linux系統(tǒng)的特點(diǎn)3Linux和Windows（xp、7、8）以及MaxOS一樣，都是一種操作系統(tǒng)。整個(gè)計(jì)算機(jī)系統(tǒng)中可以分為四塊內(nèi)容：硬件、內(nèi)核、系統(tǒng)調(diào)用、應(yīng)用程序，其中內(nèi)核和系統(tǒng)調(diào)用的組合統(tǒng)稱為操作系統(tǒng)。如圖所示。Linux是一種交互式操作系統(tǒng)。除此之外，還能使得多個(gè)用戶從不同的終端同時(shí)操作主機(jī)。也就是說(shuō)一個(gè)Linux系統(tǒng)上可以有多個(gè)用戶，這些用戶之間各自文件，各自的目錄，各自管理著自己的任務(wù)。這樣的一種系統(tǒng)稱為分時(shí)段操作系統(tǒng)。

Linux簡(jiǎn)介4Linux簡(jiǎn)介硬件內(nèi)核系統(tǒng)調(diào)用應(yīng)用程序操作系統(tǒng)51、為了研發(fā)分時(shí)操作系統(tǒng)，在各大學(xué)、科研機(jī)構(gòu)以及商業(yè)項(xiàng)目都啟動(dòng)了自己的計(jì)劃，一個(gè)項(xiàng)目叫做“Multics”是當(dāng)時(shí)最先進(jìn)的。2、在Multics開發(fā)后期，遇到了一些困難，他的一個(gè)開發(fā)團(tuán)體—貝爾實(shí)驗(yàn)室退出了這個(gè)項(xiàng)目，開發(fā)自己的操作系統(tǒng)，叫“UNIX”。3、九十年代，UNIX在服務(wù)器市場(chǎng)成為了主流的操作系統(tǒng)，開始的時(shí)候是免費(fèi)于使用者，但當(dāng)用戶逐漸增多時(shí)，UNIX系統(tǒng)開始變的商業(yè)化，并且越來(lái)越昂貴Linux發(fā)展歷史64、1991年10月，LinusTorvalds（Linux之父）開始入手，并且雇傭了一批黑客和一些對(duì)操作系統(tǒng)開發(fā)有興趣的計(jì)算機(jī)專家，開始研究自己的自己的操作系統(tǒng)，也就是Linux系統(tǒng)。

LinusTorvalds（芬蘭赫爾辛基大學(xué)）Linux發(fā)展歷史71Linux系統(tǒng)中大都為開源代碼并且免費(fèi)提供給用戶，用戶可以查看源代碼和修改定制以及再發(fā)布。2兼具圖形界面操作（需要使用帶有桌面環(huán)境的發(fā)行版）和完全的命令行操作，可以只用鍵盤完成一切操作，并且資源占用少、運(yùn)行速度快。3Linux系統(tǒng)可以實(shí)現(xiàn)多用戶共同操作同一個(gè)系統(tǒng)，且不會(huì)相互影響。Linux的特點(diǎn)8主要內(nèi)容

●Linux操作系統(tǒng)

●

shell簡(jiǎn)單介紹

●GAMIT基線處理

●算例分析

●總結(jié)9shell介紹shell語(yǔ)言的分類shell在GAMIT中的用途10

shell可以分為以下幾種

/bin/sh/bin/dash

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

/bin/csh

shell介紹11查看當(dāng)前系統(tǒng)默認(rèn)的shell編譯語(yǔ)言命令：%echo$SHELL查看當(dāng)前系統(tǒng)支持的shell編譯語(yǔ)言命令：%more/etc/shells改變默認(rèn)編譯語(yǔ)言的命令：%chsh-s/bin/*sh(bash)，重啟后即改變成功。shell語(yǔ)言分類12GAMIT是由一些程序模塊組成，并且這些程序模塊可以實(shí)現(xiàn)調(diào)用，他們?cè)贕AMIT數(shù)據(jù)處理過(guò)程中起到的作用互不相同。使用shell腳本的編寫可以將不同的程序模塊組合起來(lái)實(shí)現(xiàn)不同的作用。GAMIT自身帶了一些批處理命令，這些命令只是達(dá)到固定的結(jié)果，通過(guò)編寫shell腳本語(yǔ)言可以實(shí)現(xiàn)新的批處理過(guò)程，節(jié)約了時(shí)間，提高效率。shell在GAMIT用途13主要內(nèi)容

●Linux操作系統(tǒng)

●shell簡(jiǎn)單介紹

●

GAMIT基線處理

●算例分析

●總結(jié)14GAMIT基線處理：

參數(shù)文件更新數(shù)據(jù)處理軟件分類分步數(shù)據(jù)處理流程

自動(dòng)批處理結(jié)果評(píng)定15針對(duì)大區(qū)域GPS高精度定位解算，通常采用專業(yè)機(jī)構(gòu)開發(fā)的高精度科研軟件進(jìn)行處理。目前被廣泛使用且較為成熟的GPS數(shù)據(jù)處理軟件可以分為兩類：一、基于非差模式的定規(guī)定位軟件美國(guó)航空航天局噴氣動(dòng)力實(shí)驗(yàn)室開發(fā)的GIPSY/OASIS德國(guó)地學(xué)研究中心開發(fā)的EPOS中國(guó)武漢大學(xué)GPS研究中心開發(fā)的PANDA數(shù)據(jù)處理軟件分類：16二、基于雙差模式的定軌定位軟件瑞士伯爾尼大學(xué)天文學(xué)院開發(fā)的BERNESE美國(guó)麻省理工學(xué)院和美國(guó)斯克利普斯海洋研究所開發(fā)的GAMIT/GLOBK數(shù)據(jù)處理軟件分類：17每年更新一次：nutabl.soltal.luntab.leap.sec每月更新一次：dcb.dat(編碼文件)每日更新一次：ut1（地球自轉(zhuǎn)參數(shù)）pole.當(dāng)出現(xiàn)新的接收機(jī)：rcvant.dat新的衛(wèi)星發(fā)射時(shí)：svnav.datftp:///pub/gamit/tables（使用迅雷下載）參數(shù)文件更新：18一、創(chuàng)建相關(guān)文件夾：下面是創(chuàng)建文件夾的名字，是與控制文件process.defaults中默認(rèn)的路徑是相對(duì)應(yīng)的，如果想修改，必須也把控制文件中默認(rèn)的名字做相應(yīng)的修改。處理基線需用到的文件夾：文件夾igs：用來(lái)存放SP3精密星歷。文件夾brdc：用來(lái)存放廣播星歷，也就是n文件。文件夾rinex：用來(lái)存放觀測(cè)文件，o文件。文件夾met：用來(lái)存放氣象觀測(cè)數(shù)據(jù)文件，一般研究對(duì)象是氣象方面的時(shí)能用的到。文件夾“年紀(jì)日”：如果處理的是2012年，年積日001天的數(shù)據(jù)，就將文件改為001文件夾tables：用來(lái)存放控制文件，處理用到的參數(shù)文件，即所謂的全局行文件。sh_links.tables-yr~|links.tables數(shù)據(jù)處理流程19與網(wǎng)平差相關(guān)的文件夾hfile文件夾：存放文本形式的h文件glbf文件夾：存放二進(jìn)制形式的h文件gslon文件夾：進(jìn)行平常操作

tables：控制文件，可以直接使用基線解算的tables數(shù)據(jù)處理流程20二、IGS核心站數(shù)據(jù)的獲取廣播星歷：/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······數(shù)據(jù)處理流程21觀測(cè)文件（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)······數(shù)據(jù)處理流程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文件拷貝到對(duì)應(yīng)的創(chuàng)建文件夾中。數(shù)據(jù)處理流程23三、文件配置（1）制作測(cè)站信息文件：

拷貝GAMIT/GLOBK安裝目錄tables/文件夾下的文件至解算目錄的tables/文件夾，通過(guò)cd命令進(jìn)入該文件夾（后者），執(zhí)行：(執(zhí)行之前一定要先設(shè)置的文件的頭)

%sh_upd_stnfo-files../rinex/*o數(shù)據(jù)處理流程24（2）制作臺(tái)站先驗(yàn)坐標(biāo)lfile.執(zhí)行：

%grepPOSITION../rinex/*o>lfile.rnx（提取o文件中的測(cè)站坐標(biāo)）

%rx2aprlfile.rnx$year$doy（轉(zhuǎn)化成GAMIT識(shí)別的文件）

%gapr_to_llfile.rnx.aprlfile.$year$doy（將空間坐標(biāo)轉(zhuǎn)換成球面坐標(biāo)）這樣，就得到了測(cè)站信息文件lfile.（rinex/文件夾下），將其拷貝到tables/文件夾即可。數(shù)據(jù)處理流程25（3）文件配置數(shù)據(jù)處理流程sestbl.（策略參數(shù)文件）sistbl.（測(cè)站約束控制文件）site.default（測(cè)站控制文件）process.default（指定文件）26sestbl.27sestbl.28總結(jié)：

參數(shù)控制文件為數(shù)據(jù)處理方案的核心控制文件，其中主要的參數(shù)包括處理模式、觀測(cè)量的選取、批處理迭代方案、對(duì)流層誤差模型、天頂延遲個(gè)數(shù)、光壓模型、數(shù)據(jù)篩選、ARC參考系統(tǒng)、慣性框架、L坐標(biāo)更新。sestbl.29sistbl.30siste.default刪除“#”，輸入：chan$exptftprnxxstinfoglrepuglrepsglts31process.default32（4）數(shù)據(jù)連接：

cd命令進(jìn)入$doy（001）文件夾，執(zhí)行命令：

%ln-s../rinex/*o./ #鏈接觀測(cè)數(shù)據(jù)

%ln-s../igs/*sp3./ #鏈接精密星歷

%ln-s../brdc/*n./ #鏈接廣播星歷

%ln-s../met/*m./ #鏈接氣象數(shù)據(jù)

%sh_links.tablesJ2000$year

#鏈接表文件數(shù)據(jù)處理流程33（5）解算：

%makexp

#數(shù)據(jù)預(yù)處理

%sh_sp3fit-f*sp3-oigsf-d$year$doy#使用精密星歷生成軌道參數(shù)

%makejbrdc*njbrdc$y.$doy #生成衛(wèi)星鐘文件

%makex$expt.makex.batch #創(chuàng)建觀測(cè)數(shù)據(jù)

%fixdrvd$expt$y.$doy

#啟動(dòng)數(shù)據(jù)驅(qū)動(dòng)程序

%cshb*bat

#執(zhí)行數(shù)據(jù)解算數(shù)據(jù)處理流程34（1）創(chuàng)建文件：

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

sestbl.

解算策略配置文件

sittbl.

測(cè)站約束控制文件

process.defaults

參數(shù)指定

sites.defaults

站點(diǎn)控制文件

*.apr 固定坐標(biāo)（可以用L文件）

測(cè)站信息文件

自動(dòng)批處理35（3）開始解算：

%sh_gamit-expt$expt-d$yeardays

sh_gamit自動(dòng)批處理程序會(huì)自動(dòng)創(chuàng)建其他所需文件夾，并進(jìn)行數(shù)據(jù)解算。解算結(jié)果將保存到年積日文件夾。執(zhí)行這個(gè)命令后，會(huì)自動(dòng)創(chuàng)建一些文件夾，創(chuàng)建的這些文件夾被使用過(guò)后有的留了下來(lái)，有的被刪掉。

自動(dòng)批處理36在輸出文件中，在q$expta.$doy($expt為工程名。$doy為解算的年積日，如下面的算例用的是2012年年積日為001天，工程名為gycy因此文件名為qgycya.001。)和o$expta.$doy，o*是q*的縮略版，在autcln.post.sum中可以看控制網(wǎng)點(diǎn)的相對(duì)質(zhì)量。主要指標(biāo)：（1）查看q*文件中的，后驗(yàn)nrms值，最好在0.25附近，一般不可以超過(guò)0.3。（2）基線分量誤差的大小，參數(shù)的改正量不應(yīng)大于約束量的2倍。

（3）*.sum總結(jié)文件中的rms值，看控制網(wǎng)點(diǎn)的質(zhì)量。（4）基線向量的重復(fù)性：

結(jié)果分析37主要內(nèi)容

●Linux操作系統(tǒng)

●shell簡(jiǎn)單介紹

●GAMIT基線處理●

GLOBK網(wǎng)平差

●算例分析

●總結(jié)38文件準(zhǔn)備參數(shù)配置平差過(guò)程精度評(píng)定GLOBK網(wǎng)平差步驟：39

拷貝GAMIT/GLOBK安裝目錄下的tables/文件夾下的pmu.usno（地球自轉(zhuǎn)參數(shù)）+itrf08_comb.apr+itrf08_comb.eq（地震信息）文件到解算目錄下的tables/文件夾。

拷貝GAMIT/GLOBK安裝目錄下的tables/文件夾下的globk_comb.cmd+glorg_comb.cmd文件到解算目錄下的gsoln/文件夾。

對(duì)globk_comb.cmd文件和glorg_comb.cmd文件進(jìn)行編輯以適合你的解算任務(wù)。文件準(zhǔn)備：40

globk_comb.cmd文件對(duì)平差過(guò)程進(jìn)行控制，編輯globk_comb.cmd文件以滿足任務(wù)要求，主要有以下項(xiàng)：參數(shù)配置：參數(shù)意義說(shuō)明eq_file地震信息文件目錄（目錄全部指到文件存放處）apr_file先驗(yàn)坐標(biāo)文件目錄in_pmu地球自轉(zhuǎn)參數(shù)文件目錄apr_neu測(cè)站坐標(biāo)值約束41

glorg_comb.cmd文件對(duì)固定點(diǎn)進(jìn)行控制，以將平差結(jié)果納入到一定的坐標(biāo)框架。編輯glorg_comb.cmd文件以滿足任務(wù)要求，主要有以下項(xiàng)參數(shù)配置：參數(shù)意義說(shuō)明apr_file先驗(yàn)坐標(biāo)文件目錄stab_siteclear清除之前的固定站標(biāo)記(有這幾個(gè)字就行)stab_site固定站點(diǎn)（在“上”個(gè)后添加，在將那些固定的站寫在后面即可。“這兩行的意思是先將itrf08z8888中的固定坐標(biāo)全部刪除，然后再?gòu)南乱恍兄兄付ü潭ㄗ鴺?biāo)”）42分步處理：（1）將GAMIT的H文件（或其他解算數(shù)據(jù)）轉(zhuǎn)為GLOBK認(rèn)可的二進(jìn)制的H文件：%cdhfile/%htoglb../glbf/../tables/svnav.dath????a*將文本形式的H文件轉(zhuǎn)化為二進(jìn)制形式（擴(kuò)展名為glr或glx），并存放到glbf/文件夾。（2）創(chuàng)建H文件列表：%cd

../gsoln/%ls../glbf/h*glx>$expt_glx.gdl將所有需解算的H文件的路徑和文件名列入gdl文件。（3）創(chuàng)建待平差測(cè)站的先驗(yàn)坐標(biāo)apr文件：文件應(yīng)包括測(cè)站的三維坐標(biāo)+站速度+歷元也可使用現(xiàn)有坐標(biāo)框架的值，如itrf08_comb.apr（以該坐標(biāo)框架為例）平差過(guò)程：43（4）在gsoln/文件夾下執(zhí)行：%glred6glred.optglred.log$expt_glx.gdlglobk_comb.cmd生成測(cè)站坐標(biāo)的時(shí)間序列，并進(jìn)行測(cè)站重復(fù)性檢驗(yàn)。（5）執(zhí)行命令，進(jìn)行自由網(wǎng)平差：%globk6globk.prtglobk.log$expt_glx.gdlglobk_comb.cmd（6）執(zhí)行命令：%glorg6glorg.log$expt_glx.gdlglorg_comb.cmd進(jìn)行約束平差，將測(cè)站納入到所選的坐標(biāo)框架。平差過(guò)程：44自動(dòng)批處理：（1）配置文件：主要是配置sites.defaults

globk_comb.cmd（如上所述）

glorg_comb.cmd（如上所述）（2）執(zhí)行批處理命令：%

sh_glred-d$year$doy-expt$expt-optHGE平差過(guò)程：45結(jié)