gamit學(xué)習(xí)總結(jié)

1、. 打開LINUX系統(tǒng)，登陸。（選擇安裝語言包，關(guān)鍵更新)重啟。語言包的設(shè)置在system settings-language support- install languages 選擇chinese（simplified）點(diǎn)擊apply changes 等待系統(tǒng)自動安裝，重啟，重復(fù)步驟進(jìn)入support 選擇系統(tǒng)和菜單語言 漢語 ok（文件夾盡量保持英文狀態(tài)，以便在終端中輸入）；2.GAMIT安裝包中沒有包含海洋潮汐模型otl.grid 和mapvmf1函數(shù)模型需要自己下載連接。otl.grid在/archive/garner/gamit/tables

2、/ 中下載 （約696M）vmf1在/pub/GRIDS/ 中下載（約220M）如沒有下載otl.grid模型可以通過設(shè)置禁用該模型。模型的禁用：打開/opt/gamit/tables/sestbl. 將use otl.grid的之改為N，表示禁用海洋潮汐模型，同時(shí)將 Tides spplied 改為23(原始為31）。 獨(dú)立下載模型的連接：進(jìn)入/opt/gamit/tables 運(yùn)行終端運(yùn)行：ls -al 進(jìn)行連接 查看otl與map的連接情況rm otl.grid 移除已有的連接ln -s /opt/otlmap/otl_FES2004.grid o

3、tl.grid 進(jìn)行連接（/opt/otlmap/otl_FES2004.grid為模型的路徑）rm map.gridln -s /opt/otlmap/vmf1grd.2012_063 map.grid (map.grid 目前有vmf提供模型）ls -al 連接并查看是否已連接成功。連接結(jié)束。（連接完成后安裝更新時(shí)會自動連接）其它otl模型：􀂆 otl_FES99.grid 182.7 Mb􀂆 otl_CSR4.grid 45.6 Mb􀂆 otl_GOT00.grid 45.6 Mb􀂆 otl_NAO99b.grid

4、87.1 MbAtmosphere模型也許年更新：atml files 185 Mb/year3.sites.defaults 準(zhǔn)備介紹sites.defaults文件是站處理控制文件。格式如下：all_sites expt xstinfobrus_gps expt ftprnx ftpraw localrx xstinfo xsite glrepu glreps glts其中 ftprnx ftpraw localrx xstinfo xsite為gamit選項(xiàng)， glrepu glreps glts為GLOBK選項(xiàng)，格式具體含義如下：expt 為4個字符的項(xiàng)目名稱（工程名稱）ftprnx

5、表示該站點(diǎn)的rinex數(shù)據(jù)從FTP站點(diǎn)下載ftpraw 表示該站點(diǎn)的raw數(shù)據(jù)從FTP站點(diǎn)下載localrx 表示站點(diǎn)的rinex數(shù)據(jù)位本地?cái)?shù)據(jù)xstinfo 表示該站點(diǎn)不自動更新表中的信息xsite 該站點(diǎn)的所有天或指定日期的數(shù)據(jù)不處理glrepu 表示站點(diǎn)用于計(jì)算重復(fù)解（默認(rèn)全部）glreps 用于定義站點(diǎn)的在GLRED處理中的參考框架（默認(rèn)為IGS站點(diǎn)列表）glts 在GLRED處理中以時(shí)間序列繪制站點(diǎn)圖（默認(rèn)全部）一般通過自己編制的sites.defaults文件可以設(shè)置所有站點(diǎn)不自動更新表中的信息。示例如下：all_sites test

6、 xstinfo （test為工程名）LHAS_GPS test localrx glrepu glrepsSHAO_GPS test localrx 文件準(zhǔn)備文件可以采用sh_upd_stnfo 命令生成也可以人工編輯得到。對于部分不規(guī)范的rinex文件建議采用人編輯的方法，并參照gamit的rcvant（天線與接收）文件進(jìn)行編輯，確保后期處理的順利。在更新前應(yīng)將sites.defaults文件中的 all_sites expt xstinfo注釋掉，否則容易報(bào)錯，無法更新。采用sh_upd_stnfo 命令生成station.

