多面函数法高程拟合应用研究
GPS高程拟合中一种多面函数平滑因子求法的研究
统， 而我国采用的高程系统是正常高系统 ； 因此 ， 在
我国使用 Ｇ Ｐ Ｓ 高程数据 的时候就 要进行高程 的转 化， 从而使大地高转换为正常高。似大地水准 面至 参考椭球 面的距 离称为高程异 常＿ ２ Ｊ 。在 知道大地 高的情况下 ， 如 能准确知道该点 的高程异常 ， 就能 准确知道该点的正常高 。在实际应用 中， 通过数学 的方法 拟合 出各个 点的高 程异 常 ， 从 而得 到正 常
⑥
２ ０ １ ３ Ｓ ｃ ｉ ． Ｔ ｅ ｃ ｈ ． Ｅ ｎ ｇ ｒ ｇ ．
地球科学
ＧＰ Ｓ高程拟 合中一 种多面函数平滑 因子 求 法 的研 究
黄 继磊 罗志清
（ 昆明理工 大学 国土资 源工程学 院， 昆明 ６ ５ ０ ０ ９ ３ ）
摘
要
通 过对 多面函数正双 曲面模型 中的平滑 因子进 行分 析， 认 为平 滑 因子 的取值有 一个 最佳 区间。取值 在最佳 区间内
：
常的平均值 。 的值越大 ， 表示测区地形起伏较 大，
一
Ｉ Ｌ ｉ ‘ ． Ｉ ｏ （ ｘ ， Ｙ ； 『 ｍ ， 『 ｍ ） … Ｏ （ ｘ ， Ｙ ； ， 加 ） ｊ
当选 取 ｍ个 已测点 的 ｎ个点 为 结 点时 ， 根 据最
∑ ∑ｄ
则 的数学表面 的总和以任意精度逼 近的思想。该
方 法 自提 出以来 ， 在 测 绘 界 做 了大 量 的应 用 Ｌ ６ ］ 。该 方 法 的数学 模 型为 ：
关系。因此 ， 能否恰 当的选取平滑因子对拟合精度 有 很 大影 响 Ｊ 。马 洪 滨 等 人 提 出 了多 面 函数 平
基于遗传算法的多面函数GPS高程拟合模型_徐南
［1 ］ ［2 ］ ［3 ］ ［4 ］ ［5 ］
考
文
献
卢银宏，岳东杰，刘毅等． 线性内插法在 GPS 高程拟合中 ． 勘察科学技术，2012 ，（ 5 ） ： 52 ～ 57． 的应用 ［J］ J］． 测绘通报， 韩硕． 神经网络在 GPS 高程拟合中的应用 ［ 2006 ，（ 4 ） ： 48 ～ 50． 焦殿阳，张旭晴． 曲面拟合模型在小区域 GPS 高程拟合中 ． 测绘与空间地理信息，2012 ，（ 3 ） ： 24 ～ 26． 的应用 ［J］ 李祖锋，巨天力，张成增等． 基于重力场模型高程拟合残 差求定 GPS 正常高 ［J］． 测绘工程，2010 ，（ 4 ） ： 4 ～ 6． 刘万林，王利，赵超英． GPS 水准的有限元法与多面函数 法的加 权 综 合 模 型 ［J］． 地 球 科 学 与 环 境 学 报， 2004 ， （ 26 ） ： 48 ～ 51．
2013 年第 11 期
k
工程勘察
Geotechnical Investigation ＆ Surveying
51
y） = ξ（ x，
y； x i ， yi ） α i Q（ x， Σ i =1
（ 1）
其中，α i 为 核 函 数 加 权 系 数； Q （ x， y； x i ， y i ） 为核函数； （ x， y ） 为未测点坐标； （ x i ， y i ） 为已 知点坐标。 常用的二次核函数有正双曲面函数和倒双曲面 函数，本文选取正双曲面函数作为核函数 ： Q（ x， y； x i ， yi ） = （x 槡 － xi ）
52
工程勘察
Geotechnical Investigation ＆ Surveying
2013 年第 11 期
基于智能算法的多面函数在高程拟合中的应用
基于智能算法的多面函数在高程拟合中的应用
贾思楠;吴风华;范祺
【期刊名称】《华北理工大学学报:自然科学版》
【年(卷),期】2022(44)4
【摘要】针对多面函数在GPS高程拟合中存在中心节点和光滑因子选取困难的问题,提出采用智能算法改进多面函数中核函数的方法。
首先利用蚁群智能算法取代传统均匀格网法,选取最优中心节点;然后将核函数中光滑因子作为种群染色体,利用遗传算法搜寻染色体的最优值,通过将2种算法融合构建高精度的拟合模型。
研究结果表明,相比于传统方法,该智能算法在优化多面函数上残差波动幅度更小、拟合精度更高,可为地势波动较大的特殊地形的高程拟合提供理论依据。
【总页数】7页(P8-13)
【作者】贾思楠;吴风华;范祺
【作者单位】华北理工大学矿业工程学院;华北理工大学机械工程学院
【正文语种】中文
【中图分类】P228
【相关文献】
1.多面函数法在GPS高程拟合中的应用
2.二次多项式拟合法和多面函数拟合法在GNSS高程测量中的应用
3.基于蚁群算法的多面函数在GPS高程拟合中的应用
4.多面函数在GPS高程拟合中的应用
5.多面函数参数自适应选取方法在GPS高程拟合中的应用
因版权原因，仅展示原文概要，查看原文内容请购买。
基于蚁群算法的多面函数在GPS高程拟合中的应用
基于蚁群算法的多面函数在GPS高程拟合中的应用蒲伦;唐诗华;张紫萍;胡新凯;肖燕【摘要】针对多面函数拟合法的中心节点难以选取的问题,提出将蚁群算法引入多面函数并结合稳健估计构建高精度拟合模型的方法.利用蚁群算法在复杂地形中快速寻找特征点,与少量非特征点共同作为中心节点参与模型构建,将稳健估计加入多面函数,运用选权迭代法剔除粗差对拟合模型的影响.GPS高程拟合数据处理实例表明,基于蚁群算法的多面函数结合稳健估计的拟合方法可有效剔除粗差的影响,且拟合精度比只用均匀格网法提高26％.【期刊名称】《大地测量与地球动力学》【年(卷),期】2019(039)001【总页数】5页(P31-35)【关键词】高程异常;蚁群算法;稳健估计;多面函数;中心节点【作者】蒲伦;唐诗华;张紫萍;胡新凯;肖燕【作者单位】桂林理工大学测绘地理信息学院,桂林市雁山街319号,541006;广西空间信息与测绘重点实验室,桂林市雁山街319号,541006;桂林理工大学测绘地理信息学院,桂林市雁山街319号,541006;广西空间信息与测绘重点实验室,桂林市雁山街319号,541006;青海省生态环境遥感监测中心,西宁市南山东路116号,810007;桂林理工大学测绘地理信息学院,桂林市雁山街319号,541006;广西空间信息与测绘重点实验室,桂林市雁山街319号,541006;桂林理工大学测绘地理信息学院,桂林市雁山街319号,541006;广西空间信息与测绘重点实验室,桂林市雁山街319号,541006【正文语种】中文【中图分类】P228随着GPS技术在测绘领域的发展和完善，采用部分GPS水准重合点实现高程拟合越来越受到测绘工作者的青睐[1]。
获取高精度的高程异常是高程拟合的核心内容[2]。
对于地势起伏平缓的研究区域，常规高程拟合方法能够满足实际拟合精度的要求。
但当地势起伏较大、地形结构较复杂时，拟合结果与实际高程异常往往存在较大的偏差。
多面拟合函数法转换GPS高程
摘要传统的几何水准测量方法, 是测绘领域中测定正常高的主要方法。
这种方法虽然精度高, 但实施起来费时费力, 作业效率较低。
GPS 定位技术自问世以来, 就以其精度高, 速度快, 操作简单等优势对传统的水准测量造成了极大冲击，但是现实测绘和项目建设上需要的具有物理意义的正常高（或正高）高差，因此，如何用GPS高求定正常高，是当前测绘工作者较为热心的课题之一。
关键词：GPS水准；大地高：正常高Abstract:Traditional geometric leveling method, a major means to measure the normal height in surveying, produces high precision but can take too much time and efforts with low efficiency. GPS technology has hit hard the traditional leveling ways via its high precision, fast speed and simple operation since it came to the world. However, real leveling and project construction require a normal altitude difference with physical significance, so, how to identify the normal height with the GPS height is one of the preoccupations for surveying workers.Key words:GPS leveling;geoidheight; normal height目录摘要 (I)1 引言...................................................................................................................................................... I I2 高程系统.............................................................................................................................................. I I2.1大地高系统.................................................................................................................................. I I2.2 正高系统.................................................................................................................................... I II2.3 正常高系统................................................................................................................................ I II2.4 高程系统之间的转换关系 (IV)3 GPS高程测量的基本原理 (IV)3.1 物理大地测量方法 (IV)3.2 几何方法 (IV)4 GPS大地高转化为正常高的方法讨论 (V)4.1 平面拟合 (V)4.2 多项式曲线拟合 (V)4.3 二次多项式曲面拟合法 (V)3.4 多面函数曲面拟合 (V)4.5 加权平均值法 (VI)3.6 非参数回归法和高程异常变化梯度法 (VII)4.7 移动曲面法 (VII)4.8 固定边界3 次样条插值法 (VIII)4.9 非格网GPS 散点数据考虑地形改正法 (VIII)4.10 BP神经网络算法................................................................................................................... I X5 算法实例........................................................................................................................................... I X5.1 参与计算的数据........................................................................................................................ I X5.2平面拟合法................................................................................................................................. X I5.3二次曲面法拟合......................................................................................................................... X I5.4加权平均值法............................................................................................................................ X II5.5多面函数曲面拟合.................................................................................................................... X II6 结论.................................................................................................................................................... X II参考文献......................................................................................................................................... X III附录1.平面拟合程序： (XIV)附录2加权平均值法拟合程序： (XVI)附录3 二次多项式曲面拟合程序 (XVII)附录4多面函数曲面拟合程序 (XVIII)致谢 (XVIII)1引言我国目前采用的高程系统为正常高系统:即以似大地水准面为参考面的高程系统,确定高程通常采用的是几何水准。
几种高程拟合方法的精度分析
总结与分析
二
对于多面函数法拟合 核函数的选取，以及 平滑 因子都会对精度 产生较大影响一定要 多次尝试。
三
已知点的选取 一定要均匀， 并非已知点越多 精度就越高。
LOGO
实例分析
❖ 当核函数为锥面函数时C取1，当核函数为到双曲面时 2 取10000，以下是
这三种拟合模型的残差图。
实例分析
❖ 当选取1、7、9、10、11、13、19、20、22、26、28、33、34、36这 14个点作为已知点进行二次曲面拟合时其精度如下表：
序号
2 3 4 5 6 8 12 14 15
测量等级 三等几何水准测量 四等几何水准测量 普通几何水准测量
水准限差
允许的最大限差（mm）
12 L 20 L 30 L
注：L为已知点与检核点的距离（单位：公里）
实例分析
❖ 右图为某中型城 市的城市控制网，图 中共有37个GPS—E 级控制点。为了研究 GPS拟合原理，对以 上所有控制点都进行 了三等水准测量，并 应用稳健估计进行粗 差探测，未发现粗差。
二次曲面拟合残差（14点）
-0.050 0.002 -0.062 0.016 0.061 0.018 0.135 0.011 -0.055
二次曲面拟合残差（20点）
-0.031 -0.047 -0.140 0.052 0.080 0.028 0.002 0.061 0.037
总结
一
对于地势比较平坦， 或者高程变化比较 平缓时，二次曲面 法拟合可以满足其 精度要求
❖ 为了保证试验数 据的可靠性，其具体 数据见下表。
序号 1 2 3 4 5 6 7 8 9 10 11
X坐标 -9230.899 -10589.011 -8775.220 -7666.317 -11649.851 -8129.317 -4334.088 295.107 3802.651 -11790.336 -7892.980
多面函数法高程拟合应用研究
多 面 函 数 法 高 程 拟 合 应 用 研 究
王 世 君
摘 要： 探讨 了多面 函数算法原理， 实现 了多面 函数法高程拟合软件编制 ， 通过 实例分析 了其拟合精度及关键问题， 与 并
常 用 的二 次 曲面 函 数 法进 行 了对 比分 析 ， 出 了多 面 函 数 法 应用 时应 注 意 的 事项 ， 供 参考 。 指 以 关键 词 ： 多面 函数 ， 高程 拟 合 ， 序设 计 程 中 图分 类 号 ： Ｕ１８ Ｔ ９ 文献标识码 ： Ａ
其 中， 为待定参 数 ； ｚ， ， ， ） ｚ和 的二次核●函 Ｆ（ 为 ］● ， ●●●
数 ， 用 的 二 次 核 函数 为 ： 常
２ 光 滑 因子 的 确定 。没有 理 想 的方 法 ， 要 不 断 的试 验 改进 。 ） 需
３核心点 的选 择。这种方 法的核 心点要 求是高程 异常显 著 ） （） 点 ， ３ 即能很好地描述 该 区域 内高程 异常分 布的特 征点 ， 最好位 于
多 面 函数 法 从 几 何 观 点 出 发 ， 决 根 据 数据 点 形 成 ～ 个 平 差 高 。 解 多面函数 法高程拟合时要解决好 以下几个问题 ： 的数学 曲面 问题 。此法的基本思想是 ， 任何数学表 面和不规则 的 １核 函数 的选 取 。 由于 多 面 函数 法 基 于纯 数 学 的逼 近 理 论 ， ） 圆滑 表 面 ， 总可 以用 一 系 列 有 规 则 的 数 学 表 面 的 总 和 ， 任 意 精 以
一
在最小二乘的原则下 ， 可求得 ：
＝
（ Ｑ） Ｑ 一 Ｑ
定 的科 学 价 值 。
把求得 的系数 回代到式 （ ） 即可推 算 出任 一未知点 的 ２ 中，
基于EGM2008的多面函数高程拟合方法研究
第４１卷第１０期２０１８年１０月测绘与空间地理信息ＧＥＯＭＡＴＩＣＳ＆ＳＰＡＴＩＡＬＩＮＦＯＲＭＡＴＩＯＮＴＥＣＨＮＯＬＯＧＹＶｏｌ.４１ꎬＮｏ.１０Ｏｃｔ.ꎬ２０１８收稿日期:２０１８－０４－０８作者简介:杨沛旺(１９９２－)ꎬ男ꎬ河北保定人ꎬ助理工程师ꎬ学士ꎬ主要从事测绘㊁地理信息系统应用开发等方面的生产工作ꎮ基于ＥＧＭ２００８的多面函数高程拟合方法研究杨沛旺ꎬ徐㊀萌ꎬ罗㊀梦ꎬ李孝玲(黑龙江地理信息工程院ꎬ黑龙江哈尔滨１５００８１)摘要:随着全球定位系统ＧＰＳ高程测量技术的发展和ＥＧＭ(ＥａｒｔｈＧｒａｖｉｔａｔｉｏｎａｌＭｏｄｅｌ)模型精度以及分辨率的不断提高ꎬ基于少量水准点实现区域ＧＰＳ大地高转正常高成为可能ꎮ本文提出一种基于ＥＧＭ２００８模型ꎬ利用多面函数法来实现高程拟合的方法ꎮ结合咸阳北部某区域实验数据ꎬ采用上述方法完成区域内高程拟合ꎬ结果表明ꎬ该方法符合精度达到３ｃｍꎬ满足大比例尺地形图测图精度要求ꎮ关键词:ＥＧＭ２００８ꎻ多面函数法ꎻ高程拟合中图分类号:Ｐ２２８.４㊀㊀㊀文献标识码:Ａ㊀㊀㊀文章编号:１６７２－５８６７(２０１８)１０－０１５３－０３ＲｅｓｅａｒｃｈｏｎＭｕｌｔｉ－ＦｕｎｃｔｉｏｎＥｌｅｖａｔｉｏｎＦｉｔｔｉｎｇＭｅｔｈｏｄＢａｓｅｄｏｎＥＧＭ２００８ＹＡＮＧＰｅｉｗａｎｇꎬＸＵＭｅｎｇꎬＬＵＯＭｅｎｇꎬＬＩＸｉａｏｌｉｎｇ(ＨｅｉｌｏｎｇｊｉａｎｇＩｎｓｔｉｔｕｔｅｏｆＧｅｏｍａｔｉｃｓＥｎｇｉｎｅｅｒｉｎｇꎬＨａｒｂｉｎ１５００８１ꎬＣｈｉｎａ)Ａｂｓｔｒａｃｔ:ＷｉｔｈｔｈｅｄｅｖｅｌｏｐｍｅｎｔｏｆｔｈｅｇｌｏｂａｌｐｏｓｉｔｉｏｎｉｎｇｓｙｓｔｅｍＧＰＳｈｅｉｇｈｔｍｅａｓｕｒｅｍｅｎｔｔｅｃｈｎｏｌｏｇｙａｎｄｔｈｅｃｏｎｔｉｎｕｏｕｓｉｍｐｒｏｖｅｍｅｎｔｏｆｔｈｅａｃｃｕｒａｃｙａｎｄｒｅｓｏｌｕｔｉｏｎｏｆｔｈｅＥＧＭ(ＥａｒｔｈＧｒａｖｉｔａｔｉｏｎａｌＭｏｄｅｌ)ｍｏｄｅｌꎬｉｔｉｓｐｏｓｓｉｂｌｅｔｏａｃｈｉｅｖｅｇｅｏｄｅｔｉｃｈｅｉｇｈｔｔｕｒｎｔｏｎｏｒｍａｌｈｅｉｇｈｔｂａｓｅｄｏｎａｓｍａｌｌａｍｏｕｎｔｏｆｂｅｎｃｈｍａｒｋ.ＴｈｉｓｐａｐｅｒｐｒｏｐｏｓｅｓａｍｅｔｈｏｄｂａｓｅｄｏｎｔｈｅＥＧＭ２００８ｍｏｄｅｌｔｏａｃｈｉｅｖｅｈｅｉｇｈｔｆｉｔｔｉｎｇｕ￣ｓｉｎｇｔｈｅｐｏｌｙｈｅｄｒａｌｆｕｎｃｔｉｏｎｍｅｔｈｏｄ.ＣｏｍｂｉｎｅｄｗｉｔｈｔｈｅｅｘｐｅｒｉｍｅｎｔａｌｄａｔａｏｆａｃｅｒｔａｉｎａｒｅａｉｎｎｏｒｔｈｅｒｎＸｉａｎｙａｎｇꎬｔｈｅａｂｏｖｅｍｅｔｈｏｄｗａｓｕｓｅｄｔｏｃｏｍｐｌｅｔｅｔｈｅｒｅｇｉｏｎａｌｈｅｉｇｈｔｆｉｔｔｉｎｇ.Ｔｈｅｒｅｓｕｌｔｓｓｈｏｗｅｄｔｈａｔｔｈｅｅｘｔｅｒｎａｌａｃｃｕｒａｃｙｏｆｔｈｅｍｅｔｈｏｄｒｅａｃｈｅｄ３ｃｍꎬｗｈｉｃｈｍｅｅｔｓｔｈｅａｃｃｕｒａｃｙｒｅｑｕｉｒｅｍｅｎｔｓｆｏｒｌａｒｇｅ－ｓｃａｌｅｔｏｐｏｇｒａｐｈｉｃｍａｐｍａｐｐｉｎｇ.Ｋｅｙｗｏｒｄｓ:ＥＧＭ２００８ꎻｐｏｌｙｈｅｄｒａｌｆｕｎｃｔｉｏｎꎻｈｅｉｇｈｔｆｉｔｔｉｎｇ０㊀引㊀言众所周知ꎬ在以往的无精化水准面区域ＧＰＳ高程测量中ꎬ通常的做法是联测一定数量㊁分布合理均匀且能有效覆盖整个控制区域的ＧＰＳ水准点ꎬ然后拟合出测区的似大地水准面ꎬ进而求得ＧＰＳ网中所有待测点的正常高[１]ꎮ这种方法在日常生产实践中应用较多ꎬ在地形平坦地区理想情况下拟合精度已经可以达到四等水准测量甚至更高的精度要求[２]ꎮ但此作业过程效率低ꎬ费时费力ꎮ随着ＧＰＳ高程测量所获得的大地高精度不断提高ꎬ以及ＥＧＭ模型精度和分辨率的进一步提高ꎬ现已可以实现利用少量ＧＰＳ/水准高程点进行大地高至正常高的转换ꎮＧＰＳ直接测得的高程是相对ＷＧＳ－８４椭球面的大地高㊁ＥＧＭ模型是基于ＷＧＳ－８４椭球面的大地水准面模型ꎬ两者都与我国采用的正常高系统不一致[３]ꎮ将ＧＰＳ大地高转换为工程项目所需正常高的过程ꎬ关键在于高程异常的求解ꎮＥＧＭ２００８模型能够精确表达高程异常的长波项ꎬ而短波项主要由地形变化所引起ꎮ本文利用移去 拟合 恢复 法的思想ꎬ提出一种利用ＥＧＭ２００８重力场模型获得高程异常长波项ꎬ利用多面函数法拟合因地形起伏引起的高程异常短波项的方法来实现ＧＰＳ高程拟合ꎬ并通过咸阳北部实验数据对上述方法可行性进行验证ꎮ１㊀原理与方法１.１㊀原理与流程本文利用移去 拟合 恢复法进行高程异常求解ꎬ综合利用ＥＧＭ２００８重力场模型和多面函数拟合的方法ꎬ仅须联测少量水准点ꎬ即能高效地拟合区域内高程异常ꎮ根据物理大地测量学理论ꎬ高程异常ζ可以表示为[４]:ζ＝ζＥＧＭ＋ζｒｅｓ(１)式中ꎬζＥＧＭ为由ＥＧＭ２００８重力场模型求得的高程异常长波部分㊁ζｒｅｓ为高程异常残余部分ꎮ移去 拟合 恢复 的思想首先将ζＥＧＭ从高程异常ζ中移去ꎬ获得ζｒｅｓꎻ并利用少量水准点进行拟合ꎬ本文选用多面函数法完成ꎬ得到未知点的ζｒｅｓꎻ再将ζＥＧＭ恢复即可得到待求点高程异常ζꎮ１.２㊀ＥＧＭ２００８重力场模型计算高程异常长波项ＥＧＭ２００８是由美国国家地理空间情报局(ＮａｔｉｏｎａｌＧｅｏｓｐａｔｉａｌ－ＩｎｔｅｌｌｉｇｅｎｃｅＡｇｅｎｃｙꎬＮＧＡ)于２００８年推出的最新一代地球重力场模型ꎬ模型采用的基本空间分辨率为５ᶄˑ５ᶄꎬ最高空间分辨率达到１ᶄˑ１ᶄꎬ数据的全球地面覆盖率超过８０％[５]ꎮＥＧＭ２００８模型是以ＰＧＭ２００７Ｂ为参考ꎬ利用ＧＲＡＣＥ卫星重力数据㊁ＴＯＰＥＸ卫星测高数据㊁全球５ᶄˑ５ᶄ重力异常数据㊁地面重力数据以及地形数据ꎬ该模型的阶次完全至２１５９ꎬ相当于模型空间分辨率为９ｋｍꎮ在我国大陆ꎬＥＧＭ２００８模型高程异常的总体精度为２０ｃｍꎬ其中ꎬ华北地区可以达到９ｃｍꎬ华东华中地区为１２ｃｍꎬ西部地区为２４ｃｍꎮ该模型采用Ｂｒｕｎｓ公式计算高程异常:ζＥＧＭ＝ＧＭργð２１９０ｎ＝２(ａρ)ｎðｎｍ＝０[Ｃ－ｍｎｃｏｓ(ｍλ)＋Ｓ－ｍｎｓｉｎ(ｍλ)]Ｐ－ｍｎ(ｓｉｎφ)(２)式中ꎬＧＭ为地心引力常数ꎻρ㊁γ㊁λ和φ分别为所求点地心向径㊁正常重力值㊁经度和余纬ꎻａ为椭球长半径ꎻＣ－ｍｎ与Ｓ－ｍｎ为完全规格化位系数ꎻＰ－ｍｎ(ｓｉｎφ)为完全规格化缔合函数ꎮ目前ꎬＥＧＭ２００８模型的高程异常可通过应用软件 ＡｌｌｔｒａｎｓＥＧＭ２００８Ｃａｌｃｕｌａｔｏｒ 直接计算得到ꎮ１.３㊀多面函数法高程拟合利用已有的水准点计算得到移去长波项后的高程异常残余部分ζｒｅｓ＝ζ－ζＥＧＭꎬ将其作为已知数据ꎬ并通过多面函数法拟合残余高程异常模型ꎮ多面函数法由几何角度出发ꎬ能够解决根据数据点形成的一个平差数学曲面问题ꎮ此法核心思想是任一数学表面或不规则圆滑表面ꎬ均可用多个有规则的数学表面相加得到ꎮ在笛卡儿坐标系中ꎬ多面函数的一般式为:ζ(ｘꎬｙ)＝ｆ(ｘꎬｙ)＝ðｎｉ＝１βｉＦ(ｘꎬｙꎬｘｉꎬｙｉ)(３)式中ꎬβｉ为待求系数ꎻＦ(ｘꎬｙꎬｘｉꎬｙｉ)为ｘ和ｙ的二次核函数ꎬ其中心点在(ｘｉꎬｙｉ)处ꎬζ(ｘꎬｙ)可以表示为二次式的和ꎬ因此称为多面函数ꎮ二次核函数的选取可以是任意的ꎬ但通常为了方便计算ꎬ都采用具有对称性的距离型ꎮ常用的二次核函数为:Ｆ(ｘꎬｙꎬｘｉꎬｙｉ)＝[(ｘ－ｘｉ)２＋(ｙ－ｙｉ)２＋δ２]ｋ(４)式中ꎬδ２称为光滑因子ꎬ可以为任意正数ꎻ核函数随ｋ的取值变化而变化ꎬ当ｋ＝１/２时ꎬ核函数为正双曲面函数ꎻ当ｋ＝－１/２时ꎬ为倒双曲面函数ꎮ设有ｎ个已知点ꎬ选取其中ｍ个点(ｍ£ｎ)作为核心点ꎬ并令Ｑｉｊ＝Ｆ(ｘꎬｙꎬｘｉꎬｙｉ)ꎬ则式(３)变为:ζ(ｘꎬｙ)＝ðｎｉ＝１βｉＱｉｊ(５)则误差方程式为:ｖ１ｖ２⋮ｖｎæèççççöø÷÷÷÷＝Ｑ１１Ｑ２１Ｑ１ｍＱ２ｍ⋮⋱⋮Ｑｎ１ Ｑｍｎæèççççöø÷÷÷÷β１β２⋮βｎæèççççöø÷÷÷÷－ζ１ζ２⋮ζｎæèççççöø÷÷÷÷(６)表达成向量的形式为:ｖ＝Ｑβ－ζ(７)根据最小二乘原则ꎬ可求得:β＝(ＱＴＱ)－１ＱＴζ(８)求得β后带回(３)式中ꎬ即可得到多面函数拟合模型ꎬ并根据待求点ＧＰＳ坐标计算其高程异常残余项ꎮ由于多面函数拟合法为数学中的逼近理论ꎬ因此核函数的选取会直接影响拟合精度ꎮ此外ꎬ多面函数模型拟合精度还受核心点和光滑因子的制约ꎮ光滑因子的优劣直接作用于拟合效果ꎬ目前ꎬ对于光滑因子的确定尚无有效方法ꎬ需多次尝试ꎮ而核心点过少或精度不足时ꎬ将会为模型带来一定的系统误差ꎬ但当核心点达到一定数量时ꎬ继续增加核心点既无法提高模型精度ꎬ又费时费力ꎮ核心点的选择则需考虑区域地形特征ꎬ即能较好地描述目标区域内高程异常分布特征的位置ꎬ如最高点㊁最低点和坡度变化处ꎮ１.４㊀恢复长波项计算待求点高程异常利用拟合得到的多面函数模型计算待求点高程异常残余值ζｒｅｓꎬ并利用式(２)求解该点高程异常长波项ζＥＧＭꎬ带入式(１)即可得到待求点高程异常ꎬ进而得到该点正常高ꎮ２ 算例分析本文利用咸阳北部实测数据进行上述高程拟合方法实用性进行验证ꎮ目标区域约２１５ｋｍ２ꎬ该区域内无精化水准面ꎮ使用ＧＰＳ接收机观测和水准测量方法获得了２４个已知点ꎮ２.１㊀计算流程与数据处理由于利用Ｂｒｕｎｓ公式直接计算重力场模型高程异常比较复杂ꎬ本次采用 ＡｌｌｔｒａｎｓＥＧＭ２００８Ｃａｌｃｕｌａｔｏｒ 软件直接计算高程异常长波项ꎬ其中ꎬ重力场模型选择 Ｕｎｄ＿ｍｉｎ１ｘ１＿ｅｇｍ２００８＿ｉｓｗ＝８２＿ＷＧＳ８４＿ＴｉｄｅＦｒｅｅ＿ＳＥ ꎮ然后移去长波项ꎬ得到各已知点高程异常残余部分ꎮ经多次试验ꎬ本次拟合核函数选取正双曲面函数(ｋ＝１/２)ꎬ即:Ｆ(ｘꎬｙꎬｘｉꎬｙｉ)＝[(ｘ－ｘｉ)２＋(ｙ－ｙｉ)２＋δ２]１/２(５)多面函数高程拟合采用Ｍａｔｌａｂ编程实现ꎬ通过已知点不断尝试ꎬ选定最优光滑因子和核心点ꎬ得到最终的拟４５１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀测绘与空间地理信息㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０１８年合模型并计算待求点正常高ꎮ程序算法流程如图１所示ꎮ图１㊀多面函数拟合算法流程图Ｆｉｇ.１㊀Ｐｏｌｙｈｅｄｒａｌｆｕｎｃｔｉｏｎｆｉｔｔｉｎｇａｌｇｏｒｉｔｈｍｒｏａｄｍａｐ经过多次试算ꎬ最终得出当δ＝０.００６５时精度最高ꎬ并选定均匀分布于整个测区的１２个已知点作为核心点参与多面函数模型待求系数的解算ꎬ其余１２个已知点作为检查点ꎬ并通过计算检查点真实高程与拟合高程之差来评价模型精度ꎬ核心点和检查点分布情况如图２所示ꎮ图２㊀核心点和检查点分布图Ｆｉｇ.２㊀Ｄｉｓｔｒｉｂｕｔｉｏｎｏｆｃｏｒｅｐｏｉｎｔａｎｄｃｈｅｃｋｐｏｉｎｔ２.２㊀精度评价多面函数方法进行模型求解内符合中误差为０ꎬ外符合精度以检查点水准高程和拟合高程之差进行评价ꎮ水准高程与多面函数拟合结果对照见表１ꎮ表１㊀多面函数拟合结果Ｔａｂ.１㊀Ｒｅｓｕｌｔｏｆｐｏｌｙｈｅｄｒａｌｆｕｎｃｔｉｏｎ点号水准高程(ｍ)拟合高程(ｍ)拟合误差(ｍ)Ａ００４５２８.７４１５２８.７１１０.０３０Ａ００５５０６.１０４５０６.１４８－０.０４４Ａ００８５０５.８０９５０５.８４１－０.０３２Ａ０１０５０５.５３２５０５.５２６０.００６Ａ０１２５０９.０７０５０９.０５２０.０１８Ａ０１３４７９.４５６４７９.５０３－０.０４７Ｂ００２４９８.１４９４９８.１２０.０２９Ｂ００６４９２.８２７４９２.８２５０.００２Ｂ００７４１３.６９４４１３.６５８０.０３６Ｂ００８４０７.４８２４０７.４７１０.０１１Ｂ００９４８２.０１１４８２.００１０.０１０Ｂ０１０４１０.０５３４１０.０９２－０.０３９由表１可以看出ꎬ多面函数拟合结果与检查点水准高程的较差均小于５ｃｍꎬ拟合高程外符合中误差为３ｃｍꎬ满足大比例尺地形图测图精度要求ꎮ３㊀结束语采用 移去 拟合 恢复 法的思想ꎬ本文提出基于ＥＧＭ２００８重力场模型确定高程异常长波项ꎬ并通过多面函数法拟合高程异常残余部分的方法完成大地高至正常高的转化ꎮ实例表明ꎬ一定范围内联测少量水准点作为已知点ꎬ通过上述方法可以实现区域内大地高至正常高的转化ꎬ拟合成果可以应用于大比例尺测图中ꎮ参考文献:[１]㊀冯林刚ꎬ杨润甫ꎬ李胜.基于ＥＧＭ９６的ＧＰＳ高程转换方法[Ｊ].测绘通报ꎬ２００６(３):２２－２３ꎬ４７.[２]㊀廖中石ꎬ龚立新.基于移除 恢复法的ＧＰＳ水准高差拟合方法研究[Ｊ].北京测绘ꎬ２００８(４):４１－４２ꎬ５２.[３]㊀倪明ꎬ文加林ꎬ丁仁军.基于ＥＧＭ２００８模型在ＧＰＳ拟合高程中的应用[Ｊ].城市勘测ꎬ２０１４(４):９９－１００ꎬ１０８.[４]㊀朱亚光ꎬ高兴国ꎬ刘焱雄ꎬ等.基于ＥＧＭ２００８重力场模型的高程联测方法[Ｊ].济南大学学报:自然科学版ꎬ２０１１ꎬ２５(４):４１０－４１３.[５]㊀王建文.利用ＥＧＭ２００８模型进行区域高程异常改正的实现[Ｊ].测绘与空间地理信息ꎬ２０１６ꎬ３９(３):１８７－１８８ꎬ１９２.[编辑:任亚茹]５５１第１０期杨沛旺等:基于ＥＧＭ２００８的多面函数高程拟合方法研究。
多面函数在GPS高程拟合中的应用
多面函数在GPS高程拟合中的应用韩子清;郭杨亮;马瑞娟;张顺幸【摘要】GPS高程拟合的关键是利用少量GPS/水准重合点,实现GPS大地高与正常高的转换,函数拟合法在实际工程中应用广泛.本文对GPS高程拟合常用方法多面函数法和二次曲面法的原理进行介绍,结合实例用两种方法进行转换,对结果进行比较,证实多面函数法可以提高拟合精度,适用于高程异常较为复杂的地区.【期刊名称】《四川有色金属》【年(卷),期】2019(000)001【总页数】4页(P5-7,14)【关键词】高程异常;多面函数;二次曲面【作者】韩子清;郭杨亮;马瑞娟;张顺幸【作者单位】西安科技大学测绘科技与技术学院, 陕西西安 710054;河南省煤田地质局物探测量队,河南郑州 450009;河南省地质矿产勘查开发局第五地质勘查院, 河南郑州 450001;河南省地质矿产勘查开发局第五地质勘查院, 河南郑州 450001【正文语种】中文【中图分类】P228.41 国内外研究现状GPS技术可快速获取地面点三维坐标，在实际工作中应用广泛。
GPS可获取一点处大地高，而我国采用的是正常高，二者的起算面不同。
因此，将GPS获取的大地高转换为正常高是GPS高程应用的关键。
正常高与大地高的转换如公式（1）所示。
式中：H正－该点处的正常高；H大－该点处的大地高；ζ－该点处的高程异常；由公式（1）可知，一点处所测得的大地高减去该点处的高程异常，便可求得该点处的正常高。
因此，GPS高程异常的确定是实现GPS大地高与正常高转换的核心。
确定GPS高程异常主要有以下几种方法：联合平差法、数值逼近法、等值线图法、重力法等。
出于对拟合精度与拟合方法的实用性、易操作性考虑，数值逼近法在实际工程中应用广泛。
数值逼近法通过数学表达式拟合出一个曲面，该方法在测区规律性变化明显的地区可获得较好拟合结果。
主要方法有：多面函数拟合、多项式曲面拟合、克里格法、加权平均法、最小二乘配置、滤波与推估等。
用于GPS高程拟合的多面函数模型的应用研究_贺炳彦
3 多面函数模型与单一函数模型对比
将多面函数法应用于 GPS 网中高程异常的拟合, 拟合范围内已知点 20个 , 数据通过函数 E= 3 . 125+ 0 . 054x - 0. 254y + 0 . 016xy 模拟, 其位置如图 2 , 其 1~ 20 号点点位坐标分别为 x ( 5- 25 ), y ( 2- 8 ) 范围的 网格点坐标。其拟模拟数据如表 1 .
0 引
言
[ 1]
在传统大地测量中, 水准测量作为获取高程的最常规手段 , 具有原理简单、 误差易于检验等优点。但 是长距离水准测量的实施, 存在着劳动强度大 , 实施进度缓慢、 易产生人为误差等缺点 达到三四等水准的精度
[ 2]
。而作为现代大
地测量手段之一 , GPS 以其精度高 , 实施简便等优点, 已经得到广泛的应用, 其高程精度, 在平坦地区能够 。因此通过高程拟合的方式 , 将 GPS 测量的大地高转化为正常高 , 成为现今高 程测量的重要方法之一。
[ 4]
2 多面函数的基本方程
多面函数拟合曲面的方法是美国 H ardy 教授 1977 年提出的 , 其理论
基础是 , 任何一个圆滑的数学曲面总可以用一系列有规则的数学表面的综合, 以任意精度逼近。 GPS 高 程多面函数拟合法就是把拟合区域的高程异常 , 用多个曲面高度逼近 , 建立数学模型, 借此可以求得未知 点的高程异常, 然后根据 GPS所求的大地高来计算常规基准下的正常高。
T - 1 [ 5]
为
2 2
F ( x, y ) = a0 + a1 x + a2 y + a 3x + a4 y + a5 xy. X = [ a 0, a1, a2, a 3, a4, a5 ] B L 其中 B 为系数矩阵 , L 为已知高程异常
几种高程拟合方法的精度分析
序号 1 2 3 4 5 6 7 8 9 10 11
X坐标 -9230.899 -10589.011 -8775.220 -7666.317 -11649.851 -8129.317 -4334.088 295.107 3802.651 -11790.336 -7892.980
(x, y)
u
(x, y) a jQ(x, y, xi , yi ) j 1
aj
Q(x, y, xi , yi )
V AQ 0
A (QT Q)1QT 0
1
Q [( x xi )2 ( y yi )2 ]2
Q
[(
x
xi
)2
(
y
yi
)2
1
]2
(xi , yi )
xi , yi
大地高
12.097 11.931
11.467 11.579 13.403 12.377 11.381 11.250 11.366
正常高
4.173 3.925
3.421 3.884 5.835 4.852 3.906 2.822 3.332
高程异常
7.925 8.005
8.046 7.695 7.569 7.525 7.475 8.428 8.034
拟合值
残差值
二次曲面 锥面 倒双曲面 二次 锥面 倒双曲面
2 8.580
3 8.673
4 8.866
5 8.161
6 8.210
8 8.225
12 7.899
14
8.046
15
7.695
外符合精度
8.549 8.536 8.624 -0.031 -0.043 0.045 8.626 8.759 8.773 -0.047 0.086 0.100 8.726 9.559 8.380 -0.140 0.693 -0.487 8.213 8.048 8.198 0.052 -0.114 0.037 8.290 8.177 8.279 0.080 -0.033 0.069 8.253 8.174 8.359 0.028 -0.051 0.134 7.901 8.551 7.568 0.002 0.627 -0.357
多面函数法在GPS高程拟合中的应用
多面函数法在GPS高程拟合中的应用作者：刘晓明邸彦彬来源：《活力》2012年第04期[摘要]为了改善GPS大地高向正常高转换的精度，在局部区域内,建立多面函数模型进行高程拟合,可以达到较高的精度。
文中利用多面函数模型进行高程拟合,除选取分布均匀的GPS 水准联测点外,还对核函数形式的选取做了详细地分析，并与高程拟合中常用的二次曲面拟合法做了对比，进行了精度分析。
[关键词]GPS高程；多面函数法；核函数目前国内外应用GPS定位技术建立各类控制网时，平面控制基准的精度毋庸置疑，高程定位精度仍需进一步研究和提高。
因此如何有效利用GPS测量的高程信息把大地高转换成正常高，直接为测绘生产服务，实现GPS观测时能同时获得实用的三维坐标，是非常实际而有意义地。
对于局部区域工程控制网,可以采用数学方法建立合理的高程拟合模型,高程拟合的结果既可以验证水准测量的正确性,也可以作为精度要求较低的高程控制基准使用。
本文对多面函数方法进行了详细地研究，利用某工程控制网点的高程异常进行拟合，并与高程拟合中常用的二次曲面拟合法做了对比，进行了精度分析。
多面函数拟合法，1971年由美国哈笛(Hardy)提出。
1976年将此法应用于美国大地测量、拟合重力异常、大地水准面差距、垂线偏差等，1978年将此法用于地壳形变。
它的基本思想是任何一个规则或不规则的连续曲面均可以由若干简单面(或称单值数学面)来叠加逼近。
具体做法是在每个数据点上建立一个曲面，然后在方向上将各个旋转曲面按一定比例叠加成一张整体的连续曲面，使之严格地通过各个数据点[1][2]。
多面叠加的数学表达式为：这里Q(x,y,xi,yi)为参加插值计算的简单数学面，又称为多面函数的核函数；n为简单数学面的张数或多层叠加面的层数，它的值与分块扩充范围内参与点的个数相等；Ki(i=1,2,3,…，n)为待定参数，它代表了第i个核函数对多层叠加面的贡献。
为了计算方便，多层叠加面中的个核函数一般选用同一类型的简单函数，通常是围绕竖向轴旋转的曲面，这条竖正好通过某一参考点，例如：这里为非零参数。
多面函数在GPS高程拟合中的应用
当 ｍ ＝“时 ： ａ＝Ｑ１ ｆ （） ５
（） ６
的性质 。最小二乘配置法 中的协方差 函数是一 种统计 函数 ， 在高 程异 常资料稀少 的地 区很难确定 ， 而多 面函数的核 函数可 以按几 何关系确定 ， 它是距 离 的函数 ， 且顾 及 了待定点 和 已知点 间的相
当 ｍ＞“时 ， ： 则
Ｖ ＝Ｑ 一
标数据 中直接输 入或 从 坐标 文件 调入 也 可在 道路 数 据 中设 置。 实际水平距 离 ， ＨＲ 为对 准放样 点 尚差 的水 平 距离 ， ｄ ｄ 当 ＨＲ＝
由于有 已知的两个控制点 ， 因此采取 在坐标数 据菜单里直 接输入 ０ ０ ０表示距离正确 ， Ｆ ］ 示坐 标数 据 ， ．０ 按［ ３ 显 最后 按下 ［ ４ 就 进 Ｆ］
其中 ， 为待 定 系 数 ； ａ， ， Ｙ ） ａ Ｑ（ ．Ｙ 薯， ｉ 为核 函 数 ， 中 点在 Ｔ 其 墨 Ｙ） ｘ， （ Ｙ） ， 在实 际应 用 中， 但 需 （ ，ｉ； Ｙ均为待定点 的坐标 ；墨 ，ｉ为已知点 坐标 。核 函数 但 一般 都采 用具 有对称 性的 要 的是相对于似大地水准面 的正常高 Ｈ，两者之间的关系为 ¨ ： 在理 论上可以任意选 取 ， 为 了简便 ， ， １ Ｊ 可选择 ： ＝Ｈｒ —Ｈ８ ４ （ ） 距离型。例如 ， １
示放样点 坐标 ， ［ ３ 显 示放 样 点的距 离 和方 向角 。按 ［ １ 进 按 Ｆ］ Ｆ］ 水平 角 ， ｄ 当 ＨＲ＝００ ０时表 示方 向正确 ， ［ １ 进 入 ， 为 。 ０ ０ 按 Ｆ］ ＨＤ
入， ＨＲ为实 际测量水 平角 ， ＨＲ为 对准放样 点转 动放样 转动 的 ［ ］ 吴贵 才． ｄ ３ 全站仪 在 建 筑 工程放 样 中的应 用 ［ ］ 山西 建筑 ， Ｊ．
多面函数模型中核函数个数对GPS高程拟合精度影响分析
第33卷第6期2010年11月现 代 测 绘M odern Surveying and MappingVol.33,No.6Nov.2010多面函数模型中核函数个数对GPS高程拟合精度影响分析袁 豹,丁 旭(河海大学测绘科学与工程系,江苏南京210098)摘 要 在局部地区G PS高程拟合中,选择不同的拟合模型和方法会对G PS高程拟合的精度产生很大的影响。
本文通过对多面函数拟合模型中不同核函数中心点个数选择对G PS高程拟合精度影响的研究,通过实例计算分析得到结论:核函数中心点因个数选取的不同,对GP S高程拟合精度的影响较大,在采用多面函数模型拟合GP S 高程时,应该考虑核函数中心点个数的选择问题。
关键词 多面函数 高程拟合 核函数中心点 精度分析中图分类号:P228.4 文献标识码:A 文章编号:1672-4097(2010)06-0006-031 引 言水准测量是测定地面点高程(正常高)的主要方法之一[1],它以测量精度高、原理简单等优点广泛应用于工程及现代测绘的各个方面。
但是,大部分水准测量都具有劳动强度大以及外业施测缓慢、艰难等不足。
随着全球定位系统(GPS)技术的出现,为高程测量提供了一个崭新的技术方法。
GPS以其全球性、高效性、经济性等显著特点被广泛应用于测绘生产中。
目前利用GPS技术建立的短基线平面控制网相对精度可达10-8-10-9,垂直分量绝对精度可达毫米级,但由于GPS获取的是相对于参考椭球的大地高,而工程应用中通常采用正常高系统,因此需要对二者进行相应的转换[2]。
近些年来,很多学者对GPS高程拟合做了大量的研究[3 4],在拟合方法和拟合模型的选取方面取得了很大的进展,其中多面函数拟合模型就是一个较为新颖的方法。
在利用多面函数拟合高程时,核函数以及平滑因子的选取对GPS高程拟合的精度起着关键性的影响,有关平滑因子和多面函数的选择对GPS高程拟合精度的影响已有相关研究论文发表[5 7],而对于核函数中心点个数与GPS高程拟合精度的关系,还未有相关研究成果。
GPS高程拟合中多面函数及二次曲面函数的比较与分析
GPS高程拟合中多面函数及二次曲面函数的比较与分析摘要：结合某区域的gps和水准测量数据，用多面函数法、二次曲面函数法分别进行高程拟合，结果表明，在高程异常变化较大的区域用多面函数拟合gps高程异常时，精度较高。
关键词：gps 测量；水准测量；高程拟合；多面函数；二次曲面函数引言随着卫星定位技术的发展，gps以其独特的优势在测量领域扮演着越来越重的角色。
众所周知使用gps测量得出的点，其平面精度相当高，但是高程精度不是很高。
对于工程建设而言，保证高程控制基准的准确是必须的。
那么有没有一种处理方法使得用gps测量得出来的高程达到我们的要求呢？下面我们将结合一些实例来探究这个问题。
我们知道gps测量实在wgs-84地心坐标系中进行的，所提供的高程为相对于wgs-84椭球的大地高。
而我们要得到是海拔高。
在某一区域内，如果有一定数量的已知水准点（正常高已知），则可以在这些水准点上进行gps观测，各点的高程异常值就可根据计算得出。
根据已知点的高程异常值及其位置关系建立函数模型来模拟该区域似大地水准面的高度，再用数学内插的方法求解区域内任一点的高程异常值，进而求得该点的正常高。
本研究应用多面函数方法对某控制网gps点的高程异常进行拟合，并与水准高程进行了对比分析。
原理多面函数法是从几何观点出发，解决根据数据点形成一个平差的数学曲面问题。
其理论根据是认为“任何一个圆滑的数学表面总可以用一系列规则的数学表面总和以任意的精度逼近”。
二次曲面拟合法是认为高程异常在一定范围内变化平缓的前提下，将高程异常近似地看作是一定范围内各点坐标的曲面函数，用这一拟合函数来计算其它gps点的高程异常和正常高。
对于每一个已知点都可以列出方程，从而求出检核点的正常高。
实际运用中，如果把测区的似大地水准面假定为平面，则为平面拟合模型，要求测区面积很小且地形十分平坦，计算出来的高程异常与检核点的正常高，精度一般不高。
如果把测区的似大地水准面看成一个二次曲面，则相对符合对似大地水准面的描述。
基于EGM2008的多面函数高程拟合方法研究
基于EGM2008的多面函数高程拟合方法研究
杨沛旺;徐萌;罗梦;李孝玲
【期刊名称】《测绘与空间地理信息》
【年(卷),期】2018(041)010
【摘要】随着全球定位系统GPS高程测量技术的发展和EGM(Earth Gravitational Model)模型精度以及分辨率的不断提高,基于少量水准点实现区域GPS大地高转正常高成为可能.本文提出一种基于EGM2008模型,利用多面函数法来实现高程拟合的方法.结合咸阳北部某区域实验数据,采用上述方法完成区域内高程拟合,结果表明,该方法符合精度达到3 cm,满足大比例尺地形图测图精度要求.【总页数】3页(P153-155)
【作者】杨沛旺;徐萌;罗梦;李孝玲
【作者单位】黑龙江地理信息工程院,黑龙江哈尔滨150081;黑龙江地理信息工程院,黑龙江哈尔滨150081;黑龙江地理信息工程院,黑龙江哈尔滨150081;黑龙江地理信息工程院,黑龙江哈尔滨150081
【正文语种】中文
【中图分类】P228.4
【相关文献】
1.基于多面函数的 GPS 高程拟合模型的研究 [J], 白洪伟;吴满意
2.基于二次曲面及多面函数的GPS高程拟合研究 [J], 周明明;王山东;王璐
3.基于多面函数模型的GPS高程拟合精度分析 [J], 李秀海;韩冰
4.基于蚁群算法的多面函数在GPS高程拟合中的应用 [J], 蒲伦;唐诗华;张紫萍;胡新凯;肖燕
5.基于EGM2008模型和多面函数的似大地水准面拟合综合研究 [J], 王仲锋因版权原因，仅展示原文概要，查看原文内容请购买。
用MLS多面函数拟合法求解高程异常
用MLS多面函数拟合法求解高程异常李国鹏;刘站科;柏华岗【期刊名称】《测绘工程》【年(卷),期】2011(020)005【摘要】While obtaining the high-accuracy height-anomaly, multi-quadric function interpolation in terms of theory is not bad accuracy to approach the arbitrarily complicated curved surface arbitrarily, but the core-function and smoothing factor, the choice of weight will exceedingly affect the fitting impression. By adjusting rationally the weight function with the multi-quadric function interpolation of MLS, established on the different model of multi-quadric function interpolation to the different point of certainty, effectively improved the fitting accuracy of height-anomaly.%在求解高精度的高程异常时,多面函数在理论上可以以任意精度逼近任意复杂曲面,但其核函数与平滑因子的选取,以及定权的方式都会对拟合效果产生很大影响,而移动最小二乘(MLS)多面函数法通过对权函数合理的调整,使不同的待定点建立不同的多面函数模型,有效提高高程异常的拟合逼近精度.【总页数】3页(P33-35)【作者】李国鹏;刘站科;柏华岗【作者单位】国家测绘局第一大地测量队,陕西西安710054;国家测绘局第一大地测量队,陕西西安710054;国家测绘局第一大地测量队,陕西西安710054【正文语种】中文【中图分类】P207【相关文献】1.利用多面函数拟合法解算连续化的深度基准面模型 [J],2.二次多项式拟合法和多面函数拟合法在GNSS高程测量中的应用 [J], 魏逸飞;孙伟;黄涛;王建华3.多面函数拟合法转换GPS高程 [J], 金时华4.利用函数拟合法求解高程异常在山区的应用 [J], 高玉平5.一种基于Tikhonov正则化的改进多面函数拟合法 [J], 彭钊;陈志遥;李赫因版权原因，仅展示原文概要，查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
筑 ,2007 ,33 (12) :3482349.
Application study of polyhedral function height f itting
(上接第 297 页)
Maintenance and repair of subways with small2radius curve
WEI Xian2kun Abstract : Through the analysis of usual problems in subways with small2radius curve and their reasons , the author puts forward measures in adjusting geometric size , strengthening technical protection and treating and preventing key problems in its daily maintenance. Furthermore , the author explores the continuous improvement and technology management of those measures. Key words : subway line maintenance , track with small2radius curve , problem , maintenance and repair , technology management
[ M ] . 北京 :测绘出版社 ,1998.
The application of remote sensing technology in urban planning
L IU Ping Abstract : This article analyses t he present application of remote sensing technology in urban planning. From its self2development , it illustrates t he expanding in t he dept h and widt h of t he application of remote sensing technology in t he design and management of urban planning. It also gives advices on how to resolve t he t hings t hat hinder t he application of remote sensing technology. All of t hose provide a technical support to t he sustainable development of cities. Key words : remote sensing technology , urban planning , urban survey
0. 048 0. 011 0. 016 0. 023 - 0. 022 - 0. 012 0. 013 0. 024 0. 001 - 0. 049 0. 032 0. 046 0. 034 0. 048 - 0. 043 - 0. 035
0. 036 0. 012 - 0. 102 - 0. 031 0. 074 0. 045 - 0. 053 - 0. 067 0. 09 - 0. 028 - 0. 026 - 0. 041 - 0. 034 - 0. 065 0. 032 0. 05
差 ,必须将大地高转换为正常高 ,大地高 H 与正常高 h 之间的差
值即为高程异常 :
ζ= H - h
(1)
由式 (1) 可见 ,在 GPS 观测后求得某点的高程异常 ζ,即可精
确确定对应的正常高 。直接求解高程异常采用重力场模型方法 ,
依据重力场长波分量 、已知点大地水准面差距 、斯托克斯方程数
字积分的长波分量的球谐函数表达式和地面重力测量结果直接
点 ,即能很好地描述该区域内高程异常分布的特征点 ,最好位于
最高 、最低及坡度变化处 。
多面函数拟合法理论严密 ,但在核函数选择 、光滑因子 δ2 的
确定 、核函数节点的选取上具有一定的难度 。
2 程序设计与实现
[ 3 ] 孔祥华 ,黄泽民 ,孙建华. 遥感技术在城市综合调查中的应 用[J ] . 城市勘测 ,2003 (2) :65266.
WANG Shi2jun Abstract : Principle of polyhedral function is discussed , on this basis , programming the polyhedral function height fitting software. Fitting precision and key programs are analyzed by an example , a contrastive analysis of fitting precision is illustrated , it points out the notices when the method of polyhedral function is using , so as to provide a reference. Key words : polyhedral function , height fitting , procedure design
多面函数光滑因子及核函数的选取十分重要 ,直接影响拟合 精度 ,但是 ,它们的选择往往需要反复试算 ,因此要求使用者具有 一定的经验 。 参考文献 : [1 ] 李 鹏 ,田林亚 ,李 斌. 多面函数法 GPS 高程拟合在桥梁建
设中的应用研究[J ]. 现代测绘 ,2005 ,28 (2) :13214. [2 ] 高 伟 ,徐绍铨. GPS 高程分区拟合转换正常高的研究[J ]. 武
山 西 建 筑 第 34 卷 2008
第 年
15 5
期 月
SHANXI
ARCHI T EC TU R E
VMoaly. .34 2N0o0.815
·355 ·
文章编号 :100926825 (2008) 1520355202
多面函数法高程拟合应用研究
开发软件界面如图 2 所示 。
3 工程应用分析
为了便于检测拟合效果 ,选用了某条全线进行了四等水准测 量的普通铁路四等 GPS 控制网进行拟合测试 ,该 GPS 网为狭长的 带状网 ,全网长约 55 km ,总计 23 个控制点 (见图 3) 。
本实例核函数选为正双曲面函数 ,选用已知点共 7 个 ,分别 为: G002 , G003 , G005 , G007 , G011 , G013 , G015 ,经多次试算后选中
收稿日期 :2008202217 作者简介 :王世君 (19792 ) ,男 ,助理工程师 ,中铁第一勘察设计院集团有限公司 ,陕西 西安 710043
·356 ·
第 34 卷 2008
第 年
15 5
பைடு நூலகம்
期 月
山
西
建
筑
采用. N ET 开发环境 ,C # 语言编程实现了多面函数法高程 拟合程序 ,程序算法流程图如图 1 所示 。
G002 , G003 , G007 , G011 , G013 , G0156 点作为核心点 ,光滑因子选定 为 1。
为了多方比较验证 ,本例选用同样的 7 个点以二次曲面法进 行高程拟合 。多面函数与二次曲面两种方法的拟合结果与已知成
果对照如表 1 所示 。
表 1 不同算法拟合结果
m
点名
水准高程
[ 4 ] 张 伟 ,顾朝林. 城市规划信息技术开发及应用 [ M ] . 南京 : 东南大学出版社 ,2001.
[ 5 ] 孙东晓 ,刘 涛. 现代城市规划浅谈 [J ] . 山西建筑 ,2007 ,33
(24) :63265. [6 ] 王春峰. 用遥感和单元自动演化方法研究城市扩顺问题
(3)
其中 ,δ2 为任意常数 ,称为光滑因子 ; k 的取值有多种 , k =
1/ 2 时为正双曲面函数 , k = - 1/ 2 时为倒双曲面函数 。
设有 n 个已知点 ( xi , yi) ( i = 1 , 2 , …, n) ,选其中 m ( m ≤n)
个点 ( xj , yj ) ( j = 1 , 2 , …, m ) 为核函数的核心点 , 并令 Qij = F
( xi , yi , xj , yj) ,则式 (2) 转变为 :
m
∑ ζi =
βj Qij
(4)
j =1
由此可列出误差方程式 :
v1
Q11 Q12 … Q1 m β1
ζ1
v2
Q21 Q22 … Q2 m β2
ζ2
⁝= ⁝ ⁝ ⁝ ⁝ ⁝ - ⁝
(5)
vn
Q n1 Q n2 … Q nm βm
ζn
m
∑ ζ( x , y) = βj F ( x , y , xj , yj)
(2)
j =1
其中 ,βj 为待定参数 ; F ( x , y , xj , yj ) 为 x 和 y 的二次核函
数 ,常用的二次核函数为 :
F ( x , y , xj , yj) = [ ( x - xj) 2 + ( y - yj) 2 +δ2 ] k
因此拟合效果与核函数的选择密切相关 。核函数有很多种 ,文中
选择为正双曲面函数 ,即 :