File: jpiarea.py

package info (click to toggle)
python-geohash 0.8.3-1.1
  • links: PTS, VCS
  • area: main
  • in suites: bullseye, buster, sid, stretch
  • size: 160 kB
  • ctags: 132
  • sloc: python: 866; cpp: 715; makefile: 44; ansic: 17
file content (87 lines) | stat: -rw-r--r-- 2,176 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
# coding: UTF-8
# Coder for Japanese iarea grid code.
# NTT DoCoMo's Open iArea in Japan use a gridcode which is very similar to
# JIS X 0410, but absolutely different in detail.

def _encode_i2c(lat,lon,basebits):
	t=[]
	for i in range(basebits-3):
		t.append((lat&1)*2 + (lon&1))
		lat = lat>>1
		lon = lon>>1
	
	if basebits>=3:
		t.append(lon&7)
		t.append(lat&7)
		lat = lat>>3
		lon = lon>>3
	
	t.append(lon)
	t.append(lat)
	t.reverse()
	return ''.join([str(i) for i in t])

def encode(lat, lon):
	if lat<7 or lon<100:
		raise Exception('Unsupported location')
	
	basebits = 8
	return _encode_i2c(int(lat * (1<<basebits) * 1.5), int((lon-100.0)*(1<<basebits)), basebits)

def _decode_c2i(gridcode):
	lat = lon = 0
	base = 1
	basebits = 0
	if len(gridcode)>6:
		for i in gridcode[6:]:
			lat = (lat<<1) + int(int(i)/2)
			lon = (lon<<1) + int(i)%2
			base = base<<1
			basebits += 1
	
	if len(gridcode)>4:
		lat = int(gridcode[4:5])*base + lat
		lon = int(gridcode[5:6])*base + lon
		base = base<<3
		basebits += 3
	
	lat = int(gridcode[0:2])*base + lat
	lon = int(gridcode[2:4])*base + lon
	
	return (lat, lon, basebits)

def decode_sw(gridcode, delta=False):
	lat, lon, basebits = _decode_c2i(gridcode)
	
	if delta:
		return (float(lat)/(1.5*(1<<basebits)), float(lon)/(1<<basebits)+100.0, 1.0/(1.5*(1<<basebits)), 1.0/(1<<basebits))
	else:
		return (float(lat)/(1.5*(1<<basebits)), float(lon)/(1<<basebits)+100.0)

def decode(gridcode):
	lat, lon, basebits = _decode_c2i(gridcode)
	return ((lat<<1)+1)/float(3<<basebits), 100.0+((lon<<1)+1)/float(2<<basebits)

def bbox(gridcode):
	(a,b,c,d) = decode_sw(gridcode, True)
	return {'w':a, 's':b, 'n':b+d, 'e':a+c}

def neighbors(gridcode):
	(lat,lon,basebits)=_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<<basebits):
			continue
		if tlon<0 or tlon>(100<<basebits):
			continue
		
		ret.append(_encode_i2c(tlat,tlon,basebits))
	
	return ret

def expand(gridcode):
	ret = neighbors(gridcode)
	ret.append(gridcode)
	return ret