7、info文件方法：將tables中的文件拷貝至rinex中，在rinex文件夾下運(yùn)行sh_upd_stnfo -l sd 提取IGS站點(diǎn)的信息，得到.new然后見.new改名為，在運(yùn)行 sh_upd_stnfo -files *12o -expt test 即可。采用手工編輯的辦法：將tables中的文件拷貝至rinex中，在rinex文件夾下運(yùn)行sh_upd_stnfo -l sd 提取IGS站點(diǎn)的信息。會自動生成.new 文件。打開station.

8、info.new刪除無關(guān)條目，增加其他站點(diǎn)的信息。此時(shí)應(yīng)打開rinex文件進(jìn)行核對，確保無誤后將其保存為。應(yīng)注意將天線高參考設(shè)為DHARP（天線參考中心一般為天線底部）格式： Fixed format file with station information *TRCK SITE Station Name Ant Ht Ant N Ant E 7ODM 7ODM Seven Oaks Dam 0.0083 0.0000 0.0000 Rcvr AntCod HtCod Vers ASHZ12 ATDMRB DHPAB 9.10 Year Doy SN Start St

9、op 2000 110 0 00 00 00 24 00 00 ( = line break for this display) This file must be sorted by TRCK, then epochTRCK and SITE 4-character code match the 4-char code in RINEX file name TRCK = SITE for static analysis 在靜態(tài)分析中 Change for kinematic analysis 對于動態(tài)分析要改變 File must be sorted by this field, then

10、epoch station_name is a long (12-char) descriptionAntHt Raw antenna height as measured 原始測量高度 Critical measurement! 至關(guān)重要 HtCod Defines geometry of AntHt measurement 定義幾何類型 DHPAB is RINEX standard is vertical height to antenna reference point (ARP) （DHPAB指相對于站點(diǎn)的垂直高度） N and E offsets 基準(zhǔn)偏移 e.g., uncali

11、brated tribrach rarely used occupied point w.r.t. benchmarkRCVR 6-character codes for the receiver type not critical except for some old receivers or "buggy" firmware versions Vers version number of firmware in the receiver 接收機(jī)固件版本號 Should be correct to handle SNR coding properlyAntCod and

12、 HtCod are critical 關(guān)鍵 Sets phase center model 涉及相位中心模型 conversion of antenna height measurement 天線高度的解譯 SLBGP - Slant bottom ground plane 傾斜的底部大地平面 year and doy specify year and day of year when thevalues apply times set part of day (if it changes during UTC day) This representation will be augment

13、ed in futureversion of GAMIT5.sittbl.編輯打開sittble文件刪除無關(guān)站點(diǎn)，增加工程站點(diǎn)設(shè)置各站點(diǎn)的約束屬性。對于IGS框架點(diǎn)或你需要約束的點(diǎn)設(shè)置為0.050 0.050 0.050，其他點(diǎn)100. 100. 100.一般來說約束點(diǎn)為IGS站點(diǎn)且不少于3個。示例：SITE FIX -COORD.CONSTR.-BJFS BJFS_GPS NNN 100. 100. 100.GUAO GUAO_GPS NNN 100. 100. 100.LHAZ LHAZ_GPS NNN 100. 100. 100.SHAO SHAO_GPS NNN 100. 100. 1

14、00.URUM URUM_GPS NNN 0.050 0.050 0.05WUHN WUHN_GPS NNN 0.050 0.050 0.05KUNM KUNM_GPS NNN 0.050 0.050 0.05下載sittbl.refined 文件,里面有最近的各站的可靠性的一個評測。要想點(diǎn)位置精度高則設(shè)置IGS站的約束小一些,即認(rèn)為IGS站已經(jīng)很準(zhǔn),強(qiáng)約束,若目的是求基線后面平差則可以將IGS站的約束設(shè)置大些。6.process.defaults 文件可修改部分的解釋process.default文件中需要修改的有采樣間隔，最大歷元數(shù) ，開始時(shí)間 等，示例如下：set sint = '

