CN103136371A - 多源空间数据的剖分标识生成方法以及数据检索方法 - Google Patents

Abstract

本发明公开了多源空间数据的剖分标识生成方法以及数据检索方法，涉及地球空间信息组织、地理信息系统和遥感、测绘技术领域。其中剖分标识生成方法为：将全球划分为规则的网格，对网格进行编码；对于待处理的空间数据，找到包含其地理范围的网格，使用该网格的编码作为待处理空间数据的剖分标识；将待处理空间数据与其剖分标识相关联存储。数据检索方法为：依据上述剖分标识生成方法，对所有空间数据均进行剖分标识的生成；对于所需检索的空间数据，找到包含其地理范围的网格，使用该网格的编码作为检索号；找到与检索号一致的剖分标识，提取该剖分标识关联的空间数据即可得到所学检索的空间数据。本发明适用于空间数据的组织和检索。

Description

多源空间数据的剖分标识生成方法以及数据检索方法

技术领域

本发明涉及地球空间信息组织、地理信息系统和遥感、测绘技术领域，具体涉及一种自动监控并分析空间数据文件，进而读取空间坐标信息，生成全球统一网格标识编码，以及使用该标识对数据进行检索的方法，用于多源空间数据的关联与检索。

背景技术

地理信息系统GIS的应用，为社会经济稳定快速发展和人类进步发挥了巨大作用。

现在GIS和数据管理系统对空间数据均采用标识码ID进行唯一标识的，但由于使用单位或个人对现实世界了解程度不同、业务需求不同，对空间信息描述方法也多种多样，空间信息录入先后顺序也不一致，导致出现了同一空间数据在不同的GIS系统中描述各异、标识码ID不相同等语义不一致问题。例如，同样是北京市范围的数据，测绘系统中标识为图幅编号、遥感系统中标识为景Path/Row编号、社会统计标识为行政区划名、Google Earth系统的标识为瓦片的行列号等等。事实上，同样的问题也存在于同一信息管理系统的不同作业环节之间。正是由于纷繁而彼此独立的标识方式，造成多源空间数据关联的困难。

而且，这些空间信息标识码ID是一种“静态”的标识码，不能表达信息的多尺度的性质，也不具备表达空间位置的能力。这导致用户在进行数据共享时，必须查找各自的数据字典或者元数据才可将同一实体对应起来，势必影响数据共享的效率。

由此可见，探索一种适用于多源空间数据关联的统一标识生成方法，从而为解决上述问题提供途径，成为本项发明的应用背景与出发点。

发明内容

有鉴于此，本发明提供了一种用于多源空间数据关联的剖分标识生成方法，该方法统一了多源空间数据的标识方式，为后续提高多源空间数据的查询检索效率提供有力的支持。

为达到上述目的本发明的技术方案为：将全球划分为规则的网格，对网格进行编码；对于待处理的空间数据，确定所述空间数据标识的地理范围，找到包含所述地理范围的网格，使用该网格的编码作为所述待处理空间数据的剖分标识；将所述待处理空间数据与其剖分标识相关联存储。

进一步地，使用GeoSOT网格对全球进行划分，该方法包括如下步骤：

步骤1，对于待处理的空间数据，获取其所标识的地理空间范围的中心点的经纬度坐标P（L，B），L表示经度，B表示纬度；

若待处理空间数据所标识的是不规则的地理空间范围，则取该范围的外包矩形，找到外包矩形的中心点，获取所述中心点的经纬度坐标P（L，B）；

步骤2，设定对待处理的空间数据进行标识所使用的GeoSOT网格的层级m；

步骤3，获取所述中心点所在的m层级的GeoSOT网格的GeoSOT编码作为待处理空间数据的剖分标识，将所述待处理空间数据与其剖分标识相关联存储。

进一步地，步骤3中，中心点所在的m层级的GeoSOT网格的GeoSOT编码的获取方法具体如下：

步骤301，以步骤1中所获得的中心点的经纬度坐标P（L，B）为依据，通过设定的规则获得中心点所在的m层级的GeoSOT网格的GeoSOT编码的第1位；

所述中心点P（L，B），前缀为E或N的经纬度值表示为正数，前缀为W或S的经纬度值表示为负数；所述设定的规则为：

