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

1、GAMIT10.5數(shù)據(jù)處理圖文流程冉啟順2016.1.20目錄TOC o 1-3 h z u HYPERLINK l _Toc441083676 一、前言言 PAGEREF _Toc441083676 h 2 HYPERLINK l _Toc441083677 二、數(shù)據(jù)準(zhǔn)準(zhǔn)備 PAGEREF _Toc441083677 h 3 HYPERLINK l _Toc441083678 1.前期準(zhǔn)準(zhǔn)備 PAGEREF _Toc441083678 h 3 HYPERLINK l _Toc441083679 2.工程目目錄設(shè)置置 PAGEREF _Toc441083679 h 3 HYPERLINK l

2、_Toc441083680 3.數(shù)據(jù)下下載 PAGEREF _Toc441083680 h 3 HYPERLINK l _Toc441083681 4.更新ttablles表表， PAGEREF _Toc441083681 h 4 HYPERLINK l _Toc441083682 5.創(chuàng)立鏈鏈接 PAGEREF _Toc441083682 h 7 HYPERLINK l _Toc441083683 6.控制文文件配置置 PAGEREF _Toc441083683 h 13 HYPERLINK l _Toc441083684 6.1準(zhǔn)備備測站列列表siitess.deefauultss PAG

3、EREF _Toc441083684 h 13 HYPERLINK l _Toc441083685 6.2準(zhǔn)備備測站近近似坐標(biāo)標(biāo)文件llfille. PAGEREF _Toc441083685 h 13 HYPERLINK l _Toc441083686 6.3 sstattionn.innfo文文件的更更新 PAGEREF _Toc441083686 h 16 HYPERLINK l _Toc441083687 7 修改配配置prroceess.deffaullts、sitttbll.與sesstbll.文件件 PAGEREF _Toc441083687 h 19 HYPERLINK l _

4、Toc441083688 7.1配置置proocesss.ddefaaultts PAGEREF _Toc441083688 h 19 HYPERLINK l _Toc441083689 7.2配置置sitttbll.坐標(biāo)標(biāo)約束文文件 PAGEREF _Toc441083689 h 19 HYPERLINK l _Toc441083690 7.3配置置sesstbll.改正正模型文文件 PAGEREF _Toc441083690 h 20 HYPERLINK l _Toc441083691 三、分布處處理 PAGEREF _Toc441083691 h 22 HYPERLINK l _Toc4

5、41083692 1 .Maakexxp PAGEREF _Toc441083692 h 22 HYPERLINK l _Toc441083693 2精密星星歷處理理 PAGEREF _Toc441083693 h 22 HYPERLINK l _Toc441083694 3文件檢檢查 PAGEREF _Toc441083694 h 24 HYPERLINK l _Toc441083695 4廣播星星歷處理理 PAGEREF _Toc441083695 h 24 HYPERLINK l _Toc441083696 5文件檢檢查 PAGEREF _Toc441083696 h 25 HYPERL

6、INK l _Toc441083697 6生成中中間文件件 PAGEREF _Toc441083697 h 25 HYPERLINK l _Toc441083698 7生成bbat文文件 PAGEREF _Toc441083698 h 26 HYPERLINK l _Toc441083699 8基線解解算 PAGEREF _Toc441083699 h 27 HYPERLINK l _Toc441083700 9基線解解算精度度評定標(biāo)標(biāo)準(zhǔn) PAGEREF _Toc441083700 h 30 HYPERLINK l _Toc441083701 四、cossagpps平差差 PAGEREF _T

7、oc441083701 h 31 HYPERLINK l _Toc441083702 五、建議與與說明 PAGEREF _Toc441083702 h 34 HYPERLINK l _Toc441083703 六、相關(guān)資資料鏈接接 PAGEREF _Toc441083703 h 35前言二、數(shù)據(jù)準(zhǔn)準(zhǔn)備1.前期準(zhǔn)準(zhǔn)備首先想好要要處理什什么時段段、什么地地方的數(shù)數(shù)據(jù)本文檔以處處理bjjfs,shaao,llhazz,urrum四四個站點點20113年22月1日日的數(shù)據(jù)據(jù)為例本文檔使用用軟件版版本GAAMITT10.5處理環(huán)境：winn10上上的虛擬擬機ubbunttu155.1002.工程目目錄設(shè)