15、;30' 采樣間隔依據(jù)需要和rinex的情況可修改為15 30 60set nepc = '2880' 最大歷元數(shù) 依據(jù)觀測時(shí)長計(jì)算調(diào)整。set stime = '0 0' 開始時(shí)間（小時(shí) 分）依據(jù)rinex文件修改，一般要求時(shí)間與rinex時(shí)間差不大于0.1h 否則報(bào)錯。set use_rxc = "N" 設(shè)置為N 不依據(jù)rinex更新lfile文件。set aprf = itrf08.apr 采用08框架數(shù)據(jù)（默認(rèn)itrf2000） 如不采用GLOBK平差可不設(shè)置。set minxf = '300' Xfile文件

16、大小的最小值，小于該值的將被拋棄。依據(jù)歷元數(shù)設(shè)置。歷元數(shù)少時(shí)應(yīng)減小該值。set mailto = '*' 設(shè)置為有效郵箱以便接受gamit結(jié)算結(jié)果。7.sestbl.文件依據(jù)實(shí)際需要進(jìn)行修改。示例如下：hoice of Experiment = BASELINE （計(jì)算高精度基線）BASELINE/RELAX./ORBIT默認(rèn)為BASELINE，但一般修改為relax.(BASLINE為固定軌道）（如只需要基線文件可以選擇BASELINE）Choice of Observable = LC_AUTCLN 觀測值類型選擇LC_AUTCLN (default), LC_HELP (

17、codeless L2), L1_ONLY (L1 soln from dual freq),L2_ONLY (L2 soln from dual freq),L1,L2_INDEPENDENT (L1 + L2 from dual freq),L1&L2 (same as L1,L2_INDEPENDENT but with ion constraint);L1_RECEIVER (must add 'L1only' in autcln.cmd) 對于15km以下的短基線應(yīng)選擇L1,L2_INDEPENDENT 或L1_ONLYType of analysis =0-

18、iter 迭代多少次計(jì)算結(jié)果作為初始值再次計(jì)算一般設(shè)置1-iterZenith delay =N 天頂延遲 對于短基線選NInterval zen = 2 ;2 hrs = 13 knots/day (default is 1 ZD per day) 天頂延遲修改為1h，默認(rèn)2hUpdate T/L files = L_ONLY ;T_AND_L (default), T_ONLY, L_ONLY, NONE 人工編輯時(shí)可設(shè)為NONE# DMap = GMF ; GMF(default)/NMFH/VMF1 如果已下載VMF1 這可設(shè)為VMF1 -重新連接。# WMap = GMF ; GMF

19、(default)/NMFW/VMF1 如果已下載VMF1 這可設(shè)為VMF1Use map.grid = Y ;VMF1 grid file with mapping functions and ZHDInertial frame = J2000 ; B1950/J2000 (default = J2000)8.sh_gamit 批處理解釋> & ! sh_gamit_test.log 生成過程詳解文件sh_gamit_test.log當(dāng)采用快速星歷處理時(shí)需要設(shè)置 -orbit IGSU/IGSR#要處理最近的數(shù)據(jù)需要將gamit/gamit/fixdrv文件夾下，打開armak

20、e.f文件，將2010改為2100.然后再gamit文件夾下運(yùn)行./install_software.即重新編譯安裝。（未驗(yàn)證）-mailto lyzlrz聯(lián)網(wǎng)后發(fā)送-Basic example sh_gamit -d 1998 240 -orbit SIOF -expt emed eops usno process GPS data from day 240 of 1998 (Aug 28) use any data in the local rinex directory ftp data for sites specified by the user optionally scan mu

21、ltiple directories to find data use Final orbit estimate from Scripps Institution of Oceanography (SIO) 軌道估計(jì)數(shù)據(jù)來源 use EOP data from US Naval Observatory (USNO) IERS Bulletin A EOP數(shù)據(jù)來源 Name the run “emed" 工程名identify filesselect dataList options sh_gamit -help Select days of data to process -d yr

22、 days where yr and days are the data to be processed e.g., -d1997 153 156 179 -s yr d1 d2where d1 is start day and d2 is stop day to beprocessed e.g. 1997 153 178 -r daywhere day is number of days before current date.sh_gamit options -expt 4 char name of experiment/solution being run Default expt -o

23、rbitType of orbit to use. Default SIOF. Optional second argument gives the day-of-year for the orbit if not thesame as processing day -orbt4 char ID for orbit files. Default 4 char extracted from the orbit selected -eops Name of EOP series to use. (bull_b, bull_a, etc.) Default usno-sessinfoSession

24、- sampling interval, number of epochs, start time (default 30 2880 0 0) -rx_doy_plus nSearch RINEX file names this many days into the future for current-day data 尋找往前推遲n天的rinex數(shù)據(jù) -rx_doy_minus n Search RINEX file names this many days into the past for current-day data 尋找往后n天的rinex數(shù)據(jù) -xsite listList

25、of sites to be excluded from the processing. 解算排除的站 -netext <char>add network suffix <char> to day directories. i.e. 035r 日期文件夾添加網(wǎng)的后綴 -yrext add year prefix to day directories. Ie 1997_053 日期文件夾添加年前綴 -pres Plot gamit residuals as skyplots <Y/N/E>. YES, NO, ELEV. Default N 繪制天空殘差圖 -

26、nogifs Do NOT create skyplot GIF graphics files 不產(chǎn)生天空GIF圖片-copt file typesCompress remaining files for each dayserve when d Default x k autcln.out 壓縮指定類型的每天文件 -aopt file typesCopy the specified files for each day to the directory specified by archivepth 將指定每天文件復(fù)制到archivepth文件夾。 -dopt file typesDelet

27、e specified files after processing each day. Default c 刪除指定類型文件 -mailto Name of user and machine to send results summary. (e.g., )-snrresCalculate theoretical phase residuals from SNR.Requires raw data or RINEX with valid SNR values. -noaccDo not exclude satellites base on accura

28、cy codes in the .sp3 file. Default exclude 不排除.sp3中基于精度碼方式排除衛(wèi)星。 -xver char 1-character x-file version (6th character of x-file). Default is single-digit year -noftpUse this option if no ftp download connections are available. 不聯(lián)網(wǎng) -remakex Y/C/N Y - Remake x-files without checking. C - Remake x-file

29、if orbit and are incompatible and all RINEX files are available Default C N - Don't remake x-files under any circumstances若tables中消失，sh_gamit將會使用rinex的頭文件信息。guess_rcvant.dat 來將rinex的頭文件信息翻譯成gamit 可閱讀形式。Antenna height errors 在GLOBK可降低，但類型識別錯誤無法消除。-9.打開日期目錄查看結(jié)算結(jié)果判斷解算質(zhì)量主要有兩

30、個方面：一是數(shù)據(jù)利用率（即是否有足夠多的可用于進(jìn)行合理估計(jì)的數(shù)據(jù)參與計(jì)算）；二是數(shù)據(jù)的噪聲水平是否與選用的模型匹配。（這一項(xiàng)主要有兩個指標(biāo)，基線分量誤差的大小和 nrms（標(biāo)準(zhǔn)化的均方根誤差)，對于基線分量誤差的判斷主要依據(jù)預(yù)先設(shè)定的站點(diǎn)坐標(biāo)和軌道參數(shù)的先驗(yàn)約束，如果大于先驗(yàn)可以通過查看Q文件或者 autcln.sum 文件，是否有大量被autcln丟棄的數(shù)據(jù)。對于nrms一般應(yīng)在0.25左右。如果大于0.5就意味著周跳沒有消除，或出現(xiàn)了模型錯誤（固定站點(diǎn)坐標(biāo)錯誤或者未模型化的衛(wèi)星）等。解算結(jié)果文件autcln.sum.post中的rms殘差為顯示控制網(wǎng)站點(diǎn)的相對質(zhì)量。打開 sh_gamit_

