powerforever/py/MouthlyAnalysis.py at main · StevenLi-phoenix/powerforever · GitHub

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
#!/usr/bin/env python
# coding: utf-8

# In[2]:


import numpy as np
import matplotlib.pyplot as plt
import pandas as pd
import scipy.interpolate
import random

plot = "plot/"
npy = "npy/"

solar_rad = pd.read_csv("xunlian.csv")
solar_rad


# In[36]:


# chage value to standard
data_min = solar_rad.values[:, 2]
saved = None
saved_index = 0
hour = []
toArry = solar_rad.values
for i in range(len(solar_rad)):
    time_alpha = toArry[i,1].split(":")[0]
    if time_alpha != saved:
        saved = time_alpha
        hour.append(sum(toArry[saved_index:i,2]))
        saved_index = i
    if i % 1000 == 0:print(i)
hour.append(sum(toArry[saved_index:i,2]))


# In[37]:


hour = np.array(hour)/max(hour)*20


# In[38]:


plt.plot(hour)


# In[23]:


24*30*12


# In[24]:


len(hour)


# In[40]:


# chage value to standard
data_min = solar_rad.values[:, 2]
saved = None
saved_index = 0
week = []
toArry = solar_rad.values
for i in range(len(solar_rad)):
    time_alpha = toArry[i,0].split("/")[1]
    if time_alpha != saved:
        saved = time_alpha
        week.append(sum(toArry[saved_index:i,2]))
        saved_index = i
    if i % 1000 == 0:print(i)
week.append(sum(toArry[saved_index:i,2]))


# In[42]:


plt.plot(week)
plt.title("Mouthly photovoltaic power generation")
plt.grid()


# In[43]:


import csv

# open the file in the write mode
f = open('Mouth_data.csv', 'w')
# create the csv writer
writer = csv.writer(f)
# write a row to the csv file
writer.writerow(week)
# close the file
f.close()


# In[44]:


len(week)


# In[45]:


week[0]


# In[ ]: