GeoHash算法小记

昨天被问到一道题目，感觉挺有意思，这里记录一下:

题目:给定一个m*m的矩形区域，已知一点的(x,y)坐标,如何找到离周围xkm的所有店家?

这里需要用到一个叫做GeoHash的算法。

这个算法是什么作用呢?

GeoHash可以将地图网格化，然后把一个位置的坐标转换成一个可以排序的，可以比较的字符串编码,说到底就是将一个二维的数据映射成一个一维的数据。

核心代码:

public static String GeohashEncode(double latitude,double longitude,int precision=12){
        boolean  even=true;
        int bit=0;
        int ch=0;
        String geohash="";
        double[] lat={-90.0,90.0};
        double[] lon={-180.0,180.0};
        if(precision<1||precision>20)  precision=12;
        while(geohash.length()<precision){
            double mid;
            if(even){
                mid=(lon[0]+lon[1])/2;
                if(longitude>mid){
                    ch|=Bits[bit];
                    lon[0]=mid;

                }
                else{
                    lon[1]=mid;
                }
            }
            else{
                mid=(lon[0]+lon[1])/2;
                if(latitude>mid){
                    ch|=Bits[bit];
                    lat[0]=mid;
                }
                else{
                    lat[1]=mid;
                }
            }
            even=!even;
            if(bit<4) bit++;
            else{
                geohash+=Base32[ch];
                bit=0;
                ch=0;

            }
        }
        return geohash;
    }

结合上边的代码以及图示我们可以知道，算法的目的就是不断对(-90,90)和(-180,180)这两个区间进行划分,然后查看给定的坐标是不是在这个区间里边，如果在这个区间中,那么设置为1，如果不在这个区间中那么设置为0，直到划分的区间满足对应的精度结束循环，并且得到一串二进制数，然后再将这个二进制数变化成10进制的数，通过base32编码得到一串字符串。

base32编码表如下:

	0	1	2	3	4	5	6	7	8	9	10	11	12	13	14	15
base32	0	1	2	3	4	5	6	7	8	9	b	c	d	e	f	g
十进制	16	17	18	19	20	21	22	23	24	25	26	27	28	29	30	31
base32	h	j	k	m	n	p	q	r	s	t	u	v	w	x	y	z

下面回到这个题目,题目要求我们找到所有距离坐标Xkm的商店,我们可以采取以下的方法:

以给定的点向外拓xkm，分别对周围8个坐标进行GeoHash(如下图),然后提取字符串的公共部分，最后查找周围点的GeoHash值以xx开头的就是满足题意的商家。