GAMIT105數(shù)據(jù)處理圖文流程

1、GAMIT10.5數(shù)據(jù)處理圖文流程冉啟順2016.1.20目錄一、前言2二、數(shù)據(jù)準(zhǔn)備31.前期準(zhǔn)備32.工程目錄設(shè)置33.數(shù)據(jù)下載34.更新tables表，45.創(chuàng)立鏈接76.控制文件配置136.1準(zhǔn)備測(cè)站列表sites.defaults136.2準(zhǔn)備測(cè)站近似坐標(biāo)文件lfile.136.3 文件的更新167 修改配置process.defaults、sittbl.與sestbl.文件197.1配置process.defaults197.2配置sittbl.坐標(biāo)約束文件197.3配置sestbl.改正模型文件20三、分布處理221 .Makexp222精密星歷處理223文

2、件檢查244廣播星歷處理245文件檢查256生成中間文件257生成bat文件268基線解算279基線解算精度評(píng)定標(biāo)準(zhǔn)30四、cosagps平差31五、建議與說(shuō)明34六、相關(guān)資料鏈接35一、 前言二、數(shù)據(jù)準(zhǔn)備1.前期準(zhǔn)備首先想好要處理什么時(shí)段、什么地方的數(shù)據(jù)本文檔以處理bjfs,shao,lhaz,urum四個(gè)站點(diǎn)2013年2月1日的數(shù)據(jù)為例本文檔使用軟件版本GAMIT10.5處理環(huán)境：win10上的虛擬機(jī)ubuntu15.102.工程目錄設(shè)置在任意位置建立工程文件夾（全文以test為例），并在其目錄下建立如下的文件夾：3.數(shù)據(jù)下載在終端中打開(kāi)test，查詢2013年2月1日的年積日，和GPS周

3、等信息代碼：doy 2013 02 011 在test的目錄下再建一個(gè)文件夾以032命名的文件夾，此時(shí)，test目錄下有：2 根據(jù)步驟3的信息，到相關(guān)網(wǎng)站下載對(duì)應(yīng)的廣播星歷，精密星歷，觀測(cè)值文件分別放在文件夾brdc,igs,rinex里面。當(dāng)然還涉及到下載的數(shù)據(jù)解壓，.d文件到.o文件的轉(zhuǎn)換等操作，其中廣播星歷可以不改名字（即可以不用改成brdc0320.13n）4.更新tables表，更新準(zhǔn)則，下面的表格為準(zhǔn)：還有個(gè)dcb.dat需要每次處理都要更新表文件名表含義更新頻次內(nèi)容更新地址 pole.或pole.usno極移表每次更新/pub/gamit/

4、tables/ ut1.或ut1.usno地球自轉(zhuǎn)參數(shù)表每次更新 luntab.2012.J2000月亮表使用當(dāng)年的表 soltab.2012.J2000太陽(yáng)表使用當(dāng)年的表 nutabl.2012章動(dòng)表使用當(dāng)年的表 leap.sec跳秒表每次更新gdetic.dat,大地坐標(biāo)參數(shù)文件不常更新antmod.dat,天線相位中心改正正表有新天線時(shí)更新svnav.dat,衛(wèi)星列表有新衛(wèi)星時(shí)更新rcvant.dat,接收機(jī)天線對(duì)照表有新天線時(shí)更新svs_壞衛(wèi)星信息文件每次更新vmf1映射函數(shù)模型 使用當(dāng)年的表/pub/GRIDS/ otl

5、.grd海潮格網(wǎng)數(shù)據(jù)使用當(dāng)年的表i. 常見(jiàn)文件說(shuō)明ii. 更新的準(zhǔn)則是：“更新頻次”中的紅字部分每次處理都更新一下，同時(shí)結(jié)合本次處理與上次處理的年份是否相同來(lái)進(jìn)行年表的選擇。原本有些文件是每周更新一次，但是作為新手為了避免出錯(cuò)還是每次更新吧有些時(shí)候，有些文件找不到，或者/pub/gamit/tables/里面沒(méi)有，則可以忽略。如svs_，gdetic.dat常常找不到我在更新時(shí)，這四個(gè)文件在服務(wù)器上沒(méi)找到下圖是我更新的文件。原本還應(yīng)該有2013年的各種年表，但是我之前處理過(guò)一次2013年的數(shù)據(jù)，并且我的gamit安裝路徑/opt/

