GAMIT10.5數(shù)據(jù)處理圖文流程冉啟順2016.1.20目錄一、 前言 2二、數(shù)據(jù)準備 31.前期準備 32.工程目錄設(shè)置 33.數(shù)據(jù)下載 34.更新tables表， 45.創(chuàng)立鏈接 76.控制文件配置 136.1準備測站列表sites.defaults 136.2準備測站近似坐標(biāo)文件lfile. 136.3文件的更新 167修改配置process.defaults、sittbl.與sestbl.文件 197.1配置process.defaults 197.2配置sittbl.坐標(biāo)約束文件 197.3配置sestbl.改正模型文件 20三、分布處理 221.Makexp 222．精密星歷處理 223．文件檢查 244．廣播星歷處理 245．文件檢查 256．生成中間文件 257．生成bat文件 268．基線解算 279．基線解算精度評定標(biāo)準 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.工程目錄設(shè)置在任意位置建立工程文件夾（全文以test為例），并在其目錄下建立如下的文件夾：3.數(shù)據(jù)下載在終端中打開test，查詢2013年2月1日的年積日，和GPS周等信息代碼：doy20130201在test的目錄下再建一個文件夾以032命名的文件夾，此時，test目錄下有：根據(jù)步驟3的信息，到相關(guān)網(wǎng)站下載對應(yīng)的廣播星歷，精密星歷，觀測值文件分別放在文件夾brdc,igs,rinex里面。當(dāng)然還涉及到下載的數(shù)據(jù)解壓，.d文件到.o文件的轉(zhuǎn)換等操作，其中廣播星歷可以不改名字（即可以不用改成brdc0320.13n）4.更新tables表，更新準則，下面的表格為準：還有個dcb.dat需要每次處理都要更新表文件名表含義更新頻次內(nèi)容更新地址pole.或pole.usno極移表每次更新pub/gamit/tables/ut1.或ut1.usno地球自轉(zhuǎn)參數(shù)表每次更新luntab.2012.J2000月亮表使用當(dāng)年的表soltab.2012.J2000太陽表使用當(dāng)年的表nutabl.2012章動表使用當(dāng)年的表leap.sec跳秒表每次更新gdetic.dat,大地坐標(biāo)參數(shù)文件不常更新antmod.dat,天線相位中心改正正表有新天線時更新svnav.dat,衛(wèi)星列表有新衛(wèi)星時更新rcvant.dat,接收機天線對照表有新天線時更新svs_壞衛(wèi)星信息文件每次更新vmf1映射函數(shù)模型使用當(dāng)年的表otl.grd海潮格網(wǎng)數(shù)據(jù)使用當(dāng)年的表常見文件說明更新的準則是：“更新頻次”中的紅字部分每次處理都更新一下，同時結(jié)合本次處理與上次處理的年份是否相同來進行年表的選擇。原本有些文件是每周更新一次，但是作為新手為了避免出錯還是每次更新吧有些時候，有些文件找不到，或者里面沒有，則可以忽略。如svs_，gdetic.dat常常找不到我在更新時，這四個文件在服務(wù)器上沒找到下圖是我更新的文件。原本還應(yīng)該有2013年的各種年表，但是我之前處理過一次2013年的數(shù)據(jù)，并且我的gamit安裝路徑/opt/gamit/GAMIT10.5/tables已經(jīng)有了2013年的這些年表。將上述更新的數(shù)據(jù)復(fù)制到gamit安裝目錄下的tables文件夾內(nèi)，并替換原有文件。我的gamit安裝路徑是/opt/gamit，即我安裝路徑下的tables的路徑是：/opt/gamit/GAMIT10.5/tables然后將更新后的tables文件夾，即/opt/gamit/GAMIT10.5/tables拷貝（或者設(shè)置連接，我習(xí)慣拷貝）到test目錄下。因為我之前在/opt/gamit/GAMIT10.5/tables里面放了很多數(shù)據(jù)，所以它顯得特別大1.8G，往往在包含otl.grid的情況下tables也就7~800M吧。此后，test的結(jié)構(gòu)應(yīng)是：5.創(chuàng)立鏈接這一步相當(dāng)重要，而且常常容易出錯。將.o文件，.sp3文件，*.*n文件（觀測文件、精密星歷、廣播星歷）鏈接到032的目錄下執(zhí)行代碼：鏈接處理數(shù)據(jù)需要的各種年表等文件 此時032中有如下的文件：檢查并重新設(shè)置斷掉的鏈接，如下圖所示：如：luntab.鏈接已斷（常說的死鏈接），在032目錄下，執(zhí)行代碼：ln-s../tables/luntab.2013.J2000./luntab.（最后的點號一定要加上）系統(tǒng)報錯，即“無法創(chuàng)建符號鏈接，luntab.早已存在”，這句話的意思是，032中早已經(jīng)有了一個luntab.文件（但是我們卻看到它是個死鏈接），刪除該鏈接，重新執(zhí)行上面的代碼。新鏈接生成，并打開查看該鏈接所指向的文件是不是2013年的luntab表重點檢查6-ii的11個文件鏈接，以及soltab.,nutabl.，用類似的方法建立鏈接以下是我手動加入的部分鏈接：以下是網(wǎng)上的資料圖：其中還要注意設(shè)置：map.grid，otl.grid，atl.grid的鏈接，我在處理數(shù)據(jù)時使用的如下鏈接：在查看原始GAMIT10.5的解壓文件中的tables文件夾的鏈接atl.grid,如下圖所示：但是我在查看發(fā)現(xiàn)它原本有一個atl.grid，上傳是2010年8月28，大小也是3.1M，故我斷定，服務(wù)器上的atl.grid就是tables里面的ANU100826_grid.atl。于是我在服務(wù)器上下載了一個3.1M的atl.grid拷貝到test/tables中，然后在032里面制作了一個atl.grid的鏈接。（當(dāng)然也可以不用下載atl.grid,直接制作指向ANU100826_grid.atl的鏈接）綜上所述：更新tables常更新以下11個文件（各種年表除外）：但常常guess_rcvant.dat用不到,可以不更新制作鏈接特別注意以下幾個鏈接：在更新tables,建立完連接后執(zhí)行l(wèi)s–al檢查鏈接情況，紅色的是死鏈接：發(fā)現(xiàn)atml.grid的鏈接已斷，于是重新設(shè)置即可。至于剩下的紅色死鏈接，我全是忽略了，因為我也不會設(shè)置。（好像那個是在正式處理過程中產(chǎn)生的，可以忽略）總的來說設(shè)置鏈接是及其容易出錯，要小心又小心。6.控制文件配置6.1準備測站列表sites.defaults在test/tables下把需要處理的測站加入該表，包括IGS測站和本地測站，并確定工程名。然后將site.defaults在032下設(shè)置一個鏈接6.2準備測站近似坐標(biāo)文件lfile.制作完lfile.文件后不要忘了將lfile.拷貝到test/tables中,同時在032中設(shè)置對應(yīng)鏈接lfile.文件是測站近似點坐標(biāo)（越精確越好,10.4以前的版本里面存放的經(jīng)緯度信息，而10.4及其之后的版本直接存放x,y,z信息）GAMIT10.5的lfile和apr文件是同樣的格式，所以我們可以認為apr文件就是lfile。該文件中保存的是各個站的初始坐標(biāo)，至于該文件中的測站坐標(biāo)要達到什么樣的精度，可參考如下兩條：如果是該站作為基線解算起算點,則該測站的坐標(biāo)做保持在5cm以內(nèi)，最好不要超過10cm。如果該測站在基線解算中不作為起算點，則偽距單點定位的精度即可。對于IGS站，獲取其在特定框架下的精確坐標(biāo)方法：首先要確保test/rinex下有itrf08.apr方法1：我用的方法（可用于非IGS站）：在我的電腦上試了幾次，老是最后一步出錯，不推薦，但也有教程上是這個方法：在rinex文件夾下先生成l文件，然后生成l.apr，最后生成lfile.依次執(zhí)行代碼：grepPOSITION*.13o>lrx2aprl2013032gapr_to_lllfile.**2013032方法2：方法3： 偽距單點定位計算非IGS站的初始坐標(biāo)實例（假設(shè)shao是非IGS站）：執(zhí)行完之后會生成兩個不同格式的測站坐標(biāo)文件l和shao.apr，然后執(zhí)行雙差定位計算測站的初始坐標(biāo)實例：在此例中假設(shè)shao站是未知坐標(biāo)的測站，而bjfs是已知坐標(biāo)的測站（或起算站）綜上所述：在生成life.文件時，IGS站的坐標(biāo)可以直接從框架獲得，或者利用雙差定位求得。而非IGS站則可以根據(jù)要求精度的高低選擇：直接提取頭信息偽距單點定位雙差定位法方法1常常在偽距單點定位和雙差定位法中選擇，然后將所須的坐標(biāo)信息拷貝到lfile.中即可，并將lfile.拷貝到test/tables中,同時在032中設(shè)置對應(yīng)鏈接現(xiàn)在，我假設(shè)bjfs，lhaz兩站點是精度較高的IGS站，坐標(biāo)直接從框架中獲得，而shao,urum站是非IGS站，近似坐標(biāo)通過偽距單點定位獲得。最后按照特定的格式均復(fù)制到lfile.文件中然后將lfile.拷貝到test/tables中,同時在032中設(shè)置對應(yīng)鏈接好像不需要在032中設(shè)置鏈接，但是怕出錯，我還是設(shè)置了，希望沒有弄巧成拙6.3文件的更新最后要在032中檢查鏈接是否指向了更新后的文件在gamit全英文的官方教程（版本10.6,24頁）里面有這樣一句話：可見station.Info的鏈接可以指向下圖中選中的文件，也就是說倘若你的test/tables中沒有文件，可以自己新建并指向他們。并且這個MIT,SOPAC的文件好像可以在上更新，只是登不上去。練手的話湊合著能用。在test/tables目錄下：先建立一個sitelist的文件，并錄入測站名字（每行首字母空格），不分先后順序然后終端打開tables并執(zhí)行代碼：顯示fileshave3sites，這是錯誤的，一般來說總共有多少測站參與計算就會有多少個sites，于是回去檢查sitelist,，發(fā)現(xiàn)拉薩站寫錯了，應(yīng)該是lhaz（所以在操作過程中一定要小心又小心，切記不能將單詞寫錯，同時也要不害怕出現(xiàn)錯誤，只要根據(jù)其提示，順藤摸瓜，基本上能修改錯誤的地方，實在不行重來）還發(fā)現(xiàn)，Receviername**notfoundinrcvant.dat,網(wǎng)上百度了一下說要修改，說得又很模糊，所以不會。又去問了一下師兄，他說可以忽略，所以我也就忽略了（誰知道的話告訴我一聲）。執(zhí)行上述代碼結(jié)束后最后一行顯示：生成了一個.new，此時需要將其改名為，同時刪除原有的的鏈接。再執(zhí)行一次更新代碼：在生成的文件中會有下圖中紅框的字符串，表示更新成功：最后在032中檢查鏈接的指向，如果沒有的話，自己手動用代碼建立一個。7修改配置process.defaults、sittbl.與sestbl.文件7.1配置process.defaults它里面包含目錄信息，改動較小。需要在032中設(shè)置鏈接7.2配置sittbl.坐標(biāo)約束文件對高精度的已知坐標(biāo)強約束，待求點坐標(biāo)松約束非起算點約束量可以更大,如100m。并且如果IGS站列表中沒有你所處理數(shù)據(jù)中的IGS站，需要自己添加，格式一致。7.3配置sestbl.改正模型文件配置準則如下面的截圖：下面是我的部分設(shè)置：我沒有met.list,met.grid,也沒有map.list和map.grid,所以我設(shè)置的是N（不使用）我的atml.list和atl.list是死鏈接，我又不知道如何制作他們（好像在麻省理工的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是輸出文件當(dāng)然也可以在將auto0320.13n改名字后，執(zhí)行makejbrdc0320.13njbrdc3.0325．文件檢查6．生成中間文件7．生成bat文件8．基線解算最后一步還是出錯了，但是我百度了一下“RCnotexecutedSTOPFATALError:Stopfromreport_stat”在這篇文章的最后一頁找到解決方法，也不知道對不對，反正最后是成功解算了。在032中找到上圖的兩個文件并刪除，然后重新執(zhí)行基線解算。中間有很多省略結(jié)算結(jié)果精度還算過得去但是我往往以otesta.032為準qexpta.doy:解算記錄qexpta.doy:解算記錄oexpta.doy為qexpta.doy的簡略版hexpta.yrdoy：協(xié)方差矩陣、參數(shù)平差值9．基線解算精度評定標(biāo)準（網(wǎng)上各類教程上的標(biāo)準大同小異）標(biāo)準1：一般NRMS越小越好，一般小于0.25視為解算效果較好，本次解算的NRMS為一般NRMS越小越好，一般小于0.25視為解算效果較好，本次解算的NRMS為0.26，符合要求基線的X、Y方向誤差范圍應(yīng)控制在1cm以內(nèi)，Z方向的誤差控制在1.5cm以內(nèi)，本次解算所有基線均滿足要求本次解算x,y,z方向的誤差均有超過1cm,結(jié)果不是很理想，還過得去最后還有一個基線重復(fù)性檢驗：但是我的檢驗失敗了，是不是要進行兩天以上的數(shù)據(jù)處理才能進行重復(fù)性試驗？至此可以理解為解算完畢，然后進行基線平差，也可以使用GLBOK平差，但是我選擇了用COSAGPS進行的平差。如果想直接用GLBOK進行平差，參考安向東的處理流程。四、cosagps平差在032目錄下打開otesta.032,并在第二次出現(xiàn)基線信息的地方加入下面紅框中的一行字符串：“cosagpsforgamito-file”最后保存并復(fù)制到windows桌面下。Cosa平差打開cosagps,新建工程，設(shè)置如下，其它參數(shù)基本默認導(dǎo)入基線數(shù)據(jù)導(dǎo)入從ubuntu里面復(fù)制出來的otesta.032基線的誤差有點大（但為了走完流程，有精度強迫癥的就湊合著看吧^_^）導(dǎo)入三維已知坐標(biāo)這里就是你解算時坐標(biāo)起算數(shù)據(jù)，但是要注意框架和坐標(biāo)系的統(tǒng)一。先在點名方框內(nèi)隨便輸入數(shù)據(jù)，然后確定，如下圖在工程目錄下找到剛新建的關(guān)于起算點坐標(biāo)的文件，用記事本打開，并將已知起算數(shù)據(jù)復(fù)制粘貼到know文件中，如下圖：三維平差最后將解算出的shao,urum站的