/usr/share/pyshared/jpiarea.py is in python-geohash 0.8.3-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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
|