6、gamit/GAMIT10.5/tables已經(jīng)有了2013年的這些年表。iii. 將上述更新的數(shù)據(jù)復(fù)制到gamit安裝目錄下的tables文件夾內(nèi)，并替換原有文件。我的gamit安裝路徑是/opt/gamit，即我安裝路徑下的tables 的路徑是：/opt/gamit/GAMIT10.5/tablesiv. 然后將更新后的tables文件夾，即/opt/gamit/GAMIT10.5/tables拷貝（或者設(shè)置連接，我習(xí)慣拷貝）到test目錄下。因?yàn)槲抑霸?opt/gamit/GAMIT10.5/tables里面放了很多數(shù)據(jù)，所以它顯得特別大1.8G，往往在包含otl.grid的情況下t

7、ables也就7800M吧。此后，test的結(jié)構(gòu)應(yīng)是：5.創(chuàng)立鏈接這一步相當(dāng)重要，而且常常容易出錯(cuò)。1） 將.o文件，.sp3文件，*.*n文件（觀測(cè)文件、精密星歷、廣播星歷）鏈接到032的目錄下執(zhí)行代碼：2） 鏈接處理數(shù)據(jù)需要的各種年表等文件此時(shí)032中有如下的文件：檢查并重新設(shè)置斷掉的鏈接，如下圖所示：如：luntab.鏈接已斷（常說(shuō)的死鏈接），在032目錄下，執(zhí)行代碼：ln -s ./tables/luntab.2013.J2000 ./luntab.（最后的點(diǎn)號(hào)一定要加上）系統(tǒng)報(bào)錯(cuò)，即“無(wú)法創(chuàng)建符號(hào)鏈接，luntab.早已存在”，這句話的意思是，032中早已經(jīng)有了一個(gè)luntab.文件

8、（但是我們卻看到它是個(gè)死鏈接），刪除該鏈接，重新執(zhí)行上面的代碼。新鏈接生成，并打開(kāi)查看該鏈接所指向的文件是不是2013年的luntab表重點(diǎn)檢查6-ii的11個(gè)文件鏈接，以及soltab.,nutabl.，用類似的方法建立鏈接以下是我手動(dòng)加入的部分鏈接：以下是網(wǎng)上的資料圖：其中還要注意設(shè)置：map.grid，otl.grid，atl.grid的鏈接，我在處理數(shù)據(jù)時(shí)使用的如下鏈接：在查看原始GAMIT10.5的解壓文件中的tables文件夾的鏈接atl.grid,如下圖所示：但是我在查看/pub/GRIDS/發(fā)現(xiàn)它原本有一個(gè)atl.grid，上傳是2010

9、年8月28，大小也是3.1M，故我斷定，服務(wù)器上的atl.grid就是tables里面的ANU100826_grid.atl。于是我在服務(wù)器上下載了一個(gè)3.1M的atl.grid拷貝到test/tables中，然后在032里面制作了一個(gè)atl.grid的鏈接。（當(dāng)然也可以不用下載atl.grid,直接制作指向ANU100826_grid.atl的鏈接）綜上所述：n 更新tables常更新以下11個(gè)文件（各種年表除外）：但常常 guess_rcvant.dat 用不到,可以不更新n 制作鏈接特別注意以下幾個(gè)鏈接：n 在更新tables,建立完連接后執(zhí)行l(wèi)s al 檢查鏈接情況，紅色的是死鏈接：發(fā)