当L>0且B>0时，第1位为0；

当L<0且B>0时，第1位为1；

当L>0且B<0时，第1位为2；

当L<0且B<0时，第1位为3；

步骤302，令所述中心点的经纬度坐标P（L，B）=（L_D°L_M′L_S″，B_D°B_M′B_S″），其中L_D、L_M、L_S分别为L的度、分、秒的数值，B_D、B_M、B_S分别为B的度、分、秒的数值，使用下式将L和B分别转换为GeoSOT坐标值L_GeoSOT和B_GeoSOT：

L_GeoSOT=L_D×64×64+L_M×64+L_S；

B_GeoSOT=B_D×64×64+B_M×64+B_S；

步骤302，将步骤2中所获得的层级为m的GeoSOT网格的尺寸记为CellSize，将L_GeoSOT和B_GeoSOT分别除以CellSize，得到的结果转换为二进制数，若得到的二进制数不足m-1位，则在前面补零使其满m-1位，最终获得经纬度GeoSOT坐标值的二进制数；

步骤303，对于步骤302中所获得的经纬度GeoSOT坐标值的二进制数按位计算如下：将纬度的GeoSOT坐标值的二进制数每一位的2倍加上经度的GeoSOT坐标值的二进制数对应位，获得四进制编码；

步骤304，组合步骤301所获得的P的编码第1位以及步骤303中所获得的四进制编码，即为最终P(L，B)的GeoSOT编码，将P的GeoSOT编码作为P所属空间数据的剖分标识。

本发明还提供了一种数据检索方法，具体为：依据如权利要求1所述的多源空间数据的剖分标识生成方法，对所有的空间数据均进行剖分标识的生成；

对于所需检索的空间数据，确定其所标识的地理范围，找到包含该地理范围的网格，使用该网格的编码作为检索号；使用所述检索号与存储的剖分标识进行匹配，找到与检索号完全一致的剖分标识，提取该剖分标识关联的空间数据，即获得所需检索的空间数据。

进一步地，使用GeoSOT网格对全球进行划分，该方法包括如下步骤：

步骤一、对于空间数据，获取其所标识的地理空间范围的中心点的经纬度坐标；

若空间数据所标识的是不规则的地理空间范围，则取该范围的外包矩形，找到外包矩形的中心点，获取该中心点的经纬度坐标；

步骤二、设定对所述空间数据进行标识所使用的GeoSOT网格的层级m；

步骤三、获取上述中心点所在的m层级的GeoSOT网格的GeoSOT编码作为空间数据的剖分标识，获取所有空间数据的剖分标识，并将各空间数据与其剖分标识相关联存储；

步骤四、对于所需检索的地理空间范围内的空间数据，获取该地理空间范围的中心点的经纬度坐标；

若所需检索的空间数据所标识的是不规则的地理空间范围，则取该范围的外包矩形，找到该外包矩形的中心点，获取该中心点的经纬度坐标；

步骤五、设定对所需检索的空间数据进行剖分关联所使用的GeoSOT网格的层级n；

步骤六、获取所述步骤五中获得的中心点所在的n层级的GeoSOT网格的GeoSOT编码作为检索号；

步骤七、使用所述检索号与存储的剖分标识进行匹配，找到与检索号完全匹配一致的剖分标识，提取所述剖分标识关联的空间数据，即为所需检索的空间数据。

有益效果：

1、本发明所提供的剖分标识生成方法，基于多源空间数据所属的地理区域建立了数据标识，从而以多源空间数据的基本属性为出发点实现了标识方式的统一，能够实现对不同GIS系统中多源空间数据的关联；

2、本发明同时提出了一种数据检索方法，即使用GeoSOT全球剖分方案所获得的剖分标识，将多源空间数据及其剖分标识相关联存储，直接使用剖分标识就能够检索到该剖分标识所关联的空间数据，从而使得多源空间数据的查询检索效率大大提高。

附图说明

图1为剖分标识与空间数据的关联关系示意图；

图2为传统的同一区域不同类型空间数据的标识方法；

图3为本发明基于剖分标识的同一区域不同类型数据的标识关联示意图。

具体实施方式

