文章版权由作者李晓晖和博客园共有，若转载请于明显处标明出处：http://www.cnblogs.com/naaoveGIS/

1.背景

针对判断一个点落在面图层中哪个要素上的需求，在我之前的博客：WebGIS中一种根据网格索引判断点面关系的方法（http://www.cnblogs.com/naaoveGIS/p/5148185.html）中有详细的描述。其原理大致为：

其处理步骤为：

2.当前网格索引的几个缺点

a.网格索引方案分为了一个索引文件和一个数据文件，任何请求进入时均会先读取索引文件，再读取数据文件，那么很容易出现资源争抢情况，不利于并发支持。

b.网格的大小会严重影响到查询效率，但是如果网格建立的足够小，那么索引文件不断增大，同样会导致磁盘寻址花费的时间增多。

c.数据的读取一定要经过两次IO，一次读索引，一次读数据，会影响读取效率。

3.方案的优化，基于网格索引的索引四叉树划分

四叉树、R树等均是空间索引常用的算法，这里我选择使用四叉索引来进行进一步优化。四叉树索引原理非常简单，即将一个范围根据深度，不断平分，如图所示：

这里优化思路是：将要素首先进行四叉树平分，然后对每个叶子节点包含的范围再进行网格索引生成：

4.优化方案的详细描述

4.1索引的生成步骤

a.首先生成数据文件。

b.通过设置的四叉树深度，算出叶子节点的个数。然后通过获取到的要素四角坐标，算出叶子节点的四角范围：leafminx、leafminy、leafmaxx、leafmaxy。

c.根据要素个数和网格因子，算出整个范围内网格的个数，用整个范围的四角坐标与网格因子计算，得出一个网格的BlockXsize和BlokcYsize。

d.针对每个叶子节点，建立该节点的网格索引，索引中包含了网格与要素的对应关系。

生成文件截图：

4.2数据读取

a.读取配置获取到要素的四角范围mapminx、mapminy、mapmaxx、mapmaxy、leafgeoxsize、leafgeoysize。

b.通过mapminx、mapminy、leafgeoxsize、leafgeoysize参数算出该XY坐标所在的网格索引编号。

c.读取该网格索引，获取到该索引的leafminx、leafminy、leafmaxx、leafmaxy、blockxsize、blockysize。

d.通过leafmaxx、leafmaxy以及blockxsize、blockysize算出XY所在网格索引的字节位置pos，将磁盘指针移动至该pos处。

e.获取到索引中包含的要素信息，比如要素所在的数据文件中的datapos。

f.读取数据文件，在该文件的datapos处将详细信息读取返回。

5.方案优点总结

a.将一个大索引文件分成多个索引文件，在大量随机点并发访问时，可以将压力负载至各文件上，减少同一文件读取时的资源争抢IO瓶颈。

b.每一个索引文件大小大大减小，读取会更快，磁盘寻址也会更快。

c.为增加网格命中单个（非多个）要素的概率，可以将每个网格的大小进一步缩小，其导致的网格索引增大会平摊至每个网格索引上，从而使副作用变小。

6.进一步优化

a.在读取索引基本信息后可以将该信息缓存至内存中，减少Config文件的IO次数。

b.生成索引时，如果一个网格只包含了一个要素的信息，可以将该信息也整合至网格索引中。这样，查询时，如果查询到的网格只包含单个要素，则可以直接在索引中将要素信息获取，而不需要再对数据索引做读取操作，减少对数据索引的IO次数。

-----欢迎转载，但保留版权，请于明显处标明出处：http://www.cnblogs.com/naaoveGIS/

如果您觉得本文确实帮助了您，可以微信扫一扫，进行小额的打赏和鼓励，谢谢 ^_^