GPS高程拟合方法
GPS高程拟合方法在风场测量中的应用
在 一局 部 ＧＰ Ｓ网 中 ， 由若 干个 点 的 Ｈ舭值 作
为 已知值 ， 数值 拟 合 方 法 内插 出其 它 ＧＰ 用 Ｓ测 点
的高程 异常 ， 式 （） 求 得 各 点 的 正 常 高 。若 测 按 １可 区内测 量点 的数 量足 够 多 ， 分 布 较 为均 匀 ， 可 且 则
收稿 日期 ：０ ９ １ — ８ ２ ０ — ００ ２ １ ． ／ 球定 位 系统 ００ １全
程 常采 用正 （ ） 常 高高 程 系统 。其 正常 高 （ 常） Ｈ正 是
以似大地 水 准面起 算 的 ， 大地水 准 面差 距［ ］ 叫 ３ 。 “ 地 面上 一点 Ａ 的大 地高 Ｈ 与正 常高 Ｈ正 常的
关 系 为 ＿ ２ ］
Ｈ ， Ｈ 知 一 一 （） １
设 Ｇ Ｓ基 线 向量 网 经二 维 约束 平 差 后 ， 得 Ｐ 求 各点 的空 间 直 角 坐 标 ｚ 、 ， 大 地 高 平 差 值 Ｈ Ｙ 和 ，
摘 要 ： 风 力 发 电 场 场 址 测 量 中 ， Ｓ 定 位 测 量 技 术 可 以 得 到 高 精 度 的 中 ， 面点 的 高程 常采 用 的是 正 常 高 高程 系统 ， 但 地 因此 ， 需将 Ｇ Ｓ大 地 高 Ｐ
转 换 成 正 常 高 。 针 对 面积 较 大 或 地 形 变 化 比 较 复 杂 的 ＧＰ 测 区 ， 活 运 用 ＧＰ 高 程 拟 合 的 Ｓ 灵 Ｓ
已知 网 中有 个 点 （ 中总点 数 为 ｍ） 有 正 常高 网 具
Ｈ 。则 可确 定这 些点 的高 程异 常 ８ 为Ｌ ］ ｌ
ｉ ｆ（ — ｚ，３）＋￡ ， （） ３
或 一 Ｈｍ— Ｈ， 地 水准 面之 间的 高程差 。
（） ２
浅述GPS高程拟合的几种方法
浅述GPS高程拟合的几种方法当前我们测量中的高程系是相对于选定的某一参考面而定的，基准面有参考椭球面，大地水准面和似大地水准面，而在实际测量中，由于地球形状的不规则性，以及地球内部重力分布的不均匀性，想要得到严密的数学转换关系式是很难以实现的，高程拟合即是实现精化区域似大地水准面的一种方法，本文浅述几种高程拟合的常用方法。
标签：高程系;高程异常;GPS大地高;高程拟合;神经网络法1、高程系统1.1常见的高程系统通常应用的高程系统，主要有大地高程系统、正常高系统和正高系统。
大地高程系统是以椭球面为基准面的高程系统，由地面点沿通过该点的椭球面法线到椭球面的距离，通常以H表示。
大地高是一个几何量，它不具有物理上的意义。
利用GPS定位技术，可以直接测定观测站在WGS-84或ITRF中的大地高。
以大地水准面为基准面的高程系统，称为正高系统。
由地面点，并沿该点的铅垂线至大地水准面的距离，称为正高，通常以Hg表示。
正高实际上是无法严格确定的;正常高是指从一地面点沿过此点的正常重力线到似大地水准面的距离，似大地水准面严格说不是水准面，但接近于水准面，只是用于计算的辅助面，它与大地水准面不完全吻合，差值为正常高与正高之差。
正常高系统为我国通用的高程系统。
大地水准面与似大地水准面在海平面上是重合的，而在陆地上则既不重合也不平行。
1.2高程系统之间的关系设大地高为H，正高为Hg，正常高为Hγ，参考椭球面与大地水准面之间的差距为大地水准面差距N，参考椭球面与似大地水准面之间的差距为高程异常ξ，那么上述的3种高程系统之间存在的关系：H=Hg+N=Hγ+ξ2. GPS高程拟合原理实现方法2.1 GPS高程拟合原理由于大地水准面与椭球面一般不重合，我们把地面点P沿铅垂线投影到大地水准面P0时，P与P0间距离为正高Hg;在将点P0沿法线方向投影到椭球面上得点Q0，P0与Q0间距离称为大地水准面差距N，H=Hg+N。
似大地水准面与椭球面也不重合，它们之间的高程差称为高程异常，用ζ表示。
基于EGM2008和Kriging的GPS高程拟合算法
第３ ６卷 第 １ ２期
２ ０ １ ３年 ｌ ２月
测 绘 与 空 间地 理 信 息
Ｇ ＥＯ ＭＡＴ Ｉ Ｃ Ｓ＆ Ｓ Ｐ ＡＴ Ｉ Ａ Ｌ Ｉ ＮＦＯ ＲＭＡＴ Ｉ Ｏ Ｎ Ｔ Ｅ Ｃ ＨＮＯＬ ＯＧＹ
Ｖｏ １ ． ３６， Ｎｏ ．１ ２ Ｄｅ ｃ ．，２ ０１ ３
基于 Ｅ ＧＭ２ ０ ０ ８和 Ｋｒ ｉ ｇ ｉ ｎ ｇ的 ＧＰ Ｓ高程 拟 合算 法
文章编号 ： １ ６ ７ ２ — ５ ８ ６ ７ （ ２ ０ １ ３ ） １ ２— ０ １ ９ ６— ０ ３
ＧＰＳ Ｅｌ ｅ ｖ ａ ｔ ｉ ｏ ｎ Ｆｉ ｔ ｔ ｉ ｎｇ Ｍｅ ｔ ｈ ｏ ｄ Ｂａ ｓ ｅ ｄ ｏ ｎ Ｋｒ ｉ ｇ ｉ ｎ ｇ ａ ｎ ｄ ＥＧＭ ２ ０ ０ ８
Ｆ Ｅ ＮＧ Ｙａ ｎ—ｓ ｈ ｕ ｎ，Ｗ ＥＩ Ｌ ｅ ｉ ，Ｚ ＨＡＮＧ Ｙｕ ｅ—ｆ ｅ ｉ ，ＫＡＮＧ Ｚ ｏ ｎ ｇ—ｄ ａ ｏ
（ Ｈｅ ｎ ａ ｎ Ｐ ｒ ｏ ｖ ｉ ｎ ｃ ｉ ａ ｌ Ｇｅ ｏ ｌ ｏ ｇ ｉ ｃ ａ ｌ Ｓ ｕ ｒ ｖ ｅ ｙ ｉ ｎ ｇ ａ ｎ ｄ Ｍａ ｐ ｐ ｉ ｎ ｇ Ｉ ｎ ｓ ｉ ｔ ｔ ｕ ｔ ｅ ， Ｚ ｈ ｅ ｎ ｇ ｚ ｈ ｏ ｕ ４ ５ ０ ０ ０ ６ ， Ｃ ｈ ｉ ｎ ａ ）
提高GPS拟合高程精度的方法
由于似大地水准面不规则， 很难选取一种好的拟合方 量 和转 换方法 的关 系 ， 笔者选择 四种方 案进行 高程 约束 。 法与其相吻。二次曲面拟合对控制点 的数量要求 ７ 个以 取 网中范 围一半 的 ３ 个水准点进行 高程拟合 ； 网中 上， 取 且拟合效果不一定好。 因此， 建议首先采用平面拟合，首级 控 制采 用 Ｇ Ｓ静 态测 量 方 法 布设 上 ， Ｐ 最好三等。 的 Ｄ级 Ｇ Ｓ控 制 ， 制 面积 １８ ｋ 。高 程采 用 四等水 Ｐ 控 ８０ ｍ ，
准测量 布设 ， 联测 了多数 Ｄ级 Ｇ Ｓ控 制点 。 Ｐ Ｄ级 Ｇ Ｓ Ｐ 控 范 围要 包 含测 区范 围 ， 可 能避免 从 一端 向另一 端无 限 尽
黑江土源０． 龙国资 ２０ １ １日
５ 测 区地 形起 伏 大 小有 关 与
地形起 伏 小 ， 似大 地水 准 面较 规则 ， 拟合 精度 高 ； 则
于一 端分 布 的情 况 。方 案 ３的精 度 优于 方案 ２的精 度 ， 地 形 起伏 大 ， 似大 地水 准 面极不 规则 ， 拟合 精 度差 。 阶 二 显然是 由起 算点 的数 量不 同造 成 的 ， 即起 算点 的数量 越 多项式函数模型的几何特征为一 “ 抛物单曲面”然而一 ， 多精度 越 高 。方 案 ４的精 度 最高 ， 因为方 案 ４将 测 区划 个测 区 的似大 地水 准面 并不 会 如此 规整 ， 了使 二 阶 多 为 分为 ３块进 行 高 程 拟合 ，能够 较 好 的接 近 似大 地 水 准 项式函数与似大地水准面拟合得更好 ， 有必要为二阶多 面， 因此其 精度 最高 。 因此 ， 当测 区范 围较 大 时 ， 分 区 项 式 函数定 义 一个 区域 ， 应 即选 择 整体 拟合 和分 区拟 合 。 求 取校 正参 数才 能保证 高 程 的转换 质量 。除方 案 ４外 ，
基于GPS的高程拟合方法研究
基于GPS的高程拟合方法研究高程拟合是基于GPS数据进行地表高程估计的一种方法。
在现代测量和导航技术中，GPS被广泛应用于三维空间定位和高程测量。
由于GPS观测数据存在误差和不确定性，导致从GPS数据直接估计高程时存在一定的误差。
需要进行高程拟合来提高高程估计的精度和可靠性。
高程拟合的基本原理是通过建立GPS观测数据与地表高程之间的数学模型，利用最小二乘法等数学方法来拟合观测数据，得到地表高程的估计值。
常用的高程拟合方法包括平差法、插值法和卡尔曼滤波法等。
平差法是一种常用的高程拟合方法，主要通过将GPS观测数据与已知高程点进行权衡，利用最小二乘法来调整观测数据的权值，从而得到更精确的高程估计值。
平差法的优点是简单易行，适用于大部分高程拟合问题。
平差法的缺点是需要预先获取一定数量的已知高程点，如果没有足够的已知高程点，拟合结果可能较差。
卡尔曼滤波法是一种基于滤波理论的高程拟合方法，主要通过建立动态状态模型和观测方程来估计地表高程，利用卡尔曼滤波算法来对GPS观测数据进行滤波和优化。
卡尔曼滤波法的优点是能够考虑观测数据的权值和误差，能够在有限的观测数据中提供更精确的高程估计值。
卡尔曼滤波法的缺点是需要预先获取一定数量的已知高程点，对初始状态的选取敏感。
除了以上方法，还可以结合其他辅助数据进行高程拟合。
可以利用DEM（Digital Elevation Model）数据作为辅助数据，通过比较GPS观测数据和DEM数据的差异，来估计地表高程。
还可以利用地形特征等辅助信息，通过建立地表高程的统计模型来进行高程拟合。
城镇控制测量中GPS高程拟合方法分析
第 ２期
赵文更等 ： 城镇控制测量 中 Ｇ Ｓ Ｐ 高程拟合方法分析
１３ ７
的多层 馈 网络 ， 由输 入 层 、 出层 和 一层 或 多层 隐含 层 组 输
从上 可 以看 出人工神 经 网络 在 Ｇ Ｓ高程 拟合 方 面有 Ｐ
成 ， 之 间无反 馈连 接 ， 各层 各层 内神经 元之 间无任 何 连接 ， 仅相邻 层神经元 之间有连 接 ， 如图 １ 所示 。
析 ， 过 多次 试 算 得 到 了较 好 的结 果 。高 程 异 常 拟合 精 通 度 及运 行次数 如 图 ２所示 ， 检核 点 的拟合 残差 结果 如 表 １ 所 示 ， 用学 习集 ７个点 求 解 了二 次 曲面拟 合参 数 ， 并 对工
Ｂ Ｐ网络 的工作 方 式分 为 两 个 阶段 ： 个 阶段 为 学 习 一 阶段 ， 过 对所 选 学 习样本 的学 习 ， 断调 整各 神 经元 之 通 不 间的 连接 权值 ， 学 习结果 十 分逼 近 样本 值 ； 一 阶段 为 使 另
Ｘ１ Ｘ１
自己独特的优势 ， 研究基于人工神 经网络 的方法来转换
ＧＳ Ｐ 高程能够减少人为构建数学模型所带来的误差 ， 因此
应 能提高 Ｇ Ｓ 程转 换 的精 度 。 Ｐ高
４ 基于 Ｂ Ｐ学 习算法的高程拟合分析
Ｘ２
Ｘ２
本文以扶余城镇控制测量 中的 Ｄ， Ｅ级 Ｇ Ｓ网总共 Ｐ
合、 多面函数法 、 曲面样条拟合 、 非参数 回归法 、 移动法等
方法 。但 现在 最 常用 的仍 是 二 次 多项 式 摸 拟二 次 曲面来 完成 高程 的拟 合 。 本文 所讨 论 的 是 基 于 神 经 网络 法 进 行 Ｇ Ｓ高 程 拟 Ｐ 合 ， 是一种 自适 应 的 映射 方 法 ， 不 用 作 假 设 ， 论 上 它 它 理 也 比较合 理 ， 避免 未 知 因 素 的影 响 ， 少 模 型 误差 。这 能 减
三种GPS高程拟合计算方法
GP S 高程 拟合 精 度主 要 取决 于 GP S 大 地 高 精 度 、 重 合 点 正 常 高 的 精 度 、重 合 点 的分 布情况和高程拟合 模型的选择。
3 结语 实验结果显示 , 三种高程拟合计 算方
例介 绍了三种 实现 GP S 高程拟 合的计算方 法, 分析了三种 方曲面拟合
中 图分类 号: P 2
文 献 标识 码 : A
文章编号: 1 6 72 - 3 79 1 ( 2 00 8) 1 0( a) - 0 0 48 - 0 1
全球定 位系统( Gl obal Pos i t i o ni ng Sys t e m- GPS ) 作为 新一 代的 卫星 导航 与定 位系 统, 以 其全 球 性、 全天 候、 高精 度、 高效 益 的显著特点, 已经在测量 领域得到了广泛 的应 用[1]。
在此次实验中已 知大地高和正常高数 据 为 9 个 , 而 求 解 二 次 多项 式 曲 面 拟 合 参 数 需要 6 个 或 6 个 以 上, 当 然 了, 已 知 数据 越多求解的结果越精确, 本次实验主要是 为了验证三种方法的可行 性和正确性, 为 了 能 充 分 问题 , 我们 取 8 个 已 知 数 据 参 与 运算 , 1 6 6 号点 作 为检 核数 据。 测区 控制 点 分 布图 如 图 1 所 示。
科技资讯 2008 NO. 28 SCI ENCE & TECHNOLOGY I NF ORMATI ON 三 种 GP S 高 程 拟 合 计 算 方 法
工程技术
吴伟 盛 李胜利 ( 福建省八闽测 绘院 厦门 3 6 1 0 1 2 )
GPS高程拟合的方法及实现
（ １ ） 当前， Ｇ Ｐ Ｓ 平 面控制网 已经得到 了广泛 的运 用， 但 是Ｇ Ｐ Ｓ 用与电力基础 工业。 在 目前 推行 的电力地理信息管理系统 中，
高程却运 用得不够 ， 人们期望着 能够用Ｇ Ｐ Ｓ 高程测量代替 传统 自动化数据 的传 输需建立在各基点的准确定位上 ， 目前成熟稳 的水 准测量 。 本文对Ｇ Ｐ Ｓ 高程测 量的原理和方法 进行 了初步 的 定 的Ｇ Ｐ Ｓ 技术成为最主要的技术选择 。 设备 的状态监控 随着Ｇ Ｐ Ｓ
探 讨， 并结合 我 国Ｇ Ｐ Ｓ 高程测 量 的应 用的实 际，以数值拟 合为 地 理信息采集 技术 的应用， 已经摆脱了单纯的站 内或馈 线采集 线 路状态 的监控 从点已经推 进到线并逐渐普及 主， 建 立了高程转换 的数 学模 型。 同时用ｖ Ｃ ＋ ＋ 开发了Ｇ Ｐ Ｓ 高程控 箱采集 的模 式， 制 转换系 统， 经试 验测试 ， 在平原或 浅丘地 区 ， 在不加地 形改 到 面。 Ｇ Ｐ Ｓ 相关 技术 的研发在未 来的 电力应 用领域必将 得到广
正 的情况 下， 拟 合出的正常高 高程满足一般 工程和大 比例尺测 泛 的应 用 。 图的精度 要求 ， 在一定程 度上 降低了生产成本 。
（ ２ ） 数值拟合 的数 学模型很 多， 考虑到模 型 的通用性 ， 实用 （ ２ ） 一年里， 北斗系统发 射了四箭六星 ， 完成了系统组网， 同 时采取了一系列的措施提 高系统精度和 稳定性 ， 扩展了系统覆
性 以及计算实现 的方便性 ， 本拟合转换软件详细叙述 了四种常 盖 区域 ， 和去年试运 行相 比， 覆 盖区域 由原有 的东经８ ４ 度一 １ ６ ０ 用 的模 型： 对 加权平均值拟合和 多项式 曲面拟合详 细介绍 ， 另 度扩展 到现 在的东经５ ５ 度一 １ ８ ０ 度， 系统 定位精度 由去年发布 的 外提 到插值 拟合和多面 函数 拟合方法 。 同时， 还 考虑了利 用非 ７ ｋ － ￣２ ５ 米、 高程３ ０ 米到现在 的７ ｋ ￣ Ｚ ｌ Ｏ 米、 高程 １ ０ 米， 测速精度 由 格 网化数据进行 地形改正 的几何 方法。 当测区形状 为带状 时， 可 以采用前二种方法进行 计算， 当测 区太长 时 （ 超过 ｌ Ｏ Ｏ ｋ ｍ ） ， ０ ． ４ 米每秒到０ ． ２ 米每秒。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
GPS高程拟合方法3.1等值线图示法等值线图示法是最直接的求算高程异常的方法。
这种方法的核心思想就是内插的思想，绘制高程异常的等值线图，然后采用内插法来确定未知点的高程异常值。
具体操作十分的简单，在测区内制定分布均匀的GPS点，用水准测量的方法来测定这些点的水准高，根据公式ζ=H-Hr求出这些点的高程异常，选择适当的比例尺按照已知点的平面坐标展会在图纸内，对已知点标注出高程异常值，再确定等高距，绘制出高程异常值的等值线图。
之后就可以内插出待测点的高程异常值，进而求出待测点的正常高。
这种方法只适用地形相对平坦的地方，在此种测区内采用这种方法拟合的高程精度可达到厘米级。
测区的地形相对复杂内插出的高程异常值就不准确，而且这种内插法的精度往往取决于两个方面，分别是测区内GPS点的分布密度和已知点大地高的精确度。
首先GPS点的分布比较密集，那么内插精度就相对较高，如果比较稀疏这时候就要借助于此测区的重力测量资料，提高内插精度。
且还要注意GPS点间高程异常的非线性变化。
另外就是水准点的精度，联测时尽量选取高精度的正常高，尽可能使得出的高程异常值准确，进而才能内插出待测点高精度的高程异常值。
这种方法虽然简单易操作，但是有其弱点，就是精度不高，只有当对拟合精度要求不高的时候才使用此种方法(注：等值线法不需构造数学模型)。
3.2狭长带状区域线性拟合解析内插法作为拟合高程最常用的方法，主要思想是把似大地水准面用数学曲面近似拟合，建立所在测区内最为接近似大地水准面的数学模型，以此来计算测区内任意点的高程异常值，从而计算出正常高。
这种方法计算出的高程异常值的精度是由所采用的数学模型和似大地水准面的拟合程度所决定的。
解析内插法在选择数学模型时，首先要考虑的就是GPS点的分布情况。
GPS 点的分布情况可分为带状分布和面状分布。
若GPS点是呈线状布设，而且是以沿线似大地水准面为一条连续且光滑的曲线，这时就可以采用相对于狭长带状区域的解析内插法来内插出待定点的高程异常值，从而求出待定点的正常高。
这种线状分布的内插原理是：测区内已知水准点，用GPS测出其GPS高程，计算出已知水准点高程异常值，根据已知点的平面坐标和计算得出的高程异常值，构造出一个插值函数，这个函数是用来拟合GPS分布线上的似大地水准面的。
用这个函数内插出位置点的高程异常值。
下面是两种用来拟合线状分布的GPS高程的内插法。
3.2.1 多项式曲线拟合法多项式曲线拟合是线状分布拟合的主要方法。
多项式拟合顾名思义其插值函数是一个m次的代数多项式，若高程控制点的高程异常为ζ，坐标为x l(或y i或d i 或拟合坐标或x i-x0或y i-y0)的函数关系为下式：ζl=a0+a1x l1+a2x l2+a3x l3+……a m x l m3-1各高程控制点的已知高程异常与其拟合值之差为下式所示：r i=ζl(x)-ζi(i=0，1，2…n)3-2上式我们称之为离差。
(3-1)中x l是拟合点到参考点(x0，y0)的直线距离，x0，y0为设定的常数值。
在一般情况下都认为，x0，y0就是测区内已知点坐标的均值。
多项式曲线拟合使用起来非常方便，但是它有自身的局限性，即是使用这种方法的时候，所测路线不能太长，要限制控制点到测点的距离不能太远，通常把距离控制在300米以内。
这个要求是因为使用多项式曲线方法拟合似大地水准面，如果它拟合的范围太大，点位的高程异常变化就越复杂，削高补低的方法不能满足我们所要求的精度。
随着多项式阶数的增大，也会使拟合出的曲线振荡的更厉害，从而造成拟合的误差增大。
这些造成了多项式曲线拟合的缺陷，但是在路线较短的情况下，这种方法有足够的精度来拟合GPS点的正常高程。
在式(3-1)中用m次多项式拟合似大地水准面，这个m的值如何取定，一般情况下如果测区不是很长，地形相对平坦，那么我们通常取m取为3。
也就是说多项式为三次多项式。
若测区比较长或者是测区地形比较复杂就要依情况而定，增加多项式的次数，提高拟合精度。
依上述分析m的取值主要和测区长度以及测区的复杂程度有关。
3.2.2三次样条曲线拟合法三次样条曲线拟合法针对测线长，已知点多的测区GPS高程拟合问题。
由上述可以知道，当测线比较长已知点较多的时候，就需要构造高次的拟合多项式，当m值比较高的情况下，会出现不稳定的现象，对求解高程异常值会有比较大的影响，并且最小二乘法在求多项式系数中也会增大削高补低的误差，因此为了避免测线长、已知点多这种情况下所出现的问题，通常采用分段拟合的方法，采用三次样条函数拟合数学模型。
这种方法很好的解决了因测线长而引起的问题。
三次样条曲线的实质就是一个拼接而成的连续函数，在把测线分为多段的情况下，每段设为三次多项式函数，然后将这些多项式函数组成三次样条函数。
为了计算准确，应用中要求这种构造的曲线不仅在连接点处函数要连续。
而且还要求这个函数的一级导数还有二级导数全部要是连续的，才能保证在分段之后构造的三次样条函数后期运算中能够计算出准确的高程异常值。
设过n个已知点，ζi和x i(或y i或拟合坐标)在区间[x i,x i+1](i=1,2,…,n-1)上有三次样条函数关系:ζ(x)=ζ(x i)+(x-x i)ζ(x i,x i+1)+(x-x i)(x− x i+1)ζ(x,x i,x i+1)3-3式中，x为待定点坐标,ζ(x i,x i+1)为一阶差商,ζ(x i,x i+1)=(ζi+1−ζi)/(x i+1-x i)；ζ(x，x i,x i+1)为二阶差商,ζ(x,x i,x i+1)=1[ζ,,(x i)+ζ,,(x)+ζ,,(x i+1),而6ζ,,(x)(i=1,2,…,n-1)，满足系数矩阵为对称三角阵的线性方程组(x i-x i+1)ζ,,(x i−1)+2(x i+1,x i−1)ζ,,(x i)+(x i+1,x i)ζ,,(x i+1)=6[ζ,,(x,x i+1)-ζ(x,x i)]3-4ζ(x0)=ζ,,(x n)=0用追赶法解上面方程组，可求出ζ,,(x i)和ζ(x i,x i+1)，而ζ,,(x)=ζ,,(x i)+(x-x i)ζ,,(x i,x i+1)3-5这种做法有诸多好处，其中优点有三点：其一计算简便，其二保留了多项式的优点，其三克服了多项式的缺点。
多项式的缺点是单个多项式会有不灵活不稳定的现象。
由于三次样条曲线的种种优点，往往在实际中当遇到测线长已知点多的情况下采用此方法拟合高程。
3.3曲面拟合法曲面拟合法是用于GPS点的分布在一定区域的时候，且可以选择数学曲面拟合该区域的似大地水准面，构造适当的数学模型，计算该区域内的高程异常值，然后求出正常高。
这种拟合法的主体思想和曲线拟合法异曲同工的。
具体思想是：已知测区的若干已知水准点，并且用GPS测定这些点的高程，利用公式求得这些点的高程异常，有了已知点的高程异常，已知点的平面坐标是已知的，所以利用其平面坐标(x，y)和高程异常值ζ构造出来的数学模型拟合最为接近于该测区的似大地水准面，然后内插出未知点的高程异常值ζ，进而求出正常高。
3.3.1多项式曲面拟合法测站点的大地高H与正常高h之间有如下关系：h=H-ζ3-6多项式函数拟合法的基本思想是在小区域GPS网内，将似大地水准面看成曲面(或平面)，将高程异常表示为平面坐标(x，y)的函数，通过网中起算点(既进行了GPS测量又进行了几何水准联测的点)已知的高程异常确定测区的似大地水准面形状，求出其余各点的高程异常，然后根据式(3-6)求出其他点的正常高，其数学模型为：ζ=f(x，y)+ε3-7 式中f(x，y)是拟合的似大地水准面；ε是拟合误差，f(x，y)=a0+a1x +a2y+a3x2+a4xy+a5y2+…3-8x=B-B0y=L-L0B0=1n ∑B L0=1n∑L其中：n为GPS网中点的数量，(B,L)为已知点的大地坐标，a0，a1，a2，a3，a4，a5……为拟合待定参数；x，y为各GPS点的平面坐标的近似值，一般取起算点的平面坐标减去网中全部点平面坐标的均值。
〔1〕二次曲面拟合取(3-8)式中的一、二次项后将大地水准面拟合为：f(x，y)=a0+a1x+a2y+a3x2+a4xy+a5y23-9 即得二次曲面拟合模型：ζ=[a0a1a2a3a4a5][ 1 x y x2 x y y2 ] T+ε3-10每一个起算点可以组成一个上式，若共存在m个这样的起算点，则可列出m个方程：ζ1=a0+a1x1+a2y1+a3x12+a4x1y1+a5y11+ε1ζ2=a0+a1x2+a2y2+a3x22+a4x2y2+a5y21+ε2……ζm=a0+a1x m+a2y m+a3x m2+a4x m y m+a5y m1+εm3-11从而组成误差方程：V=-Bx+L3-12上式中，B=[1 x1 y1x12y12 x1y11 x2y2x22 y22 x2y2…1 x m y m x m2 y m2x m y m]，x=[a0 a1 a2 a3 a4 a5]T，L=[ζ1 ζ2 …ζm]T，解得x=(B T PB)−1B T PL3-13解算出a i即可求出网中其余点的高程异常，并利用式(3-6)求出各未知点的正常高h。
〔2〕多项式平面拟合在小范围或平原地区，可以认为大地水准面趋近于平面。
此时，可选用公式（3-8）的前三项，将大地水准面拟合为：f(x，y)=a0+a1x+a2y3-14 拟合模型为：ζ=[a0a1a2][ 1 x y ] T+ε3-15其中，a i(i=0,1,2)为未知参数，此时要求公共点至少3个。
〔3〕多项式相关平面拟合：也叫做四参数曲面拟合，若选用公式(3-8)的前三项和第五项进行拟合，则拟合曲面的表达式变为：f(x，y)=a0+a1x+a2y+a3xy 3-16拟合模型为：ζ=[a0a1a2a5][ 1 x y xy] T+ε3-17 其中，a i(i=0,1,2,3)为未知参数，此时需要公共点至少4个。
3.3.2移动曲面拟合法移动曲面拟合法是一种局部逼近法，其基本思想是以每一个内差点为中心，利用内差点周围数据点的值，建立一个拟合曲面，使其到各个数据点的距离之加权平方和为极小，而这个曲面在内插点上的值就是所求的内插值。
设P为内插的点，下面对P构造相应的曲面。
本文取如下的二次多项式曲面为例：f(x，y)=a0+a1x +a2y+a3x2+a4xy+a5y23-18设选取数据点的坐标为(x i,y i)，i=1,2,…,n；n≥6且设内插点P的坐标为(x p,y p)。
将(x i,y i)改化到以P为原点的局部坐标系中，即：x i,=x i-x py i,=y i-y p3-19形成新的坐标( x i,，y i,)，为移动坐标。