GAMIT10.5數(shù)據(jù)處理圖文流程

GAMIT10.5數(shù)據(jù)處理圖文流程冉啟順目錄一、前言2二、數(shù)據(jù)準備31.前期準備32.工程目錄設(shè)置33.數(shù)據(jù)下載34.更新tables表，35.創(chuàng)立36.控制文件配置36.1準備測站列表sites.defaults36.2準備測站近似坐標文件lfile.36.3文件的更新37修改配置process.defaults、sittbl.與sestbl.文件37.1配置process.defaults37.2配置sittbl.坐標約束文件37.3配置sestbl.改正模型文件3三、分布處理31.Make*p32．精細星歷處理33．文件檢查34．播送星歷處理35．文件檢查36．生成中間文件37．生成bat文件38．基線解算39．基線解算精度評定標準3四、cosagps平差3五、建議與說明3六、相關(guān)資料3前言二、數(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)下載對應(yīng)的播送星歷，精細星歷，觀測值文件分別放在文件夾brdc,igs,rine*里面。當(dāng)然還涉及到下載的數(shù)據(jù)解壓，.d文件到.o文件的轉(zhuǎn)換等操作，其中播送星歷可以不改名字〔即可以不用改成brdc0320.13n〕4.更新tables表，更新準則，下面的表格為準：還有個dcb.dat需要每次處理都要更新表文件名表含義更新頻次容更新地址pole.或pole.usno極移表每次更新pub/gamit/tables/ut1.或ut1.usno地球自轉(zhuǎn)參數(shù)表每次更新月亮表使用當(dāng)年的表太陽表使用當(dāng)年的表nutabl.2012章動表使用當(dāng)年的表leap.sec跳秒表每次更新gdetic.dat,坐標參數(shù)文件不常更新antmod.dat,天線相位中心改正正表有新天線時更新svnav.dat,衛(wèi)星列表有新衛(wèi)星時更新rcvant.dat,接收機天線對照表有新天線時更新svs_e*壞衛(wèi)星信息文件每次更新vmf1映射函數(shù)模型使用當(dāng)年的表otl.grd海潮格網(wǎng)數(shù)據(jù)使用當(dāng)年的表常見文件說明更新的準則是：“更新頻次〞中的紅字局部每次處理都更新一下，同時結(jié)合本次處理與上次處理的年份是否一樣來進展年表的選擇。原本有些文件是每周更新一次，但是作為新手為了防止出錯還是每次更新吧有些時候，有些文件找不到，或者里面沒有，則可以忽略。如svs_e*，gdetic.dat常常找不到我在更新時，這四個文件在效勞器上沒找到下列圖是我更新的文件。原本還應(yīng)該有2013年的各種年表，但是我之前處理過一次2013年的數(shù)據(jù)，并且我的gamit安裝路徑/opt/gamit/GAMIT10.5/tables已經(jīng)有了2013年的這些年表。將上述更新的數(shù)據(jù)復(fù)制到gamit安裝目錄下的tables文件夾，并替換原有文件。我的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的構(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，故我斷定，效勞器上的atl.grid就是tables里面的ANU100826_grid.atl。于是我在效勞器上下載了一個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準備測站近似坐標文件lfile.制作完lfile.文件后不要忘了將lfile.拷貝到test/tables中,同時在032中設(shè)置對應(yīng)lfile.文件是測站近似點坐標〔越準確越好,10.4以前的版本里面存放的經(jīng)緯度信息，而10.4及其之后的版本直接存放*,y,z信息〕GAMIT10.5的lfile和apr文件是同樣的格式，所以我們可以認為apr文件就是lfile。該文件中保存的是各個站的初始坐標，至于該文件中的測站坐標要到達什么樣的精度，可參考如下兩條：如果是該站作為基線解算起算點,則該測站的坐標做保持在5cm以，最好不要超過10cm。如果該測站在基線解算中不作為起算點，則偽距單點定位的精度即可。對于IGS站，獲取其在特定框架下的準確坐標方法：首先要確保test/rine*下有itrf08.apr方法1：我用的方法〔可用于非IGS站〕：在我的電腦上試了幾次，老是最后一步出錯，不推薦，但也有教程上是這個方法：在rine*文件夾下先生成lfile.rn*文件，然后生成lfile.rn*.apr，最后生成lfile.依次執(zhí)行代碼：grepPOSITION*.13o>lfile.rn*r*2aprlfile.rn*2013032gapr_to_llfile.rn*.aprlfile.**2013032方法2：方法3：偽距單點定位計算非IGS站的初始坐標實例〔假設(shè)shao是非IGS站〕：執(zhí)行完之后會生成兩個不同格式的測站坐標文件lfile.shao和shao.apr，然后執(zhí)行雙差定位計算測站的初始坐標實例：在此例中假設(shè)shao站是未知坐標的測站，而bjfs是坐標的測站〔或起算站〕綜上所述：在生成life.文件時，IGS站的坐標可以直接從框架獲得，或者利用雙差定位求得。而非IGS站則可以根據(jù)要求精度的上下選擇：直接提取頭信息偽距單點定位雙差定位法方法1常常在偽距單點定位和雙差定位法中選擇，然后將所須的坐標信息拷貝到lfile.中即可，并將lfile.拷貝到test/tables中,同時在032中設(shè)置對應(yīng)現(xiàn)在，我假設(shè)bjfs，lhaz兩站點是精度較高的IGS站，坐標直接從框架中獲得，而shao,urum站是非IGS站，近似坐標通過偽距單點定位獲得。最后按照特定的格式均復(fù)制到lfile.文件中然后將lfile.拷貝到test/tables中,同時在032中設(shè)置對應(yīng)好似不需要在032中設(shè)置，但是怕出錯，我還是設(shè)置了，希望沒有弄巧成拙6.3文件的更新最后要在032中檢查是否指向了更新后的文件在gamit全英文的官方教程〔版本10.6,24頁〕里面有這樣一句話：可見station.Info的可以指向下列圖中選中的文件，也就是說倘假設(shè)你的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í)行上述代碼完畢后最后一行顯示：生成了一個，此時需要將其改名為，同時刪除原有的的。再執(zhí)行一次更新代碼：在生成的文件中會有下列圖中紅框的字符串，表示更新成功：最后在032中檢查的指向，如果沒有的話，自己手動用代碼建立一個。7修改配置process.defaults、sittbl.與sestbl.文件7.1配置process.defaults它里面包含目錄信息，改動較小。需要在032中設(shè)置7.2配置sittbl.坐標約束文件對高精度的坐標強約束，待求點坐標松約束非起算點約束量可以更大,如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.Make*p在032目錄下執(zhí)行代碼：make*p2．精細星歷處理如下列圖執(zhí)行代碼：3．文件檢查執(zhí)行下列圖的代碼：對于代碼中的gigsf3.032,其中g(shù)表示文件類型，中間igsf表示軌道，第一步輸入的：，3表示年的最后一位，032是年積日。4．播送星歷處理其中auto0320.13n沒有改名字作為輸入文件，jauto3.032是輸出文件當(dāng)然也可以在將auto0320.13n改名字后，執(zhí)行makejbrdc0320.13nrdc3.0325．文件檢查6．生成中間文件7．生成bat文件8．基線解算最后一步還是出錯了，但是我百度了一下“RCnote*ecutedSTOPFATALError:Stopfromreport_stat〞在這篇文章的最后一頁找到解決方法，也不知道對不對，反正最后是成功解算了。.doc88./p-2.html在032中找到上圖的兩個文件并刪除，然后重新執(zhí)行基線解算。中間有很多省略結(jié)算結(jié)果精度還算過得去但是我往往以otesta.032為準qe*pta.doy:解算記錄qe*pta.doy:解算記錄oe*pta.doy為qe*pta.doy的簡單版he*pta.yrdoy：協(xié)方差矩陣、參數(shù)平差值9．基線解算精度評定標準〔網(wǎng)上各類教程上的標準小異〕標準1：一般NRMS一般NRMS越小越好，一般小于0.25視為解算效果較好，本次解算的NRMS為0.26，符合要求基線的*、Y方向誤差圍應(yīng)控制在1cm以，Z方向的誤差控制在1.5cm以，本次解算所有基線均滿足要求本次解算*,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ù)根本默