SQL Server 2008空间数据应用系列二：空间索引（Spatial Index）基础

　　在前一篇博文中我们学习到了一些关于地理信息的基础知识，也学习了空间参照系统，既地球椭球体、基准、本初子午线、计量单位、投影等相关理论知识，我们可以使用这些空间参照系统组件来定义一系列应用于地球空间上的几何图像来表示地理空间中的特定功能，表示着地球上一个一个特定的位置点。本篇主要介绍地理空间索引的概念以及微软SQL Server 2008 R2中的空间索引的应用。

一、空间索引

　　空间索引是指依据空间对象的位置和形状或空间对象之间的某种空间关系按一定的顺序排列的一种数据结构，其中包含空间对象的概要信息。作为一种辅助性的空间数据结构，空间索引介于空间操作算法和空间对象之间，它通过筛选作用，大量与特定空间操作无关的空间对象被排除，从而提高空间操作的速度和效率。

　　来至百科的定义：为便于空间目标的定位及各种空间数据操作，按要素或目标的位置和形状或空间对象之间的某种空间关系来组织和存储数据的结构。

　　关于更多空间索引的知识点需要深入了解EPSG和SRID，本文主要介绍基于微软SQL Server 2008的空间数据特性实现介绍，故不详细介绍这些基础理论知识点，有兴趣的可以访问文末的推荐资料链接，了解更多关于EPSG核SRID的相关知识。

　　SQL Server 2008中的空间数据类型作为CLR系统类型来执行。SQL Server 2008增加了数据库中的CLR类型的最大规模，提高了原来在SQL Server 2005中的8000字节的限制，这使得它可以存储非常复杂的空间数据元素，例如通过许多点定义的多边形。通过在关系表中存储空间数据，SQL Server 2008 使得可以结合空间数据到其他任何商业数据类型中去;这消除了对维护一个单独的只用于空间数据存储的维护要求，并使得可以做高性能查询，它不需要结合从多个外部来源获得的数据。

　　在SQL Server 2008中对空间索引的支持进一步增强了对空间数据的查询操作。你可以用一个集成在SQL Server数据库引擎中的适合的多级网格索引来检索空间数据。空间索引包含一个基于网格的层级，在其中每一级索引又细分为由上一级所定义的网格区域。下图为官方给出的一个空间索引的概率模型图：

　　SQL Server查询优化器会作出基于成本的决策，决定对给定的查询使用哪种索引，并且因为空间索引是数据库引擎的一个完整部分，可以作出关于是否使用特殊的空间索引的基于成本的决策，就像其它索引一样。

二、SQL Server 2008的空间参考系统

　　微软SQL Server 2008的空间系统表(实际上是一个系统视图)"sys.spatial_reference_systems"中存储了详细的支持所有大地测量空间参考系统标准，该表中的每一行数据对应于一个唯一的空间参考系统，可以基于这些空间参考系统在SQL Server 2008中来定义任何形式的空间数据。

select * from sys.spatial_reference_systems

　　系统表“sys.spatial_reference_systems”中实际上仅仅只是存储了应用于SQL Server 2008中地理参照坐标的参考系统，除了空间参考系统可以使用这张表中列出的参考系统标准，也可以基于这些标准来定义基于任何地理投影的空间数据。

三、空间参照文本格式（WKT）

　　Well-Known Text (WKT) Representation of Spatial Reference Systems：顾名思义用众所周知的文字形式来描述看见坐标参考系统，它与EPSG（2005年已经换名）的坐标文字表述模型一致。WKT字符串应用于确定空间参考系统的组成部分的描述，其定义主要包含以下10项内容：

　　1、一个总体的坐标系名

　　2、一个地理图形坐标系统名
　　3、一个基准面定义
　　4、一个椭球体的名字。长半轴(semi-major axis)和反扁率(inverse flattening)
　　5、本初子午线(prime meridian)名和其与格林威治子午线的偏移值
　　6、投影方法类型（如横轴莫卡托）
　　7、投影参数列表（如中央经线等）
　　8、一个单位的名称和其到米和弧度单位的转换参数
　　9、轴线的名称和顺序
　　10、在预定义的权威坐标系中的编码（如EPSG）

select well_known_text
    from sys.spatial_reference_systems
        where authorized_spatial_reference_id=4326
            and authority_name='EPSG'