8、置置在任意位置置建立工工程文件件夾（全全文以ttestt為例），并在在其目錄錄下建立立如下的的文件夾夾：3.數(shù)據(jù)下下載在終端中打打開teest，查查詢20013年年2月11日的年年積日，和和GPSS周等信信息代碼：dooy 220133 022 011在testt的目錄錄下再建建一個文文件夾以以0322命名的的文件夾夾，此時時，teest目目錄下有有：根據(jù)步驟33的信息息，到相相關(guān)網(wǎng)站站下載對對應(yīng)的廣廣播星歷歷，精密密星歷，觀觀測值文文件分別別放在文文件夾bbrdcc,iggs,rrineex里面面。當(dāng)然還涉及及到下載載的數(shù)據(jù)據(jù)解壓，.d文件件到.oo文件的的轉(zhuǎn)換等等操作，其其中廣播播星歷可可

9、以不改改名字（即即可以不不用改成成brddc03320.13nn）4.更新ttablles表表，更新準(zhǔn)則，下面的的表格為為準(zhǔn)：還有個dccb.ddat需需要每次次處理都都要更新新表文件名表含義更新頻次內(nèi)容更新地地址 poole.或polle.uusnoo極移表每次更新 HYPERLINK pub/ggamiit/ttablles/ utt1.或或ut11.ussno地球自轉(zhuǎn)參參數(shù)表每次更新 lunntabb.20012.J20000月亮表使用當(dāng)年的的表 solltabb.20012.J20000太陽表使用當(dāng)年的的表 nuutabbl.220122章動表使用當(dāng)年的的表 leaap.ssec跳秒表

10、每次更新gdetiic.ddat,大地坐標(biāo)參參數(shù)文件件不常更新antmood.ddat,天線相位中中心改正正正表有新天線時時更新svnavv.daat,衛(wèi)星列表有新衛(wèi)星時時更新rcvannt.ddat,接收機天線線對照表表有新天線時時更新svs_eexclludee HYPERLINK .innfo壞衛(wèi)星信息息文件每次更新vmf1映射函數(shù)模模型 使用當(dāng)年年的表 HYPERLINK otl.ggrd海潮格網(wǎng)數(shù)數(shù)據(jù)使用當(dāng)年的的表常見文件說說明更新的準(zhǔn)則則是：“更新頻頻次”中的紅紅字部分分每次處處理都更更新一下下，同時時結(jié)合本本次處理理與上次次處理的的年份是是否相同同來進(jìn)行行年表的的選擇。原原本有些

11、些文件是是每周更更新一次次，但是是作為新新手為了了避免出出錯還是是每次更更新吧有些時候，有些文文件找不不到，或者 HYPERLINK 里面沒有，則可以以忽略。如svss_exxcluude HYPERLINK .inffo，gdeeticc.daat常常常找不到到我在更新時時，這四個個文件在在服務(wù)器器上沒找找到下圖是我更更新的文文件。原本還還應(yīng)該有有20113年的的各種年年表，但但是我之之前處理理過一次次20113年的的數(shù)據(jù)，并并且我的的gammit安安裝路徑徑/oppt/ggamiit/GGAMIIT100.5/tabbless已經(jīng)有有了20013年年的這些些年表。將上述更新新的數(shù)據(jù)據(jù)復(fù)制到

12、到gammit安安裝目錄錄下的ttablles文文件夾內(nèi)內(nèi)，并替換換原有文文件。我的gammit安安裝路徑徑是/oopt/gammit，即我安安裝路徑徑下的ttablles 的路徑徑是：/oppt/ggamiit/GGAMIIT100.5/tabbless然后將更新新后的ttablles文文件夾，即/oppt/ggamiit/GGAMIIT100.5/tabbless拷貝（或或者設(shè)置置連接，我我習(xí)慣拷拷貝）到到tesst目錄錄下。因為我之前前在/oopt/gammit/GAMMIT110.55/taablees里面面放了很很多數(shù)據(jù)據(jù)，所以以它顯得得特別大大1.88G，往往往在包包含ottl.g

13、gridd的情況況下taablees也就就78800MM吧。此后，teest的的結(jié)構(gòu)應(yīng)應(yīng)是：5.創(chuàng)立鏈鏈接這一步相當(dāng)當(dāng)重要，而且常常常容易易出錯。將.o文件件，.ssp3文文件，*.*nn文件（觀觀測文件件、精密密星歷、廣廣播星歷歷）鏈接接到0332的目目錄下執(zhí)行代碼：鏈接處理數(shù)數(shù)據(jù)需要要的各種種年表等等文件此時0332中有有如下的的文件：檢查并重新新設(shè)置斷斷掉的鏈鏈接，如如下圖所所示：如：lunntabb.鏈接接已斷（常常說的死死鏈接），在032目錄下，執(zhí)行代碼：ln -s ./tabbless/luuntaab.220133.J220000 ./lunntabb.（最后的點點號一定定要加上

