经纬度搜索(1)-Geohash算法原理

geohash作为Solr的位置信息搜索算法，有必要了解其基本的实现原理。geohash算法的wiki链接为[url]http://en.wikipedia.org/wiki/Geohash[/url]，这里再结合自己的理解，重新复述一下。

[b]由经纬度变成base32标识[/b]
geohash的思想，是将所有的经纬度坐标，通过geohash，变成一个唯一的base32标识。他将世界上的所有区域进行分块，每个维度都是32块，进而将范围逐渐变小、变小，最后的一堆数字，就成了这个base32的唯一标识。比如说“DRT2Y”这个标识，首先确定D：
[img]http://dl.iteye.com/upload/attachment/0062/1008/26a4c3f4-6430-313a-a558-1520ef6fc6b3.jpg[/img]
OK,找到D区域后，再将D进行分块，精确到DR：
[img]http://dl.iteye.com/upload/attachment/0062/1010/7d082ee0-0f7c-3cc8-9b72-2be5e2cffc8f.jpg[/img]
类似的，继续找到DRT：
[img]http://dl.iteye.com/upload/attachment/0062/1006/0b1f7999-df0c-34ab-b1a0-d9eb9ee39c16.jpg[/img]
以此类推，慢慢的精确到你指定的点，这就变成了唯一的base32标识。

base 32的对应表如下：[img]http://dl.iteye.com/upload/attachment/0062/1028/5db16d31-acfc-30e3-8e88-609b84beef85.jpg[/img]

[b]encode：由经纬度到base32[/b]
使用wiki上的例子，纬度为42.6，经度为-5.6的点，转化为base32的话要如何转呢？
拿纬度来进行说明，纬度的范围为-90到90，将这个范围划为两段，则为[-90,0]、[0,90]，然后看给定的纬度在哪个范围，在前面的范围的话，就设当前位为0，后面的话值便为1.然后继续将确定的范围1分为2，继续以确定值在前段还是后段来确定bit的值。就这样慢慢的缩小范围，一般最多缩小13次就可以了(经纬度的二进制位相加最多25位，经度13位，纬度12位)。这时的中间值，将跟给定的值最相近。
结合图，看的更明白些：
[img]http://dl.iteye.com/upload/attachment/0062/1060/c820ae7f-e656-3739-a575-c7494b2953ef.jpg[/img]
其中val可以认为是下一次的中间值，而err是误差，通过mid+err=val，err当前值是上一次误差的一半。

这样，取出图中的bit位：1011 1100 1001，同样的方法，将经度(范围-180到180)算出来为
0111 1100 0000 0。

得到了经纬度的二进制位后，下面需要将两者进行结合：从经度、纬度的循环，每次取其二进制的一位(不足位取0),合并为新的二进制数：01101 11111 11000 00100 00010。每5位为一个十进制数，结合base32对应表映射为base32值为：ezs42。这样就完成了encode的过程 :)

[b]decode：由base32到经纬度[/b]
相信知道了encode的过程后，decode的过程很容易推导了。不清楚的可以参考wiki。这里就不细说了。

[b]after[/b]
接下来，我会写一下Lucene/Solr是如何应用geohash的，包括创建索引、查询索引，以及距离的计算等。