-
Notifications
You must be signed in to change notification settings - Fork 22
/
Copy pathlos_loop.py
114 lines (98 loc) · 3.7 KB
/
los_loop.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
# link: https://github.com/lehaifeng/T-GCN/tree/master/data
import numpy as np
import pandas as pd
import json
import sys
import util
import h5py
outputdir = 'output/LOS_LOOP'
util.ensure_dir(outputdir)
dataurl = 'input/LOS-LOOP/'
dataname = outputdir+'/LOS_LOOP'
los_dataset = h5py.File(dataurl+'Los_traffic.h5', 'r')
los_speed = los_dataset.get('df')['axis0']
geo_id_list = []
for i in range(len(los_speed)):
geo_id_list.append(int(los_speed[i]))
# -------------------- Create .geo table--------------------------------------
geo = []
for i in range(len(geo_id_list)):
# id type
geo.append([geo_id_list[i], 'Point', '[]'])
geo = pd.DataFrame(geo, columns=['geo_id', 'type', 'coordinates'])
geo.to_csv(dataname+'.geo', index=False)
# -------------------- Create .rel table--------------------------------------
los_adj = pd.read_csv(dataurl+"los_adj.csv", header=None)
adj = np.mat(los_adj)
rel = []
rel_id_counter = 0
for i in range(adj.shape[0]):
for j in range(adj.shape[1]):
rel.append([rel_id_counter, 'geo', geo_id_list[i], geo_id_list[j], adj[i, j]])
rel_id_counter += 1
rel = pd.DataFrame(rel, columns=['rel_id', 'type', 'origin_id', 'destination_id', 'link_weight'])
rel.to_csv(dataname + '.rel', index=False)
# -----------------------------------------------------------------------------------------
# --------------------------Create .dyna table---------------------------------------------
los_speed = los_dataset.get('df')['block0_values']
speed = np.mat(los_speed)
dyna = []
dyna_id_counter = 0
def num2time(num):
if num < 288*31:
month = 3
day = num // 288 + 1
elif num < 288*(31+30):
month = 4
day = (num - 288*31) // 288 + 1
elif num < 288*(31+30+31):
month = 5
day = (num - 288*(31+30)) // 288 + 1
else:
month = 6
day = (num - 288*(31+30+31)) // 288 + 1
hour = (num % 288) // 12
minute = num % 12
if day > 31:
print(num)
sys.exit()
month = '0' + str(month)
day = str(day) if day > 9 else '0' + str(day)
hour = str(hour) if hour > 9 else '0' + str(hour)
minute = str(5*minute) if 5*minute > 9 else '0' + str(5*minute)
if minute == '01':
print(num)
sys.exit()
time = '2012-' + month + '-' + day + 'T' + hour + ':' + minute + ':' + '00Z'
return time
for j in range(speed.shape[1]):
for i in range(speed.shape[0]):
time = num2time(i)
dyna.append([dyna_id_counter, 'state', time, geo_id_list[j], speed[i, j]])
dyna_id_counter += 1
dyna = pd.DataFrame(dyna, columns=['dyna_id', 'type', 'time', 'entity_id', 'traffic_speed'])
dyna.to_csv(dataname + '.dyna', index=False)
# --------------------------Create .json table---------------------------------------------
config = dict()
config['geo'] = dict()
config['geo']['including_types'] = ['Point']
config['geo']['Point'] = {}
config['rel'] = dict()
config['rel']['including_types'] = ['geo']
config['rel']['geo'] = {'link_weight': 'num'}
config['dyna'] = dict()
config['dyna']['including_types'] = ['state']
config['dyna']['state'] = {'entity_id': 'geo_id', 'traffic_speed': 'num'}
config['info'] = dict()
config['info']['data_col'] = ['traffic_speed']
config['info']['weight_col'] = 'link_weight'
config['info']['data_files'] = ['LOS_LOOP']
config['info']['geo_file'] = 'LOS_LOOP'
config['info']['rel_file'] = 'LOS_LOOP'
config['info']['output_dim'] = 1
config['info']['time_intervals'] = 300
config['info']['init_weight_inf_or_zero'] = 'inf'
config['info']['set_weight_link_or_dist'] = 'dist'
config['info']['calculate_weight_adj'] = False
config['info']['weight_adj_epsilon'] = 0.1
json.dump(config, open(outputdir+'/config.json', 'w', encoding='utf-8'), ensure_ascii=False)