下面结合附图并举实施例，对本发明进行详细描述。

本发明一种用于多源空间数据关联的剖分标识生成方法的基本思想是：将全球划分为规则的网格，对网格进行编码；对于待处理的空间数据，确定该空间数据标识的地理范围，找到包含上述地理范围的网格，使用该网格的编码即为上述待处理空间数据的剖分标识，将上述待处理空间数据与剖分标识相关联存储。

本发明通过将全球划分为规则的网格，并对网格进行唯一编码，则空间数据相对于其所属地理区域就可以使用该编码进行统一标识，该方法让原本相互独立的多源空间数据具有了统一的标识基准，使多源空间数据建立基于空间区域的关联关系，从而实现了跨部门同一区域的多源空间数据的关联搜索。

结合以上基本思想，本实施例使用GeoSOT全球剖分方案中对应的GeoSOT网格对全球进行规则划分，其中GeoSOT全球剖分方案参见北京大学提出的专利申请：“一种统一现有经纬度剖分网格的方法”（公开号为CN102609525，申请日为2012年2月10日），该专利申请公开了一种GeoSOT地理网格设计方案，用于解决全球地理空间剖分和标识问题。该方案采用全四叉树递归剖分，将地球表面空间从全球至厘米级共进行了32级剖分，每个GeoSOT剖分层级均有其对应大小的GeoSOT网格，GeoSOT网格上下层级之间的面积之比是1/4。该方案对GeoSOT网格进行编码所产生的GeoSOT编码有四进制、一维二进制、二维二进制和十进制等四种形式。

使用GeoSOT网格对全球进行划分时，本方法包括如下步骤：

步骤1，获取待处理空间数据所标识的地理空间范围的中心点的经纬度坐标P（L，B），L表示经度，B表示纬度；

若待处理空间数据所标识的是不规则的地理空间范围，则取该范围的外包矩形，找到外包矩形的中心点，获取所述中心点的经纬度坐标P（L，B）。

例如，空间数据data的空间范围为（E115.1477°-E118.2757°，N39.4118°-N41.2402°），空间分辨率r=240m。该数据中心点P（L，B）的经纬度坐标为（E116.71167°，N40.32682°），其中L=116.71167，B=40.32682。则L和B的度分秒D°M′S″格式分别为（E116°42′42″，N40°19′37″）。在经纬度的表示中，通常将前缀为E或N的经纬度值表示为正数，前缀为W或S的经纬度值表示为负数。

步骤2，设定对待处理的空间数据进行剖分关联所使用的GeoSOT网格的层级m。

对于空间数据，其具有包含地理空间范围信息的属性文件或者元数据信息，在本实施例中可以根据待处理空间数据的属性文件或者元数据信息，获得空间数据所标识的的地理空间范围，同时获取空间分辨率或者比例尺。

根据空间分辨率或者比例尺信息设定待关联的空间数据对应GeoSOT网格的层级数m，并获得层级m的GeoSOT网格大小CellSize。

当待关联空间数据为影像数据时，获得该影像数据的分辨率r，其中r以度为单位，在实际的应用中可将空间分辨率的其他单位转换为度。

在这种情况下，用户会首先根据应用业务的需求指定GeoSOT网格对应的影像其像素大小为2ⁿ×2ⁿ，使用r以及上述像素大小，计算理想GeoSOT网格的尺寸，根据已知的GeoSOT网格的经纬度范围（如表1），找到最接近理想GeoSOT网格尺寸的GeoSOT网格，其所属层级即为层级m。

假设欲获得空间分辨率r为240米的空间数据Data，将其单位米转换为单位度即r=0.0002925°。假设用户指定GeoSOT网格对应的影像像素大小为512×512像素，则理想GeoSOT网格的边长尺寸应为r×512=0.14976°=8.9856′，根据已知的每个层级GeoSOT网格大小，如表1，可获知8.9856′在第11层级的16′以及第12层级的8′之间，而且更接近8′，因此判断第12层级的GeoSOT网格更接近理想GeoSOT网格，因此本实施例中选择第12层级的GeoSOT网格。那么，网格的大小CellSize=8×64″=512″。

表1各层级GeoSOT网格的经纬度范围

