chaoshu的個人博客分享 http://www.ueservicedoffices.com/u/chaoshu

博文

[轉載]globk中的控制文件

已有 1062 次閱讀 2017-9-6 21:59 |個人分類:GAMIT/GLOBK|系統分類:科研筆記|文章來源:轉載

1、globk_comb.cmd

* This group must be first
eq_file ../tables/eq_renames
make_svs ../tables/sat1.apr
com_file comb.com
srt_file comb.srt
sol_file comb.sol
* Earth-rotation values (not needed if pmu free in final combination)
in_pmu ../tables/pmu.bull_b
* File(s) for a priori coordinates
apr_file ../tables/itrf00.apr
* Sites to include in the solution (default is all)
use_site clear
use_site mate ankr_gps kit3 nico sofi zwen kant
use_site –zwen
* (1) Max chi**2, (2) Max prefit diff, (3) Max rotation; defaults are 100 10000 10000
max_chi 30 50 2000.0
* Apply the pole tide whenever not applied in GAMIT
app_ptid ALL
* Invoke glorg for stabilization
org_cmd glorg_comb.cmd
org_opt PSUM GDLF CMDS
* Set minimal globk print options since using glorg output
prt_opt NOPR
* Allow the network to be loose since using glorg for stabilization
apr_neu all 10 10 10 0 0 0
* Satellites are loose if using global data -- ignored if GAMIT BASELINE mode
x apr_svs all 100 100 100 10 10 10 1R
* but tight if not combining with global data
x apr_svs all .05 .05 .05 .005 .005 .005 .01 .01 FR
* EOP loose if estimating rotation in glorg
apr_wob 100 100 10 10
apr_ut1 100 10
* Write out an h-file if needed for future combinations
x out_glb H------.GLX

上面第一組命令定義了glred/globk中使用的文件,這些命令必須放在控制文件的前面,這樣程序才能很好的運行。對于時間序列計算指令,除了 eq_file所有的都可以認為是臨時文件。下面將要介紹的eq_file使得我們可以排除測站或重命名測站以解釋由于地震或儀器變化造成的地形突變。
緊接著的兩個記錄具體指定了用于地球旋轉和坐標的先驗文件。如果在此忽略的話,將會采用h文件中的,這樣在很多情況也是可行的。
use_site命令允許我們選擇解算中包含的測站。這里我們展示了多余的命令以舉例說明它們的用法:默認的是all;use_site clear命令將刪除所有測站,在它的后面應該跟著你希望解算中包含的測站清單;在列出的測站名稱前加負號,將會在清單中刪除它。對于從全球處理的h文件 中選擇測站,我們可以執行use_pos以在一定的地理區域內使用或排除測站(查看globk的幫助文件或者globk參考手冊的3.1節)
當數據是“壞”的時候 max_chi命令允許自動刪除h文件,在globk中檢查高卡方值探測,坐標相對于先驗值調整較大,或者測站網較大的旋轉。
包含glorg命令文件的名字會告訴globk運行完globk后執行glorg for stabilization(框架定義)。org_opt選項規定了glorg打印總結中包含的內容(org文件)。PSUM表示位置平差的總 結;GDLF表示列出輸入的h文件的名字;CMDS表示顯示globk 和glorg命令文件。如果你正在執行glorg,然后globk打印將不起作用,且通過prt_opt NOPR壓縮。
apr_形式的命令告訴globk估計參數和指定使用的先驗值。如果你通過在glorg中估計旋轉和轉換固定框架(通常是穩健網的情況下),然后坐標 (apr_neu)和EOP(apr_ut1 apr_wob)約束應該是松弛的,如上所示。如果只在glorg中估計轉換,在小區域中可選的是少于4個站,然后應該嚴密約束EOP。為了應用有限的約 束,要為每個站設置合適的apr_neu,且不能執行glorg。
out_glb命令將會在執行時告訴globk將解算輸出以h文件的形式輸出。這對于聚集按周或按年的平均值和將自己的處理與全球或局部h文件聯合解算很有用。


2、glorg_comb.cmd

* Stabilization using IGS stations
apr_file ../tables/itrf00.apr
* Regional stabilization using coordinates from a prior solution
x apr_file ../tables/globk_vel_990306.apr
* Define the stabilization frame
stab_site clear
stab_site mate kit3 nico sofi zwen
* Set parameters to estimate in stabilization
pos_org xtran ytran ztran xrot yrot zrot
* Default height ratio is 10; increase if heights questionable
cnd_hgtv 1000 1000 2.0 10.0
* Default iterations is 4, relative weight ratio 50%, horizontal sigma cuoff 4.0
x stab_ite 4 0.8 3.