31、doy.summary，otesta.doy, qtesta.doy 文件查看解算總結(jié)和基線解算情況，解算成功與否判斷：Postfit nrms: *E+00 是否在0.25左右， 最差不能大于0.5，基線分量的誤差大小是否小于預(yù)設(shè)值。至此基線解算全部結(jié)束?；€解算結(jié)果：sh_gamit_doy.summary基線解算結(jié)果： qexpta.doy,oexpta.doy,hexpta.yrdoy10.若軟件限制了最大測站數(shù)據(jù)，最大衛(wèi)星數(shù)目，最大對流層延遲參數(shù)，最大解算歷元數(shù)目。更改于Makefile.config文件。然后重新編。 11.tables表介紹。12.格式問題o文件名需要小寫，并且內(nèi)

32、容格式符合UINX文件名改小寫：sh_casefold -dir d -files *格式轉(zhuǎn)UNIX: sh_dos2unix *13.gamit模塊以及結(jié)構(gòu)14. Data Types L1 and L2 phase used for parameter estimation P1 and P2 range measurements used to help clean phase data used to estimate receiver clock offset with tolerance of 1microsecond (= 300 m) LC phase proportional

33、 to geometric range sensitive insensitive to ionospheric delay linear combination of L1 and L2 LG phase measures ionospheric delay only insensitive to geomery Widelane WL linear combination of L1, L2, P1 and P2 Melbourne-Wubbena Should be constant (+/- noise)15. Aim of GAMITBasic aim of GAMIT is to

