基于地形改正的GPS高程异常曲面的拟合.pdfVIP

第 1期 矿 山 测 量 NO．1 2008年 3月 MINE SURVEYING Mar．2008 基于地形改正的GPS高程异常曲面的拟合 汤凯煌 (1．岳阳市房地产大厦测绘中心，湖南岳阳 414000) 摘要：分析了在地形起伏较大地区传统高程拟合方法效果不好的原因，提出了基于地形改正的GPS 曲面拟合模型。针对工程测量实践中用于地形改正的 DEM数据获取难的问题，文中尝试使用 NASA对外公布的90米分辨率的 DEM数据，工程实践结果表明文中提出的新方法和思路是可行 的。 关键词：高程异常；曲面拟合；地形改正；DEM 中图分类号：P228．4 文献标识码：B 文章编号：1001—358X(2008)0l～0035—03 GPS几何水准是近年来发展起来的测定水准的 表示似大地水准面长波特征，在100～200kin内 新方法，在GPS工程应用上特别常见，GPS几何水准 变化均匀； 表示波长一般在20～100kin之间的局 方法确定高程异常实质是用一个光滑的平面或曲面 部似大地水准面特征； 。表示 20kin以下波长特 来代替似大地水准面，这在平坦地区是适用的，然而 征 一 。 对于山区，地形起伏较大，似大地水准面变化剧烈， 用光滑的平面或曲面来代替似大地水准面，必然会 带来较大的偏差 。近年的研究表明：对于某一 地形起伏较大的测区而言，由地形起伏引起的短波 部分严重影响着GPS高程拟合的精度，有时甚至达 到cm级，因此必须地形改正。目前，虽然它的计算 公式已经比较成熟，但是应用到工程实践中的并不 多，其中的主要原因是地面DEM数据的获取和计算 的复杂。文献[3]中曾尝试用顾及非格网的数据进 行地形改正，但是计算量大并且比较复杂，同时参与 内插大地高的点并不能真实的体现测区的DEM模 图l lOO km内似大地水准面波长特征 型。本文尝试用 NASA公开的美国航天飞机全球测 图任务(SRTM)所获取的DEM数据进行地形改正， 在小区域 GPS工程网中，一般情况下，在中长波 然后用多项式进行曲面拟合。 项，大地水准面与似大地水准面变化趋势基本相同 并且变化的趋势很平缓，可以用多项式曲面拟合来 1 地形改正的几何方法 实现。而对于高程异常的短波项主要是由地形起伏 精确地确定似大地水准面(即高程异常)所用的 所引起的，地形起伏越大，似大地水准面越粗糙；反 数据类型是地球重力场模型或平均重力异常及数字 之越光滑。因此本文提出基于地形改正的GPS多项 高程模型。将这些数据进行综合处理，可给出表示 式曲面拟合模型。 高程异常至关重要的各种波长的最优解，这一解可 地形起伏可看作重力场的噪声，设地球上有一 表达如下： 个参考面(平滑面)，通过计算此面以外、以内质量的 = GM+ + rc (1) 影响，即得到地形起伏的影响。若要计算周围地形 式中， 为地球重力场模型； 为重力异常； 表 起伏对某点高程异常的影响，只需将周围地形起伏 示地形起伏对高程异常的影响，如图 1所示。其中， 点的扰动位进行叠加，即可求得 计算方法 。