GAMIT基線解算 批處理

分步處理2023年001,002天（shao hkkt hkks，後面兩個站點為香港cors站）

資料準備工作：

1）mkdir demo

2）cd demo

3）mkdir rinex brdc igs

4）cd brdc

wget -c

wget -c

解壓:uncompress *

5）cd igs

wget -nd

wget -nd

解壓:uncompress *

6）cd rinex

sh_get_rinex -archive sopac cddis  -yr 2017 -doy 001 -ndays 1 -sites shao

sh_get_rinex -archive sopac cddis  -yr 2017 -doy 002 -ndays 1 -sites shao

7）鏈結表檔案

sh_setup -yr 2017

（在 sestbl. 檔案中將海潮模型設定為：

use otl.list = n; ocean tidal loading list file from oso

use otl.grid = y; ocean tidal loading grid file, gamit-format converted from oso）

8）更新station.info

在rinex目錄下：

sed -n '2,6p' ../tables/station.info > station.info

sh_upd_stnfo -files *.17o

cp station.info ../tables

9）更新l檔案

在table表目錄下執行

# 提取觀測檔案中的概略座標

，完成此步驟後繼續看下面

rx2apr lfile.rnx 2017 001# 生成 lfile.rnx.apr 檔案

gapr_to_l lfile.rnx.apr lfile. 2017 001# 生成 lfile. 檔案

sh_gamit -expt text -d 2017 001 -orbit igsf -pres elev -noftp -dopt d ao x c > sh_gamit.log

如：xdg-open hkks.002_elev_res.png

xdg-open hkks.002.png

空間面積解算

通常情況下，在地球表面的乙個小範圍區域，投影變形較小，我們用平面面積代替球面面積，兩點間的弧長等效為兩點間的距離，兩條相交的曲線間的夾角等效為兩線段間的夾角。經證明 當邊長小於200km時，橢球面上三角形的解算完全可以在平面上進行。此時，橢球面三角形與球面三角形各對應角的差異小於0.001秒，邊長差...

LPMS IMU姿態解算

參考文章 ahrs姿態解算說明 加速度 陀螺儀 磁力計原理及原始資料分析 ahrs俗稱航姿參考系統，ahrs由加速度計，磁場計，陀螺儀構成，ahrs的真正參考來自於地球的重力場和地球的磁場 他的靜態終精度取決於對磁場的測量精度和對重力的測量精度 而則陀螺決定了他的動態效能。如何用加速度計和陀螺儀的資...

姿態解算解析

用四元素法進行姿態解算，其步驟如下 1.四元素初始化 第一次解算時計算 靜止狀態下由加速度，磁力計值求取初始 gamma theta psi gamma frac theta frac psi frac 初始化四元素 left begin end right left begin cos frac ...