14、上）系統(tǒng)報錯，即“無法創(chuàng)創(chuàng)建符號號鏈接，lunntabb.早已已存在”，這句話話的意思思是，0032中中早已經(jīng)經(jīng)有了一一個luuntaab.文文件（但但是我們們卻看到到它是個個死鏈接接），刪刪除該鏈鏈接，重重新執(zhí)行行上面的的代碼。新鏈接生成成，并打開開查看該該鏈接所所指向的的文件是是不是220133年的llunttab表表重點檢查66-ii的11個文件件鏈接，以以及sooltaab.,nuttabll.，用類似似的方法法建立鏈鏈接以下是我手手動加入入的部分分鏈接：以下是網(wǎng)上上的資料料圖：其中還要注注意設(shè)置置：mapp.grrid，ootl.griid，aatl.griid的鏈鏈接，我在處處理數(shù)

15、據(jù)據(jù)時使用用的如下下鏈接：在查看原始始GAMMIT110.55的解壓壓文件中中的taablees文件件夾的鏈鏈接attl.ggridd,如下下圖所示示：但是我在查查看 HYPERLINK 發(fā)現(xiàn)現(xiàn)它原本本有一個個atll.grrid，上傳是是20110年88月288，大小小也是33.1MM，故我我斷定，服服務(wù)器上上的attl.ggridd就是ttablles里里面的AANU11008826_griid.aatl。于是我在服服務(wù)器上上下載了了一個33.1MM的attl.ggridd拷貝到到tesst/ttablles中中，然后在在0322里面制制作了一一個attl.ggridd的鏈接接。（當(dāng)然然也

16、可以以不用下下載attl.ggridd,直接接制作指指向ANNU10008226_ggridd.attl的鏈鏈接）綜上所述：更新tabbless常更新新以下111個文文件（各各種年表表除外）：但常常 gguesss_rrcvaant.datt 用不不到,可可以不更更新制作鏈接特特別注意意以下幾幾個鏈接接：在更新taablees,建建立完連連接后執(zhí)執(zhí)行l(wèi)ss aal 檢檢查鏈接接情況，紅色的的是死鏈鏈接：發(fā)現(xiàn)atmml.ggridd的鏈接接已斷，于于是重新新設(shè)置即即可。至至于剩下下的紅色色死鏈接接，我全全是忽略略了，因因為我也也不會設(shè)設(shè)置。（好像像那個ssesssionn.innfo是是在正式式

17、處理過過程中產(chǎn)產(chǎn)生的，可可以忽略略）總的來說設(shè)設(shè)置鏈接接是及其其容易出出錯，要要小心又又小心。6.控制文文件配置置6.1準(zhǔn)備備測站列列表siitess.deefauultss在testt/tabbless下把需要要處理的的測站加加入該表表，包括括IGSS測站和和本地測測站，并并確定工工程名。然然后將ssitee.deefauultss在0332下設(shè)設(shè)置一個個鏈接 6.2準(zhǔn)備備測站近近似坐標(biāo)標(biāo)文件llfille.制作完lffilee.文件件后不要要忘了將將lfiile.拷貝到到tesst/ttablles中中,同時時在0332中設(shè)設(shè)置對應(yīng)應(yīng)鏈接lfilee.文件件是測站站近似點點坐標(biāo)（越越精確越

18、越好,110.44以前的的版本里里面存放放的經(jīng)緯緯度信息息，而110.44及其之之后的版版本直接接存放xx,y,z信息息）GAMITT10.5 的的lfiile 和appr 文文件是同同樣的格式，所以以我們可可以認(rèn)為為aprr 文件件就是llfille。該該文件中中保存的的是各個站的初始始坐標(biāo)，至至于該文文件中的的測站坐坐標(biāo)要達(dá)達(dá)到什么么樣的精精度，可可參考如下兩兩條：如果是該站站作為基基線解算算起算點點,則該測測站的坐坐標(biāo)做保保持在5cm 以以內(nèi)，最最好不要要超過110cmm。如果該測站站在基線線解算中中不作為為起算點點，則偽偽距單點點定位的的精度即可。對于IGSS站，獲獲取其在在特定框框架

19、下的的精確坐坐標(biāo)方法法：首先要確保保tesst/rrineex下有有itrrf088.appr方法1：我用的方法法（可用用于非IIGS站站）：在在我的電電腦上試試了幾次次，老是是最后一一步出錯錯，不推推薦，但但也有教教程上是是這個方方法：在rineex文件件夾下先先生成ll文件，然后生生成l.aprr，最后后生成llfille.依次執(zhí)行代代碼：grep POSSITIION *.113o llrx2appr ll 20013 0322gapr_to_l ll lffilee. * 220133 0332方法2：方法3：偽距單點定定位計算算非IGGS站的的初始坐坐標(biāo)實例例（假設(shè)設(shè)shaao是非非

