EGM重力场模型在工程中的应用

---

关于EGM基本介绍

EGM：Earth Gravitational Model，地球重力模型。

地球表面是极其复杂的表面，但是海水面占据全球表面的大部分，且比较规则，在某种假设下，可认为海水面是重力等位面，并把它延伸到大陆下，可以得到一个遍及全球的等位面。德国的J.B.Listing于1872年把它命名为大地水准面。
大地水准面是最接近平均海水面的等位面。

使用不同的椭球（ellipsoid）近似不同的大地水准面（Geoid），我们可以得到WGS84椭球、GRS80椭球等等。因为椭球只是近似表达大地水准面，在地球不同位置上，椭球和大地水准面是有差异的。

EGM可以用来计算两个面之间的差异：

NGA发布了三种重力场模型：

• EGM84
• EGM96
• EGM2008
• EGM2020（官网暂未提供）

EGM96提供360阶的球谐系数
EGM2008提供2190阶且2159次的球谐系数

---

以EGM2008的举例说明EGM2008的作用：
图.EGM2008图谱化展示

图.WGS84椭球大地高与EGM2008的关系

工程应用

1. 利用EGM模型补充改正网上公开的dem。
   网上公开的DEM，比如SRTM、AW3D30、GDEM等等，有的采用WGS84椭球大地高，有的加入了EGM96模型的高程异常值，需要利用EGM模型进行补充改正。

常见DEM高程基准信息如下图（参考SRTM、ASTER、ALOS、TanDEM-X等常用免费数字高程模型（DEM）基本信息汇总）

---

2. RTK使用CORS测量的2000大地高转85高
   常用的方法：

   • 找测区周边85高的控制点来计算高程拟合参数转换
   • 找CORS主管部门转换
   • 基于EGM2008大地水准面模型转换

   关于第三种方法【基于EGM2008大地水准面模型转换】的应用：
   (1) 基本原理：$85高=2000大地高-高程异常$
   (2) 应用情况：
   a. 测区周边无85高的控制点情况下：

    	a.1. 提取测区测点，利用高程异常计算工具计算EGM2008高程异常
    	a.2. 利用基本原理和EGM2008高程异常，转换2000大地高为85高
   

   b. 测区周边存在85高的控制点情况下：

    	b.1. 提取测区测点，利用高程异常计算工具计算异常
    	b.2. 提取控制点的85高和2000大地高，计算控制点的高程异常（= 控制点2000大地高 - 控制点85高）
    	b.3. 提取控制点位置，利用高程异常计算工具计算EGM2008高程异常
    	b.4. 计算控制点高程异常综合改正数（= 控制点的高程异常 - 控制点的EGM2008高程异常）
    	b.5. 如果存在多个控制点，需要按照一定的方法进行插值求解高程异常综合改正数
    	b.6. 根据上述步骤得到各个测点的高程异常综合改正数，计算每个测点的高程异常拟合参数（= 测点的EGM2008高程异常 + 测点的高程异常综合改正数）
    	b.7. 利用基本原理和测点的高程异常综合改正数，转换2000大地高为85高
   

   附.高程异常计算工具:
   a. Alltrans EGM2008 Calculator
   b. GAGE-Geoid Height Calculator
   c. Online geoid calculations using the GeoidEval utility
   d. GeographicLib开发库
   关于GeographicLib的CMAKE编译：

   cd @geographiclib的根目录
   @cmake安装目录\bin\cmake.exe -G “Visual Studio 16” -A x64 -D CMAKE_INSTALL_PREFIX=“@geographiclib的运行环境编译结果安装目录/cmake_install_vs2019_x64” -B “@geographiclib的cmake编译结果目录\cmake_build_vs2019_x64”

   关于GeographicLib的使用：

   api调用
   应用调用：
   Gravity应用计算高程异常：echo 26:28:22.168N 100:6:23.513E 0 | Gravity -n egm2008 -d “@GeographicLib的gravity数据安装目录” -H
   –
   关于Gravity应用的一些选项说明：
   -n：指定使用哪个EGM模型。例如EGM96、EGM2008
   -d：指定gravity数据的所在目录。"Gravity -h"可以列出默认指定的gravity数据所在目录，默认目录可以通过环境变量GEOGRAPHICLIB_GRAVITY_PATH或者GEOGRAPHICLIB_DATA设置
   -H：根据EGM模型计算高程异常。Gravity应用计算的高程异常比GeoidEval精确，因为GeoidEval是根据栅格数据来采样

   注.使用EGM2008有两点需考虑:

    a. EGM2008与85高两者似大地水准面高程基准有细微偏差；
    b. EGM2008参考椭球是WGS84，WGS84与CGCS2000坐标之间大地高的存在细微偏差
   

---

参考引用

1. NOAA-A tutorial on Datums
2. eye4software-Using the EGM2008 geoid model
3. 没有控制点如何把CORS测的2000大地高转为85高
4. SRTM、ASTER、ALOS、TanDEM-X等常用免费数字高程模型（DEM）基本信息汇总