-
Notifications
You must be signed in to change notification settings - Fork 1
/
jmalib.py
162 lines (131 loc) · 4.69 KB
/
jmalib.py
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
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
# -*- coding: utf-8 -*-
import datetime
import re
import imglib
color_reduction_table = {
(255, 0, 0): (255, 0, 0), # 雨雲 80mm/h 以上
(255, 0, 255): (255, 0, 255), # 雨雲 50-80mm/h
(255, 153, 0): (255, 153, 0), # 雨雲 30-50mm/h
(255, 255, 0): (255, 255, 0), # 雨雲 20-30mm/h
( 0, 255, 0): ( 0, 255, 0), # 雨雲 10-20mm/h
( 0, 0, 255): ( 0, 0, 255), # 雨雲 5-10mm/h
( 51, 102, 255): ( 51, 102, 255), # 雨雲 1- 5mm/h
(153, 204, 255): (153, 204, 255), # 雨雲 0- 1mm/h
( 96, 57, 19): ( 96, 128, 96), # 観測点
(230, 230, 230): (255, 255, 255), # 都道府県境界
(102, 102, 102): (102, 102, 102), # 海岸境界/グリッド
(255, 255, 255): (255, 255, 255), # 海岸境界
(116, 123, 114): (116, 123, 114), # 海岸境界
(160, 160, 160): (160, 160, 160), # 海岸境界
}
railfall_table = {
(255, 0, 0): (80, 100),
(255, 0, 255): (50, 80),
(255, 153, 0): (30, 50),
(255, 255, 0): (20, 30),
( 0, 255, 0): (10, 20),
( 0, 0, 255): ( 5, 10),
( 51, 102, 255): ( 1, 5),
(153, 204, 255): ( 0, 1),
}
def parse_radar_js(src):
regexp = re.compile(r"(\d{4})(\d{2})(\d{2})(\d{2})(\d{2})-(\d{2})\.png")
result = []
for (year, month, day, hour, minute, ordinal) in regexp.findall(src):
time = datetime.datetime(
int(year), int(month), int(day),
int(hour), int(minute))
result.append((time, int(ordinal)))
return result
def parse_nowcast_js(src):
regexp = re.compile(r"(\d{4})(\d{2})(\d{2})(\d{2})(\d{2})-(\d{2})\.png")
result = []
for (year, month, day, hour, minute, ordinal) in regexp.findall(src):
time = datetime.datetime(
int(year), int(month), int(day),
int(hour), int(minute))
result.append((time, int(ordinal)))
return result
def get_latest_time(time_and_ordinals):
times = [time for (time, ordinal) in time_and_ordinals]
times.sort()
return times[-1]
def get_current_radar_time(fetcher):
url = "http://www.jma.go.jp/jp/radnowc/hisjs/radar.js"
source = fetcher(url)
parsed = parse_radar_js(source)
return get_latest_time(parsed)
def get_current_nowcast_time(fetcher):
url = "http://www.jma.go.jp/jp/radnowc/hisjs/nowcast.js"
source = fetcher(url)
parsed = parse_nowcast_js(source)
return get_latest_time(parsed)
def create_radar_image_url(area, time):
url = "http://www.jma.go.jp/jp/radnowc/imgs/radar"
url += "/" + ("%03i" % area)
url += "/" + time.strftime("%Y%m%d%H%M")
url += "-00.png"
return url
def create_nowcast_image_url(area, time, ordinal):
url = "http://www.jma.go.jp/jp/radnowc/imgs/nowcast"
url += "/" + ("%03i" % area)
url += "/" + time.strftime("%Y%m%d%H%M")
url += "-" + ("%02i" % ordinal)
url += ".png"
return url
def create_image_url(area, time, ordinal):
if ordinal == 0:
return create_radar_image_url(area, time)
else:
return create_nowcast_image_url(area, time, ordinal)
def get_image(area, time, ordinal, fetcher):
url = create_image_url(area, time, ordinal)
binary = fetcher(url)
return binary
def is_water_color(rgb):
if rgb == (192, 192, 192): return True
elif rgb == (193, 193, 193): return True
else:
r, g, b = rgb
return (b >= g)
def is_ground_color(rgb):
r, g, b = rgb
return (g > b)
def color_reduction(rgb):
ret = color_reduction_table.get(rgb)
if ret == None:
if is_ground_color(rgb):
return ( 96, 128, 96)
elif is_water_color(rgb):
return ( 64, 96, 128)
else:
return (192, 0, 192)
else:
return ret
def get_minimum_rainfall_from_rgb(rgb):
return railfall_table.get(rgb, (None, None))[0]
def get_maximum_rainfall_from_rgb(rgb):
return railfall_table.get(rgb, (None, None))[1]
def get_rainfall_from_bitmap(bitmap):
pixels = bitmap.get_pixels()
minimums1 = [get_minimum_rainfall_from_rgb(pixel) for pixel in pixels]
maximums1 = [get_maximum_rainfall_from_rgb(pixel) for pixel in pixels]
minimums2 = [value for value in minimums1 if value != None]
maximums2 = [value for value in maximums1 if value != None]
minimum = min(minimums2) if len(minimums2) > 0 else 0
maximum = max(maximums2) if len(maximums2) > 0 else 0
return (minimum, maximum)
def get_full_time_table(radar_time, nowcast_time):
time_table = [((radar_time, 0), radar_time)]
for i in range(6):
ordinal = i + 1
present_time = nowcast_time + datetime.timedelta(minutes = ordinal * 10)
time_table.append(((nowcast_time, ordinal), present_time))
return time_table
def get_time_table(radar_time, nowcast_time):
full_table = get_full_time_table(radar_time, nowcast_time)
time_table = []
for (time, ordinal), present_time in full_table:
if ordinal > 0 and present_time <= radar_time: continue
time_table.append(((time, ordinal), present_time))
return time_table