定義框架通過在glorg中具體指定一系列先驗坐標已知且滿足框架定義要求的測站來實現(通常水平式1-2mm,垂直2-5mm),想要估計的參數是為了 使這些測站的坐標調整量最小。如果在glorg控制文件中包含apr文件,它的坐標將會取代gamit或globk中使用的坐標,所以對于三個程序沒有必 要使用相同的先驗坐標。對于glorg,著重關心的是穩定測站的坐標。
pos_org命令指定了估計的參數,只進行平移(忽略后三個記錄) or 轉換和旋轉.stab_site的使用規則同globk控制文件中的use_site一樣。
最后兩個記錄控制固定過程中的(in the stabilization)準觀測值(測站坐標)的權。對于cnd_hgtv,前兩個值具體指定相對于水平坐標的downweighting of height(兩個值分別對應位置和速度)。默認的是10,假設高度的不確定性事水平坐標的三倍;但是如果你想定義一個水平框架或知道高度估計不可靠,你 可以提高這個比率。最后兩個值限制高度不確定性,允許包含一個測站和刪除不穩定的測站????。對于stab_ite,第一個值是 stabilization過程中迭代的次數,第二個值是迭代過程中重新加權的數量,第三個是在刪除測站前允許的相對于不確定性的殘差比率。你可以在定義 框架的過程中改變最后一個值以調整包含或排除的測站。這些命令的詳細描述可看glorg的幫助文件和globk參考手冊的3.2節。
在固定過程中允許平移和旋轉是absorb空間相關噪聲最普遍的方法,但是如果你只有少量的測站或異常的幾何體時要小心。例如,如果你只有三個站,包含旋 轉,你從九個觀測值(每個站三個坐標)中估計六個參數,將會導致只有三度的冗余度。而且即使有很多測站,如果其中一個離其他的站很遠的話,長基線會使 absorb所有的旋轉自由度。如果測站可靠,不關心測站殘差的話這并不是一個致命的缺點,但是你需要知道this chararistic of the stabilization.
如果你需要提供坐標或速度的step changes,或在某個時期排除測站,可以在globk中使用地震文件(eq_file)。這個文件有兩種記錄。第一個重命名測站(為globk創建新的邏輯上的測站)基于它們接近地震的程度:
# November 3, 2002, UT=22:12 Lat = 63.52 Lon -147.53 Depth=5 km M 7.9
eq_def DN 63.52 -147.53 700 5. 2002 11 3 22 12
eq_cosei DN 1.0 1.0 1.0 1.8 1.8 0.7
eq_post DN 100 0.1 0.1 0.1 1.8 1.8 0.7 Postseismic: 0.1 mm**2/day and spatial about same
eq_rename DN
eq_file文件中包含以上各行,開始和結束時間晚于2002.11.3的22:12的距離Denali地震700公里內的測站將會被重命名。 eq_cosei和eq_post命令描述了co-seismic和post-seismic變化應用的限制(查看globk幫助文件和globk參考手 冊的3.1節)。第二種類型重命名測站的命令按照具體的說明解釋數據質量和儀器的變化:
rename mate_gps mate_1ps 1999 6 18 0 0
rename kit3_gps kit3_xcl 1997 5 14 0 0 1997 5 18 0 0
rename madr_gps madr_xhi emed 1996 9 26 0 0 1996 10 27 0 0
第一行說明在Matera使用新接收機導致的位置變化;第二行是刪除Kitab四天中的壞數據;第三行說明在h文件中名字包含字節emed的使用了不正確 的天線高度?梢允褂胓lorg中的equate特征強制madr_gps和madr_xhi使用相同的水平調整(查看第四章)。也可以使用rename 命令改變測站位置而不用引入一個新的測站,應用已知的當地測量的偏差或地震模型中假定的偏移:
rename tigr_gni tigr_gni 2001 2 28 0 0 2100 1 1 0 0 0.0032 0.0033 0.0000 NEU
該記錄使用了測站坐標北和東方向基于先驗模型3mm的偏差。



http://www.ueservicedoffices.com/blog-858128-1074604.html

上一篇:[轉載]globk中的衛星軌道約束
下一篇:[轉載]Antenna Phase Centers

0

該博文允許注冊用戶評論 請點擊登錄 評論 (0 個評論)

數據加載中...

Archiver|手機版|科學網 ( 京ICP備14006957 )

GMT+8, 2019-5-26 01:16

Powered by ScienceNet.cn

Copyright © 2007- 中國科學報社

返回頂部
时时彩平台