CN105953772A - 一种高程异常的测定方法及测定装置 - Google Patents

Abstract

本发明公开了一种高程异常的测定技术，属于测量测绘领域。一种高程异常的测定方法，包括：获取设备当前经纬度信息；根据获得的经纬度信息在预设的高程异常网点查找表确定对应的第一目标网点；从高程异常网点查找表获取第一目标网点的高程异常，以及距离第一目标网点预设间隔距离的第二目标网点的高程异常及对应经纬度信息；根据获得的高程异常、经纬度信息及预置的多项式曲面公式计算得到多项式曲面的拟合参数值；根据所述经纬度信息、多项式曲面公式及拟合参数值计算得到所述经纬度信息的高程异常。本发明还公开一种高程异常的测定装置。通过本发明提供的高程异常的测定技术，可以解决GPS测绘工程中高程异常测量困难且精度不高的问题。

Description

一种高程异常的测定方法及测定装置

技术领域

本发明涉及测量测绘领域，尤其涉及一种应用于GPS测绘工程中高程异常的测定方法及测定装置。

背景技术

GPS测绘工程中，高程异常是似大地水准面至地球椭球面的高度。

目前，在GPS测绘工程中确定高程异常的传统方法主要有以下几种：1、采用GPS水准测量高程异常；2、采用重力场模型求解高程异常；3、根据不同区域的地理形状，通过测量一定数量GPS-水准点，采用曲线、曲面等拟合方法进行高程异常拟合。

在我国的丘陵以及山区地区，采用GPS-水准测量高程异常工程量大且在山区地区测量是比较困难；并且由于缺乏我国丘陵以及山区等地区精确的重力场数据导致求解高程异常算法复杂并且精度不高；而采用曲线曲面等拟合方法只能局限于局部区域，并不能适合包括所有地形区域。

所有，有必要提供一种新的测量点高程异常的测定方法，以解决GPS测绘工程上高程异常测量困难且精度不高的问题。

发明内容

有鉴于此，本发明要解决的技术问题是提供一种测量点高程异常的测定方法及测定装置，以解决GPS测绘工程中高程异常测量困难且精度不高的问题。

为解决上述技术问题，本发明提供一种高程异常的测定方法，所述测定方法包括：

获取设备当前的经纬度信息；

根据上述获得的经纬度信息在预设的高程异常网点查找表中确定对应的第一目标网点；所述第一目标网点为所述设备当前所处位置对应的网点；

从高程异常网点查找表获取所述第一目标网点的高程异常，以及距离所述第一目标网点预设间隔距离的N个第二目标网点的高程异常及对应经纬度信息；

根据上述获取到的高程异常、经纬度信息及预置的多项式曲面公式计算得到所述多项式曲面的拟合参数值；

根据所述经纬度信息、所述多项式曲面公式及所述拟合参数值计算得到所述经纬度信息的高程异常。

进一步地，所述预设的高程异常网点查找表是以经纬度值的整数为网点坐标建立的，各网点间隔距离根据具体的精度要求来预设。

进一步地，所述根据上述获得的经纬度信息在预设的高程异常网点查找表中确定对应的第一目标网点包括：

将所述经纬度信息对应的经纬度值按预设的调整规则调整为整数；

将调整后的经纬度值所对应的网点设定为所述第一目标网点。

进一步地，所述预设的调整规则包括以下取整规则之一：四舍五入、向上取整、向下取整。

进一步地，所述根据上述获得的经纬度信息在预设的高程异常网点查找表中确定对应的第一目标网点包括：

将所述经纬度信息，与预设的高程异常网点查找表中的各网点的经纬度信息进行比对，将与所述经纬度信息最接近的对应网点确定为第一目标网点。

进一步地，所述N个第二目标网点是以所述第一目标网点为中心，根据预设间隔距离按照预置的取点规则从高程异常网点查找表获取与所述第一目标网点相邻的网点来确定。

进一步地，所述预置的取点规则能满足特定形状的边上覆盖至少N个第二目标网点。

进一步地，所述预置的多项式曲面公式为：

f(x,y)＝a₀+a₁x+a₂y+a₃x²+a₄xy+a₅y²

