基于全球重力模型的区域分割方法、系统、设备及介质

摘要

本发明提供的基于全球重力模型的区域分割方法，包括获取分割区域数据，读取全球重力模型的参数数据，根据分割区域经差、经度参数、分割区域纬差、纬度参数计算分割列数和分割行数，根据西南角大地坐标、分割区域经差、分割区域纬差计算出分割区域的西北角大地坐标、东北角大地坐标以及东南角大地坐标；读取全球重力模型文件中的重力异常数据；根据分割区域各角的大地坐标、分割列数、分割行数、纬度参数以及经度参数读取全球重力模型文件得到头数据，将头数据与重力异常数据均写入预设区域水准文件中，得到区域重力模型文件。本发明提供的基于全球重力模型的区域分割方法，可以得到较为精准的区域重力模型文件。

说明书

技术领域

本发明涉及测绘计算领域，尤其涉及基于全球重力模型的区域分割方法、系统、设备及介质。

背景技术

国内GNSS高程测量方法一般采用单点校正、平面拟合、曲面拟合等方式，对于点位分布不均匀的测区，高程测量精度很难保证。全球重力模型文件，由于文件过大，很难被APP直接使用，因此需要对全球重力模型文件进行区域分割才能快速的供GNSS高程测量使用，但是目前对于全球重力模型文件的区域分割均是采用人为预估性的对全球重力模型文件进行分割，导致分割的结果并不精准。

发明内容

为了克服现有技术的不足，本发明的目的之一在于提供基于全球重力模型的区域分割方法，其能解决目前对全球重力模型文件进行分割的结果并不精准的问题。

本发明的目的之二在于提供基于全球重力模型的区域分割系统，其能解决目前对全球重力模型文件进行分割的结果并不精准的问题。

本发明的目的之三在于提供一种电子设备，其能解决目前对全球重力模型文件进行分割的结果并不精准的问题。

本发明的目的之四在于提供一种计算机可读存储介质，其能解决目前对全球重力模型文件进行分割的结果并不精准的问题。

本发明的目的之一采用以下技术方案实现：

基于全球重力模型的区域分割方法，包括以下步骤：

获取区域分割范围数据，获取从全球重力模型文件中的进行区域分割的分割区域数据，所述分割区域数据包括分割区域的西南角大地坐标、分割区域经差以及分割区域纬差；

读取参数数据，读取全球重力模型的参数数据，所述参数数据包括全球重力模型总行数、全球重力模型总列数、纬度参数以及经度参数；

计算分割区域参数，根据分割区域经差以及经度参数计算出分割区域的分割列数，根据分割区域纬差以及纬度参数计算出分割区域的分割行数，根据西南角大地坐标、分割区域经差、分割区域纬差计算出分割区域的西北角大地坐标、东北角大地坐标以及东南角大地坐标；

读取重力异常数据，读取全球重力模型文件中的重力异常数据；

生成区域重力模型文件，根据西南角大地坐标、西北角大地坐标、东北角大地坐标、东南角大地坐标、分割列数、分割行数、纬度参数以及经度参数并以预设插值方式读取全球重力模型文件得到头数据，将所述头数据与所述重力异常数据均写入预设区域水准文件中，得到区域重力模型文件。

进一步地，所述经度参数为经度间隔，所述根据分割区域经差、经度参数计算出分割区域的分割列数具体为：根据分割区域经差、经度间隔计算得到分割区域的分割列数。

进一步地，所述纬度参数为纬度间隔，所述根据分割区域纬差、纬度参数计算出分割区域的分割行数具体为：根据分割区域纬差、纬度间隔计算得到分割区域的分割行数。

进一步地，所述读取参数数据具体为：从全球重力模型文件中位于分割区域的西北角大地坐标的位置循环读取全球重力模型文件中的重力异常数据。

进一步地，根据西南角大地坐标、分割区域经差、分割区域纬差计算出分割区域的西北角大地坐标、东北角大地坐标以及东南角大地坐标具体为：根据西南角大地坐标、分割区域纬差计算得到西北角大地坐标，根据西南角大地坐标、分割区域经差计算得到东南角大地坐标，根据西南角大地坐标、分割区域经差、分割区域纬差计算得到东北角大地坐标。

进一步地，所述全球重力模型文件为EMG2008全球重力模型文件中的1分格网文件。

本发明的目的之二采用以下技术方案实现：

