一种利用工厂模式实现坐标转换模型的参数解算方法与流程

1.本发明涉及坐标转换技术领域，具体涉及一种利用工厂模式实现坐标转换模型的参数解算方法。


背景技术：

2.近年来，全省国土资源地理空间数据开始全面采用2000国家大地坐标系,省级、市级、县级国土资源主管部门所属的存量数据，需要完成向2000国家大地坐标系的转换工作；采用相对独立平面坐标系的，也需要和2000国家大地坐标系建立有效联系；各级国土资源主管部门之间的数据上传、下发等全面采用2000国家大地坐标系。全面采用2000国家大地坐标系，就需要对现存数据进行坐标转换，这就离不开转换模型参数的解算。
3.对现有数据进行坐标转换时，根据不同的图幅范围以及转换需求，常用的转换方法有坐标模型转换与图幅改正量转换。如何在拥有足够的有效同名点(同名点就是同一个点在不同坐标系上的点)，并且根据实际情况选取了正确的转换模型后，如何进行转换模型的参数解算成为进行坐标转换时亟待解决的问题。


技术实现要素：

4.鉴于以上技术问题，本发明的目的在于提供及一种利用工厂模式实现坐标转换模型的参数解算方法，可实现转换模型的参数解算。
5.本发明采用以下技术方案：
6.一种利用工厂模式实现坐标转换模型的参数解算方法，包括以下步骤：
7.通过工厂类定义创建对象的公共接口；所述公共接口用于获取给定信息，所述给定信息包括待转换坐标系向大地坐标系进行转换采用的转换模型信息；所述创建对象的过程延迟到至少一个实体类进行；所述至少一个实体类存储有法方程；
8.根据所述采用的转换模型信息确定与所述转换模型信息关联的第一实体类；
9.所述第一实体类根据存储的法方程和所述转换模型信息计算所述转换模型的参数。
10.可选的，所述法方程满足以下公式：
11.x＝(b
t
b)-1
l
12.其中，b为同名点在待转换坐标系下的大地坐标值按照转换数学模型计算出来的系数矩阵，x为转换数学模型中多个参数组成的多行 1列的转换参数矩阵，l为转换数学模型计算出来的常数项。
13.可选的，所述法方程通过以下步骤得到：
14.获取转换模型的误差方程；利用最小二乘法的原理，通过间接平差法，对所述误差方程进行求解，得到转换模型的参数满足的法方程；所述误差方程满足以下公式：
15.v＝b
×
x-l；
16.其中，b为同名点在待转换坐标系下的大地坐标值按照转换数学模型计算出来的
系数矩阵，x为转换数学模型中多个参数组成的多行 1列的转换参数矩阵，l为转换数学模型计算出来的常数项。
17.可选的，所述实体类包括工具类，所述工具类用于实现矩阵的存储、矩阵的乘积、矩阵的转置与矩阵的逆的计算。
18.可选的，所述转换模型包括布尔莎七参数模型、二维七参数模型、三维七参数模型、三维四参数模型和二维四参数模型中的一种或多种。
19.可选的，当所述转换模型为布尔莎七参数模型时，所述转换模型的参数满足以下公式：
[0020][0021]
其中，b为同名点在待转换坐标系下的大地坐标值按照转换数学模型计算出来的系数矩阵，x为转换数学模型中多个参数组成的多行 1列的转换参数矩阵，l为转换数学模型计算出来的常数项； x1、y1、z1分别为坐标点转换前在原坐标系的坐标值， x2、y2、z2分别为坐标点转换后在目标坐标系的坐标值，t
x
、ty、tz、r
x
、ry、rz和d分别为转换模型的七参数。
[0022]
可选的，当所述转换模型为二维七参数模型时，所述转换模型的参数满足以下公式：
[0023][0024][0025][0026]
其中，b为同名点在待转换坐标系下的大地坐标值按照转换数学模型计算出来的系数矩阵，x为转换数学模型中多个参数组成的多行 1列的转换参数矩阵，l为转换数学模型计算出来的常数项；e2为第一偏心率平方，无量纲；m、n分别为地球基本元素，即m、n分别为子午线曲率和卯酉圈曲率半径，单位为米；b、l、δe、δl分别为点位纬度、经度、坐标点在两个坐标系下的纬度差、经度差，经纬度单位为弧度，经纬度差值单位为弧度秒；ρ
″‑
180
×
3600
/π，单位弧度秒；a、δa 分别为椭球长半轴和长半轴差，单位为米；f、δf分别为椭球扁率和扁率差，无量纲；t
x
、ty、tz、r
x
、ry、rb、d分别为转换模型的七参数。
[0027]
可选的，当所述转换模型为三维七参数模型时，所述转换模型的参数满足以下公式：
[0028][0029][0030][0031]
其中，b为同名点在待转换坐标系下的大地坐标值按照转换数学模型计算出来的系数矩阵，x为转换数学模型中多个参数组成的多行1 列的转换参数矩阵，l为转换数学模型计算出来的常数项；e2＝2f-f2为第一偏心率的平方，无量纲；为地球椭球子午圈曲率半径，单位为米；地球椭球卯酉圈曲率半径，单位为米；b、l、h、δb、δl、δh分别为点位的纬度、经度、大地高，以及其在两个坐标系下纬度差、经度差、大地高差，其中经纬度单位为弧度，其差值单位为弧度秒，大地高及其差值单位均为米；ρ
″
＝180
×
3600/π，单位为弧度秒；a、δa分别为椭球长半轴和长半轴差，单位为米；f、δf分别为椭球扁率和扁率差，无量纲；t
x
、ty、t
zrx
、ry、ra、d分别为转换模型的七参数。
[0032]
一种电子设备，其特征在于，包括：至少一个处理器，以及与所述至少一个处理器通信连接的存储器，其中，所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行所述的利用工厂模式实现坐标转换模型的参数解算方法。
[0033]
一种计算机存储介质，其上存储有计算机程序，所述计算机程序在被处理器执行时，实现所述的利用工厂模式实现坐标转换模型的参数解算方法。
[0034]
相比现有技术，本发明的有益效果在于：
[0035]
本发明通过工厂类定义创建对象的公共接口，所述创建对象的过程延迟到至少一个实体类进行；所述至少一个实体类存储有法方程；根据所述采用的转换模型信息确定与所述转换模型信息关联的第一实体类；所述第一实体类根据存储的法方程和所述转换模型信息计算所述转换模型的参数，可实现转换模型的参数解算，并且利用了工厂模式来管理不同转换模型的参数解算，可屏蔽模型解算代码的具体实现，保证了模型实现代码的安全性，降低了不同功能模块代码间的耦合性。
附图说明
[0036]
图1为本发明一实施例提供的一种利用工厂模式实现坐标转换模型的参数解算方法的流程示意图；
[0037]
图2为本发明一实施例提供的一种转换模型参数求解的工厂模式设计的示意图；
[0038]
图3为本发明一实施例提供的一种电子设备的结构示意图。
具体实施方式
[0039]
下面，结合附图以及具体实施方式，对本发明做进一步描述，需要说明的是，在不相冲突的前提下，以下描述的各实施例之间或各技术特征之间可以任意组合形成新的实施例：
[0040]
实施例一：
[0041]
请参照图1所示，图1示出了本发明的一种利用工厂模式实现坐标转换模型的参数解算方法，包括以下步骤：
[0042]
步骤s1:通过工厂类定义创建对象的公共接口；所述公共接口用于获取给定信息，所述给定信息包括待转换坐标系向大地坐标系进行转换采用的转换模型信息；所述创建对象的过程延迟到至少一个实体类进行；所述至少一个实体类存储有法方程；
[0043]
在坐标转换过程中，根据不同的图幅范围与需求，最常用的坐标转换模型分别有二维四参数模型、三维四参数模型、二维七参数模型、三维七参数模型以及布尔莎模型。而这些模型就是改正数法，理论基础是大地坐标微积分公式，采用广义大地坐标微积分公式直接求出大地坐标改正数，其基本思想就是要将基准转换和坐标系转换融合到一个模型中，从而建立起不同基准下的大地坐标间的直接关系。因此，这类转换模型的参数的求解，都可采用间接平差法，即利用最小二乘法进行求解。由于在利用最小二乘法进行求解的过程中，有相似的求解法方程，因此在程序的实现过程中，可利用工厂模式来管理不同转换模型的参数解算。
[0044]
需要说明的是，法方程是测量平差计算中的一个重要的方程，它是由平差函数模型线性化后的方程式，以及在vtpv＝min的条件下所导得的方程式合并而成的联立方程组。
[0045]
具体的，所述法方程满足以下公式：
[0046]
x＝(b
t
b)-1
l
[0047]
其中，b为同名点在待转换坐标系下的大地坐标值按照转换数学模型计算出来的系数矩阵，x为转换数学模型中多个参数组成的多行 1列的转换参数矩阵，l为转换数学模型计算出来的常数项。
[0048]
在本实施例中，所述法方程可通过以下步骤得到：
[0049]
获取转换模型的误差方程；利用最小二乘法的原理，通过问接平差法，对所述误差方程进行求解，得到转换模型的参数满足的法方程；所述误差方程满足以下公式：
[0050]
v＝b
×
x-l；
[0051]
其中，b为同名点在待转换坐标系下的大地坐标值按照转换数学模型计算出来的系数矩阵，x为转换数学模型中多个参数组成的多行 1列的转换参数矩阵，l为转换数学模型计算出来的常数项。
[0052]
具体的，所述法方程的计算原理如下：
[0053]
利用最小二乘平差法得出求解模型参数的通用法方程：
[0054]
需要说明的是，最小二乘法平差是在残差向量v和权矩阵p满足 v
t
pv＝min的条件下，求取测量值和参数的最佳估值，并进行精度估计的的理论和方法。
[0055]
用最小二乘法进行测量平差时，未知量估值的数学期望等于未知量的数学期望(估值无偏)，且估值的方差为最小，所获得的估值是最佳估值。
[0056]
在求取坐标转换参数时，已知量包括同名点坐标值、椭球参数，在采用间接平差法时，上述已知量可视为常数，在转换数学模型已知的前提下，其误差方程可视为：
[0057]
v＝b
×
x-l
[0058]
在该误差方程中，b为同名点在待转换坐标系下的大地坐标值按照转换数学模型计算出来的系数矩阵，x为转换数学模型中多个参数组成的多行1列的转换参数矩阵，l则为转换数学模型计算出来的常数项。
[0059]
按照最小二乘法的原理，上述误差方程中x必须满足vtpv＝min 的要求，又因为转换数学模型中每个参数互为独立量，故可用求函数自由极值的方式进行计算，将其转置后得到：
[0060]bt
pv＝0
[0061]
把上述式子代入误差方程，以便消除v，得到下列方程：
[0062]bt
pbx-b
t
pt＝0
[0063]
该式为转换数学模型进行转换参数求解时的间接平差方法的法方程。对其进行转换后可得到下列方程式：
[0064]
x＝(b
t
pb)-1
p
[0065]
利用上述方程式，即可解出待转换坐标系向目标坐标系进行转换时，转换数学模型中的转换参数的最小二乘解，其中在上述式子中， p为同名点坐标观测值的权阵，因为同名点的大地坐标值时采用等精度观测方法获得的，故p可视为单位矩阵，不对平差结果产生影响，即可得到下列方程式，即求解模型参数的最终通用法方程：
[0066]
x＝(b
t
b)-1
l
[0067]
上述得到的最终通用法方程即为转化模型的法方程：
[0068]
其中，b为同名点在待转换坐标系下的大地坐标值按照转换数学模型计算出来的系数矩阵，x为转换数学模型中多个参数组成的多行1列的转换参数矩阵，1则为转换数学模型计算出来的常数项。
[0069]
步骤s2：根据所述采用的转换模型信息确定与所述转换模型信息关联的第一实体类；
[0070]
在具体实施中传统的坐标转换方法，通常是不同的转换模型都实现一个类，但在无论在代码管理、程序拓展、调用便利上，都存在诸多弊端。
[0071]
本技术考虑到尽管坐标转换中转换模型不尽相同，但通过建立法方程x＝(b^tb)^(-1)l，作为通用的解算过程，可使用设计模式中的工厂模式来进行管理。
[0072]
应理解，工厂模式这种类型的设计模式属于创建型模式，该模式提供了一种创建对象的方法。在工厂模式种，创建对象时往往不会对调用者暴露创建逻辑，并且时通过一个共同的接口来指向新创建的对象，其意图在于定义一个创建对象的接口，让其子类自己决定实例化哪一个种类的工厂类，工厂模式使得创建的过程延迟到子类进行。其优点主要有
以下几点：1)一个调用者想创建一个对象，只要知道其名称就可以了；2)扩展性高，如果想增加一个产品，只要扩展一个实体类就可以；通过工厂类管理不同的实体类，增加新的产品就增加一个工厂类管理的实体类；3)屏蔽产品的具体实现，调用者只关心产品的接口。
[0073]
以上优点解决了关于实现、管理坐标转换参数解算的问题。
[0074]
请参照图2所示，图2示出了本发明的转换模型参数求解的工厂模式设计的示意图；首先创建一个公共接口和实现接口的实体类；然后需要定义一个工厂类。假设factorypatterndemo需要获取通过 solvingmodelfactory获取solvingmodel的对象，它将需要向 solvingmodel传递信息(转换模型信息，例如para7_2d、para7_3d、para7_bursa转换模型信息)，以便获取它所需要的对象类型。
[0075]
在具体实施中，具体实现转换模型参数求解代码的主要难点主要在于solvingmodel接口以及工厂类solvingmodelfactory的设计。解算时所使用的法方程x＝(b
t
b-1
l涉及矩阵及其几个基本运算，因此实现一个对应的工具类matrix也是必不可少的。之后，只需要继承 solvingmodel接口，实现不同模型的实体类。
[0076]
一、工具类matrix
[0077]
工具类matrix的实现，主要在于对矩阵的存储、矩阵的乘积、转置与矩阵的逆的实现，即可满足转换模型参数求解的需求。
[0078]
具体的，接口设计如下：
[0079]
1、solvingmodel接口
[0080]
求解坐标转换模型的参数，即使模型不同，都离不开通用法方程 x＝(b
t
b)-1
l的解算，因此无论是哪个坐标转换模型的参数解算，都需要实现对获取矩阵b、l的实现，以及根据矩阵
x
获取具体的参数集合。因此，solvingmodel接口的设计如下：
[0081]
[0082][0083]
2)工厂类solvingmodelfactory
[0084]
工厂类solvingmodelfactory主要作用为外部调用者提供获取 solvingmodel以及解算模型参数的公开方法。因此，工厂类solvingmodelfactory的代码如下：
[0085]
[0086]
[0087]
[0088][0089]
3)继承solvingmodel接口的实体类
[0090]
继承solvingmodel接口的实体类是针对某一模型实现解算算法的真正完成者，此时，只需要结合化简后的方程式，实现solvingmodel 接口中的方法即可。
[0091]
以布尔莎七参数模型为例，其实体类的代码设计如下：
[0092][0093]
具体的，布尔莎七参数模型实体类getlset实现代码如下：
[0094]
在具体实施中，后续若想对程序进行拓展，实现更多的转换模型的要求时，只需要创建新的继承solvingmodel接口的实体类，并在工厂类solvingmodelfactory中的getsolvingmodel方法中插入新的分支，无需对原有实现的转换模型代码进行任何改动，保证了原有实现模型代码的稳定性，同时使代码具有优秀的可拓展性。而对于调用者而言，他们无需关注模型解算的具体实现过程，只需要关心准备调用模型的名称、参与计算的同名点集合以及调用的接口，屏蔽了模型解算代码的具体实现，保证了模型实现代码的安全性，降低了不同功能模块代码间的耦合性。
[0095]
步骤s3:所述第一实体类根据存储的法方程和所述转换模型信息计算所述转换模型的参数。
[0096]
需要说明的是，实体类主要是作为数据管理和业务逻辑处理层面上存在的类别；它们主要在分析阶段区分实体类的主要职责是存储和管理系统内部的信息，它也可以有行为，甚至很复杂的行为，但这些行为必须与它所代表的实体对象密切相关。实体类是用于对必须存储的信息和相关行为建模的类。实体对象(实体类的实例)用于保存和更新一些现象的有关信息，例如：事件、人员或者一些现实生活中的对象。
[0097]
具体的，所述实体类包括工具类，所述工具类用于实现矩阵的存储、矩阵的乘积、矩阵的转置与矩阵的逆的计算。
[0098]
在具体实施中，所述转换模型包括布尔莎七参数模型、二维七参数模型、三维七参数模型、三维四参数模型和二维四参数模型中的一种或多种。
[0099]
具体的，1)当所述转换模型为布尔莎七参数模型时，根据法方程得到所述转换模型的参数的过程如下：
[0100]
已知转换公式为：
[0101][0102]
式中，x1、y1、z1分别为原坐标系的坐标值，x2、y2、z2分别为目 标坐标系的坐标值，
t
x
、ty、tz、r
x
、ry、rz和d为转换模型的七参数。
[0103]
为方便观察，我们可对上述公式进行转换，得出观测方程为：
[0104][0105]
由上述观测方式，结合3.1所述转换过程可解得对应 x＝(b
t
b)-1
l中的参数分别为：
[0106][0107]
即所述转换模型的参数满足以下公式：
[0108][0109]
其中，b为同名点在待转换坐标系下的大地坐标值按照转换数学 模型计算出来的系数矩阵，x为转换数学模型中多个参数组成的多行 1列的转换参数矩阵，l为转换数学模型计算出来的常数项； x1、y1、z1分别为坐标点转换前在原坐标系的坐标值， x2、y2、z2分别为坐标点转换后在目标坐标系的坐标值， t
x
、ty、tz、r
x
、ry、rz和d分别为转换模型的七参数。
[0110]
2)当所述转换模型为二维七参数模型时，根据法方程得到所述 转换模型的参数的过程如下：
[0111]
已知转换公式为：
[0112][0113]
式中，e2为第一偏心率平方，无量纲；m、n分别为地球基本元素， 即m、n分别为子午线曲率和卯酉圈曲率半径，单位为米； b、l、δb、δl分别为点位纬度、经度、坐标点在两个坐标系下的纬度 差、经度差，经纬度单位为弧度，经纬度差值单位为弧度秒； ρ
″
＝180
×
3600/π，单位弧度秒；a、δa分别为椭球长半轴和长半轴差， 单位为米；f、δf分别为椭球扁
率和扁率差，无量纲； t
x
、ty、tz、r
x
、ry、rz、d分别为转换模型的七参数。
[0114]
为方便观察，我们可对上述公式进行转换，得出观测方程为：
[0115][0116]
由上述观测方式，结合所述转换过程可解得对应x＝(b
t
b)-1
l中的 参数分别为：
[0117][0118][0119][0120]
即所述转换模型的参数满足以下公式：
[0121][0122][0123][0124]
其中，b为同名点在待转换坐标系下的大地坐标值按照转换数学 模型计算出来的系数矩阵，x为转换数学模型中多个参数组成的多行 1列的转换参数矩阵，l为转换数学模型计算出来的常数项；e2为第 一偏心率平方，无量纲；m、n分别为地球基本元素，即m、n分别为子 午线曲率和卯酉圈曲率半径，单位为米；b、l、δb、δl分别为点位纬 度、经度、坐标点在两个坐标系下的纬度差、经度差，经纬度单位为 弧度，经纬度差值单位为弧度秒；ρ
″
＝180
×
3600/π，单位弧度秒；a、δa 分别为椭球长半轴和长半轴差，单位为米；f、δf分别为椭球扁率和 扁率差，无量纲；t
x
、ty、tz、r
x
、ry、rz、d分别为转换模型的七参 数。
[0125]
3)当所述转换模型为三维七参数模型时，根据法方程得到所述 转换模型的参数的过程如下：
[0126]
已知转换公式为：
[0127][0128]
式中，e2＝2f-f2为第一偏心率的平方，无量纲； 为地球椭球子午圈曲率半径，单位为米； 地球椭球卯酉圈曲率半径，单位为米；b、l、h、δb、δl、δh 分别为点位的纬度、经度、大地高，以及其在两个坐标系下纬度差、 经度差、大地高差，其中经纬度单位为弧度，其差值单位为弧度秒， 大地高及其差值单位均为米；ρ
″
＝180
×
3600/π，单位为弧度秒；a、δa 分别为椭球长半轴和长半轴差，单位为米；f、δf分别为椭球扁率和 扁率差，无量纲；t
x
，ty，tz，r
x
，ry，rz，d分别为转换模型的七参数。
[0129]
为方便观察，我们可对上述公式进行转换，得出观测方程为：
[0130][0131]
由上述观测方式，结合所述转换过程可解得对应x＝(b
t
b)-1
l中的 参数满足以下公式：
[0132][0133][0134][0135]
其中，b为同名点在待转换坐标系下的大地坐标值按照转换数学模 型计算出来的系数矩阵，x为转换数学模型中多个参数组成的多行1 列的转换参数矩阵，l为转换数学模型计算出来的常数项；e2＝2f-f2为第一偏心率的平方，无量纲；为地球椭球子午圈曲 率半径，单位为米；地球椭球卯酉圈曲率半径，单位为 米；b、l、h、δb、δl、δh分别为点位的纬度、经度、大地高，以及其在两 个坐标系下纬度差、经度差、大地高差，其中经纬度单位为弧度，其 差值单位为弧度秒，大地高及其差值单位均为米；ρ
″
＝180
×
3600/π， 单位为弧度秒；a、δa分别为椭球长半轴和长半轴差，单位为米；f、δf 分别为椭球扁率和扁率差，无量纲；t
x
，ty，tz，r
x
，ry，rz，d分别为 转换模型的七参数。
[0136]
4)当所述转换模型为三维四参数模型时，根据法方程得到所述 转换模型的参数的过程如下：
[0137]
已知转换公式为：
[0138][0139]
式中，b0、l0为区域中心p0点的大地经、纬度，单位为弧度； t
x
、ty、tz、a为转换模型的四参数。
[0140]
为方便观察，我们可对上述公式进行转换，得出观测方程为：
[0141][0142]
由上述观测方式，结合所述转换过程可解得对应x＝(b
t
b)-1
l中的 参数分别为：
[0143][0144]
其中，b为同名点在待转换坐标系下的大地坐标值按照转换数学模 型计算出来的系数矩阵，x为转换数学模型中多个参数组成的多行1 列的转换参数矩阵，l为转换数学模型计算出来的常数项；b0、l0为区 域中心p0点的大地经、纬度，单位为弧度；t
x
、ty、tz、a为转换模型 的四参数。
[0145]
5)当所述转换模型为二维四参数模型时，根据法方程得到所述 转换模型的参数的过程如下：
[0146]
已知转换公式为：
[0147][0148]
为方便化简为观测方程，设定得到新方程： [0149]
对上述公式进行转换，得出观测方程为：
[0150][0151]
由上述观测方式，结合3.1所述转换过程可解得对应x＝(b
t
b)-1
l 中的参数分别为：
[0152][0153]
其中，b为同名点在待转换坐标系下的大地坐标值按照转换数学模 型计算出来的系数矩阵，x为转换数学模型中多个参数组成的多行1 列的转换参数矩阵，l为转换数学模型计算出来的常数项； t
x
、ty、a、m为转换模型的四参数。
[0154]
在上述实现过程中，通过工厂类定义创建对象的公共接口，所述 创建对象的过程延迟到至少一个实体类进行；所述至少一个实体类存 储有法方程；根据所述采用的转换模型信息确定与所述转换模型信息 关联的第一实体类；所述第一实体类根据存储的法方程和所述转换模 型信息计算所述转换模型的参数，可实现一种或多种转换模型的参数 解算，并且利用了工厂模式来管理不同转换模型的参数解算，屏蔽了 模型解算代码的具体实现，保证了模型实现代码的安全性，降低了不 同功能模块代码间的耦合性，从而可使实现过程的代码更加便于阅 读、管理、调用与拓展。
[0155]
实施例三：
[0156]
图3为本技术实施例提供的一种电子设备的结构示意图，在本申 请中可以通过图3所示的示意图来描述用于实现本技术实施例的本 发明一种利用工厂模式实现坐标转换模型的参数解算方法的电子设 备100。
[0157]
如图3所示的一种电子设备的结构示意图，电子设备100包括一 个或多个处理器102、一个或多个存储装置104，这些组件通过总线 系统和/或其它形式的连接机构(未示出)互连。应当注意，图3所示 的电子设备100的组件和结构只是示例性的，而非限制性的，根据需 要，所述电子设备可以具有图3示出的部分组件，也可以具有图3未 示出的其他组件和结构。
[0158]
所述处理器102可以是中央处理单元(cpu)或者具有数据处理能 力和/或指令执行能力的其它形式的处理单元，并且可以控制所述电 子设备100中的其它组件以执行期望的功能。
[0159]
所述存储装置104可以包括一个或多个计算机程序产品，所述计 算机程序产品可以包括各种形式的计算机可读存储介质，例如易失性 存储器和/或非易失性存储器。所述易失性存储器例如可以包括随机 存取存储器(ram)和/或高速缓冲存储器(cache)等。所述非易失性存 储器例如可以包括只读存储器(rom)、硬盘、闪存等。在所述计算机 可读存储介质上可以存储一个或多个计算机程序指令，处理器102可 以运行所述程序指令，以实现下文所述的本技术实施例中(由处理器 实现)的功能以及/或者其它期望的功能。在所述计算机可读存储介质 中还可以存储各种应用程序和各种数据，例如所述应用程序使用和/ 或产生的各种数据等。
[0160]
本发明还提供一种计算机存储介质，其上存储有计算机程序，本 发明的方法如果以软件功能单元的形式实现并作为独立的产品销售 或使用时，可以存储在该计算机存储介质中。基于这样的理解，本发 明实现上述实施例方法中的全部或部分流程，也可以通过计算机程序 来指令相关的硬件来完成，所述的计算机程序可存储于一计算机存储 介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施 例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机 程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间 形式等。所述计算机存储介质可以包括：能够携带所述计算机程序代 码的任何实体或装置、记录介质、u盘、移动硬盘、磁碟、光盘、计 算机存储器、只读存储器(rom，read-only memory)、随机存取 存储器(ram，random access memory)、电载波信号、电信信号 以及软件分发介质等。需要说明的是，所述计算机存储介质包含的内 容可以根据司法管辖区内立法和专利实践的要求进行适当的增减，例 如在某些司法管辖区，根据立法和专利实践，计算机存储介质不包括 电载波信号和电信信号。
[0161]
对本领域的技术人员来说，可根据以上描述的技术方案以及构 思，做出其它各种相应的改变以及形变，而所有的这些改变以及形变 都应该属于本发明权利要求的保护范围之内。