上式中，f(x,y)表示高程异常，x表示纬度，y表示经度，a₀、a₁、a₂、a₃、a₄和a₅为所述多项式曲面的拟合参数。

为解决上述技术问题，本发明还提供一种高程异常的测定装置，所述测定装置包括：

第一获取模块，用于获取设备当前的经纬度信息；

确定模块，用于根据上述获得的经纬度信息在预设的高程异常网点查找表中确定对应的第一目标网点；所述第一目标网点为所述设备当前所处位置对应的网点；

第二获取模块，用于从高程异常网点查找表获取所述第一目标网点的高程异常，以及距离所述第一目标网点预设间隔距离的N个第二目标网点的高程异常及对应经纬度信息；

第一计算模块，用于根据上述获取到的高程异常、经纬度信息及预置的多项式曲面公式计算得到所述多项式曲面的拟合参数值；

第二计算模块，用于根据所述经纬度信息、所述多项式曲面公式及所述拟合参数值计算得到所述经纬度信息的高程异常。

与现有技术相比，本发明提供的一种高程异常的测定方法及测定装置，通过获取设备当前的经纬度信息，根据获得的经纬度信息在预设的高程异常网点查找表中确定对应的第一目标网点，并将所述第一目标网点为作所述设备当前所处位置对应的网点；再从高程异常网点查找表获取所述第一目标网点的高程异常，以及距离所述第一目标网点预设间隔距离的N个第二目标网点的高程异常及对应经纬度信息；然后根据上述获取到的高程异常、经纬度信息及预置的多项式曲面公式计算得到所述多项式曲面的拟合参数值；最后根据所述经纬度信息、所述多项式曲面公式及所述拟合参数值计算得到所述经纬度信息的高程异常。通过本发明以上的技术方案，就可以在GPS测量过程中根据卫星定位到的测量点的经度纬度自动进行高程异常拟合，既解决了采用传统GPS-水准测量工程量大且在山区地区测量困难的问题，也避免了我国由于丘陵以及山区等地区重力场数据不足而导致求解高程异常算法复杂并且精度不高的难题，也能够简化网点查找表的网点精度。

附图说明

图1为本发明实施例提供的一种高程异常的测定方法的流程图；

图2为本发明实施例提供的一种高程异常的测定方法中选取网点的示意图；

图3为本发明实施例提供的一种高程异常的测定装置的结构示意图。

具体实施方式

为了使本发明所要解决的技术问题、技术方案及有益效果更加清楚、明白，以下结合附图和实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

如图1所示，本发明提供一种高程异常的测定方法，所述测定方法包括：

1、获取设备当前的经纬度信息。

2、根据上述获得的经纬度信息在预设的高程异常网点查找表中确定对应的第一目标网点；所述第一目标网点为所述设备当前所处位置对应的网点。

其中：

预设的高程异常网点查找表是以经纬度值的整数为网点坐标建立的，该高程异常网点查找表中的各网点间隔距离可以根据具体的精度要求来预设，例如预设间隔距离可以为1°、2°、3°等；高程异常网点查找表除了包括经度纬度的信息之外，还包括有高程异常的信息。

上述根据所述经纬度信息在预设的高程异常网点查找表中确定对应的第一目标网点，包括：

将所述经纬度信息对应的经纬度值按预设的调整规则调整为整数；其中，预设的调整规则可以根据具体的精度要求来进行设定，例如可以采用四舍五入、向上取整、向下取整等规则来进行取整；

将调整后的经纬度值所对应的网点设定为所述第一目标网点。

另一方面，上述根据所述经纬度信息在预设的高程异常网点查找表中确定对应的第一目标网点，也可以包括：

将设备当前的所述经纬度信息，与预设的高程异常网点查找表中的各网点的经纬度信息进行比对，将与设备当前的所述经纬度信息最接近的对应网点确定为第一目标网点。

3、从高程异常网点查找表获取所述第一目标网点的高程异常，以及距离所述第一目标网点预设间隔距离的N个第二目标网点的高程异常及对应经纬度信息。

其中：