34、generate positions/orbits forone day of data. Standard product of GAMIT is aminimally constrained solution16. gg/tables/rcvant.datFile contains all receivers and antennas known to GAMIT Translates user-specified 6-character code to IGS SINEX standard. The file in templates directory provides mapping

35、 from RINEX header information to 6-char codes. Example linesrcv ASHTECHZ-XII ASHZ12rcv Z-XII3 ASHZ12swv 5F 1.00 Explanation given at top of file used a represent space Sometimes need to be careful that pattern correctly gets correct type. Should not be a problem for IGS rinex files default lines al

36、lowed, but be careful17. 模型說明：1. otl 潮汐改正2. vmfl GMF 投影函數(shù)3. atml大氣荷載模型，對高程影響較大，可消除周跳波動，可靠性需要進(jìn)一步證實(shí)4. atl大氣抄襲荷載模型和met氣象模型下載地址：1)海洋潮汐。例如otl_FES2004.grid放在軟件talbels目錄下。鏈接到otl.grid。2)大氣負(fù)荷。例如atmldisp_cm.2006，每年更新一次。連接到atmldisp_grid.20063)vmfl投影函數(shù)柵格，例如vmflgrd2006，連接到map.grid.2006。每年更新一次

37、。18. 解算所需文件說明類型文件名說明：數(shù)據(jù)文件bjfs0900.06o bjfs 為站點(diǎn)代碼， 090 為年積日，0為時(shí)段號， 06 為2006 年，o 為觀測文件auto0900 06n 廣播星歷文件，前4 個字節(jié)改為站點(diǎn)代碼則為站點(diǎn)廣播星歷文件igs13685 sp3 igs 為事后精密星歷，1368 為GPS周，5為星期五igr13685 sp3 igr 代表快速精密星歷igu13685 sp3 igu 代表預(yù)報(bào)精密星歷配置文件lgigg6 090 l 代表L 文件， gigg 代表工程名，6代表2006 年， 090 為年積日station Oct測站信息文件，記錄各測站接收機(jī)和天

38、線類型session Info 時(shí)段信息文件，可以在數(shù)據(jù)處理過程中自動創(chuàng)建sestbl解算參數(shù)控制文件，設(shè)置與解算有關(guān)的參數(shù)，達(dá)到不同的解算目的sittbl設(shè)置測站坐標(biāo)約束、大氣模型等參數(shù)表文件ut1 國際時(shí)間系統(tǒng)表pole 極移表leap sec 跳秒表nutabl 章動表luntal 月亮表soltab 太陽表gdetic dat 地球形狀參數(shù)表rcvant dat 接收機(jī)及其天線型號對照表svnav dat 星號對照表svs_excludeDat壞衛(wèi)星信息表hi dat 儀器高參數(shù)表dcb dat 差分碼偏差表autcln cmd 周跳自動探測和修改命令表antmod dat 天線相位中

39、心改正表ot1 list /grid全球格網(wǎng)或各測站的海洋潮汐參數(shù)表19 .crx2rnx 壓縮的d文件轉(zhuǎn)換成o文件Usage: crx2rnx file - -f -s -h stdin and stdout are used if input file name is not given. - : output to stdout -f : force overwrite of output file -s : skip strange epochs (default:stop with error) This option may be used for salvaging usable

40、 data when middle of the Compact RINEX file is missing. The data after the missing part, are, however, useless until the compression operation of all data are initialized at some epoch. Combination with use of -e option of RNX2CRX (ver.4.0 or after) may be effective. Caution : It is assumed that no