10、現(xiàn)atml.grid的鏈接已斷，于是重新設(shè)置即可。至于剩下的紅色死鏈接，我全是忽略了，因?yàn)槲乙膊粫?huì)設(shè)置。（好像那個(gè)是在正式處理過(guò)程中產(chǎn)生的，可以忽略）總的來(lái)說(shuō)設(shè)置鏈接是及其容易出錯(cuò)，要小心又小心。6.控制文件配置6.1準(zhǔn)備測(cè)站列表sites.defaults在test/tables下把需要處理的測(cè)站加入該表，包括IGS測(cè)站和本地測(cè)站，并確定工程名。然后將site.defaults在032下設(shè)置一個(gè)鏈接 6.2準(zhǔn)備測(cè)站近似坐標(biāo)文件lfile.制作完lfile.文件后不要忘了將lfile.拷貝到test/tables中,同時(shí)在032中設(shè)置對(duì)應(yīng)鏈接lfile.文件是測(cè)站近似點(diǎn)

11、坐標(biāo)（越精確越好,10.4以前的版本里面存放的經(jīng)緯度信息，而10.4及其之后的版本直接存放x,y,z信息）GAMIT10.5 的lfile 和apr 文件是同樣的格式，所以我們可以認(rèn)為apr 文件就是lfile。該文件中保存的是各個(gè)站的初始坐標(biāo)，至于該文件中的測(cè)站坐標(biāo)要達(dá)到什么樣的精度，可參考如下兩條：l 如果是該站作為基線解算起算點(diǎn),則該測(cè)站的坐標(biāo)做保持在5cm 以內(nèi)，最好不要超過(guò)10cm。l 如果該測(cè)站在基線解算中不作為起算點(diǎn)，則偽距單點(diǎn)定位的精度即可。對(duì)于IGS站，獲取其在特定框架下的精確坐標(biāo)方法：首先要確保test/rinex下有itrf08.apr方法1：我用的方法（可用于非IGS站

12、）：在我的電腦上試了幾次，老是最后一步出錯(cuò)，不推薦，但也有教程上是這個(gè)方法：l 在rinex文件夾下先生成lfile.rnx文件，然后生成lfile.rnx.apr，最后生成lfile.依次執(zhí)行代碼：grep POSITION *.13o > lfile.rnxrx2apr lfile.rnx 2013 032gapr_to_l lfile.rnx.apr lfile. * 2013 032方法2：方法3：l 偽距單點(diǎn)定位計(jì)算非IGS站的初始坐標(biāo)實(shí)例（假設(shè)shao是非IGS站）：執(zhí)行完之后會(huì)生成兩個(gè)不同格式的測(cè)站坐標(biāo)文件lfile.shao 和shao.apr，然后執(zhí)行l(wèi) 雙差定位計(jì)算測(cè)

13、站的初始坐標(biāo)實(shí)例：在此例中假設(shè)shao站是未知坐標(biāo)的測(cè)站，而bjfs是已知坐標(biāo)的測(cè)站（或起算站）綜上所述：在生成life.文件時(shí)，IGS站的坐標(biāo)可以直接從框架獲得，或者利用雙差定位求得。而非IGS站則可以根據(jù)要求精度的高低選擇：ü 直接提取頭信息ü 偽距單點(diǎn)定位ü 雙差定位法ü 方法1常常在偽距單點(diǎn)定位和雙差定位法中選擇，然后將所須的坐標(biāo)信息拷貝到lfile.中即可，并將lfile.拷貝到test/tables中,同時(shí)在032中設(shè)置對(duì)應(yīng)鏈接現(xiàn)在，我假設(shè)bjfs，lhaz兩站點(diǎn)是精度較高的IGS站，坐標(biāo)直接從框架中獲得，而shao,urum站是非IGS站，

14、近似坐標(biāo)通過(guò)偽距單點(diǎn)定位獲得。最后按照特定的格式均復(fù)制到lfile.文件中然后將lfile.拷貝到test/tables中,同時(shí)在032中設(shè)置對(duì)應(yīng)鏈接好像不需要在032中設(shè)置鏈接，但是怕出錯(cuò)，我還是設(shè)置了，希望沒(méi)有弄巧成拙6.3 文件的更新最后要在032中檢查鏈接是否指向了更新后的文件在gamit全英文的官方教程（版本10.6,24頁(yè)）里面有這樣一句話：可見(jiàn)station.Info的鏈接可以指向下圖中選中的文件，也就是說(shuō)倘若你的test/tables中沒(méi)有文件，可以自己新建并指向他們。并且這個(gè)MIT