基于全球重力模型的区域分割系统，包括获取模块、读取模块、计算模块以及生成模块，所述获取模块用于获取从全球重力模型文件中的进行区域分割的分割区域数据，所述分割区域数据包括分割区域的西南角大地坐标、分割区域经差以及分割区域纬差；所述读取模块用于读取全球重力模型的参数数据，所述参数数据包括全球重力模型总行数、全球重力模型总列数、纬度参数以及经度参数；所述读取模块还用于读取全球重力模型文件中的重力异常数据；所述计算模块用于根据分割区域经差、经度参数计算出分割区域的分割列数，根据分割区域纬差以及纬度参数计算出分割区域的分割行数，根据西南角大地坐标、分割区域经差、分割区域纬差计算出分割区域的西北角大地坐标、东北角大地坐标以及东南角大地坐标；所述生成模块用于根据西南角大地坐标、西北角大地坐标、东北角大地坐标、东南角大地坐标、分割列数、分割行数、纬度参数以及经度参数并以预设插值方式读取全球重力模型文件得到头数据，将所述头数据与所述重力异常数据均写入预设区域水准文件中，得到区域重力模型文件。

进一步地，所述读取模块包括读取参数数据单元以及读取重力异常数据单元，所述读取参数数据单元用于读取全球重力模型的参数数据，所述参数数据包括全球重力模型总行数、全球重力模型总列数、纬度参数以及经度参数；所述读取重力异常数据单元用于读取全球重力模型文件中的重力异常数据。

本发明的目的之三采用以下技术方案实现：

一种电子设备，包括：处理器；

存储器；以及程序，其中所述程序被存储在所述存储器中，并且被配置成由处理器执行，所述程序包括用于执行本申请中基于全球重力模型的区域分割方法。

本发明的目的之四采用以下技术方案实现：

一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行本申请中基于全球重力模型的区域分割方法。

相比现有技术，本发明的有益效果在于：本申请中的基于全球重力模型的区域分割方法，包括获取从全球重力模型文件中的进行区域分割的分割区域数据，分割区域数据包括分割区域的西南角大地坐标、分割区域经差以及分割区域纬差；读取全球重力模型的参数数据，参数数据包括全球重力模型总行数、全球重力模型总列数、纬度参数以及经度参数；根据分割区域经差以及经度参数计算出分割区域的分割列数，根据分割区域纬差以及纬度参数计算出分割区域的分割行数，根据西南角大地坐标、分割区域经差、分割区域纬差计算出分割区域的西北角大地坐标、东北角大地坐标以及东南角大地坐标；读取全球重力模型文件中的重力异常数据；根据西南角大地坐标、西北角大地坐标、东北角大地坐标、东南角大地坐标、分割列数、分割行数、纬度参数以及经度参数并以预设插值方式读取全球重力模型文件得到头数据，将头数据与重力异常数据均写入预设区域水准文件中，得到区域重力模型文件。通过计算得到分割区域的各个角的坐标、分割行数以及分割列数等参数精准读取全球重力模型文件中的头数据，从而得到较为精准的区域重力模型文件。

上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，并可依照说明书的内容予以实施，以下以本发明的较佳实施例并配合附图详细说明如后。本发明的具体实施方式由以下实施例及其附图详细给出。

附图说明

此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：

图1为本发明的基于全球重力模型的区域分割方法流程示意图。

具体实施方式

下面，结合附图以及具体实施方式，对本发明做进一步描述，需要说明的是，在不相冲突的前提下，以下描述的各实施例之间或各技术特征之间可以任意组合形成新的实施例。

如图1所示，本申请中的基于全球重力模型的区域分割方法，包括以下步骤：

获取区域分割范围数据，获取从全球重力模型文件中的进行区域分割的分割区域数据，所述分割区域数据包括分割区域的西南角大地坐标、分割区域经差以及分割区域纬差。本实施例中的全球重力模型文件为1分格网文件。

读取参数数据，读取全球重力模型的参数数据，所述参数数据包括全球重力模型总行数、全球重力模型总列数、纬度参数以及经度参数。在本实施例中，所述经度参数为经度间隔，所述纬度参数为纬度间隔。