层级	网格大小	层级	网格大小	层级	网格大小
						0	512°	11	16′	22	0.5″
1	256°	12	8′	23	0.25″
						2	128°	13	4′	24	0.125″
3	64°	14	2′	25	0.0625″
						4	32°	15	1′	26	0.03125″
5	16°	16	32″	27	0.015625″
						6	8°	17	16″	28	0.0078125″
7	4°	18	8″	29	0.00390625″
						8	2°	19	4″	30	0.001953125″
9	1°	20	2″	31	0.0009765625″
						10	32′	21	1″	32	0.00048828125″

当待关联空间数据为矢量时，则获取矢量的比例尺信息，按照如表2所示的地图比例尺与GeoSOT网格对应关系获取层级level，同时根据表1获取层级level的GeoSOT网格大小CellSize。

表2地图比例尺与GeoSOT网格对应关系

地图比例尺	网格大小	层级
			1:100万	2°	8
1:50万	1°	9
			1:25万	2′	14
1:10万	2′	14
			1:5万	1′	15
1:2.5万	2″	20
			1:1万	1″	21
1:5000	0.5″	22

步骤3，获取所述中心点所在的m层级的GeoSOT网格的GeoSOT编码即为P所属空间数据的剖分标识，将上述待处理空间数据与剖分标识相关联存储。

在该步骤中，可以根据前述“一种统一现有经纬度剖分网格的方法”的专利中所公开的方法对中心点所在的m层级的GeoSOT网格进行编码，从而获得GeoSOT编码作为P所属空间数据的剖分标识，同时也可以使用如下具体方法获得编码：

步骤301，以步骤1中所获得的中心点的经纬度坐标P（L，B）中，L和B的所指代的经纬度信息为依据，通过设定的规则获得中心点所在的m层级的GeoSOT网格的GeoSOT编码的第1位剖分编码值；

经纬度坐标空间的原点为本初子午线与赤道的交点，因此对于中心点的经纬度，前缀为E或N的经纬度值表示为正数，前缀为W或S的经纬度值表示为负数。

针对本实施例所使用GeoSOT全球地理剖分编码方案，其GeoSOT坐标是基于经纬度坐标空间定义的。故对于GeoSOT坐标来说，其是以本初子午线与赤道的交点为原点，因此纬度坐标P（L，B）中L与B的前缀所表示的正负关系对于GeoSOT坐标依旧适用。

GeoSOT全球地理剖分方案是使用GeoSOT网格对全球地理进行全四叉树递归剖分，其剖分体系由32级网格构成，网格上下级别之间的面积之比大致都为4:1，是均匀变化的。在使用GeoSOT网格对全球地理进行全四叉树递归剖分时，首先应当是将全球地理剖分为4个一级GeoSOT网格，对这4个一级GeoSOT网格进行编码，则对于本实施例来说，确定了P点的经纬度坐标之后，即可根据设定的规则确定其GeoSOT编码的第1位，其中设定的规则为：

当L>0且B>0时，第1位剖分编码值为0；

当L<0且B>0时，第1位剖分编码值为1；

当L>0且B<0时，第1位剖分编码值为2；

当L<0且B<0时，第1位剖分编码值为3。

步骤302，令所述中心点的经纬度坐标P（L，B）=（L_D°L_M′L_S″，B_D°B_M′B_S″），其中L_D、L_M、L_S分别为L的度、分、秒的数值，B_D、B_M、B_S分别为B的度、分、秒的数值，使用下式将L和B分别转换为秒级数值：

L_GeoSOT=L_D×64×64+L_M×64+L_S；

B_GeoSOT=B_D×64×64+B_M×64+B_S。

其中L_GeoSOT和B_GeoSOT记为GeoSOT坐标值。

本实施例所使用GeoSOT坐标，是基于经纬度坐标空间定义的，但是将经纬度坐标空间做了三次扩展：将180°×360°空间扩展到512°×512°；将每度的60'空间扩展到64'；将每分的60"空间扩展到64"，仍旧保持整数“度、分、秒”量纲。因此将经纬度坐标P（L，B）的度、分、秒转换为秒级数值时，度乘以64×64，分乘以64从而实现转换。

