-
Notifications
You must be signed in to change notification settings - Fork 0
/
FlowAnalysis.py
75 lines (58 loc) · 2.05 KB
/
FlowAnalysis.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
"""
Created on Thu Sep 22 16:37:29 2016
@author: utente
Flow Analysis
"""
import pandas as pd
import numpy as np
from collections import OrderedDict
import matplotlib.pyplot as plt
flow = pd.read_excel('C:/Users/d_floriello/Documents/Flows.xlsx', sheetname = 'flow')
pun = pd.read_excel('C:/Users/d_floriello/Documents/Flows.xlsx', sheetname = 'pun', header = None)
flow = flow.set_index(flow['Date'])
rng = pd.date_range('2016-09-01', '2016-09-23', freq = 'D')
diz = OrderedDict()
CF = []
CN = []
F = []
#D = []
for d in rng:
CF.append(flow[flow.columns[2]].ix[flow.index == d].mean())
CN.append(flow[flow.columns[3]].ix[flow.index == d].mean())
F.append(flow[flow.columns[4]].ix[flow.index == d].mean())
# D.append(np.mean(flow[flow.columns[2]].ix[flow.index == d] - flow[flow.columns[4]].ix[flow.index == d]))
diz['CF'] = np.array(CF)
diz['CN'] = np.array(CN)
diz['F'] = np.array(F)
#diz['diff'] = np.array(D)
diz['anom'] = np.array(list([1,1,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,1]))
diz['pun'] = pun[pun.columns[0]].ix[0:22]
df= pd.DataFrame.from_dict(diz).set_index(rng)
df[df.columns[0:3]].plot()
plt.figure()
df[df.columns[-1]].plot()
################################################################
dem = pd.read_excel('C:/Users/d_floriello/Documents/Demand.xlsx')
rng = pd.date_range('2016-09-01', '2016-09-24', freq = 'D')
diz = OrderedDict()
sard = []
sici = []
sud = []
csud = []
cnor = []
nord = []
for d in rng:
sard.append(dem['SARD'].ix[dem.index == d].mean())
sici.append(dem['SICI'].ix[dem.index == d].mean())
sud.append(dem['SUD'].ix[dem.index == d].mean())
csud.append(dem['CSUD'].ix[dem.index == d].mean())
cnor.append(dem['CNOR'].ix[dem.index == d].mean())
nord.append(dem['NORD'].ix[dem.index == d].mean())
diz['SARD'] = np.array(sard)
diz['SICI'] = np.array(sici)
diz['SUD'] = np.array(sud)
diz['CSUD'] = np.array(csud)
diz['CNOR'] = np.array(cnor)
diz['NORD'] = np.array(nord)
de = pd.DataFrame.from_dict(diz).set_index(rng)
de.plot()