高程拟合探讨 目录1、什么是TGO2、TGO固定差高程拟合差距3、TGO和七参数4、选择合适的高程拟合方法5、实例南方工程之星5.1加权平均固定差5.2平面拟合4-6个点5.3曲面拟合≥75.4垂直平差TGO5.5七参数6、高程转换、重力场模型、RTK 外业校正完整技术文档1、什么是TGOTGO方法是天宝Trimble公司的Trimble Geomatics Office软件中用于高程拟合的一种算法。1.1原理TGO算法是在几个控制点的中心拟合出一个相对准确的原点假设测区内高差拟合值均匀变化通过北、东方向的高程拟合值变化量来确定待转换坐标的高程拟合值。1.2计算所需参数TGO软件的高程转换模型包括五个参数分别是常数平差、北斜坡、东斜坡、原点北、原点东。通过这些参数来建立高程拟合模型实现从大地高到正常高的转换。1.3适用场景TGO方法一般适用于控制点数量为3 - 5个且分布能形成一个面同时高差拟合值相对均匀的情况。例如在一些地形相对平坦、测区范围适中且控制点分布较为合理的区域使用TGO方法进行高程拟合能够获得较好的精度。1.4操作步骤在TGO软件中首先要进行二维网约束平差然后选择合适的大地水准面模型如EGM96模型加载该模型后固定必要数量的点的正常高通过平差计算得出其余点的正常高。与其他高程拟合方法相比TGO方法在满足其适用条件时能够利用多个控制点的信息考虑到测区的平面位置关系对高程进行拟合通常比简单的固定差改正精度要高又不像曲面拟合那样对控制点数量和分布要求过于严格在实际应用中具有一定的灵活性和实用性。2、TGO固定差高程拟合差距在GPS测量中TGO固定差高程拟合差距受多种因素影响具体如下- 测区范围一般来说测区范围越大TGO固定差高程拟合的差距可能越大。例如在1公里范围内固定差改正数变化可能在5cm以内能满足1500测图的5cm高程精度要求而在20公里范围内固定差改正数可能会变化一米多。- 控制点分布TGO算法需要至少3个控制点并且要求控制点分布能形成一个面且高差拟合值相对均匀。如果控制点呈狭长型的线性分布或者分布不合理会导致高程拟合精度下降差距增大。- 地形起伏TGO方法适用于山地等地形但对于高程起伏很大、地形复杂的测区TGO固定差高程拟合可能无法很好地适应会产生较大差距此时可能更适合采用曲面拟合等方法。- 高程异常变化高程异常值受多种因素影响在小范围内小于5km一般可当做一个固定值处理但测区范围大于5km时高程异常变化带来的影响就需要考虑否则会使TGO固定差高程拟合产生较大差距。3、TGO和七参数是在测量领域中用于不同目的的两种方法或参数集以下是它们的比较原理方面- TGO在几个控制点的中心拟合出一个相对准确的原点假设测区内高差拟合值均匀变化通过北、东方向的高程拟合值变化量来确定待转换坐标的高程拟合值。- 七参数用于两个不同坐标系之间的转换包含三个平移参数X、Y、Z方向的平移量、三个旋转参数绕X、Y、Z轴的旋转量以及一个尺度因子考虑了空间位置、方向和尺度的变化。适用场景方面- TGO适用于控制点数量为3 - 5个且分布能形成一个面同时高差拟合值相对均匀的情况。在地形相对平坦、测区范围适中的区域能获得较好精度。- 七参数适用于测区范围在几十公里以上的大区域且需要至少3个以上的控制点。在进行大范围的空间坐标转换如不同椭球之间的坐标系统转换时使用能保证较高的平面精度但高程精度可能需分区转换来保证。精度方面- TGO在满足适用条件时能获得一定精度但如果控制点分布不合理或测区地形复杂、高差变化大精度会受影响。- 七参数由于考虑了更多的转换因素且适用于大区域当控制点选择合理、分布均匀时平面转换精度较高但高程精度在远距离时可能会出现一定误差需采取措施保证精度。操作复杂度方面- TGO在TGO软件中操作相对简单通过二维网约束平差等步骤固定必要数量点的正常高即可进行计算。- 七参数求解七参数需要至少3个控制点且对控制点的分布要求较高不能近似在一条直线上。计算过程相对复杂需要使用专业的测量软件或工具进行求解。4、选择合适的高程拟合方法需综合考虑测区条件、精度要求等因素具体可参考以下要点1. 测区范围与地形特点- 小范围5km且地形平坦可选用固定差改正如TGO固定差假设测区内高程异常为固定值操作简单且效率高。- 中等范围5-20km或地形起伏较小优先考虑平面拟合或多项式曲面拟合如二次曲面通过控制点拟合平面或曲面模型适应一定范围内的高程异常变化。- 大范围20km或地形复杂如山地、丘陵推荐多面函数拟合或BP神经网络拟合能更灵活地适应复杂地形下的高程异常非线性变化提升拟合精度。2. 控制点数量与分布- 控制点少2-3个且分布均匀可尝试固定差改正或平面拟合但需注意控制点需覆盖测区主要范围避免线性分布。- 控制点多≥4个且分布成面优先选择多项式曲面拟合如二次、三次曲面利用更多点的信息拟合曲面提高模型精度若控制点分布极不规则可考虑多面函数拟合通过局部插值增强适应性。3. 精度要求与应用场景- 低精度需求如厘米级小范围测区可直接用固定差改正中等范围可用平面拟合或简单的曲面拟合。- 高精度需求如毫米级需结合高精度大地水准面模型如EGM2008或采用GPS水准联合拟合通过实测水准点与GPS点结合提升高程转换精度。- 工程测量如道路、建筑若测区狭长如道路沿线可采用沿线分段拟合避免大范围拟合误差若为区域性工程优先用曲面拟合并保证控制点均匀分布。4. 已有数据与模型支持- 若存在区域大地水准面模型如省级或市级精化模型可直接调用模型进行高程异常改正效率高且精度有保障。- 无现成模型但有实测水准点利用水准点与GPS点构建拟合模型如多项式拟合实测数据越密集拟合效果越好。5. 其他方法对比与选择建议方法优点缺点适用场景固定差改正操作简单计算快仅适用于极小范围5km平坦区域低精度需求平面拟合模型简单需控制点少无法适应地形起伏中等范围、地形平坦区域多项式曲面拟合适应性强精度较高高阶多项式可能过拟合5km、地形略有起伏的区域多面函数拟合局部插值精度高计算量大需较多控制点地形复杂、控制点分布不均区域GPS水准联合拟合精度最高依赖实测水准成本高需大量水准数据高精度工程测量、大地水准面精化总结步骤1. 先确定测区范围、地形及控制点数量2. 根据精度要求筛选可行方法如低精度选固定差高精度选曲面或联合拟合3. 结合已有数据如大地水准面模型、水准点选择效率与精度平衡的方案4. 可通过对比不同方法的拟合残差如均方根误差RMSE验证效果选择残差最小的方法。5、实例南方工程之星5.1加权平均固定差5.2平面拟合4-6个点5.3曲面拟合≥75.4垂直平差TGO5.5七参数6、高程转换、重力场模型、RTK 外业校正完整技术文档本文档整合大地测量高程全套核心知识点包含高程基础概念、重力场模型阶次解析、EGM 全球模型与省级精化大地水准面对比、RTK 外业校正各类算法TGO 平面拟合、二次曲面、四 / 七参数、一体化七参数弊端、坐标数据文件说明、误差来源与分场景作业方案适配地籍、矿山、沉降监测、土方勘测等工程场景。第一章 高程基础核心定义1.1 两类核心高程大地高椭球高 H地面点沿法线到参考椭球面的垂直距离GNSS 接收机直接输出coorblh.dat 文件中 H 字段无工程使用价值仅作为转换中间量。正常高85 国家高程海拔 h地面点到 1985 青岛平均海水面似大地水准面高度施工放样、规划、国土报批唯一法定高程。1.2 高程异常 ζ转换核心公式ζ高程异常椭球面与似大地水准面垂直差值已知大地高减去高程异常即可得到法定海拔所有重力场模型、高程拟合算法本质都是求解测区各点 ζ 值。1.3 两套基准面关键区别似大地水准面国内工程基准依托青岛验潮站平均海平面省级精化模型以此约束直接输出 85 高程EGM/EIGEN 全球模型基准全球海洋平均海面和我国 85 基准存在全域固定常数偏差山西地区约 0.10~0.15m。第二章 全球重力场模型阶次、分辨率与精度2.1 “阶” 的物理含义重力场采用球谐函数拆分最大阶数决定模型最小分辨尺度最小可分辨半波长公式赤道 1°≈111km。阶数越高可捕捉越小范围的局部重力起伏山区精度越高阶数仅代表重力细节分辨能力无法消除全球模型与国内 85 高程的零点固定差。2.2 主流模型参数、精度分级2.2.1 老旧淘汰模型工程禁止正式使用EIGEN‑CG01C360 阶最小分辨尺度 55km仅能刻画省级大范围重力趋势山西平地误差 5~13cm山区最大 0.7m2004 年发布仅用于历史数据复盘。EGM96360 阶缺陷同 CG01C软件基本淘汰。2.2 通用工程临时模型踏勘、DEM 粗修正EGM20082159 阶理论原生分辨率 5 弧分≈9km市面 1′格网仅插值加密不提升真实精度山西汾河平地高程异常误差 ±10cm东山丘陵单点最大偏差 0.7m缺陷零点存在 0.1m 固定偏移差值随地形非线性变化不能仅用单一常数全局修正。2.2.3 科研高阶全球模型仅学术计算XGM2019e、EIGEN‑6C41949~2159 阶精度略优于 EGM2008但仍为分米级无法满足国内工程报批规范。2.3 EGM2008 与省级精化似大地水准面差值四大来源差值≠单一固定常数由三部分叠加全域固定零点差10~15cm全球平均海面 vs 青岛 85 基准全省统一常数仅能消除一小部分偏差。分辨率与重力测点密度差主要误差EGM 原生 9km 网格内陆地面重力数据稀疏省级模型 1.5′~2.5′网格融合航空重力、全省数千 GPS 水准点可还原乡镇级重力起伏。地形改正精度差异EGM 使用全球低精度地形省级模型搭配 1:1 万本地 DEM 精细地形均衡改正山区误差大幅拉开。拟合约束条件不同省级模型强制贴合实测水准点消除区域坡度渐变误差EGM 无国内控制点约束存在东西、南北缓慢渐变偏差。2.4 精度层级排序从高至低测区现场二次曲面 / TGO 高程拟合厘米级外业首选市级城区精化似大地水准面涉密±1~1.5cm省级新版似大地水准面自然资源厅申领±2~5cmXGM2019e / EIGEN‑6C4科研分米级EGM2008踏勘粗算专用EIGEN‑CG01C、EGM96老旧淘汰高程异常值计算文件解读这些不是 DEM 高程模型是全球地球重力位模型GGM后缀.gfc是 ICGEM 标准球谐系数文件核心用途计算高程异常 N实现 GNSS 大地高 h → 正常高海拔 H转换\(\boldsymbol{H h-N}\)测绘外业 RTK、似大地水准面精化、重力勘探都要用。 全部都是2190 阶超高阶模型SGG-UGM-1 为 2159 阶分辨率约 5′约 9km。逐个介绍1. EGM2008.gfc最经典工程最常用发布2008 年美国 NGA特点国内测绘、CASS、工程之星、RTK 里默认标配模型没有 GOCE 卫星重力数据。优缺点全球通用软件兼容性最好山区精度弱于后面几款新一代模型。✅适用普通工程测绘、GNSS 高程粗转换、教学国内大部分中小项目首选。2. EIGEN-6C 系列德国 GFZEIGEN-6C.gfc / EIGEN-6C2.gfc / EIGEN-6C3stat.gfc / EIGEN-6C4.gfc迭代版本EIGEN‑6C4 是最终稳定版C/C2 是早期测试版C3stat 是静态重力场版本。数据源GRACEGOCE 卫星重力 地面重力 海洋重力2014 年定型。✅适用大地水准面精化、山区、我国西部 / 高原区域精度优于 EGM2008国内很多科研论文基准模型。3. GECO.gfc2017 年融合 GOCE 卫星重力 EGM2008 短波信息定位折中款兼顾 GOCE 长波重力与 EGM2008 的短波细节。✅适用陆海过渡带、海岸区域重力场计算。4. SGG‑UGM 系列武汉大学出品国产超高阶重力场模型SGG‑UGM‑1.gfc20192159 阶SGG‑UGM‑2.gfc20222190 阶升级版亮点专门优化了中国区域重力场采用国内大量地面重力数据在我国陆域高程异常精度普遍优于 EGM2008、EIGEN-6C4。✅适用国内国土、测绘项目国内似大地水准面精化优先推荐海南、山西项目都很适合。5. XGM2019e_2159.gfc2019 年德国 GFZ2190 阶超高阶模型特点全球海洋、近海区域精度突出GOCE 卫星重力信息利用充分近海岛礁、海洋测绘优势明显。✅适用沿海、海岛海南海域、海洋重力、海岸带项目。快速选型建议测绘工程场景普通外业 RTK 高程转换、软件兼容性优先EGM2008国内陆地项目山西、琼中陆地测绘SGG‑UGM‑2国产国内精度最好高原、山区大地水准面精化EIGEN‑6C4海南沿海、海岛、海域项目XGM2019e补充提示.gfc不能直接当 DEM 打开要在重力场计算软件ICGEM 网页计算器、GrafLab、SHM、ArcPy 自编球谐代码读取球谐系数计算每个点位高程异常。工程上全球重力模型只是基础国内正式测绘项目还需要用本地 GNSS 水准点做改正得到区域似大地水准面。第三章 RTK 外业校正全套方法平面转换 高程解耦3.1 平面坐标转换仅处理 X、Y不干扰高程3.1.1 四参数适用100km² 内中小测区如应县、榆次城区控制点≥3 个模型二维平面转换平移、旋转、尺度仅同椭球高斯平面转换误差规律远离控制点误差线性缓慢增大无非线性爆炸问题。3.1.2 布尔莎七参数标准合规用法适用大范围、跨区域测绘、沉降监测最优 6 个均匀包围控制点硬性规则两组输入坐标必须全部为椭球 BLH 大地坐标coorblh.dat目标 H 只能是 CGCS2000 椭球高作用WGS84 地心三维坐标转换至 CGCS2000 三维坐标仅完成基准切换输出依旧是椭球高必须额外搭配高程拟合才能得到 85 海拔。3.1.3 一体化七参数违规仅临时踏勘将 85 正常高强行替代目标椭球高联立求解 XYZ 七参数致命缺陷高程异常非线性曲面被线性几何模型强行拟合误差规律控制点内残差人为抹平超出范围误差非线性加速放大山区数公里高程偏差超 1m平面坐标同步偏移限制禁止用于放样、勘测定界、入库成果。3.2 高程拟合算法独立修正大地高→85 高程不改动平面 XY3.2.1 TGO 平面拟合一阶线性模型数学模型天宝 TGO 原生算法适用纯平原、高差100m 城区应县平地控制点≥3 个优势计算稳定、不易过拟合外推误差匀速线性增长短距离外推可用全域残差≤5cm。3.2.2 二次曲面拟合二阶非线性模型模型增加曲率、交叉项可模拟高程异常曲面弯曲适用测区包含丘陵、山地太原东山必须≥6 个全域均匀控制点优势山区全域厘米级精度短板远距离外推误差非线性爆发控制点分布不均极易失真。3.3 标准合规作业组合方案城镇平地正式项目七参数平面基准转换 TGO 平面高程拟合含丘陵山地、沉降监测七参数 二次曲面高程拟合小型工地短期作业四参数 TGO 高程拟合仅看图、土方概算单用 EGM2008不可用于正式施工第四章 分场景实操作业规范4.1 正式项目供地报批、地籍、建筑放样、InSAR 地表沉降监测最优方案向自然资源主管部门申领省级最新似大地水准面涉密格网直接插值高程异常无涉密数据替代方案测区均匀布设 6 个 GPS 水准控制点采用「标准七参数 二次曲面高程拟合」平面、高程完全解耦成果合规可入库严禁一体化七参数、仅依靠 EGM2008 批量换算高程。4.2 临时外业踏勘、土方粗略测算、DEM 批量出图、奥维地图看图高程转换采用 EGM2008 模型纯城区平地可统一增加 0.10m 抵消全球零点偏差包含山地丘陵区域禁止全局固定修正高程偏差无法控制。4.3 历史数据复盘、科研大地测量计算使用 XGM2019e、EIGEN‑6C4 等高阶全球重力场模型不用于工程施工。第五章 核心总结要点椭球高≠海拔二者差值为高程异常 ζ所有高程转换核心是求解 ζ全球重力场模型存在固定零点偏差 空间渐变非线性误差无法用单一常数完全修正阶数仅决定重力细节分辨能力工程高程校正必须遵循平面转换、高程拟合完全解耦原则一体化七参数外推误差爆炸仅能临时粗略参考TGO 本质是一阶平面高程拟合适合平原二次曲面适配山地但不可远距离外推coorblh.dat 存储 WGS84/CGCS2000 椭球经纬度与大地高是基准转换、高程拟合的标准输入文件国内省级精化似大地水准面精度全面碾压所有公开全球 EGM/EIGEN 模型是正式测绘成果首选高程基准数据。