在本实施例中，其值分别为L_GeoSOT=477866、B_GeoSOT=165093。

步骤302，将步骤2中所获得的层级为m的GeoSOT网格的尺寸记为CellSize，将L_GeoSOT和B_GeoSOT分别除以CellSize，得到的结果分别转换为二进制数，若得到的二进制数不足m-1位，则在前面补零使其满m-1位，最终获得经纬度GeoSOT坐标值的二进制数。

本步骤中，已知中心点P的经纬度坐标的秒级数值，CellSize为m层级GeoSOT网格的尺寸，也是以秒作为单位，因此二者相除就可以得到中心点P在m层级的哪个GeoSOT网格中，同时除法运算的结果也是GeoSOT网格的一种编码，是十进制编码。通过该除法运算得到的是二维编码。

本实施例中，希望采用一维四进制数来表达中心点P的GeoSOT编码。由于在GeoSOT全球剖分方案中，对GeoSOT网格进行编码时，按照递归层级一级一级地进行编码，该GeoSOT网格所处的GeoSOT层级越高、GeoSOT网格越细，获得的编码越长。对于m层级的GeoSOT网格，其四进制编码为m位，而在步骤301中已经获得了第1位的编码数值，还需要m-1位四进制编码。为得到m-1位的四进制编码，需将上述得到的二维十进制编码（L_GeoSOT，B_GeoSOT）分别转换为m-1位二进制编码，不足m-1补位0。

在本实施例中，对L_GeoSOT=477866和B_GeoSOT=165093除以步骤3中所获得的12层级的GeoSOT网格大小CellSize=512″，将结果转换为二进制，即

(L_GeoSOT/Cellsize)₂=（477866/512）₁₀=（1110100101）₂

(B_GeoSOT/Cellsize)₂=（165093/512）₁₀=（101000010）₂

由于12层级的GeoSOT网格的四进制编码应当为12位，而在步骤301中已经获得了第1位的编码数值，还需要11位。为得到11位的四进制编码，需将上述得到的二维二进制编码补位为11位，本实施例中对于上述得到的二进制结果前面补0，得到：

(L_GeoSOT/Cellsize)₂=（1110100101）₂=（01110100101）₂

(B_GeoSOT/Cellsize)₂=（101000010）₂=（00101000010）₂

步骤303，将步骤302中所获得的经纬度GeoSOT坐标值的二进制数按位计算如下：将纬度的GeoSOT坐标值的二进制数的每一位的2倍加上经度的GeoSOT坐标值的二进制数对应位，获得四进制编码。

假设步骤302中B_GeoSOT/Cellsize的二进制数值表示为b₁b₂...b_m-1，L_GeoSOT/Cellsize的二进制数值表示为l₁l₂...l_m-1，则采用下述方式四进制GeoSOT坐标值为(P_GeoSOT)₄=p₁p₂...p_m-1。二进制与四进制之间的关系为，p_i=b_i×2+l_i，1≤i≤m-1。

在本实施例中，该步骤转换结果为：(P_GeoSOT)₄=01312100121。

步骤304，组合步骤301所获得的P的第1级和第2级剖分编码值以及将步骤6中所获得的四进制编码，即为最终P(L，B)的GeoSOT编码，将P的GeoSOT编码作为P所属空间数据的剖分标识。在实际的应用中可以将最终P(L，B)的GeoSOT编码进行转换，如可以转换为十进制、二进制一维等编码格式。

在本实施例中，将步骤303的结果01312100121追加在步骤301所获得的第一位编码0的后面，得到001312100121——即为空间数据关联的剖分标识。

将上述待处理空间数据与其剖分标识相关联存入数据库中，即将步骤3中计算出的剖分编码001312100121和空间数据相关联，使多源空间数据按照地理位置和区域建立关联关系，如图1所示。以上方案实现了空间数据关联的剖分标识生成，从技术方案层面建立了剖分标识与空间数据之间良好的关联关系，适合进行关联检索，这是现有技术所不具备的。

基于上述剖分标识的生成方案，本发明同时提供了一种空间数据的检索方法，具体为：依据上述方案所提供的多源空间数据的剖分标识生成方法，对所有的空间数据均进行剖分标识的生成；

