Imported Upstream version 0.2

Author: ieure

PKG-INFO

@@ -0,0 +1,10 @@
+Metadata-Version: 1.0
+Name: python-geohash
+Version: 0.2
+Summary: UNKNOWN
+Home-page: UNKNOWN
+Author: UNKNOWN
+Author-email: UNKNOWN
+License: UNKNOWN
+Description: UNKNOWN
+Platform: UNKNOWN

README

@@ -0,0 +1,6 @@
+* python-geohash
+python-geohash is a fast, accurate python geohashing library.
+
+** LICENSE
+Code is licensed under Apache License 2.0, MIT Licence and NEW BSD License.
+You can choose one of the licences above.

geohash.py

@@ -0,0 +1,153 @@
+# coding: UTF-8
+"""
+Copyright (C) 2009 Hiroaki Kawai <[email protected]>
+"""
+__all__ = ['encode','decode','decode_exactly','bbox', 'neighbors', 'expand']
+
+_base32 = '0123456789bcdefghjkmnpqrstuvwxyz'
+_base32_map = {}
+for i in range(len(_base32)):
+	_base32_map[_base32[i]] = i
+del i
+
+# def _binstr(b, length=32):
+# 	t = []
+# 	while length>0:
+# 		if b&1:
+# 			t.append('1')
+# 		else:
+# 			t.append('0')
+# 		b=b>>1
+# 		length-=1
+# 	
+# 	t.reverse()
+# 	return ''.join(t)
+
+def _encode_i2c(lat,lon,lat_length,lon_length):
+	precision=(lat_length+lon_length)/5
+	if lat_length < lon_length:
+		a = lon
+		b = lat
+	else:
+		a = lat
+		b = lon
+	
+	ret = ''
+	while precision>0:
+		c = ((a&4)<<2) + ((b&2)<<2) + ((a&2)<<1) + ((b&1)<<1) + (a&1)
+		ret += _base32[c]
+		t = a>>3
+		a = b>>2
+		b = t
+		precision-=1
+	
+	return ret[::-1]
+
+def encode(latitude, longitude, precision=12):
+	if latitude > 90.0 or latitude < -90.0:
+		raise Exception("invalid latitude.")
+	while longitude < -180.0:
+		longitude += 360.0
+	while longitude >= 180.0:
+		longitude -= 360.0
+	
+	lat = (latitude+90.0)/180.0
+	lon = (longitude+180.0)/360.0
+	
+	lat_length=lon_length=precision*5/2
+	if precision%2==1:
+		lon_length+=1
+	
+	lat = int((1<<lat_length)*lat)
+	lon = int((1<<lon_length)*lon)
+	
+	return _encode_i2c(lat,lon,lat_length,lon_length)
+
+def _decode_c2i(hashcode):
+	lon = 0
+	lat = 0
+	bit_length = 0
+	lat_length = 0
+	lon_length = 0
+	for i in hashcode:
+		t = _base32_map[i]
+		if bit_length%2==0:
+			lon = lon<<3
+			lat = lat<<2
+			lon += (t>>2)&4
+			lat += (t>>2)&2
+			lon += (t>>1)&2
+			lat += (t>>1)&1
+			lon += t&1
+			lon_length+=3
+			lat_length+=2
+		else:
+			lon = lon<<2
+			lat = lat<<3
+			lat += (t>>2)&4
+			lon += (t>>2)&2
+			lat += (t>>1)&2
+			lon += (t>>1)&1
+			lat += t&1
+			lon_length+=2
+			lat_length+=3
+		
+		bit_length+=5
+	
+	return (lat,lon,lat_length,lon_length)
+
+def decode(hashcode, delta=False):
+	(lat,lon,lat_length,lon_length) = _decode_c2i(hashcode)
+	
+	lat = (lat<<1) + 1
+	lon = (lon<<1) + 1
+	lat_length += 1
+	lon_length += 1
+	
+	latitude  = 180.0*lat/(1<<lat_length) - 90.0
+	longitude = 360.0*lon/(1<<lon_length) - 180.0
+	if delta:
+		latitude_delta  = 180.0/(1<<lat_length)
+		longitude_delta = 360.0/(1<<lon_length)
+		return latitude,longitude,latitude_delta,longitude_delta
+	
+	return latitude,longitude
+
+def decode_exactly(hashcode):
+	return decode(hashcode, True)
+
+## hashcode operations below
+
+def bbox(hashcode):
+	(lat,lon,lat_length,lon_length) = _decode_c2i(hashcode)
+	
+	ret={}
+	ret['n'] = 180.0*(lat+1)/(1<<lat_length) - 90.0
+	ret['s'] = 180.0*lat/(1<<lat_length) - 90.0
+	ret['e'] = 360.0*(lon+1)/(1<<lon_length) - 180.0
+	ret['w'] = 360.0*lon/(1<<lon_length) - 180.0
+	return ret
+
+def neighbors(hashcode):
+	(lat,lon,lat_length,lon_length) = _decode_c2i(hashcode)
+	ret = []
+	tlat = lat
+	for tlon in (lon-1, lon+1):
+		ret.append(_encode_i2c(tlat,tlon,lat_length,lon_length))
+	
+	tlat = lat+1
+	if not tlat >> lat_length:
+		for tlon in (lon-1, lon, lon+1):
+			ret.append(_encode_i2c(tlat,tlon,lat_length,lon_length))
+	
+	tlat = lat-1
+	if tlat >= 0:
+		for tlon in (lon-1, lon, lon+1):
+			ret.append(_encode_i2c(tlat,tlon,lat_length,lon_length))
+	
+	return ret
+
+def expand(hashcode):
+	ret = neighbors(hashcode)
+	ret.append(hashcode)
+	return ret