所述N个第二目标网点是以所述第一目标网点为中心，根据预设间隔距离按照预置的取点规则从高程异常网点查找表获取与所述第一目标网点相邻的网点来确定。N的取值个数至少要满足计算预置的多项式曲面公式的拟合参数的要求，具体来说，N的取值个数至少要大于或等于预置的多项式曲面公式的拟合参数个数减1。

预置的取点规则主要能满足特定形状的边上覆盖至少N个第二目标网点即可，例如：预置的取点规则可以按“田”字形状来取点，即以田字中心点为第一目标网点，取田字四周与第一目标网点相邻的8个网点为第二目标网点；或者，预置的取点规则可以按四边形状取点，即以四边形的中心点为第一目标网点，取四边形四周与第一目标网点相邻的4个网点为第二目标网点；或者，预置的取点规则可以按五角形状取点，即以五角形的中心点为第一目标网点，取五角形各个顶点与第一目标网点相邻的至少4个网点为第二目标网点。

从高程异常网点查找表获取所述第一目标网点的高程异常以及上述N个第二目标网点的高程异常及对应经纬度信息。

4、根据上述获取到的高程异常、经纬度信息及预置的多项式曲面公式计算得到所述多项式曲面的拟合参数值。

其中，预置的多项式曲面公式为：

f(x,y)＝a₀+a₁x+a₂y+a₃x²+a₄xy+a₅y²

上式中，f(x,y)表示高程异常，x表示纬度，y表示经度，a₀、a₁、a₂、a₃、a₄和a₅为所述多项式曲面的拟合参数。

根据以上获取到的各个网点的高程异常、经纬度信息的相关信息，利用以上预置的多项式曲面公式进行拟合，得到拟合参数a₀、a₁、a₂、a₃、a₄和a₅的数值。

将上述得到的拟合参数a₀、a₁、a₂、a₃、a₄和a₅的数值返回到上述的多项式曲面公式，即可得到带有具体拟合参数值的多项式曲面公式如下：

f(x,y)＝a₀+a₁x+a₂y+a₃x²+a₄xy+a₅y²

上式中，f(x,y)表示高程异常，x表示纬度，y表示经度，a₀、a₁、a₂、a₃、a₄和a₅为所述多项式曲面的拟合参数。

另一方面，本发明除了可以采用上述预置的这种二次多项式曲面公式之外，还可以采用其他的多项式曲面公式，例如，三次曲面拟合法：

f(x,y)＝a₀+a₁x+a₂y+a₃x²+a₄xy+a₅y²+a₆x³+a₇x²y+a₈xy²+a₉y³

但需要说明的是，不管是采用何种多项式曲面公式，对于上述第二目标网点的取值是根据确定选用的多项式曲线公式来进行确定的。

5、根据所述经纬度信息、所述多项式曲面公式及所述拟合参数值计算得到所述经纬度信息的高程异常。

以下结合一个具体的应用实施例来对以上提到的本发明技术方案作一步的说明。

如图2所示。根据本发明的技术方案，提出一种高程异常的测定方法的一个具体实施例。

本实施例中，设备当前的经纬度信息为纬度22.5°和经度113.9°，现在需要得到该经纬度的精确高程异常。

根据本发明的技术方案，一种高程异常的测定方法，该方法包括：

1、获取设备当前的经纬度信息；

本实施例中，获取到设备当前的经纬度信息为纬度22.5°和经度113.9°。

2、根据上述获得的经纬度信息在预设的高程异常网点查找表中确定对应的第一目标网点；所述第一目标网点为所述设备当前所处位置对应的网点。

其中，预设的高程异常网点查找表是以经纬度值的整数为网点坐标建立的，在本实施例中，网点预设间隔距离为1°。

根据本实施例设备所在的经纬度信息(纬度22.5°和经度113.9°)，预设的高程异常网点查找表的部分相关信息如下：

高程异常网点查找表

序号	经度	纬度	高程异常
				1	113	26	-12.729
2	114	26	-7.969
				3	115	26	-5.426
4	113	25	-10.163
				5	114	25	-7.573
6	115	25	-3.542
				7	113	24	-9.759
8	114	24	-5.642
				9	115	24	-1.984
10	113	23	-7.525
				11	114	23	-4.1
12	115	23	-0.102
				13	113	22	-6.096
