地球物理大地测量学计算系列之二十三过指定点重力等位面构造-尧图手机网站定制

1、过指定点重力等位面构造

利用高阶地球重力场位系数模型，计算过指定点(经度，纬度，大地高)重力等位面相对于指定点的大地高差和重力差格网。利用局部重力场数据，能进一步精化过指定点重力等位面的大地高差和重力差格网。

[功能]利用参考地球重力场位系数模型，计算过指定点(B,L,H)重力等位面的重力(mGal)和大地高(m)模型值格网。
[输入]重力等位面范围格网文件，参考重力场位系数模型文件。
等位面范围格网文件，只用于确定等位面格网的经纬度范围和分辨率，程序忽略格网的格值数据。
位系数模型文件第1行约定为位系数尺度参数:地心引力常数GM与椭球长半轴a。尺度参数(GM、a)是位系数模型的表示基准。
[参数设置]输入指定点大地坐标(经纬度大地高)，与位系数模型最大计算阶数。
程序要求指定点落在等位面格网范围内。位系数模型阶数不宜过大，程序需要迭代逼近计算，耗费时间较长。
[输出]过指定点重力等位面的重力和大地高模型值格网文件。

过指定点重力等位面的重力模型值格网文件

过指定点重力等位面的大地高模型值格网文件

2、程序编制

由于无需考虑高程系统，参考文献1的代码需要修改，可根据上面的公式进行修改，也可通过deepseek修改，程序主界面如下图：

计算的大地高、大地高差和重力差，见下图。

大地高

大地高差，大地高-指定点的大地高

重力差

从上述文件可以看出：大地高差、指定点的重力位和重力与章老师程序一致，但是重力差不一致。指定点的重力和重力位为：978718.0644 ； 62623481.0494

重力差的第一点经纬度和大地高为：27.0208333333333，106.020833333333，1335.9711

根据以前的程序计算第一点的正常位和扰动位、正常重力和扰动重力，见下图。

正常位+扰动位=62623780.97+ （-299.9164 ）= 62623481.05，与指定点重力位相等。