对于所需检索的地理空间范围内的空间数据，找到包含该地理范围的网格，使用该网格的编码作为检索号；使用检索号与存储的剖分标识进行匹配，找到与检索号完全匹配一致的剖分标识，提取所述剖分标识关联的空间数据，即为所需检索的空间数据。

使用GeoSOT网格对全球进行划分时，该方法包括如下步骤：

步骤一、对于空间数据，获取其所标识的地理空间范围的中心点的经纬度坐标；

若空间数据所标识的是不规则的地理空间范围，则取该范围的外包矩形，找到外包矩形的中心点，获取所述中心点的经纬度坐标；

步骤二、设定对空间数据进行标识所使用的GeoSOT网格的层级m；

步骤三、获取所述中心点所在的m层级的GeoSOT网格的GeoSOT编码作为待处理空间数据的剖分标识，具体做法可参照剖分标识生成方法方案中的步骤301～304。通过上述步骤一至步骤三的操作，获取所有空间数据的剖分标识，并将各空间数据与其剖分标识相关联存储；

步骤四、对于所需检索的地理空间范围内的空间数据，获取该地理空间范围的中心点的经纬度坐标；

若所需检索的空间数据所标识的是不规则的地理空间范围，则取该范围的外包矩形，找到该外包矩形的中心点，获取该中心点的经纬度坐标；

步骤五、设定对所需检索的空间数据进行剖分关联所使用的GeoSOT网格的层级n；

步骤六、对于步骤四中获得的中心点，找到其所在的n层级的GeoSOT网格的GeoSOT编码作为检索号；

步骤七、使用所述检索号与存储的剖分标识进行匹配，找到与检索号完全匹配一致的剖分标识，提取所述剖分标识关联的空间数据，即为所需检索的空间数据。

本发明所阐述的基于剖分标识的关联存储与传统方法之间的差别可以从图2和图3对比中说明：

图2表示传统的同一区域不同类型空间数据的标识方法，尽管这些数据表达的是同一个区域，但由于数据针对不同部门不同应用，因此它的编码多样且彼此独立，互不相通。

图3表示本发明基于剖分标识的同一区域不同类型数据的标识关联示意图，本方法使得多源空间数据具有了统一的数据标识基础，让同一空间区域的多源空间数据具有了相同区位编码。因此，基于GeoSOT网格的统一剖分标识实现对多源空间数据的关联。

综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明所阐释原理的前提下，所做的任何修改、等同替换、局部应用等，均应包含在本发明的保护范围之内。

Claims (7)

1.多源空间数据的剖分标识生成方法，其特征在于：将全球划分为规则的网格，对网格进行编码；

对于待处理的空间数据，确定所述空间数据标识的地理范围，找到包含所述地理范围的网格，使用该网格的编码作为所述待处理空间数据的剖分标识；

将所述待处理空间数据与其剖分标识相关联存储。

2.如权利要求1所述的多源空间数据的剖分标识生成方法，其特征在于，使用GeoSOT网格对全球进行划分，该方法包括如下步骤：

步骤1，对于待处理的空间数据，获取其所标识的地理空间范围的中心点的经纬度坐标P（L，B），L表示经度，B表示纬度；

若待处理空间数据所标识的是不规则的地理空间范围，则取该范围的外包矩形，找到外包矩形的中心点，获取所述中心点的经纬度坐标P（L，B）；

步骤2，设定对待处理的空间数据进行标识所使用的GeoSOT网格的层级m；

步骤3，获取所述中心点所在的m层级的GeoSOT网格的GeoSOT编码作为待处理空间数据的剖分标识，将所述待处理空间数据与其剖分标识相关联存储。

3.如权利要求2所述的多源空间数据的剖分标识生成方法，其特征在于，所述待处理空间数据为影像，则所述步骤2中获取待处理影像的空间分辨率信息；层级m的选择方法为：用户指定理想网格所对应的影像的像素大小，根据待处理影像的分辨率计算理想网格的尺寸，找到最接近理想网格尺寸的GeoSOT网格，该GeoSOT网格所处层级即为层级m。

