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参数进行星间差分分析。