计算分割区域参数，根据分割区域经差以及经度参数计算出分割区域的分割列数，根据分割区域纬差以及纬度参数计算出分割区域的分割行数，根据西南角大地坐标、分割区域经差、分割区域纬差计算出分割区域的西北角大地坐标、东北角大地坐标以及东南角大地坐标。具体为：根据分割区域经差、经度间隔计算得到分割区域的分割列数，分割列数＝分割区域经差/经度间隔；根据分割区域纬差、纬度间隔计算得到分割区域的分割行数，分割行数＝分割区域纬差/纬度间隔。根据西南角大地坐标、分割区域纬差计算得到西北角大地坐标，根据西南角大地坐标、分割区域经差计算得到东南角大地坐标，根据西南角大地坐标、分割区域经差、分割区域纬差计算得到东北角大地坐标。本实时例中，全球重力模型文件为1分格网文件，即一个格网点表示一个坐标点，假设西南角大地坐标为(AL

读取重力异常数据，读取全球重力模型文件中的重力异常数据；具体为：从全球重力模型文件中位于分割区域的西北角大地坐标的位置循环读取全球重力模型文件中的重力异常数据。

生成区域重力模型文件，根据西南角大地坐标、西北角大地坐标、东北角大地坐标、东南角大地坐标、分割列数、分割行数、纬度间隔以及经度间隔并以预设插值方式读取全球重力模型文件得到头数据，将所述头数据与所述重力异常数据均写入预设区域水准文件中，得到区域重力模型文件。

本申请中提供了基于全球重力模型的区域分割系统，包括获取模块、读取模块、计算模块以及生成模块，所述获取模块用于获取从全球重力模型文件中的进行区域分割的分割区域数据，所述分割区域数据包括分割区域的西南角大地坐标、分割区域经差以及分割区域纬差；所述读取模块用于读取全球重力模型的参数数据，所述参数数据包括全球重力模型总行数、全球重力模型总列数、纬度参数以及经度参数；所述读取模块还用于读取全球重力模型文件中的重力异常数据；所述计算模块用于根据分割区域经差、经度参数计算出分割区域的分割列数，根据分割区域纬差以及纬度参数计算出分割区域的分割行数，根据西南角大地坐标、分割区域经差、分割区域纬差计算出分割区域的西北角大地坐标、东北角大地坐标以及东南角大地坐标；所述生成模块用于根据西南角大地坐标、西北角大地坐标、东北角大地坐标、东南角大地坐标、分割列数、分割行数、纬度参数以及经度参数并以预设插值方式读取全球重力模型文件得到头数据，将所述头数据与所述重力异常数据均写入预设区域水准文件中，得到区域重力模型文件。本实施例中，所述读取模块包括读取参数数据单元以及读取重力异常数据单元，所述读取参数数据单元用于读取全球重力模型的参数数据，所述参数数据包括全球重力模型总行数、全球重力模型总列数、纬度参数以及经度参数；所述读取重力异常数据单元用于读取全球重力模型文件中的重力异常数据。

本申请中提供了一种电子设备，包括：处理器；

存储器；以及程序，其中所述程序被存储在所述存储器中，并且被配置成由处理器执行，所述程序包括用于执行本申请中基于全球重力模型的区域分割方法。

本申请中提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行本申请中基于全球重力模型的区域分割方法。

本申请中的基于全球重力模型的区域分割方法，包括获取从全球重力模型文件中的进行区域分割的分割区域数据，分割区域数据包括分割区域的西南角大地坐标、分割区域经差以及分割区域纬差；读取全球重力模型的参数数据，参数数据包括全球重力模型总行数、全球重力模型总列数、纬度参数以及经度参数；根据分割区域经差以及经度参数计算出分割区域的分割列数，根据分割区域纬差以及纬度参数计算出分割区域的分割行数，根据西南角大地坐标、分割区域经差、分割区域纬差计算出分割区域的西北角大地坐标、东北角大地坐标以及东南角大地坐标；读取全球重力模型文件中的重力异常数据；根据西南角大地坐标、西北角大地坐标、东北角大地坐标、东南角大地坐标、分割列数、分割行数、纬度参数以及经度参数并以预设插值方式读取全球重力模型文件得到头数据，将头数据与重力异常数据均写入预设区域水准文件中，得到区域重力模型文件。通过计算得到分割区域的各个角的坐标、分割行数以及分割列数等参数精准读取全球重力模型文件中的头数据，从而得到较为精准的区域重力模型文件。

以上，仅为本发明的较佳实施例而已，并非对本发明作任何形式上的限制；凡本行业的普通技术人员均可按说明书附图所示和以上而顺畅地实施本发明；但是,凡熟悉本专业的技术人员在不脱离本发明技术方案范围内，利用以上所揭示的技术内容而做出的些许更动、修饰与演变的等同变化，均为本发明的等效实施例；同时,凡依据本发明的实质技术对以上实施例所作的任何等同变化的更动、修饰与演变等，均仍属于本发明的技术方案的保护范围之内。