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

1、數(shù)據(jù)處理圖文流程冉啟順目錄一、 前言二、數(shù)據(jù)準(zhǔn)備1.前期準(zhǔn)備首先想好要處理什么時(shí)段、什么地方的數(shù)據(jù)本文檔以處理bjfs,shao,lhaz,urum四個(gè)站點(diǎn)2013年2月1日的數(shù)據(jù)為例本文檔使用軟件版本GAMIT處理環(huán)境：win10上的虛擬機(jī)2.工程目錄設(shè)置在任意位置建立工程文件夾（全文以test為例），并在其目錄下建立如下的文件夾：3.數(shù)據(jù)下載在終端中打開test，查詢2013年2月1日的年積日，和GPS周等信息代碼：doy 2013 02 011 在test的目錄下再建一個(gè)文件夾以032命名的文件夾，此時(shí)，test目錄下有：2 根據(jù)步驟3的信息，到相關(guān)網(wǎng)站下載對(duì)應(yīng)的廣播星歷，精密星歷，觀測(cè)

2、值文件分別放在文件夾brdc,igs,rinex里面。當(dāng)然還涉及到下載的數(shù)據(jù)解壓，.d文件到.o文件的轉(zhuǎn)換等操作，其中廣播星歷可以不改名字（即可以不用改成）4.更新tables表，更新準(zhǔn)則，下面的表格為準(zhǔn)：還有個(gè)需要每次處理都要更新表文件名表含義更新頻次內(nèi)容更新地址 pole.或極移表每次更新pub/gamit/tables/ ut1.或地球自轉(zhuǎn)參數(shù)表每次更新月亮表使用當(dāng)年的表太陽表使用當(dāng)年的表章動(dòng)表使用當(dāng)年的表跳秒表每次更新gdetic.dat,大地坐標(biāo)參數(shù)文件不常更新antmod.dat,天線相位中心改正正表有新天線時(shí)更新svnav.dat,衛(wèi)星列表有新衛(wèi)星時(shí)更新rcvant.dat,接收

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

4、經(jīng)有了2013年的這些年表。iii. 將上述更新的數(shù)據(jù)復(fù)制到gamit安裝目錄下的tables文件夾內(nèi)，并替換原有文件。我的gamit安裝路徑是/opt/gamit，即我安裝路徑下的tables 的路徑是：/tablesiv. 然后將更新后的tables文件夾，即/tables拷貝（或者設(shè)置連接，我習(xí)慣拷貝）到test目錄下。因?yàn)槲抑霸?tables里面放了很多數(shù)據(jù)，所以它顯得特別大，往往在包含otl.grid的情況下tables也就7800M吧。此后，test的結(jié)構(gòu)應(yīng)是：5.創(chuàng)立鏈接這一步相當(dāng)重要，而且常常容易出錯(cuò)。1） 將.o文件，.sp3文件，*.*n文件（觀測(cè)文件、精密星歷、廣播星歷

5、）鏈接到032的目錄下執(zhí)行代碼：2） 鏈接處理數(shù)據(jù)需要的各種年表等文件此時(shí)032中有如下的文件：檢查并重新設(shè)置斷掉的鏈接，如下圖所示：如：luntab.鏈接已斷（常說的死鏈接），在032目錄下，執(zhí)行代碼：ln -s ./tables/luntab.2013.J2000 ./luntab.（最后的點(diǎn)號(hào)一定要加上）系統(tǒng)報(bào)錯(cuò)，即“無法創(chuàng)建符號(hào)鏈接，luntab.早已存在”，這句話的意思是，032中早已經(jīng)有了一個(gè)luntab.文件（但是我們卻看到它是個(gè)死鏈接），刪除該鏈接，重新執(zhí)行上面的代碼。新鏈接生成，并打開查看該鏈接所指向的文件是不是2013年的luntab表重點(diǎn)檢查6-ii的11個(gè)文件鏈接，以及

6、soltab.,nutabl.，用類似的方法建立鏈接以下是我手動(dòng)加入的部分鏈接：以下是網(wǎng)上的資料圖：其中還要注意設(shè)置：，atl.grid的鏈接，我在處理數(shù)據(jù)時(shí)使用的如下鏈接：在查看原始的解壓文件中的tables文件夾的鏈接atl.grid,如下圖所示：但是我在查看發(fā)現(xiàn)它原本有一個(gè)，上傳是2010年8月28，大小也是，故我斷定，服務(wù)器上的atl.grid就是tables里面的。于是我在服務(wù)器上下載了一個(gè)的拷貝到test/tables中，然后在032里面制作了一個(gè)atl.grid的鏈接。（當(dāng)然也可以不用下載atl.grid,直接制作指向的鏈接）綜上所述：n 更新tables常更新以下11個(gè)文件（各

7、種年表除外）：但常常 guess_rcvant.dat 用不到,可以不更新n 制作鏈接特別注意以下幾個(gè)鏈接：n 在更新tables,建立完連接后執(zhí)行l(wèi)s al 檢查鏈接情況，紅色的是死鏈接：發(fā)現(xiàn)的鏈接已斷，于是重新設(shè)置即可。至于剩下的紅色死鏈接，我全是忽略了，因?yàn)槲乙膊粫?huì)設(shè)置。（好像那個(gè)是在正式處理過程中產(chǎn)生的，可以忽略）總的來說設(shè)置鏈接是及其容易出錯(cuò)，要小心又小心。6.控制文件配置準(zhǔn)備測(cè)站列表在test/tables下把需要處理的測(cè)站加入該表，包括IGS測(cè)站和本地測(cè)站，并確定工程名。然后將在032下設(shè)置一個(gè)鏈接 準(zhǔn)備測(cè)站近似坐標(biāo)文件lfile.制作完lfile.文件后不要忘了將lfile.拷