　　基于sys.spatial_reference_systems表定义了多种不同的空间参照系统的WKT描述，如上SQL语句查询出的是我们平时在地图应用开发中使用最多的EPSG:4326空间参照系统标准，则对于的WKT如下：

GEOGCS["WGS 84",
        DATUM["World Geodetic System 1984",
        ELLIPSOID["WGS 84", 6378137, 298.257223563]],
        PRIMEM["Greenwich", 0],
        UNIT["Degree", 0.0174532925199433]
      ]

　　以上KWT字符串确定了EPSG:4326空间参照系统的组成部分元素描述，其主要使用了四种空间配置，分别为：基准（Datum）、椭圆（Ellipsoid）、本初子午线（Primem）和Unit（计量单位）。关于这些空间参考系统参数在上一篇文中中介绍过，这里就不在做重复介绍，以上KWT字符串表述的意思大致是这样的：“WGS 84坐标系，采用1984年订立的世界地理坐标系统标准，地理标准半径为：6378137........”。

四、SQL Server 2008的空间索引

　　SQL Server 2008 引入了对空间数据和空间索引的支持。“空间索引”是一种扩展索引，允许您对空间列编制索引。空间列是包含空间数据类型（如 geometry 或 geography）数据的表列。本节中的主题介绍了空间索引。

　　在 SQL Server 2008 中，空间索引(存储在:sys.spatial_indexes表中)使用 B 树构建而成，也就是说，这些索引必须按 B 树的线性顺序表示二维空间数据。因此，将数据读入空间索引之前，SQL Server 2008 先实现对空间的分层均匀分解。索引创建过程会将空间分解成一个四级“网格层次结构”。这些级别指的是“第 1 级”（顶级）、“第 2 级”、“第 3 级”和“第 4 级”。

　　每个后续级别都会进一步分解其上一级，因此上一级别的每个单元都包含下一级别的整个网格。在给定级别上，所有网格沿两个轴都有相同数目的单元（例如 4x4 或 8x8），并且单元的大小都相同。下图显示了网格层次结构每个级别的右上角单元被分解成 4x4 网格的情况。事实上，所有单元都是以这种方式分解的。因此，以此为例，将一个空间分解成四个级别的 4x4 网格际上会总共产生 65,536 个第四级单元。针对空间索引进行的空间分解与应用程序数据使用的度量单位无关。

　　网格层次结构的单元是利用多种 Hilbert 空间填充曲线以线性方式编号的。然而，出于演示目的，这里使用的是简单的按行编号，而不是由 Hilbert 曲线实际产生的编号。在下图中，几个表示建筑物的多边形和表示街道的线已经放进了一个 4x4 的 1 级网格中。第 1 级单元的编号为 1 到 16，编号从左上角的单元开始。

　　沿网格轴的单元数目确定了网格的“密度”：单元数目越大，网格的密度越大。例如，8x8 网格（产生 64 个单元）的密度就大于 4x4 网格（产生 16 个单元）的密度。网格密度是以每个级别为基础定义的。网格配置单元数目低：4X4 =16，中8X8 = 64，高16X16 =256，默认设置所有级别都为中。

　　您可以通过指定非默认的网格密度控制分解过程。例如，在不同级别指定不同网格密度对于基于索引空间的大小和空间列中的对象来优化索引可能非常有用。空间索引的网格密度显示在 sys.spatial_index_tessellations 目录视图的 level_1_grid、level_2_grid、level_3_grid 和 level_4_grid 列中。

　　将索引空间分解成网格层次结构后，空间索引将逐行读取空间列中的数据。读取空间对象（或实例）的数据后，空间索引将为该对象执行“分割过程”。分割过程通过将对象与其接触的网格单元集（“接触单元”）相关联使该对象适合网格层次结构。从网格层次结构的第 1 级开始，分割过程以“广度优先”方式对整个级别进行处理。在可能的情况下，此过程可以连续处理所有四个级别，一次处理一个级别。

　　注：以上内容截取SQL Server 2008帮助文档，详细请查看联机帮助文档。

　　[2]、SRID：http://en.wikipedia.org/wiki/SRID

　　[3]、WKT：http://en.wikipedia.org/wiki/Well-known_text

　　[4]、空间索引概论：http://www.cnblogs.com/terryfeng/archive/2009/05/27/1490803.html