歡迎來到Linux教程網
Linux教程網
Linux教程網
Linux教程網
Linux教程網 >> Linux基礎 >> Linux教程 >> 使用GLOBK批處理命令進行網平差

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

日期:2017/2/28 13:54:27   编辑:Linux教程

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)。

GAMIT相關閱讀:http://www.linuxidc.com/GAMIT/

GAMIT10.4及最新版本的GAMIT下載在Linux公社5號FTP服務器,具體下載見 http://www.linuxidc.net/thread-1186-1-1.html

Copyright © Linux教程網 All Rights Reserved