8、貝到test/tables中,同時(shí)在032中設(shè)置對(duì)應(yīng)鏈接lfile.文件是測(cè)站近似點(diǎn)坐標(biāo)（越精確越好以前的版本里面存放的經(jīng)緯度信息，而及其之后的版本直接存放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)，最好不要超過10cm。l 如果該測(cè)站在基線解算中不作為起算點(diǎn)，則偽距單點(diǎn)定位的精度即可。對(duì)于IGS站，獲取其在特定框架下的精確坐標(biāo)方法：首先要確保test/ri

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

10、標(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站，近似坐標(biāo)通過偽距單點(diǎn)定位獲得。最后按照特定的格式均復(fù)制到lfile.文件中然

11、后將lfile.拷貝到test/tables中,同時(shí)在032中設(shè)置對(duì)應(yīng)鏈接好像不需要在032中設(shè)置鏈接，但是怕出錯(cuò)，我還是設(shè)置了，希望沒有弄巧成拙 文件的更新最后要在032中檢查鏈接是否指向了更新后的文件在gamit全英文的官方教程（版本10.6,24頁）里面有這樣一句話：可見station.Info的鏈接可以指向下圖中選中的文件，也就是說倘若你的test/tables中沒有文件，可以自己新建并指向他們。并且這個(gè)MIT,SOPAC的文件好像可以在上更新，只是登不上去。練手的話湊合著能用。在test/tables目錄下：先建立一個(gè)sitelist的

12、文件，并錄入測(cè)站名字（每行首字母空格），不分先后順序然后終端打開tables并執(zhí)行代碼：顯示 files have 3 sites，這是錯(cuò)誤的，一般來說總共有多少測(cè)站參與計(jì)算就會(huì)有多少個(gè)sites，于是回去檢查sitelist, ，發(fā)現(xiàn)拉薩站寫錯(cuò)了，應(yīng)該是lhaz（所以在操作過程中一定要小心又小心，切記不能將單詞寫錯(cuò)，同時(shí)也要不害怕出現(xiàn)錯(cuò)誤，只要根據(jù)其提示，順藤摸瓜，基本上能修改錯(cuò)誤的地方，實(shí)在不行重來）還發(fā)現(xiàn)，Recevier name * not found in rcvant.dat,網(wǎng)上百度了一下說要修改，說得又很模糊，所以不會(huì)。又去問了一下師兄，他說可以忽略，

13、所以我也就忽略了（誰知道的話告訴我一聲）。執(zhí)行上述代碼結(jié)束后最后一行顯示：生成了一個(gè)，此時(shí)需要將其改名為，同時(shí)刪除原有的的鏈接。再執(zhí)行一次更新代碼：在生成的文件中會(huì)有下圖中紅框的字符串，表示更新成功：最后在032中檢查鏈接的指向，如果沒有的話，自己手動(dòng)用代碼建立一個(gè)。7 修改配置、sittbl.與sestbl.文件7.1配置它里面包含目錄信息，改動(dòng)較小。需要在032中設(shè)置鏈接配置sittbl.坐標(biāo)約束文件 對(duì)高精度的已知坐標(biāo)強(qiáng)約束，待求點(diǎn)坐標(biāo)松約束非起算點(diǎn)約束量可以更大,如100m。并且如果IGS站列表中沒有你所處理數(shù)據(jù)中的IGS站，需要自己添加，格式一致。配置sestbl.改正模型文件配置準(zhǔn)

14、則如下面的截圖：下面是我的部分設(shè)置：我沒有met.list,met.grid,也沒有和map.grid,所以我設(shè)置的是N（不使用）我的和是死鏈接，我又不知道如何制作他們（好像在麻省理工的FTP上可以下載，但是這個(gè)FTP好久都登不上去了），于是全選N，而我有atl.grid和 ，于是選Y修改完了，別忘了保存，同時(shí)為了保險(xiǎn)起見到032里面再檢查一次幾個(gè)配置文件是否修改、鏈接好了至此，數(shù)據(jù)準(zhǔn)備階段完成。下面我們開始計(jì)算。三、分布處理1 .Makexp在032目錄下執(zhí)行代碼：makexp2精密星歷處理如下圖執(zhí)行代碼：3文件檢查執(zhí)行下圖的代碼：對(duì)于代碼中的gigsf3.032,其中g(shù)表示文件類型，中間i

15、gsf表示軌道，第一步輸入的：，3表示年的最后一位，032是年積日。4廣播星歷處理其中沒有改名字作為輸入文件，jauto3.032是輸出文件當(dāng)然也可以在將改名字后，執(zhí)行 makej brdc0320.13n jbrdc5文件檢查6生成中間文件7生成bat文件8基線解算最后一步還是出錯(cuò)了，但是我百度了一下“RC not executed STOP FATAL Error: Stop from report_stat”在這篇文章的最后一頁找到解決方法，也不知道對(duì)不對(duì)，反正最后是成功解算了。在032中找到上圖的兩個(gè)文件并刪除，然后重新執(zhí)行基線解算。中間有很多省略結(jié)算結(jié)果精度還算過得去但是我往往以為準(zhǔn)

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

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