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

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

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

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

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

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

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

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

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

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

10、現(xiàn)atml.grid的鏈接已斷，于是重新設置即可。至于剩下的紅色死鏈接，我全是忽略了，因為我也不會設置。（好像那個是在正式處理過程中產(chǎn)生的，可以忽略）總的來說設置鏈接是及其容易出錯，要小心又小心。6.控制文件配置6.1準備測站列表sites.defaults在test/tables下把需要處理的測站加入該表，包括IGS測站和本地測站，并確定工程名。然后將site.defaults在032下設置一個鏈接 6.2準備測站近似坐標文件lfile.制作完lfile.文件后不要忘了將lfile.拷貝到test/tables中,同時在032中設置對應鏈接lfile.文件是測站近似點

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

12、）：在我的電腦上試了幾次，老是最后一步出錯，不推薦，但也有教程上是這個方法：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 偽距單點定位計算非IGS站的初始坐標實例（假設shao是非IGS站）：執(zhí)行完之后會生成兩個不同格式的測站坐標文件lfile.shao 和shao.apr，然后執(zhí)行l(wèi) 雙差定位計算測

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

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

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

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

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

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

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

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

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

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