41、change in # of data types happens in the lost part of the data. -h : display help message exit code = 0 (success) = 1 (error) = 2 (warning)20某個人的參數(shù)設(shè)置：表2 解算控制參數(shù)的選擇控制選項(xiàng)控制參數(shù)處理模式RELAX迭代方案1-ITER觀測量LC_AUTCLN天頂延遲參數(shù)25 Number Zen = 25 ; number of zenith-delay parameters (default 1)梯度參數(shù)3-Number gradients海潮模型o

42、tl.list光壓模型BERNE慣性框架J2000 0對流層誤差模型SAAS截止高度角 /( °) 10 - Elevation Cutoff測站約束( BJFS) /m 0 05 /0 05 /0 1測站約束( WUHN) /m 0 025 /0 01 /0 1測站約束( KUNM) /m 9 999 /9 999 /9 999 9其他默認(rèn)某篇論文寫到;不同截止高度角對基線精度的影響表4 不同截止高度角對基線(BJFS_KUNM)精度的影響精度指標(biāo)高度角/( °)5 10 15 20 25nrms ± 0 241 04 ± 0 241 04 ±

43、; 0 231 09 ± 0 222 52 ± 0 213 98X /m ± 0 041 9 ± 0 041 9± 0 045 0± 0 048 5 ± 0 051 2Y /m ± 0 040 1 ± 0 040 1 ± 0 043 4 ± 0 049 1 ± 0 056 3Z /m ± 0 044 5± 0 044 5 ± 0 047 8 ± 0 051 1 ± 0 055 1L /L 4.90 × 10 E-9

44、4 90 × 10 E-95 49 × 10 E-9 6 53 × 10 E-9 7 67 × 10 E-9不過高度角太低時(shí)路徑延遲對水汽變化非常敏感，因此在一般的解算情況下選擇10°的截止高度角比較適宜。不同天頂延遲參數(shù)對基線精度的影響表5 不同天頂延遲參數(shù)對基線(BJFS_KUNM)精度的影響精度指標(biāo)天頂延遲參數(shù)7 13 19 25Nrms ± 0 239 47 ± 0 242 77 ± 0 241 66 ± 0 241 04X /m ± 0 042 5 ± 0 041 9

45、77; 0 041 8 ± 0 041 9Y /m ± 0 040 8 ± 0 040 2 ± 0 039 8 ± 0 040 1Z /m ± 0 045 9 ± 0 044 7 ± 0 044 2 ± 0 044 5L /L 4 95 × 10 E-9 5 05 × 10 E-9 4 90 × 10 E-9 4 90 × 10 E-9一般情況下，天頂延遲參數(shù)取為13，即每2 h 產(chǎn)生一個天頂延遲估值是比較合適的。天頂延遲參數(shù)變化對長基線的精度影響不是非常明顯。選取

46、不同觀測量對基線精度的影響表6 選取不同觀測量對基線(BJFS_KUNM)精度的影響精度指標(biāo)觀測量LC_AUTCLN L1L2 LC_HELP LC_ONLYnrms ± 0 241 04 ± 0 239 28 ± 0 239 27 ± 0 239 27X /m ± 0 041 9 ± 0 047 3 ± 0 047 9 ± 0 047 9Y /m ± 0 040 1 ± 0 044 1 ± 0 044 6 ± 0 044 6Z /m ± 0 044 5 ±

47、; 0 053 8 ± 0 055 3 ± 0 055 3L /L 4 90 × 10 E-9 6 53 × 10 E-9 6 73 × 10 E-9 6 73 × 10 E-9LC_AUTCLN( 觀測值是無電離層的線性組合) 不同測站坐標(biāo)約束對基線(BJFS_KUNM)精度的影響精度指標(biāo)測站約束/m0.01/0.01/0.02 0.05/0.05/0.06 0.1/0.1/0.1 0.5/0.5/0.5 1/1/1nrms ± 0 239 38 ± 0 241 18 ± 0 241 07 ± 0 240 64 ± 0 240 59X /m ± 0 013 3 ± 0 037 2 ± 0 047 4 ± 0 119 4 ± 0 223 9Y /m ± 0 019 9 ± 0 034 8 ± 0 040 1 ± 0 081 1 ± 0 146 2Z /m ± 0 01