14	114	22	-2.172
				15	115	22	1.468
16	113	21	-4.024
				17	114	21	-0.385
18	115	21	3.713
				19	113	20	-1.701
20	114	20	1.977

将所述经纬度信息对应的经纬度值按照预设调整规则调整为整数；在本实施例中，预设调整规则采用四舍五入规则来进行取整。

本实施例中，设备所在经纬度信息为[纬度22.5°、经度113.9°]，经过四舍五入取整后经纬度信息为[纬度，经度]＝[23°，114°]。

将调整后的经纬度[23°，114°]对应的网点设定为所述第一目标网点。

3、从高程异常网点查找表获取所述第一目标网点的高程异常，以及距离所述第一目标网点预设间隔距离的N个第二目标网点的高程异常及对应经纬度信息。其中，所述N个第二目标网点是以所述第一目标网点为中心，根据预设间隔距离按照预置的取点规则从高程异常网点查找表获取与所述第一目标网点相邻的网点来确定。

由于N的取值个数至少要满足计算预置的多项式曲面公式的拟合参数的要求，具体来说，N的取值个数至少要大于或等于预置的多项式曲面公式的拟合参数个数减1。

所以，在本实施例中，由于预置的多项式曲面公式的有a₀、a₁、a₂、a₃、a₄和a₅共6个拟合参数，所以N的取值个数至少要大于或等于5个(6-1＝5)，本实施例中，N取8。

预置的取点规则按“田”字形状来取点，即以田字中心点为第一目标网点，取田字四周与第一目标网点相邻的8个网点为第二目标网点。

如图2所示。在本实施中，以[纬度，经度]＝[23°，114°]对应的网点为第一目标网点，从高程异常网点查找表获取所述第一目标网点的高程异常为-4.1。以该第一目标网点为中心，按田字形状在高程异常网点查找表中选取离该第一目标网点最近的8个第二目标网点分别是：[22，114]、[24，114]、[22，113]、[23，113]、[24，113]、[22，115]、[23，115]、[24，115]，从高程异常网点查找表获取对应的高程异常分别为[-4.1]、[-2.172]、[-5.642]、[-6.096]、[-7.525]、[-9.759]、[1.468]、[-0.102]、[-1.984]。

4、根据上述获取到的高程异常、经纬度信息及预置的多项式曲面公式计算得到所述多项式曲面的拟合参数值。

其中，预置的多项式曲面公式为：

f(x,y)＝a₀+a₁x+a₂y+a₃x²+a₄xy+a₅y²

上式中，f(x,y)表示高程异常，x表示纬度，y表示经度，a₀、a₁、a₂、a₃、a₄和a₅为所述多项式曲面的拟合参数。

在本实施例中，将以上获到的9个网点的经纬度、高程异常的具体数值，分别代入到上述多项式曲面公式中，分别如下：

f(x₁,y₁)＝a₀+a₁*23+a₂*114+a₃*23²+a₄*23*114+a₅*114²＝-4.1

f(x₂,y₂)＝a₀+a₁*22+a₂*114+a₃*22²+a₄*22*114+a₅*114²＝-2.172

f(x₃,y₃)＝a₀+a₁*24+a₂*114+a₃*24²+a₄*24*114+a₅*114²＝-5.642

f(x₄,y₄)＝a₀+a₁*22+a₂*113+a₃*22²+a₄*22*113+a₅*113²＝-6.096

f(x₅,y₅)＝a₀+a₁*23+a₂*113+a₃*23²+a₄*23*113+a₅*113²＝-7.525

f(x₆,y₆)＝a₀+a₁*24+a₂*113+a₃*24²+a₄*24*113+a₅*113²＝-9.759

f(x₇,y₇)＝a₀+a₁*22+a₂*115+a₃*22²+a₄*22*115+a₅*115²＝1.468

f(x₈,y₈)＝a₀+a₁*23+a₂*115+a₃*23²+a₄*23*115+a₅*115²＝-0.102

f(x₉,y₉)＝a₀+a₁*24+a₂*115+a₃*24²+a₄*24*115+a₅*115²＝-1.984

以上公式中，相关的经纬度信息如下：