15、,SOPAC的文件好像可以在上更新，只是登不上去。練手的話湊合著能用。在test/tables目錄下：先建立一個(gè)sitelist的文件，并錄入測(cè)站名字（每行首字母空格），不分先后順序然后終端打開(kāi)tables并執(zhí)行代碼：顯示 files have 3 sites，這是錯(cuò)誤的，一般來(lái)說(shuō)總共有多少測(cè)站參與計(jì)算就會(huì)有多少個(gè)sites，于是回去檢查sitelist, ，發(fā)現(xiàn)拉薩站寫(xiě)錯(cuò)了，應(yīng)該是lhaz（所以在操作過(guò)程中一定要小心又小心，切記不能將單詞寫(xiě)錯(cuò)，同時(shí)也要不害怕出現(xiàn)錯(cuò)誤，只要根據(jù)其提示，順藤摸瓜，基本上能修改錯(cuò)誤的地方，實(shí)在不行重

16、來(lái)）還發(fā)現(xiàn)，Recevier name * not found in rcvant.dat,網(wǎng)上百度了一下說(shuō)要修改，說(shuō)得又很模糊，所以不會(huì)。又去問(wèn)了一下師兄，他說(shuō)可以忽略，所以我也就忽略了（誰(shuí)知道的話告訴我一聲）。執(zhí)行上述代碼結(jié)束后最后一行顯示：生成了一個(gè).new，此時(shí)需要將其改名為，同時(shí)刪除原有的的鏈接。再執(zhí)行一次更新代碼：在生成的文件中會(huì)有下圖中紅框的字符串，表示更新成功：最后在032中檢查鏈接的指向，如果沒(méi)有的話，自己手動(dòng)用代碼建立一個(gè)。7 修改配置process.de

17、faults、sittbl.與sestbl.文件7.1配置process.defaults它里面包含目錄信息，改動(dòng)較小。需要在032中設(shè)置鏈接7.2配置sittbl.坐標(biāo)約束文件 對(duì)高精度的已知坐標(biāo)強(qiáng)約束，待求點(diǎn)坐標(biāo)松約束非起算點(diǎn)約束量可以更大,如100m。并且如果IGS站列表中沒(méi)有你所處理數(shù)據(jù)中的IGS站，需要自己添加，格式一致。7.3配置sestbl.改正模型文件配置準(zhǔn)則如下面的截圖：下面是我的部分設(shè)置：我沒(méi)有met.list,met.grid,也沒(méi)有map.list和map.grid,所以我設(shè)置的是N（不使用）我的atml.list和atl.list是死鏈接，我又不知道如何制作他們（好像

18、在麻省理工的FTP上可以下載，但是這個(gè)FTP好久都登不上去了），于是全選N，而我有atl.grid和atml.grid ，于是選Y修改完了，別忘了保存，同時(shí)為了保險(xiǎn)起見(jiàn)到032里面再檢查一次幾個(gè)配置文件是否修改、鏈接好了至此，數(shù)據(jù)準(zhǔn)備階段完成。下面我們開(kāi)始計(jì)算。三、分布處理1 .Makexp在032目錄下執(zhí)行代碼：makexp2精密星歷處理如下圖執(zhí)行代碼：3文件檢查執(zhí)行下圖的代碼：對(duì)于代碼中的gigsf3.032,其中g(shù)表示文件類型，中間igsf表示軌道，第一步輸入的：，3表示年的最后一位，032是年積日。4廣播星歷處理其中auto0320.13n沒(méi)有改名字作為輸入文件，jauto3.032是