jpgrid.py

@@ -0,0 +1,160 @@
+# coding: UTF-8
+# Coder for Japanese grid square code. (JIS C 6304 / JIS X 0410)
+# 行政管理庁告示第143号 http://www.stat.go.jp/data/mesh/
+
+def _encode_i2c(lat, lon, base1):
+	t=[]
+	while base1>80:
+		t.append(1 + (lat&1)*2 + (lon&1))
+		lat = lat>>1
+		lon = lon>>1
+		base1 = base1>>1
+	
+	if base1==80:
+		t.append(lon%10)
+		t.append(lat%10)
+		lat = lat/10
+		lon = lon/10
+		base1 = base1/10
+	elif base1==16: # Uni5
+		t.append(1 + (lat&1)*2 + (lon&1))
+		lat = lat>>1
+		lon = lon>>1
+		base1 = base1>>1
+	elif base1==40: # Uni2
+		t.append(5)
+		t.append(lon%5*2)
+		t.append(lat%5*2)
+		lat = lat/5
+		lon = lon/5
+		base1 = base1/5
+	
+	if base1==8:
+		t.append(lon%8)
+		t.append(lat%8)
+		lat = lat>>3
+		lon = lon>>3
+		base1 = base1>>3
+	
+	t.append(lon)
+	t.append(lat)
+	t.reverse()
+	return ''.join([str(i) for i in t])
+
+def encode(latitude, longitude, base1=80):
+	return _encode_i2c(int(latitude*base1*1.5), int((longitude-100)*base1), base1)
+
+#def _encode_i2c(lat, lon, base1):
+def _decode_c2i(gridcode):
+	base1 = 1
+	lat = lon = 0
+	codelen = len(gridcode)
+	if codelen>0:
+		lat = int(gridcode[0:2])
+		lon = int(gridcode[2:4])
+	
+	if codelen>4:
+		lat = (lat<<3) + int(gridcode[4:5])
+		lon = (lon<<3) + int(gridcode[5:6])
+		base1 = base1<<3
+	
+	if codelen>6:
+		if codelen==7:
+			i = int(gridcode[6:7])-1
+			lat = (lat<<1) + i/2
+			lon = (lon<<1) + i%2
+			base1 = base1<<1
+		else:
+			lat = lat*10 + int(gridcode[6:7])
+			lon = lon*10 + int(gridcode[7:8])
+			base1 = base1*10
+	
+	if codelen>8:
+		if gridcode[8:]=='5':
+			lat = lat>>1
+			lon = lon>>1
+			base1 = base1>>1
+		else:
+			for i in gridcode[8:]:
+				i = int(i)-1
+				lat = (lat<<1) + i/2
+				lon = (lon<<1) + i%2
+				base1 = base1<<1
+	
+	return (lat, lon, base1)
+
+def decode_sw(gridcode, delta=False):
+	(lat, lon, base1) = _decode_c2i(gridcode)
+	
+	lat = lat/(base1*1.5)
+	lon = lon/float(base1) + 100.0
+	
+	if delta:
+		return (lat, lon, 1.0/(base1*1.5), 1.0/base1)
+	else:
+		return (lat, lon)
+
+def decode(gridcode):
+	(lat, lon, base1) = _decode_c2i(gridcode)
+	
+	lat = (lat<<1) + 1
+	lon = (lon<<1) + 1
+	base1 = base1<<1
+	return (lat/(base1*1.5), lon/float(base1) + 100.0)
+
+def bbox(gridcode):
+	(a,b,c,d) = decode_sw(gridcode, True)
+	return {'w':a, 's':b, 'n':b+d, 'e':a+c}
+
+
+## short-cut methods
+def encodeLv1(lat, lon):
+	return encode(lat,lon,1)
+
+def encodeLv2(lat, lon):
+	return encode(lat,lon,8)
+
+def encodeLv3(lat, lon):
+	return encode(lat,lon,80)
+
+def encodeBase(lat,lon):
+	return encodeLv3(lat,lon)
+
+def encodeHalf(lat,lon):
+	return encode(lat,lon,160)
+
+def encodeQuarter(lat,lon):
+	return encode(lat,lon,320)
+
+def encodeEighth(lat,lon):
+	return encode(lat,lon,640)
+
+def encodeUni10(lat,lon):
+	return encodeLv2(lat,lon)
+
+def encodeUni5(lat, lon):
+	return encode(lat,lon,16)
+
+def encodeUni2(lat, lon):
+	return encode(lat,lon,40)
+
+
+def neighbors(gridcode):
+	(lat,lon,base1)=_decode_c2i(gridcode)
+	ret = []
+	for i in ((0,-1),(0,1),(1,-1),(1,0),(1,1),(-1,-1),(-1,0),(-1,1)):
+		tlat=lat+i[0]
+		tlon=lon+i[1]
+		if tlat<0 or tlat>(90*base1):
+			continue
+		if tlon<0 or tlon>(100*base1):
+			continue
+		
+		ret.append(_encode_i2c(tlat,tlon,base1))
+	
+	return ret
+
+def expand(gridcode):
+	ret = neighbors(gridcode)
+	ret.append(gridcode)
+	return ret