4.如权利要求2所述的多源空间数据的剖分标识生成方法，其特征在于，所述待处理空间数据为矢量，则在步骤2中获取待处理空间数据的比例尺信息，层级m的选择方法为：根据设定的矢量比例尺与GeoSOT网格尺寸的关系，通过比例尺信息获取GeoSOT网格的尺寸，该尺寸的GeoSOT网格所处层级即为层级m。

5.如权利要求2所述的多源空间数据的剖分标识生成方法，其特征在于，所述步骤3中，中心点所在的m层级的GeoSOT网格的GeoSOT编码的获取方法具体如下：

步骤301，以步骤1中所获得的中心点的经纬度坐标P（L，B）为依据，通过设定的规则获得中心点所在的m层级的GeoSOT网格的GeoSOT编码的第1位；

所述中心点P（L，B），前缀为E或N的经纬度值表示为正数，前缀为W或S的经纬度值表示为负数；所述设定的规则为：

当L>0且B>0时，第1位为0；

当L<0且B>0时，第1位为1；

当L>0且B<0时，第1位为2；

当L<0且B<0时，第1位为3；

步骤302，令所述中心点的经纬度坐标P（L，B）=（L_D°L_M′L_S″，B_D°B_M′B_S″），其中L_D、L_M、L_S分别为L的度、分、秒的数值，B_D、B_M、B_S分别为B的度、分、秒的数值，使用下式将L和B分别转换为GeoSOT坐标值L_GeoSOT和B_GeoSOT：

L_GeoSOT=L_D×64×64+L_M×64+L_S；

B_GeoSOT=B_D×64×64+B_M×64+B_S；

步骤302，将步骤2中所获得的层级为m的GeoSOT网格的尺寸记为CellSize，将L_GeoSOT和B_GeoSOT分别除以CellSize，得到的结果转换为二进制数，若得到的二进制数不足m-1位，则在前面补零使其满m-1位，最终获得经纬度GeoSOT坐标值的二进制数；

步骤303，对于步骤302中所获得的经纬度GeoSOT坐标值的二进制数按位计算如下：将纬度的GeoSOT坐标值的二进制数每一位的2倍加上经度的GeoSOT坐标值的二进制数对应位，获得四进制编码；

步骤304，组合步骤301所获得的P的编码第1位以及步骤303中所获得的四进制编码，即为最终P(L，B)的GeoSOT编码，将P的GeoSOT编码作为P所属空间数据的剖分标识。

6.一种数据检索方法，其特征在于，具体为：依据如权利要求1所述的多源空间数据的剖分标识生成方法，对所有的空间数据均进行剖分标识的生成；

对于所需检索的地理范围内的空间数据，包含该地理范围的网格，使用该网格的编码作为检索号；使用所述检索号与存储的剖分标识进行匹配，找到与检索号完全一致的剖分标识，提取该剖分标识关联的空间数据，即获得所需检索的空间数据。

7.如权利要求6所述的一种数据检索方法，其特征在于，使用GeoSOT网格对全球进行划分，该方法包括如下步骤：

步骤一、对于空间数据，获取其所标识的地理空间范围的中心点的经纬度坐标；

若空间数据所标识的是不规则的地理空间范围，则取该范围的外包矩形，找到外包矩形的中心点，获取该中心点的经纬度坐标；

步骤二、设定对所述空间数据进行标识所使用的GeoSOT网格的层级m；

步骤三、获取上述中心点所在的m层级的GeoSOT网格的GeoSOT编码作为空间数据的剖分标识，获取所有空间数据的剖分标识，并将各空间数据与其剖分标识相关联存储；

步骤四、对于所需检索的地理范围内的空间数据，，获取该地理范围的中心点的经纬度坐标；

若所需检索的空间数据所标识的是不规则的地理空间范围，则取该范围的外包矩形，找到该外包矩形的中心点，获取该中心点的经纬度坐标；

步骤五、设定对所需检索的空间数据进行剖分关联所使用的GeoSOT网格的层级n；

步骤六、获取所述步骤五中获得的中心点所在的n层级的GeoSOT网格的GeoSOT编码作为检索号；

步骤七、使用所述检索号与存储的剖分标识进行匹配，找到与检索号完全匹配一致的剖分标识，提取所述剖分标识关联的空间数据，即为所需检索的空间数据。

Images