【正文】 GAMIT/CosaGPS軟件的使用③ ? 數(shù)據(jù)處理階段 ? 7) 生成 bat文件 ? 輸入命令： ? fixdrv ? 需要準(zhǔn)備的文件： ，以及該文件所列出的所有文件 ? 輸出的文件： ，以及 ? 8）基線解算 ? 輸入命令： ? csh ? 需要準(zhǔn)備文件： ，以及 ? 輸出文件： q、 o、 h等所有解算的結(jié)果文件。 2022/6/2 武 漢 大學(xué) 測(cè)繪 學(xué)院 61 GAMIT/CosaGPS軟件的使用③ ? 數(shù)據(jù)處理階段 ? B.―batch‖處理 ?使用 GAMIT軟件進(jìn)行基線解算有一定程度的了解之后，可以不用每次只執(zhí)行一條命令，可利用UNIX/LINX系統(tǒng)自帶的 shell編程語言，定制適合自己的批處理腳本 Shell命令。 ?值得注意的是： Shell腳本命令只對(duì)數(shù)據(jù)處理進(jìn)行批處理，數(shù)據(jù)準(zhǔn)備及文件準(zhǔn)備建議手動(dòng)操作。 ?以某公路隧道 GPS二等網(wǎng) 2022年 doy=086第 1時(shí)段基線解算為例，定制的批處理腳本文件 容如下所示。 2022/6/2 武 漢 大學(xué) 測(cè)繪 學(xué)院 62 GAMIT/CosaGPS軟件的使用③ ? 數(shù)據(jù)處理階段 ? “ batch‖處理 ? !/bin/csh f ? Step1： customized your experiment ? set expt = 39。shsh39。 ? set yr_4 = 39。202239。 ? set yr_2 = 39。0939。 ? set yr_1 = 39。939。 ? set doy = 39。08639。 ? set sp3type = 39。igsr39。 ? set sp3file = 39。39。 ? set brdcfile = 39。39。 ? set session = 39。039。 2022/6/2 武 漢 大學(xué) 測(cè)繪 學(xué)院 63 GAMIT/CosaGPS軟件的使用③ ? 數(shù)據(jù)處理階段 ? “ batch‖處理 ? step2： link the tables into current directory $yr_4 $doy $expt ? step3： remove broken files and restartlink files into current directory rm luntab. rm soltab. rm nutabl. rm pole. rm ut1. rm pmu. 2022/6/2 武 漢 大學(xué) 測(cè)繪 學(xué)院 64 GAMIT/CosaGPS軟件的使用③ ? 數(shù)據(jù)處理階段 ? “ batch‖處理 ln s ../tables/ ./luntab. ln s ../tables/ ./soltab. ln s ../tables/ ./nutabl. ln s ../tables/ ./pole. ln s ../tables/ ./ut1. ln s ../tables/ ./pmu. 2022/6/2 武 漢 大學(xué) 測(cè)繪 學(xué)院 65 GAMIT/CosaGPS軟件的使用③ ? 數(shù)據(jù)處理階段 ? “ batch‖處理 ? step4： remove rinex files and restrartlink rinex files into current directory rm *.09o rm *.sp3 rm *.09n ln ../rinex/${doy}/${doy}${session}/*.09o ./ ln ../brdc/$brdcfile ./ ln ../igs/$sp3file ./ 2022/6/2 武 漢 大學(xué) 測(cè)繪 學(xué)院 66 GAMIT/CosaGPS軟件的使用③ ? 數(shù)據(jù)處理階段 ? “ batch‖處理 ? step5： preparation for the experiment ? rm *.fatal ? rm ? sh_makexp expt $expt orbt $sp3type yr $yr_4 doy $doy sess 99 srin nav $brdcfile apr l${expt}${yr_1}.$doy sinfo 30 00 00 2880 ? sh_sp3fit f $sp3file o $sp3type d $yr_4 $doy r BERNE t u ? sh_check_sess sess $doy type gfile file g${sp3type}${yr_1}.$doy ? makej $brdcfile jbrdc${yr_1}.$doy ? sh_check_sess sess $doy type jfile file jbrdc${yr_1}.$doy ? makex $ ? fixdrv d${expt}${yr_1}.$doy 2022/6/2 武 漢 大學(xué) 測(cè)繪 學(xué)院 67 GAMIT/CosaGPS軟件的使用③ ? 數(shù)據(jù)處理階段 ? ： ? A)檢查所有測(cè)站點(diǎn)是否都參與計(jì)算； ? B）標(biāo)準(zhǔn)均方根 nrms解算結(jié)果為 。 ?若是 nrms大于 ，表明處理結(jié)果是有問題的（例如：周跳沒有修復(fù)；測(cè)站的起算點(diǎn)坐標(biāo)有問題等）。 ?若是 nrms小于 ，則認(rèn)為成功解算，否則需要檢查原因，重新處理； ? C) 檢查各基線文件處理精度是否符合 gamit的實(shí)際精度。 2022/6/2 武 漢 大學(xué) 測(cè)繪 學(xué)院 68 GAMIT/CosaGPS軟件的使用③ ? CosaGPS軟件的使用 ? 簡(jiǎn)要操作步驟 ? 文件 —新建工程 /打開工程 ? GPS數(shù)據(jù)處理 —基準(zhǔn)數(shù)據(jù) —三維已經(jīng)坐標(biāo) ? 二維已知坐標(biāo) ? 一維高程點(diǎn) ? 輸入地面邊長(zhǎng) ? 輸入地面方位 ? 三維已知道坐標(biāo)是指地心空間直角坐標(biāo)（ X,Y,Z）或大地坐標(biāo)（ B,L,H）；二維已經(jīng)坐標(biāo)是指平面直角坐標(biāo)（ x， y），可以是指任意坐標(biāo)系中的坐標(biāo)；一維高程是指正常高 h。 2022/6/2 武 漢 大學(xué) 測(cè)繪 學(xué)院 69 GAMIT/CosaGPS軟件的使用③ ? CosaGPS軟件的使用 ? 簡(jiǎn)要操作步驟 ? GPS數(shù)據(jù)處理 —基線數(shù)據(jù)（選取獨(dú)立基線） ? GPS數(shù)據(jù)處理 —設(shè)置（查看或設(shè)置參數(shù)值） ? GPS數(shù)據(jù)處理 —GPS三維向量網(wǎng)（三維平差） ? GPS數(shù)據(jù)處理 —設(shè)置（查看或設(shè)置參數(shù)值） ? GPS數(shù)據(jù)處理 —GPS二維聯(lián)合平差 ? GPS數(shù)據(jù)處理 —GPS高程擬合 2022/6/2 武 漢 大學(xué) 測(cè)繪 學(xué)院 70 主要內(nèi)容 ? GPS數(shù)據(jù)處理軟件介紹 ? GAMIT/GLOBK/CosaGPS軟件 ? GAMIT/GLOBK軟件的安裝 ? GAMIT/CosaGPS軟件的使用 ? 工程數(shù)據(jù)處理實(shí)例 2022/6/2 武 漢 大學(xué) 測(cè)繪 學(xué)院 71 工程數(shù)據(jù)處理實(shí)例 ? 數(shù)據(jù)來源：某國土測(cè)繪院在 2022年 5月 2日在某公路隧洞進(jìn)行了 GPS控制測(cè)量。我們?nèi)∵@天的一個(gè)時(shí)段為例，來說明利用 GMAIT進(jìn)行高精度數(shù)據(jù)處理。 ? 數(shù)據(jù)準(zhǔn)備階段 ? 文件準(zhǔn)備階段 ? 基線解算階段 ? 網(wǎng)平差階段 ? 結(jié)束語