基于北斗的微形变预测方法、装置、设备及可读存储介质与流程

1.本发明涉及卫星监测技术领域，具体而言，涉及基于北斗的微形变预测方法、装置、设备及可读存储介质。


背景技术：

2.目前，基于位移形变建立的预测模型通常采用非线性预测方法，常用的经典模型主要包括时间序列分析、灰色gm(1,1)预测模型、卡尔曼滤波、回归分析模型、支持向量机、人工神经网络等，另外，关于混合模型的预测研究也是目前较为流行一种方法。
3.但是目前的位移形变监测预警系统往往仅针对某一个对象或某一类场景，系统的通用性不强，应对突发性变化能力较弱，一些新的监测方案虽然精度较高，但相应成本较高。
4.由于各应用环境下位移形变的复杂性，很难仅通过单个模型准确预测监测体的位移和演变过程，并且在实际的数据采集过程中，往往收集得到的数据会受到系统中噪声和干扰的影响，数据也可能受特殊情况的影响而导致出现缺失的情况。


技术实现要素：

5.本发明的目的在于提供一种基于北斗的微形变预测方法、装置、设备及可读存储介质，以改善上述问题。为了实现上述目的，本发明采取的技术方案如下：现有技术中采用常用的经典模型包括时间序列分析、灰色gm(1,1)预测模型、卡尔曼滤波、回归分析模型、支持向量机、人工神经网络等，然而只是仅针对某一个对象或某一类场景，系统的通用性不强，应对突发性变化能力较弱，一些新的监测方案虽然精度较高，但相应成本较高。
6.而本发明提供了一套通用性高、成本低、准确度高的基于北斗gnss微形变数据的最佳权重预测方法、设备及系统。
7.第一方面，本技术提供了一种基于北斗的微形变预测方法，包括：获取位移数据集，所述位移数据集包括至少两组位移数据，所述位移数据是通过卫星测量得到的所述目标的三维位移量，所述三维位移量为第二时刻相对于第一时刻的位移量；将所述位移数据集中的每个所述位移数据分别进行计算，分别得到每个所述位移数据的位移形变数据，所述位移形变数据为在所述位移数据集中通过预测算法模块提取出的数据;根据所述位移形变数据，建立预测模型,并通过所述预测模型计算出所述目标在第三时刻相对于第二时刻的位移量预测结果。
8.可选地，所述将位移数据集中的每个位移数据分别进行计算，分别得到位移形变数据,包括:通过箱型图分析方法对每个所述位移数据进行异常值判断；获取每个所述位移数据的异常值判断结果，若判断结果正常，则记做第一信息，若
判断结果不正常，则剔除；对判断结果正常的所述第一信息的空缺部分进行填补;获取填补后的所述第一信息,将所述第一信息记做位移形变数据。
9.可选地，所述根据所述位移形变的数据，建立预测模型，包括：将所述位移形变数据经过预设公式进行转化得到转化数据集，转化公式如下：其中，式中为所有位移原始数据的均值，为所有位移数据的标准差,为均值方差归一化的结果，为需要处理的数据；更新所述位移形变数据为所述转化数据集。
10.可选地，所述根据所述位移形变的数据，建立预测模型，包括：更新所述位移形变数据为所述转化数据集,根据所述转化数据集，建立第一模型，所述第一模型包括更新门单元（rt）和重置门单元（zt），所述第一模型的每个门单元和隐藏层的计算公式如下：层的计算公式如下：层的计算公式如下：层的计算公式如下：式中，为sigmoid激活函数，为双曲正切激活函数，为当前的输入，为当前输出，表示候选激活向量，为之前的输出，和是更新门单元的权重，和是重置门单元的权重，和是输出候选值，、和分别是更新门单元、重置门单元和输出候选值的偏置向量。
11.可选地,所述根据所述位移形变的数据，建立预测模型，还包括：更新所述位移形变数据为所述转化数据集,根据所述转化数据集，建立第二模型；通过所述第二模型，求解自回归过程（ar），计算公式如下：；式中，，，，表示ar的系数，为当前输出，，，，为之前的输入，为第t个时间序列的白噪声;根据所述自回归过程（ar）获取到移动平均过程（ma），计算如下：式中，式中，、，表示对应输入的误差，为随机干扰项序列，，，，表示对应的权重系数,为白噪声；根据所述自回归过程（ar）和所述移动平均过程（ma），求解得到自回归移动平均
（arma），计算公式如下：（arma），计算公式如下：表示误差，式中，，，，表示ar的系数，为当前输出，，，，为之前的输入，、，表示对应输入的误差，，，，表示对应的权重系数。
12.可选地,所述根据所述位移形变的数据，建立预测模型,并通过所述预测模型计算出第三时刻最终位移预测结果，包括：更新所述位移形变数据为所述转化数据集,根据所述转化数据集，建立第三模型；通过第三模型,求解回归系数，计算公式如下：式中，为样本点x变量的自变量矩阵，ki为角矩阵。
13.可选地,所述根据所述位移形变的数据，建立预测模型,并通过所述预测模型计算出第三时刻最终位移预测结果，还包括：根据所述位移形变的数据作为预测模型的输入信息,得到每个预测模型的位移数据预测结果,将所述预测结果和真实情况作差,得到每个基本模型的拟合误差;根据每个基本模型的拟合误差,计算得到每个基本模型对应的平均平方误差,其计算公式为:式中，表示每个基本模型的拟合误差；根据每个基本模型的平均平方误差,计算出每个基本模型的权值，其计算公式为：式中，表示每个基本模型的平均平方误差；根据每个基本模型的权值和所有基本模型,建立最佳权重线性组合模型，其最佳权重线性组合模型公式为：
式中，表示每个基本模型的预测数据，每个模型的权重表示为，的和为1；求解所述最佳权重线性组合模型，得到最终位移预测结果。
14.所述最佳权重预测模块，包括数据预处理阶段、构建四种基本模型和权值计算三个流程，可预测三维空间位移数据中的指定方向的值，其中数据预处理阶段包括异常数据处理、缺失数据处理和卡尔曼滤波降噪；构建基本模型包括岭回归、svm、gru和arima四种，并分别对近n期数据进行预测。
15.第二方面，本技术还提供了一种基于北斗的微形变预测装置，包括：获取模块、提取模块和计算模块，其中：获取模块:用于获取位移数据集，所述位移数据集包括至少两组位移数据，所述位移数据是通过卫星测量得到的所述目标的三维位移量，所述三维位 移量为第二时刻相对于第一时刻的位移量；提取模块:用于将所述位移数据集中的每个所述位移数据分别进行计算，分别得到每个所述位移数据的位移形变数据，所述位移形变数据为在所述位移数据集中通过预测算法模块提取出的数据;计算模块:用于根据所述位移形变数据，建立预测模型,并通过所述预测模型计算出所述目标在第三时刻相对于第二时刻的位移量预测结果。
16.第三方面，本技术还提供了一种基于北斗的微形变预测设备，包括：存储器，用于存储计算机程序；处理器，用于执行所述计算机程序时实现所述基于北斗的微形变预测方法的步骤。
17.第四方面，本技术还提供了一种可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现上述基于北斗的微形变预测方法的步骤。
18.本发明的有益效果为：由于微形变组合模型基于最小方差准则利用了最大价值的信息，该组合模型集成了所有单个模型的有用信息，平衡了各模型间的优缺点，因此误差比任何单个模型都低，避免了单个模型偶尔出现误差较大的情况，大大提高了单个模型的预测精度，拥有更强的泛化能力；并且通用性高、成本低、准确度高；解决了传统单一位移形变预测模型无法准确描述变形和演化规律的局限性。
19.本发明的其他特征和优点将在随后的说明书阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明实施例了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
附图说明
20.为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
21.图1为本发明实施例中所述的北斗的微形变预测方法流程示意图;图2为本发明实施例中所述的北斗的微形变预测装置结构示意图；图3为本发明实施例中所述的北斗的微形变预测设备结构示意图。
22.1、获取模块；11、第一计算单元；2、提取模块；3、计算模块；33、第一预测单元；34、第二预测单元；35、第三预测单元；36、第四预测单元；37、第五预测单元；801、处理器；802、存储器；803、多媒体组件；804、输入/输出(i/o)接口；805、通信组件。
具体实施方式
23.为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
24.应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。同时，在本发明的描述中，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
25.实施例1：现有技术中仅针对某一个对象或某一类场景，系统的通用性不强，应对突发性变化能力较弱，一些新的监测方案虽然精度较高，但相应成本较高。
26.由于各应用环境下位移形变的复杂性，很难仅通过单个模型准确预测监测体的位移和演变过程，并且在实际的数据采集过程中，往往收集得到的数据会受到系统中噪声和干扰的影响，数据也可能受特殊情况的影响而导致出现缺失的情况。
27.此发明可应用于地质灾害、人工构筑物等领域的监测及预警，主要包括数据预处理阶段和组合模型最佳权重计算两大流程：在数据预处理阶段的流程包括异常数据处理、缺失数据插值、使用卡尔曼滤波方法降噪。该最佳权重预测方法基本模型由岭回归、支持向量机、门控循环单元以及自回归移动平均构成，根据各基本模型与实际测量值的残差计算得到各模型的最佳权重而建立多个模型的组合模型，提高了预测准确度，解决了传统单一位移形变预测模型无法准确描述变形和演化规律的局限性。
28.本发明中主要包括北斗gnss监测网与最佳权重预测两大模块，卫星测量实则为北斗gnss监测网，包括1个gnss基准站和至少1个gnss监测站，其中gnss监测站包括gnss接收机、gnss接收天线、数据处理与控制系统、gprs模块等；数据通讯系统使用北斗通讯方式；供电系统使用太阳能供电，监测设备配置单晶硅太阳能板和铅酸蓄电池，可保证监测设备连续工作时间不少于20自然日。
29.所谓北斗的微形变为:在相邻两个间隔时间相等的时刻下,某一时刻相对于所述某一时刻的相邻时刻在三维空间中的位移数据的差叫做微形变,微形变就是细微的位移变化,可应用于地质灾害、建筑、桥梁等领域，表示为自然灾害或人工构筑物在两个不同时间点所相对应的两个位移之差,而北斗gnss高精度监测站终端主机与控制中心通信采用有线
或无线通信的方式，支持4g或5g网络，最终将采集的数据发送至数据中心的数据库中。
30.本实施例提供了一种基于北斗的微形变预测方法。
31.参见图1，图中示出了本方法包括步骤s100、步骤s200、步骤s300。
32.s100、获取位移数据集，位移数据集包括至少两组位移数据，位移数据是通过卫星测量得到的所述目标的三维位移量，三维位移量为第二时刻相对于第一时刻的位移量；可以理解的是，在本步骤中，目标可以是人工构筑物或者是自然灾害场景下的边坡，卫星测量为基于北斗gnss监测网而获得的高精度三维空间的位移数据，通过监测获取到人工构筑物或者是自然灾害发生的位移数据，将位移数据存入数据库中，而位移数据就是第二时刻相对于第一时刻的位移量。
33.第一时刻和第二时刻的间隔时间相等，相邻的两个时刻的间隔时间相等。
34.s200、将位移数据集中的每个位移数据分别进行计算，分别得到每个位移数据的位移形变数据，位移形变数据为在位移数据集中通过预测算法模块提取出的数据;可以理解的是，在本步骤中，将位移数据存入数据库中，通过预测算法模块从数据库中读取多组位移数据，并且对位移数据进行预处理，处理步骤依次为：异常数据处理、缺失数据处理和卡尔曼滤波降噪，具体方法如下，共三个步骤：s201、通过箱型图分析方法对每个位移数据进行异常值判断；获取每个位移数据的异常值判断结果，若判断结果正常，则记做第一信息，若判断结果不正常，则剔除；将将数据序列按照从小到大的顺序排列得到新的数据序列，其中位数记为m，异常值判定标准表示为：式中，u为上四分位数，是区间的中位数；l为下四分位数，是区间的中位数；为四分位距，其值通过u
ꢀ‑ꢀ
l计算得到；k为步长系数，通常取k为1.5；异常值的判断标准即为超过上界和下界的数据，将这些异常数据进行剔除。
35.s202、对判断结果正常的第一信息的空缺部分进行填补：根据缺失数据量大小选择合适的方法对缺失数据进行填补，当某段连续的缺失数据量超过10个时，采用所有n期数据的平均值进行填补；当某段连续的缺失数据量小于等于10个时，利用插值方法进行填补，其计算公式为：式中，和分别代表缺失数据端前后两点的数据，求区间内某一点x对应的y值；
s203、使用卡尔曼滤波方法对数据进行滤波降噪，其具体步骤如下：步骤s2031、指定初始条件，并初始化系统状态和误差协方差矩阵；步骤s2032、时间更新过程：计算先验状态估计和先验误差协方差；步骤s2033、状态更新过程：计算卡尔曼增益矩阵，并使用当前观测值更新系统状态并调整误差协方差；步骤s2034、如果算法没有结束，则返回步骤3.2，直到过滤算法结束，即到达最后一个观测值；执行完三个步骤后，获取填补后的第一信息,将第一信息记做位移形变数据。
36.s300、根据位移形变数据，建立预测模型,并通过预测模型计算出人工构筑物或者是自然灾害场景下的边坡在第三时刻相对于第二时刻的位移量预测结果；可以理解的是，在本步骤中，紧固预处理过后的多组位移形变数据进行均值方差归一化处理后，分别建立四种预测模型，包括机器学习模型与时间序列方法，具体四种预测模型分别为岭回归、svm、gru和arima；并且利用建立好的基本模型对近期多组数据进行预测，也就是通过预测模型计算出人工构筑物或者是自然灾害场景下的边坡在第三时刻相对于第二时刻的位移量预测结果。
37.将位移形变数据经过预设公式进行转化得到转化数据集，转化公式如下：其中，式中为所有位移原始数据的均值，为所有位移数据的标准差,为均值方差归一化的结果，为需要处理的数据。
38.更新位移形变数据为转化数据集,四种预测模型分别为：1.建立基本岭回归模型，根据转化数据集，建立模型；求解回归系数，计算公式如下：式中，为样本点x变量的自变量矩阵，ki为角矩阵。
39.2.建立基本svm模型，其公式表示为：式中，为拉格朗日乘子，为处理后的拉格朗日乘子，为偏执，核函数选择为高斯核函数，其公式为：
式中，代表内核宽度参数gamma，一定程度上决定了数据在特征空间的分布范围，svm还有一个重要的参数为正则项的惩罚系数c，用于平衡模型拟合的风险程度。
40.3.根据转化数据集，建立模型包括更新门单元（rt）和重置门单元（zt），计算公式如下：如下：如下：如下：式中，为sigmoid激活函数，为双曲正切激活函数，为当前的输入，为当前输出，表示候选激活向量，为之前的输出，和是更新门单元的权重，和是重置门单元的权重，和是输出候选值，、和分别是更新门单元、重置门单元和输出候选值的偏置向量。
41.4.根据转化数据集，建立模型，求解自回归过程（ar），计算公式如下：；式中，，，，表示ar的系数，为当前输出，，，，为之前的输入，为第t个时间序列的白噪声;根据自回归过程（ar）获取到移动平均过程（ma），计算如下：式中，式中，、，表示对应输入的误差，为随机干扰项序列，，，，表示对应的权重系数,为白噪声；根据自回归过程（ar）和移动平均过程（ma），求解得到自回归移动平均（arma），计算公式如下：算公式如下：表示误差，式中，，，，表示ar的系数，为当前输出，，，，为之前的输入，、，表示对应输入的误差，，，，表示对应的权重系数。
42.进而，根据位移形变的数据作为预测模型的输入信息,得到每个预测模型的位移数据预测结果,将预测结果和真实情况作差,得到每个基本模型的拟合误差;根据每个基本模型的拟合误差,计算得到每个基本模型对应的平均平方误差,其计算公式为:式中，表示每个基本模型的拟合误差；根据每个基本模型的平均平方误差,计算出每个基本模型的权值，其计算公式为：式中，表示每个基本模型的平均平方误差；每个基本模型的拟合误差表示为：式中，表示时序位移数据，表示每个基本模型的预测数据。
43.每个基本模型的平均平方误差可表示为：其中，m表示基本模型的个数，表示每个基本模型的拟合误差。
44.对于每个基本模型的权值，其计算公式为：式中，表示每个基本模型的平均平方误差。
45.根据每个基本模型的权值和所有基本模型,建立最佳权重线性组合模型，其最佳权重线性组合模型公式为：
式中，表示每个基本模型的预测数据，每个模型的权重表示为，的和为1；求解最佳权重线性组合模型，得到最终位移预测结果；当需要对新的位移数据进行预测时，则需要重新执行以上所有步骤。
46.实施例2：如图2所示，本实施例提供了一种基于北斗的微形变预测装置，参见图2装置包括：获取模块1、提取模块2和计算模块3，其中：获取模块1:用于获取位移数据集，位移数据集包括至少两组位移数据，位移数据是通过卫星测量得到的人工构筑物或者是自然灾害场景下的边坡的三维位移量，三维位移量为第二时刻相对于第一时刻的位移量；第一计算单元11：将位移数据集中的每个位移数据分别进行计算，分别得到位移形变数据,包括:通过箱型图分析方法对每个位移数据进行异常值判断；获取每个位移数据的异常值判断结果，若判断结果正常，则记做第一信息，若判断结果不正常，则剔除；对判断结果正常的第一信息的空缺部分进行填补;获取填补后的第一信息,将第一信息记做位移形变数据。
47.提取模块2:用于将位移数据集中的每个位移数据分别进行计算，分别得到每个位移数据的位移形变数据，位移形变数据为在位移数据集中通过预测算法模块提取出的数据;计算模块3:用于根据位移形变数据，建立预测模型,并通过预测模型计算出人工构筑物或者是自然灾害场景下的边坡在第三时刻相对于第二时刻的位移量预测结果；将更新位移形变数据为转化数据集，其中第一预测单元33：转化公式如下：其中，式中为所有位移原始数据的均值，为所有位移数据的标准差,为均值方差归一化的结果，为需要处理的数据；将更新位移形变数据为转化数据集，其中第二预测单元34：更新位移形变数据为转化数据集,根据转化数据集，建立第一模型，第一模型包括更新门单元（rt）和重置门单元（zt），第一模型的每个门单元和隐藏层的计算公式如下：（zt），第一模型的每个门单元和隐藏层的计算公式如下：（zt），第一模型的每个门单元和隐藏层的计算公式如下：（zt），第一模型的每个门单元和隐藏层的计算公式如下：
式中，为sigmoid激活函数，为双曲正切激活函数，为当前的输入，为当前输出，表示候选激活向量，为之前的输出，和是更新门单元的权重，和是重置门单元的权重，和是输出候选值，、和分别是更新门单元、重置门单元和输出候选值的偏置向量。
48.将更新位移形变数据为转化数据集，其中第三预测单元35：通过第二模型，求解自回归过程（ar），计算公式如下：；式中，，，，表示ar的系数，为当前输出，，，，为之前的输入，为第t个时间序列的白噪声;根据自回归过程（ar）获取到移动平均过程（ma），计算如下：式中，式中，、，表示对应输入的误差，为随机干扰项序列，，，，表示对应的权重系数,为白噪声；根据自回归过程（ar）和移动平均过程（ma），求解得到自回归移动平均（arma），计算公式如下：（ma），求解得到自回归移动平均（arma），计算公式如下：表示误差，式中，，，，表示ar的系数，为当前输出，，，，为之前的输入，、，表示对应输入的误差，，，，表示对应的权重系数。
49.将更新位移形变数据为转化数据集，其中第四预测单元36：通过第三模型,求解回归系数，计算公式如下：式中，为样本点x变量的自变量矩阵，ki为角矩阵。
50.将更新位移形变数据为转化数据集，其中第五预测单元37：根据位移形变的数据作为预测模型的输入信息,得到每个预测模型的位移数据预测结果,将预测结果和真实情况作差,得到每个基本模型的拟合误差;根据每个基本模型的拟合误差,计算得到每个基本模型对应的平均平方误差,其计算公式为:
式中，表示每个基本模型的拟合误差；根据每个基本模型的平均平方误差,计算出每个基本模型的权值，其计算公式为：式中，表示每个基本模型的平均平方误差；根据每个基本模型的权值和所有基本模型,建立最佳权重线性组合模型，其最佳权重线性组合模型公式为：式中，表示每个基本模型的预测数据，每个模型的权重表示为，的和为1；求解最佳权重线性组合模型，得到最终位移预测结果。
51.需要说明的是，关于上述实施例中的装置，其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。
52.相应于上面的方法实施例，本实施例中还提供了一种基于北斗的微形变预测设备，下文描述的一种基于北斗的微形变预测设备与上文描述的一种基于北斗的微形变预测方法可相互对应参照。
53.图3是根据示例性实施例示出的一种基于北斗的微形变预测设备800的框图。如图3所示，该基于北斗的微形变预测设备800可以包括：处理器801，存储器802。该基于北斗的微形变预测设备800还可以包括多媒体组件803，i/o接口804，以及通信组件805中的一者或多者。
54.其中，处理器801用于控制该基于北斗的微形变预测设备800的整体操作，以完成上述的基于北斗的微形变预测方法中的全部或部分步骤。存储器802用于存储各种类型的数据以支持在该基于北斗的微形变预测设备800的操作，这些数据例如可以包括用于在该基于北斗的微形变预测设备800上操作的任何应用程序或方法的指令，以及应用程序相关的数据，例如联系人数据、收发的消息、图片、音频、视频等等。该存储器802可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，例如静态随机存取存储器(static random access memory，简称sram)，电可擦除可编程只读存储器(electrically erasable programmable read
‑
only memory，简称eeprom)，可擦除可编程只读存储器(erasable programmable read
‑
only memory，简称eprom)，可编程只读存储器(programmable read
‑
only memory，简称prom)，只读存储器(read
‑
only memory，简称rom)，磁存储器，快闪存储器，磁盘或光盘。多媒体组件803可以包括屏幕和音频组件。其中屏幕例如可以是触摸屏，音频组件用于输出和/或输入音频信号。例如，音频组件可以包括一个麦克风，麦克风用于接收外部音频信号。所接收的音频信号可以被进一步存储在存储器802或通过通信组件805发
送。音频组件还包括至少一个扬声器，用于输出音频信号。i/o接口804为处理器801和其他接口模块之间提供接口，上述其他接口模块可以是键盘，鼠标，按钮等。这些按钮可以是虚拟按钮或者实体按钮。通信组件805用于该基于北斗的微形变预测设备800与其他设备之间进行有线或无线通信。无线通信，例如wi
‑
fi，蓝牙，近场通信(near fieldcommunication，简称nfc)，2g、3g或4g，或它们中的一种或几种的组合，因此相应的该通信组件805可以包括：wi
‑
fi模块，蓝牙模块，nfc模块。
55.在一示例性实施例中，基于北斗的微形变预测设备800可以被一个或多个应用专用集成电路(application specific integrated circuit，简称asic)、数字信号处理器(digitalsignal processor，简称dsp)、数字信号处理设备(digital signal processing device，简称dspd)、可编程逻辑器件(programmable logic device，简称pld)、现场可编程门阵列(field programmable gate array，简称fpga)、控制器、微控制器、微处理器或其他电子元件实现，用于执行上述的基于北斗的微形变预测方法。
56.在另一示例性实施例中，还提供了一种包括程序指令的计算机可读存储介质，该程序指令被处理器执行时实现上述的基于北斗的微形变预测方法的步骤。例如，该计算机可读存储介质可以为上述包括程序指令的存储器802，上述程序指令可由基于北斗的微形变预测设备800的处理器801执行以完成上述的基于北斗的微形变预测方法。
57.实施例4：相应于上面的方法实施例，本实施例中还提供了一种可读存储介质，下文描述的一种可读存储介质与上文描述的一种基于北斗的微形变预测方法可相互对应参照。
58.一种可读存储介质，可读存储介质上存储有计算机程序，计算机程序被处理器执行时实现上述方法实施例的基于北斗的微形变预测方法的步骤。
59.该可读存储介质具体可以为u盘、移动硬盘、只读存储器(read
‑
only memory，rom)、随机存取存储器(random access memory，ram)、磁碟或者光盘等各种可存储程序代码的可读存储介质。
60.以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
61.以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。