geohash_encode_decode

def goehash_encode(longitude, latitude, precision):

    base32 = '0123456789bcdefghjkmnpqrstuvwxyz'

    lon_interval, lat_interval =  (-180.0, 180.0), (-90.0, 90.0)

    geohash = []

    bits = [ 16, 8, 4, 2, 1 ]

    bit = 0

    ch = 0

    even = True

    while len(geohash) < precision:

        if even:

            mid = (lon_interval[0] + lon_interval[1]) / 2

            if longitude > mid:

                ch |= bits[bit]

                lon_interval = (mid, lon_interval[1])

            else:

                lon_interval = (lon_interval[0], mid)

        else:

            mid = (lat_interval[0] + lat_interval[1]) / 2

            if latitude > mid:

                ch |= bits[bit]

                lat_interval = (mid, lat_interval[1])

            else:

                lat_interval = (lat_interval[0], mid)

        even = not even

        if bit < 4:

            bit += 1

        else:

            geohash += base32[ch]

            bit = 0

            ch = 0

    return ''.join(geohash)

def goehash_decode(geohash):

    base32 = '0123456789bcdefghjkmnpqrstuvwxyz'

    decodemap = {}

    for i in range(len(base32)):

        decodemap[base32[i]] = i 

    del i

    lon_interval, lat_interval = (-180.0, 180.0), (-90.0, 90.0)

    lat_err, lon_err = 90.0, 180.0

    is_even = True

    print decodemap

    for c in geohash:

        cd = decodemap[c]

        for mask in [16, 8, 4, 2, 1]: 

            if is_even: # adds longitude info

                lon_err /= 2

                if cd & mask:

                    lon_interval = ((lon_interval[0]+lon_interval[1])/2, lon_interval[1])

                else:

                    lon_interval = (lon_interval[0], (lon_interval[0]+lon_interval[1])/2)

            else:     # adds latitude info

                lat_err /= 2

                if cd & mask:

                    lat_interval = ((lat_interval[0]+lat_interval[1])/2, lat_interval[1])

                else:

                    lat_interval = (lat_interval[0], (lat_interval[0]+lat_interval[1])/2)

            is_even = not is_even

    lat = (lat_interval[0] + lat_interval[1]) / 2 

    lon = (lon_interval[0] + lon_interval[1]) / 2 

    print lon_err, lat_err

    return ','.join([str(lon), str(lat), str(lon_err), str(lat_err)])

print goehash_encode(109.22222601, 34.01730767, 8)

print goehash_decode("wnrz4rjxt")