2023年2月第1期城市勘测UrbanGeotechnicalInvestigation＆SurveyingFeb．2023No．1引文格式:胡玉祥，刘宝华，张洪德等．基于EGM2008和移动曲面模型的GNSS高程拟合方法研究及应用［J］．城市勘测，2023(1):79－82．文章编号:1672－8262(2023)01－79－04中图分类号:P226文献标识码:A基于EGM2008和移动曲面模型的GNSS高程拟合方法研究及应用胡玉祥1，2，刘宝华1，2，张洪德1，2，尹相宝1，2，孟庆年1，2*收稿日期:2022—02—17作者简介:胡玉祥(1988—)，男，硕士，工程师，注册测绘师，主要从事GNSS数据处理与应用、城市轨道交通测量等应用研究工作。E－mail:huyuxiang1223@163．com(1.青岛市勘察测绘研究院，山东青岛266032;2.山东省海岸带调查监测工程技术协同创新中心，山东青岛266032)摘要:随着GNSS高程测量的广泛使用，如何高精度求取高程异常成为作业的关键。鉴于沿海岛礁缺乏地形资料，提出一种基于EGM2008地球重力场模型提取高程异常长波项，移动曲面函数拟合残余高程异常的方法。结合青岛沿海岛礁实测数据，通过与常规直接拟合方法进行对比分析，结果表明，该方法能够提高GNSS高程拟合的精度，于沿海岛礁可以达到厘米级，可以替代四等水准测量。关键词:高程异常;EGM2008;移动曲面函数;高程拟合;四等水准1引言GNSS高程测量以其效率高、灵活性强等优势在工程测量中得到了广泛的应用，尤其对于实测常规水准测量困难的地区。而实际应用中需要将GNSS测得的大地高转化为正常高，关键在于精确求取高程异常［1］。高程异常的常规获取方法是采用曲面函数拟合方式获取，然而此方法只能表示出测区的大致变化趋势，对于细节方面表现较差，一些工程中不能满足其对高程测量精度的要求。随着卫星重力计划的实施和逐步完善，尤其是EGM2008高阶重力场模型［2］的公布，使得GNSS高程转换的精度得到了很大提升［3，4］。EGM2008的空间分辨率为5'×5'(约为9km×9km)，能够精确表达高程异常中的长波项［3～5］。对于地形起伏较大的地区，仅考虑高程异常的长波项不足以反映其变化，这时往往需要考虑由地形引起的短波项的影响。高程异常可以分为3部分:重力场模型高程异常长波项、地形起伏变化高程异常短波项和残余高程异常。实际应用中，常利用数值逼近算法考虑地形起伏变化的影响［6，7］，而忽略高程异常的长波项;或只考虑高程异常的长波项，不考虑地形起伏的影响［7～9］;抑或综合考虑高程异常的长波项和短波项的影响。对于高程异常残差的拟合方法很多，比如二次...