20、IGSS站）：執(zhí)行完之后后會生成成兩個不不同格式式的測站站坐標(biāo)文文件l 和shao.aprr，然后后執(zhí)行雙差定位計計算測站的初初始坐標(biāo)標(biāo)實例：在此例中假假設(shè)shhao站站是未知知坐標(biāo)的的測站，而bjjfs是是已知坐坐標(biāo)的測測站（或或起算站站）綜上所述：在生成liife.文件時時，IGGS站的的坐標(biāo)可可以直接接從框架架獲得，或或者利用用雙差定定位求得得。而非IGSS站則可可以根據(jù)據(jù)要求精精度的高高低選擇擇：直接提取頭頭信息偽距單點定定位雙差定位法法方法1常常在偽距距單點定定位和雙雙差定位位法中選選擇，然然后將所所須的坐坐標(biāo)信息息拷貝到到lfiile.中即可可，并將lffilee.拷貝貝到tees

21、t/tabbless中,同同時在0032中中設(shè)置對應(yīng)應(yīng)鏈接現(xiàn)在，我假假設(shè)bjjfs，lhaaz兩站站點是精精度較高高的IGGS站，坐坐標(biāo)直接接從框架架中獲得得，而sshaoo,urrum站站是非IIGS站站，近似坐坐標(biāo)通過過偽距單單點定位位獲得。最后按按照特定定的格式式均復(fù)制制到lffilee.文件件中然后將lffilee.拷貝貝到teest/tabbless中,同同時在0032中中設(shè)置對應(yīng)應(yīng)鏈接好像不需要要在0332中設(shè)設(shè)置鏈接接，但是是怕出錯錯，我還還是設(shè)置置了，希希望沒有有弄巧成成拙6.3 sstattionn.innfo文文件的更更新最后要在0032中中檢查sstattionn.inn

22、fo鏈鏈接是否否指向了了更新后后的sttatiion.inffo文件件在gamiit全英英文的官官方教程程（版本本10.6,224頁）里面有這樣一句話：可見staatioon.Inffo的鏈鏈接可以以指向下下圖中選選中的文文件，也也就是說說倘若你你的teest/tabbless中沒有有staatioon.iinfoo文件，可以自自己新建建并指向向他們。并且這這個MIIT,SSOPAAC的文文件好像像可以在在 HYPERLINK 上更新新，只是登登不上去去。練手的的話湊合合著能用用。在testt/taablees目錄錄下：先建立一個個sitteliist的的文件，并錄入入測站名名字（每每行首字字

23、母空格格），不不分先后后順序然后終端打打開taablees并執(zhí)執(zhí)行代碼碼：顯示staatioon.iinfoo fiiless haave 3 ssitees，這是錯錯誤的，一般來來說總共共有多少少測站參參與計算算就會有有多少個個sittes，于是回去檢查sitelist, ，發(fā)現(xiàn)拉薩站寫錯了，應(yīng)該是lhaz（所以在操作過程中一定要小心又小心，切記不能將單詞寫錯，同時也要不害怕出現(xiàn)錯誤，只要根據(jù)其提示，順藤摸瓜，基本上能修改錯誤的地方，實在不行重來）還發(fā)現(xiàn)，RReceevieer nnamee * noot ffounnd iin rrcvaant.datt,網(wǎng)上上百度了了一下說說要修改改，說

24、得得又很模模糊，所所以不會會。又去去問了一一下師兄兄，他說說可以忽忽略，所所以我也也就忽略略了（誰誰知道的的話告訴訴我一聲聲）。執(zhí)行上述代代碼結(jié)束束后最后后一行顯顯示：生成了了一個sstattionn.innfo.neww，此時需需要將其其改名為為staatioon.iinfoo，同時刪刪除原有有的sttatiion.inffo的鏈鏈接。再執(zhí)行一次次更新代代碼：在生成的sstattionn.innfo文文件中會會有下圖圖中紅框框的字符符串，表表示更新新成功：最后在0332中檢檢查sttatiion.inffo鏈接接的指向向，如果果沒有的的話，自自己手動動用代碼碼建立一一個。7 修改配配置prr

