[发明专利]一种基于Matlab的区域似大地水准面精化方法

权利要求书

1.一种基于Matlab的区域似大地水准面精化方法，其特征在于，包括如下步骤：

S1.采集待测区域的重力数据G；对所述重力数据G进行粗差剔除，然后将经过粗差剔除后得到的重力数据进行预处理，得到格网化法耶重力异常G_pre；

S2.利用移去恢复法对所述格网化法耶重力异常G_pre进行处理,得到待测区域的大地水准面N；利用KTH法对所述格网化法耶重力异常G_pre进行处理，得到待测区域的大地水准面N；

S3.根据水准数据提供的高程异常和所述大地水准面N，完成与国家高程基准的融合，得到融合结果R₃；根据水准数据提供的高程异常和所述大地水准面N，完成与国家高程基准的融合，得到融合结果R₄；

S4.利用水准点对所述融合结果R₃和所述融合结果R₄进行精度验证，选择精度高的融合结果作为最终的待测区域似大地水准面。

2.根据权利要求1所述的基于Matlab的区域似大地水准面精化方法，其特征在于，步骤S1的具体内容为：

利用重力仪测量待测区域的重力数据G，所述重力数据G包括待测区域的绝对重力值G_abs和正常重力值G_normal，利用公式ξ^G＝G_abs-G_normal求取待测区域的重力异常ξ^G；根据空间改正方法对所述重力异常ξ^G进行处理得到待测区域的空间重力异常β^G；对所述空间重力异常β^G经过层间改正、局部地形改正和均衡改正后得到待测区域的均衡重力异常ω^G，将所述均衡重力异常ω^G减去层间改正、局部地形改正和均衡改正，得到格网化法耶重力异常G_pre。

3.根据权利要求1所述的基于Matlab的区域似大地水准面精化方法，其特征在于，步骤S2中利用RCR法计算待测量区域的重力大地水准面的具体内容为：根据重力场模型计算大地水准面差距N_GM，计算地形间接影响N_T，计算残差大地水准面N_RES，则最终的大地水准面N＝N_GM+N_T+N_RES。

4.根据权利要求1所述的基于Matlab的区域似大地水准面精化方法，其特征在于，步骤S2中利用KTH法计算待测量区域的重力大地水准面的具体内容为：对大地水准面进行预估，得到大地水准面预估值经过垂线延拓校正δN_DWC、综合地形校正综合大气校正和大地水准面球近似的椭球校正δN_e后得到如下式所示的重力大地水准面