(x₁,y₁)＝(23,114)、(x₂,y₂)＝(22,114)、(x₃,y₃)＝(24,114)、(x₄,y₄)＝(22,113)、(x₅,y₅)＝(23,113)、(x₆,y₆)＝(24,113)、(x₇,y₇)＝(22,115)、(x₈,y₈)＝(23,115)、(x₉,y₉)＝(24,115)。

高程异常信息如下：

f(x₁,y₁)＝-4.1、f(x₂,y₂)＝-2.172、f(x₃,y₃)＝-5.642、f(x₄,y₄)＝-6.096、

f(x₅,y₅)＝-7.525、f(x₆,y₆)＝-9.759、f(x₇,y₇)＝1.468、f(x₈,y₈)＝-0.102、

f(x₉,y₉)＝-1.984。

采用最小二乘法[At＝(B*B')^-1*B*S]，计算得出拟合参数值为：

[a₀,a₁,a₂,a₃,a₄,a₅]＝[-690.15157,-2.17333,9.04041,-0.12183,0.05272,-0.02833]

其中：

At＝[a₀,a₁,a₂,a₃,a₄,a₅]

S＝[f(x₁,y₁),f(x₂,y₂),f(x₃,y₃),f(x₄,y₄),f(x₅,y₅),f(x₆,y₆),f(x₇,y₇),f(x₈,y₈),f(x₉,y₉)]

$B=\left[\begin{array}{c}1,x_1,y_1,x_1^2,x_1{y}_1,y_1^2\\ 1,x_2,y_2,x_2^2,x_2{y}_2,y_2^2\\ 1,x_3,y_3,x_3^2,x_3{y}_3,y_3^2\\ 1,x_4,y_4,x_4^2,x_4{y}_4,y_4^2\\ 1,x_5,y_5,x_5^5,x_5{y}_5,y_5^2\\ 1,x_6,y_6,x_6^2,x_6{y}_6,y_6^2\\ 1,x_7,y_7,x_7^2,x_7{y}_7,y_7^2\\ 1,x_8,y_8,x_8^2,x_8{y}_8,y_8^2\\ 1,x_9,y_9,x_9^2,x_9{y}_9,y_9^2\end{array}\right]$

将上述得到的拟合参数值[a₀,a₁,a₂,a₃,a₄,a₅]返回到上述多项式曲面公式后，即可得到多项式曲面公式如下：

f(x,y)＝a₀+a₁x+a₂y+a₃x²+a₄xy+a₅y²

上式中，f(x,y)表示高程异常，x表示纬度，y表示经度，拟合参数[a₀,a₁,a₂,a₃,a₄,a₅]＝[-690.15157,-2.17333,9.04041,-0.12183,0.05272,-0.02833]。

5、根据所述经纬度信息、所述多项式曲面公式及所述拟合参数值计算得到所述经纬度信息的高程异常。

在本实施例中，设备所在经纬度信息为纬度22.5°和经度113.9°，将以上经纬度信息代入以上多项式曲面公式中，即可得到设备所在经纬度的高程异常为-3.4155。

f(x,y)＝a₀+a₁*22.5+a₂*113.9+a₃*22.5²+a₄*22.5*113.9+a₅*113.9²＝-3.4155

上式中，f(x,y)表示高程异常，纬度x＝22.5，经度y＝113.9，拟合参数[a₀,a₁,a₂,a₃,a₄,a₅]＝[-690.15157,-2.17333,9.04041,-0.12183,0.05272,-0.02833]。

如图3所示，本发明还提供一种高程异常的测定装置100，该测定装置100包括：第一获取模块10、确定模块20、第二获取模块30、第一计算模块40、第二计算模块50；

第一获取模块10，用于获取设备当前的经纬度信息；

确定模块20，用于根据上述获得的经纬度信息在预设的高程异常网点查找表中确定对应的第一目标网点；所述第一目标网点为所述设备当前所处位置对应的网点；

第二获取模块30，用于从高程异常网点查找表获取所述第一目标网点的高程异常，以及距离所述第一目标网点预设间隔距离的N个第二目标网点的高程异常及对应经纬度信息；

第一计算模块40，用于根据上述获取到的高程异常、经纬度信息及预置的多项式曲面公式计算得到所述多项式曲面的拟合参数值；

第二计算模块50，用于根据所述经纬度信息、所述多项式曲面公式及所述拟合参数值计算得到所述经纬度信息的高程异常。

该第一获取模块10，用于获取设备当前的经纬度信息；

该确定模块20，用于根据上述获得的经纬度信息在预设的高程异常网点查找表中确定对应的第一目标网点；所述第一目标网点为所述设备当前所处位置对应的网点；

其中：

预设的高程异常网点查找表是以经纬度值的整数为网点坐标建立的，该高程异常网点查找表中的各网点间隔距离可以根据具体的精度要求来预设，例如预设间隔距离可以为1°、2°、3°等；高程异常网点查找表除包括经度纬度的信息之外，还包括有高程异常的信息。

上述根据所述经纬度信息在预设的高程异常网点查找表中确定对应的第一目标网点包括：

将所述经纬度信息对应的经纬度值按预设的调整规则调整为整数；其中，预设的调整规则可以根据具体的精度要求来进行设定，例如可以采用四舍五入、向上取整、向下取整等规则来进行取整。

将调整后的经纬度值所对应的网点设定为所述第一目标网点。

该第二获取模块30，用于从高程异常网点查找表获取所述第一目标网点的高程异常，以及距离所述第一目标网点预设间隔距离的N个第二目标网点的高程异常及对应经纬度信息；

其中：

所述N个第二目标网点是以所述第一目标网点为中心，根据预设间隔距离按照预置的取点规则从高程异常网点查找表获取与所述第一目标网点相邻的网点来确定。

预置的取点规则主要能满足特定形状的边上覆盖至少N个点即可，例如：预置的取点规则可以按“田”字形状来取点，即以田字中心点为第一目标网点，取田字四周与第一目标网点相邻的8个网点为第二目标网点；或者，预置的取点规则可以按四边形状取点，即以四边形的中心点为第一目标网点，取四边形四周与第一目标网点相邻的4个网点为第二目标网点；或者，预置的取点规则可以按五角形状取点，即以五角形的中心点为第一目标网点，取五角形各个顶点与第一目标网点相邻的至少4个网点为第二目标网点。

从高程异常网点查找表获取所述第一目标网点的高程异常以及上述N个第二目标网点的高程异常及对应经纬度信息。

该第一计算模块40，用于根据上述获取到的高程异常、经纬度信息及预置的多项式曲面公式计算得到所述多项式曲面的拟合参数值；

其中，预置的多项式曲面公式为：

f(x,y)＝a₀+a₁x+a₂y+a₃x²+a₄xy+a₅y²

上式中，f(x,y)表示高程异常，x表示纬度，y表示经度，a₀、a₁、a₂、a₃、a₄和a₅为所述多项式曲面的拟合参数。

根据以上获取到的各个网点的高程异常、经纬度信息的相关信息，利用以上预置的多项式曲面公式进行拟合，得到拟合参数a₀、a₁、a₂、a₃、a₄和a₅的数值。

将上述得到的拟合参数a₀、a₁、a₂、a₃、a₄和a₅的数值返回到上述的多项式曲面公式，即可得到带有具体拟合参数值的多项式曲面公式如下：

f(x,y)＝a₀+a₁x+a₂y+a₃x²+a₄xy+a₅y²

上式中，f(x,y)表示高程异常，x表示纬度，y表示经度，a₀、a₁、a₂、a₃、a₄和a₅为所述多项式曲面的拟合参数。

该第二计算模块50，用于根据所述经纬度信息、所述多项式曲面公式及所述拟合参数值计算得到所述经纬度信息的高程异常。

本发明提供的一种高程异常的测定方法及测定装置，通过获取设备当前的经纬度信息，根据获得的经纬度信息在预设的高程异常网点查找表中确定对应的第一目标网点，并将所述第一目标网点为作所述设备当前所处位置对应的网点；再从高程异常网点查找表获取所述第一目标网点的高程异常，以及距离所述第一目标网点预设间隔距离的N个第二目标网点的高程异常及对应经纬度信息；然后根据上述获取到的高程异常、经纬度信息及预置的多项式曲面公式计算得到所述多项式曲面的拟合参数值；最后根据所述经纬度信息、所述多项式曲面公式及所述拟合参数值计算得到所述经纬度信息的高程异常。通过本发明以上的技术方案，就可以在GPS测量过程中根据卫星定位到的测量点的经度纬度自动进行高程异常拟合，既解决了采用传统GPS-水准测量工程量大且在山区地区测量困难的问题，也避免了我国由于丘陵以及山区等地区重力场数据不足而导致求解高程异常算法复杂并且精度不高的难题，也能够简化网点查找表的网点精度。

以上参照附图说明了本发明的优选实施例，并非因此局限本发明的权利范围。本领域技术人员不脱离本发明的范围和实质内所作的任何修改、等同替换和改进，均应在本发明的权利范围之内。

Claims (9)

1.一种高程异常的测定方法，其特征在于，所述测定方法包括：

获取设备当前的经纬度信息；

根据上述获得的经纬度信息在预设的高程异常网点查找表中确定对应的第一目标网点；所述第一目标网点为所述设备当前所处位置对应的网点；

从所述预设的高程异常网点查找表获取所述第一目标网点的高程异常，以及距离所述第一目标网点预设间隔距离的N个第二目标网点的高程异常及对应经纬度信息；

根据上述获取到的高程异常、经纬度信息及预置的多项式曲面公式计算得到所述多项式曲面的拟合参数值；

根据所述经纬度信息、所述多项式曲面公式及所述拟合参数值计算得到所述经纬度信息的高程异常。

2.如权利要求1所述的测定方法，其特征在于，所述预设的高程异常网点查找表是以经纬度值的整数为网点坐标建立，各网点间隔距离根据具体的精度要求来预设。

3.如权利要求1所述的测定方法，其特征在于，所述根据上述获得的经纬度信息在预设的高程异常网点查找表中确定对应的第一目标网点包括：

将所述经纬度信息对应的经纬度值按预设的调整规则调整为整数；

将调整后的经纬度值所对应的网点设定为所述第一目标网点。

4.如权利要求3所述的测定方法，其特征在于，所述预设的调整规则包括以下取整规则之一：四舍五入、向上取整、向下取整。

5.如权利要求1所述的测定方法，其特征在于，所述根据上述获得的经纬度信息在预设的高程异常网点查找表中确定对应的第一目标网点包括：

将所述经纬度信息，与预设的高程异常网点查找表中的各网点的经纬度信息进行比对，将与所述经纬度信息最接近的对应网点确定为第一目标网点。

6.如权利要求3或5所述的测定方法，其特征在于，所述N个第二目标网点是以所述第一目标网点为中心，根据预设间隔距离按照预置的取点规则从高程异常网点查找表获取与所述第一目标网点相邻的网点来确定。

7.如权利要求6所述的测定方法，其特征在于，所述预置的取点规则能满足特定形状的边上覆盖至少N个第二目标网点。

8.如权利要求6所述的测定方法，其特征在于，所述预置的多项式曲面公式为：

f(x,y)＝a₀+a₁x+a₂y+a₃x²+a₄xy+a₅y²

上式中，f(x,y)表示高程异常，x表示纬度，y表示经度，a₀、a₁、a₂、a₃、a₄和a₅为所述多项式曲面的拟合参数。

9.一种高程异常的测定装置，其特征在于，所述测定装置包括：

第一获取模块，用于获取设备当前的经纬度信息；

确定模块，用于根据上述获得的经纬度信息在预设的高程异常网点查找表中确定对应的第一目标网点；所述第一目标网点为所述设备当前所处位置对应的网点；

第二获取模块，用于从高程异常网点查找表获取所述第一目标网点的高程异常，以及距离所述第一目标网点预设间隔距离的N个第二目标网点的高程异常及对应经纬度信息；

第一计算模块，用于根据上述获取到的高程异常、经纬度信息及预置的多项式曲面公式计算得到所述多项式曲面的拟合参数值；

第二计算模块，用于根据所述经纬度信息、所述多项式曲面公式及所述拟合参数值计算得到所述经纬度信息的高程异常。