25、oceess.deffaullts、sitttbll.與sesstbll.文件件7.1配置置proocesss.ddefaaultts它里面包含含目錄信信息，改改動較小小。需要要在0332中設(shè)設(shè)置鏈接接7.2配置置sitttbll.坐標(biāo)標(biāo)約束文文件 對高精精度的已已知坐標(biāo)標(biāo)強約束束，待求求點坐標(biāo)標(biāo)松約束束非起算點約約束量可可以更大大,如1100mm。并且且如果IIGS站站列表中中沒有你你所處理理數(shù)據(jù)中中的IGGS站，需需要自己己添加，格格式一致致。7.3配置置sesstbll.改正正模型文文件配置準(zhǔn)則如如下面的的截圖：下面是我的的部分設(shè)設(shè)置：我沒有meet.llistt,meet.ggridd

26、,也沒沒有maap.llistt和maap.ggridd,所以以我設(shè)置置的是NN（不使使用）我的atmml.llistt和attl.llistt是死鏈鏈接，我我又不知知道如何何制作他他們（好好像在麻麻省理工工的FTTP上可可以下載載，但是這這個FTTP好久久都登不不上去了了），于于是全選選N，而而我有attl.ggridd和attml.griid ，于于是選YY修改完了，別別忘了保保存，同同時為了了保險起起見到0032里里面再檢檢查一次次幾個配配置文件件是否修修改、鏈接好了了至此，數(shù)據(jù)據(jù)準(zhǔn)備階階段完成成。下面我我們開始始計算。三、分布處處理1 .Maakexxp在032目目錄下執(zhí)執(zhí)行代碼碼：m

27、aakexxp2精密星星歷處理理如下圖執(zhí)行行代碼：3文件檢檢查執(zhí)行下圖的的代碼：對于代碼中中的giggsf33.0332,其其中g(shù)表表示文件件類型，中中間iggsf表表示軌道道，第一一步輸入入的：，33表示年年的最后后一位，0032是是年積日日。4廣播星星歷處理理其中autto03320.13nn沒有改改名字作作為輸入入文件，jjautto3.0332是輸輸出文件件當(dāng)然也可以以在將aautoo03220.113n改改名字后后，執(zhí)行行 maakejj brrdc003200.133n jjbrddc3.03225文件檢檢查6生成中中間文件件7生成bbat文文件8基線解解算最后一步還還是出錯錯了，

28、但是我我百度了了一下“RC nott exxecuutedd STOOP FFATAAL EErroor: Stoop ffromm reeporrt_sstatt”在這篇文文章的最最后一頁頁找到解解決方法法，也不知知道對不不對，反正最最后是成成功解算算了。 HYPERLINK 在032中中找到上上圖的兩兩個文件件并刪除除，然后后重新執(zhí)執(zhí)行基線線解算。中間有很多多省略結(jié)算結(jié)果精精度還算算過得去去但是我往往往以ottestta.0032為為準(zhǔn)qexpta.doy:解算記錄qexpta.doy:解算記錄oexpta.doy為qexpta.doy的簡略版hexpta.yrdoy：協(xié)方差矩陣、參數(shù)平差

29、值9基線解解算精度度評定標(biāo)標(biāo)準(zhǔn)（網(wǎng)上各類類教程上上的標(biāo)準(zhǔn)準(zhǔn)大同小小異）標(biāo)準(zhǔn)1：一般NRMS越小越好，一般小于0.25視為解算效果較好，本次解算的NRMS為一般NRMS越小越好，一般小于0.25視為解算效果較好，本次解算的NRMS為0.26，符合要求基線的X、YY方向誤誤差范圍圍應(yīng)控制制在1ccm以內(nèi)內(nèi)，Z方方向的誤誤差控制制在1.5cmm以內(nèi)，本本次解算算所有基基線均滿滿足要求求本次解算xx,y,zz方向的的誤差均均有超過過1cmm,結(jié)果果不是很很理想，還過得得去最后還有一一個基線線重復(fù)性性檢驗：但是我的檢檢驗失敗敗了，是不是是要進(jìn)行行兩天以以上的數(shù)數(shù)據(jù)處理理才能進(jìn)進(jìn)行重復(fù)復(fù)性試驗驗？至此可以理理解為解解算完畢畢，然后進(jìn)進(jìn)行基線線平差，也可以以使用GLLBOKK平差，但但是我選選擇了用用COSSAGPPS進(jìn)行行的平差差。如果想直接接用GLLBOKK進(jìn)行平平差，參參考安向向東的處處理流程程。四、cossagpps平差差在032目目錄下打打開ottestta.0032,并在第第二次出出現(xiàn)基線線信息的的地方加加入下面面紅框中中的一行行字符串串：“cossagpps ffor gammit o-ffilee”最后保保存并復(fù)復(fù)制到wwinddowss桌面下下。Cosa平平差打開cossagpps,新新建工程程，設(shè)置如如下，其它