使用GLOBK批次處理命令進行網平差

1. 說明

首先，你要先完成基線解算。

sh_glred 在工程級目錄執行，但在解目錄（gsoln）中執行，該目錄也可以由 process.defaults 檔案中 glbpth 引數指定。

2. 引數設定

2.1 globk_comb.cmd 和 glorg_comb.cmd

若這兩個檔案不存在，sh_glred 儘管可以生成他們，但推薦手動建立。

執行命令“sh_glred -cmd”，將在gsoln資料夾下生成上述設定的模板檔案以供修改。

需在其中設定先驗坐標檔案（.apr）、解算中要用到的測站、定義參考框架時用到的測站和約束等（若使用IGS站點作為固定點，則可以使用自帶的itrf 檔案作為先驗坐標）。

2.2 sites.dafaults

每步處理中所使用的測站都由sites.dafaults 中標記指定：在 glred 組合中包含的測站由 glrepu 標記，在glorg 中用來定義參考框架的測站由 glreps 標記，要生成時間序列圖的測站由 glts 標記。預設使用所有測站並繪圖。

注意：若在sites.defaults檔案中使用glreps標記平差時的固定點，則在glorg_combNaNd中需將“stab_site clear”選項註釋掉（前邊加“x”）。

3. 命令執行

完成上述操作之後，你就可以進行批次處理網平差了。在工程目錄執行：

sh_glred -s<yr1 doy1 yr2doy2>  -expt <expt>  -net <networks>  -local  -nettext <char>  -yrext<year> -ncomb <num>  -stnfo<station.info>  -cmd  -opt <A F H L U G E K C R>

其中-s 指定開始時間（yr1 doy1）和結束時間（yr2 doy2），expt 是 4 個字元的工程名，network 表示包含SOPAC 的 H-檔案，SOPAC預設選項是all（所有），或者僅組合某些解，假如-local 被指定了，指令碼會僅僅處理可以獲得資料範圍內的天；否則它會自動去 SOPAC 查詢 H-檔案，處理所有指定天數內的資料。要處理的天也可用 -d yr days或者採用-r days 詳細指定。可選引數-ncomb用來指定每一個組合包含多少天，預設為 1 天，但也可以用指令碼產生每周或每月的本地或全球檔案的平均值。

下面是opt選項可選引數的詳細釋義（摘自GAMIT 10.35中文文件）

•<R> 在開始之前把舊的 h*.gl?檔案從 glfpth 目錄下刪除掉。

•<F> 從 SOPAC 網站上下載全球 h 檔案，根據<networks>引數項提供的名稱搜尋檔案。

•<H> 對所有現有的或者連結到 glfpth(通常是 procdir/glbf)的 ascii 檔案執行 htoglb。

•<LA> 連結所有本地存檔的 ascii H 檔案(即所有的 h[net|expt]?.yyddd)到組合中去。

•<LB> 連結所有本地存檔的 ascii H 檔案(即所有的 h[net|expt]?.yyddd)到組合中去。從 htnd 開始向下搜尋，其中 htnd 是系統預設的。

•<LC> 連結本地生成的組合二進位制 H 檔案(名為*.GLX) 。

•<U> 對二進位制 h 檔案執行 hfupd(尚未試)。

•<G> 執行 glred 以進行組合或重新組合。

•<E> 執行 ensum 和 sh_baseline 進行繪圖。

•<C> 後壓縮 ascii H 檔案，刪除所有的連結，分別複製 ascii H 檔案和二進位制 H 檔案為 glbpth/ascii_yyyy    和 glbpth/bin_yyyy, 其中 glbpth 是由 process.defaults 指定的。

•<A> 完成所有選項

例：  % sh_glred -s $year1 $sdoy $year2 $edoy -expt $expt  -local -opt H G E

$year1和$sdoy為要進行網平差的開始年與年積日，$year2和$edoy為結束年與年積日。

4. 解算結果

命令執行完成後，在gsoln/資料夾下即可找到平差結果檔案*.prt和*.org。

固定點有“*”號標記。

PS：雖然這種方法可以快速得到網平差結果和站速度，但是在GAMIT/GLOBK軟體的官方指南中對這種平差方法是不推薦的。推薦方法是只使用該方法合併GPS網得到站坐標，然後使用globk命令呼叫globk_vel.cmd和glorg_vel.cmd得到站速度（因整理時間較晚，有些可能記得不對，請讀者自行查閱 GLOBK_Ref）。