19、輸出文件當(dāng)然也可以在將auto0320.13n改名字后，執(zhí)行 makej brdc0320.13n jbrdc3.0325文件檢查6生成中間文件7生成bat文件8基線解算最后一步還是出錯(cuò)了，但是我百度了一下“RC not executed STOP FATAL Error: Stop from report_stat”在這篇文章的最后一頁(yè)找到解決方法，也不知道對(duì)不對(duì)，反正最后是成功解算了。在032中找到上圖的兩個(gè)文件并刪除，然后重新執(zhí)行基線解算。中間有很多省略結(jié)算結(jié)果精度還算過(guò)得去但是我往往以otesta.032為準(zhǔn)qexpta.doy:解算記錄oexpta.doy為qexpta.doy的簡(jiǎn)略

20、版hexpta.yrdoy：協(xié)方差矩陣、參數(shù)平差值9基線解算精度評(píng)定標(biāo)準(zhǔn)（網(wǎng)上各類教程上的標(biāo)準(zhǔn)大同小異）標(biāo)準(zhǔn)1：一般NRMS越小越好，一般小于0.25視為解算效果較好，本次解算的NRMS為0.26，符合要求基線的X、Y方向誤差范圍應(yīng)控制在1cm以內(nèi)，Z方向的誤差控制在1.5cm以內(nèi)，本次解算所有基線均滿足要求本次解算x,y,z方向的誤差均有超過(guò)1cm,結(jié)果不是很理想，還過(guò)得去最后還有一個(gè)基線重復(fù)性檢驗(yàn)：但是我的檢驗(yàn)失敗了，是不是要進(jìn)行兩天以上的數(shù)據(jù)處理才能進(jìn)行重復(fù)性試驗(yàn)？至此可以理解為解算完畢，然后進(jìn)行基線平差，也可以使用GLBOK平差，但是我選擇了用COSAGPS進(jìn)行的平差。如果想直接用GL

21、BOK進(jìn)行平差，參考安向東的處理流程。四、cosagps平差1. 在032目錄下打開(kāi)otesta.032,并在第二次出現(xiàn)基線信息的地方加入下面紅框中的一行字符串：“cosagps for gamit o-file”最后保存并復(fù)制到windows桌面下。2. Cosa平差打開(kāi)cosagps,新建工程，設(shè)置如下，其它參數(shù)基本默認(rèn)3. 導(dǎo)入基線數(shù)據(jù)導(dǎo)入從ubuntu里面復(fù)制出來(lái)的otesta.032基線的誤差有點(diǎn)大（但為了走完流程，有精度強(qiáng)迫癥的就湊合著看吧_）4. 導(dǎo)入三維已知坐標(biāo)這里就是你解算時(shí)坐標(biāo)起算數(shù)據(jù)，但是要注意框架和坐標(biāo)系的統(tǒng)一。先在點(diǎn)名方框內(nèi)隨便輸入數(shù)據(jù)，然后確定，如下圖在工程目錄下找

22、到剛新建的關(guān)于起算點(diǎn)坐標(biāo)的文件，用記事本打開(kāi)，并將已知起算數(shù)據(jù)復(fù)制粘貼到know文件中，如下圖：5. 三維平差最后將解算出的shao,urum站的坐標(biāo)與原本查詢到的坐標(biāo)進(jìn)行對(duì)比即可。以下圖片里面的坐標(biāo)就是我在ubuntu里面采用：獲得的，對(duì)比后發(fā)現(xiàn)相差不是太大，但是還是不太理想（不管怎么樣將軟件走通了）五、建議與說(shuō)明六、相關(guān)資料鏈接天線查詢/ANTCAL/TEQC/software/teqc/廣播星歷下載/pub/rinex/YYYY/DOY/ (以auto開(kāi)頭的文件，下載后改為brdc開(kāi)頭)/pub/gps/data/daily/YYYY/brdc/  （以brdc開(kāi)頭）ftp:/igs.ensg.ign.fr/pub/igs/data/YYYY/DOY/  （以brdc開(kāi)頭）精密星歷下載/pub/ftp