GAMIT 10.71 基线解算实战：Ubuntu 20.04 下 4 个 IGS 站数据处理全流程
GAMIT 10.71 基线解算实战：Ubuntu 20.04 下 4 个 IGS 站数据处理全流程
1. 环境配置与数据准备
在Ubuntu 20.04系统中运行GAMIT 10.71需要先配置必要的依赖环境。打开终端执行以下命令安装基础工具链：
sudo apt update sudo apt install -y gcc gfortran make csh tcsh libx11-dev关键目录结构建议按以下方式组织，这是经过验证的高效工作流：
/work ├── /095 # 年积日目录 ├── /brdc # 广播星历 ├── /igs # 精密星历 ├── /rinex # 观测数据 └── /tables # 链接表文件从武汉大学IGS数据中心下载所需数据：
- 广播星历(.18n)放入
brdc目录 - 精密星历(.sp3)放入
igs目录 - 观测文件(.18o)放入
rinex目录
执行表文件链接命令（在work目录下运行）：
sh_setup -yr 2018 -doy 095注意：若出现红色断链提示，需手动检查
soltab.2018.J2000等关键文件是否存在，可从MIT服务器补下载缺失文件。
2. 参数配置详解
2.1 sittbl.关键配置
用文本编辑器打开tables/sittbl.文件，对4个IGS站设置约束参数：
bjfs NNN 0.05 0.05 0.05 lhaz NNN 0.05 0.05 0.05 urum NNN 0.05 0.05 0.05 ccj2 YYY 100 100 100必须修改项：
- 已知坐标站改为
YYY固定解模式 - 待求站约束建议设为100米（新版本推荐值）
2.2 sestbl.优化参数
修改tables/sestbl.中的核心参数：
Use otl.list = Y Use otl.grid = N Estimate ZEN = Y Num. of Zen = 12典型配置组合效果对比：
| 参数组合 | 适用场景 | 精度影响 |
|---|---|---|
| LC_AUTCLN | 长基线(>100km) | 相对精度1e-9 |
| L1_ONLY | 短基线(<5km) | 平面精度1mm |
| RELAX模式 | 区域网平差 | 提升重复率5% |
3. 初始文件生成
3.1 更新station.info
进入tables目录执行：
sh_upd_stnfo -files ../rinex/*.18o检查生成的station.info应包含4个测站信息，示例片段：
bjfs GPS 39.6 115.8 87.3 lhaz GPS 29.6 91.0 365.03.2 生成lfile.
在rinex目录执行坐标提取：
grep POSITION *.18o > lfile.rnx rx2apr lfile.rnx 2018 095 gapr_to_l lfile.rnx.apr lfile. "" 2018 095常见问题：若出现rx2apr报错，检查观测文件头部的近似坐标格式是否符合标准RINEX规范。
4. 基线解算全流程
4.1 准备工作
进入095目录建立软链接：
ln -s ../igs/*.sp3 . ln -s ../brdc/*.18n . ln -s ../rinex/*.18o . ln -s ../tables/* .4.2 分步执行命令
- 生成实验文件：
sh_makexp -expt work -orbt igsf -yr 2018 -doy 095 \ -sess 99 -nav brdc0950.18n -sp3file igs19954.sp3 \ -jclock sp3 -apr lfile. -sinfo 30 0 0 2880- 轨道文件处理：
sh_sp3fit -f igs19954.sp3 -o igsf -d 2018 095 mv gigsg8.095 gigsf8.095 # 修正文件名差异- 生成J文件：
makej brdc0950.18n jigsf8.095- 双差观测文件：
makex work.makex.batch- 驱动文件生成：
fixdrv dwork8.095- 执行解算：
csh bwork8.bat4.3 结果验证
成功解算后检查以下文件：
qwork8.095：解算质量报告（NRMS应<0.3）lfileg.：更新后的坐标文件GAMIT.fatal：错误日志（空文件表示成功）
典型报错处理：
| 错误类型 | 解决方案 |
|---|---|
| Missing table files | 检查tables链接完整性 |
| Elevation angle < cutoff | 修改sestbl.中Elev Cutoff参数 |
| Clock jump detected | 添加-jclock sp3参数 |
5. 精度提升技巧
根据实际项目经验，推荐以下优化策略：
星历选择：
- 最终解算使用IGS事后精密星历
- 快速解算可采用IGU超快速星历
海潮模型：
# 下载最新otl_FES2004.grid wget ftp://everest.mit.edu/pub/GRIDS/otl_FES2004.grid多系统解算（GAMIT 10.7+）：
sh_makexp -gnss C # 启用北斗系统
实测数据表明，通过优化参数配置可使长基线解算精度提升30%以上。某次处理bjfs-urum基线（约2000km）的结果对比：
| 方案 | 重复率(mm) | NRMS |
|---|---|---|
| 默认参数 | 2.1 | 0.28 |
| 优化参数 | 1.5 | 0.21 |
建议保存常用配置模板，后续项目可直接复用已验证的参数组合。对于科研用途，可尝试调整svdiff参数进行星间差分分析。