#!/usr/bin/python
# -*- coding: utf-8 -*-
#Falta
#Posar el nom del colectiu
#L1Ard2 Taller #L2Ard2 Taller #L3Ard2 Taller
#L1Ard3 Coworking #L2Ard3 Fugart #L3Ard3 Comú
#L1Ard4 Comú #L2Ard4 Comú #L3Ard4 Cooperatives
#L1Ard5 General #L2Ard5 General #L3Ard5 General
#L1Ard6 Cronopis #L2Ard6 Cronopis #L3Ard6 Cronopis
#L1Ard7 Estudi i Mig #L2Ard7 Taller de Cosir #L3Ard7 LAC
#Arduino 2,5 i 6 SUMEM
#Arduino 3,4,7 Separat
#Falta controlar si falta un dia les dades!!!
import numpy as np
import matplotlib.pyplot as plt
import os, sys, bz2
from pylatex import Document, Section, Figure, SubFigure, NoEscape, Command
from pylatex.base_classes.command import Options, Arguments
from pylatex.utils import verbatim, bold
from pylatex.package import Package
if __name__ == '__main__':
#Inicialitzem les dades
PL1=4
SL1=5
PL2=8
SL2=9
PL3=12
SL3=13
MESOS=12
Mesos=['Gen', 'Feb', 'Marc', 'Abril', 'Maig', 'Juny', 'Juliol', 'Agost', 'Set', 'Oct', 'Nov', 'Des']
MaxDiari=0
MaxMensual=0
consum_mensualL1=[[0]*12 for x in range(12)]
consum_mensualL2=[[0]*12 for x in range(12)]
consum_mensualL3=[[0]*12 for x in range(12)]
DadesMesos=[0]*8
#Volem saber el mes de la ultima lectura
for root, dirs, files in os.walk('Factura/Dades'):
file=sorted(files)[-1]
mesactual=file.split('-')[1].split('.')[0][:-2]
MesActual=Mesos[int(mesactual[2:])-1]
if (MesActual == 'Gen') or (MesActual == 'Marc') or (MesActual == 'Maig') or (MesActual == 'Juliol') or (MesActual == 'Agost') or (MesActual == 'Oct') or (MesActual == 'Des'):
numdies=31
elif MesActual == ('Abril') or MesActual == ('Juny') or MesActual == ('Set') or MesActual == ('Nov'):
numdies=30
elif (mesactual[0:2] == 16) or (mesactual[0:2] == 20) or (mesactual[0:2] == 24) or (mesactual[0:2] == 28) or (mesactual[0:2] == 32):
numdies=29
else:
numdies=28
DIES=numdies+1
valorDiariL1=[[0]*32 for i in range(8)]
valorDiariL2=[[0]*32 for i in range(8)]
valorDiariL3=[[0]*32 for i in range(8)]
print "Ultim mes que tenim dades: "+Mesos[int(mesactual[2:])-1]+" de l\'any 20"+ mesactual[0:2]
#Creem els directoris on guardarem els grafics si no existeixen
if not os.path.exists('Grafics'):
os.makedirs('Grafics')
#for arduino in range(2,8):
for arduino in 5,7:
print "Analitzant dades arduino ",arduino
CAMI="Factura/Dades/"+str(arduino)+"/"
if not os.path.exists('Grafics/'+str(arduino)):
os.makedirs('Grafics/'+str(arduino))
if (arduino=="2") or (arduino=="4") or (arduino=="7"):
CampL1=SL1
CampL2=SL2
CampL3=SL3
else:
CampL1=PL1
CampL2=PL2
CampL3=PL3
#Ens movem per dins del directori
for root, dirs, files in os.walk(CAMI):
Mes_aux=sorted(files)[0].split('-')[1].split('.')[0][:-2]
aux_mes=0
for file in sorted(files):
if file[0]=="D":
consum=[[],[],[]]
#Obrim el fitxers per llegir
f=bz2.BZ2File('Factura/Dades/'+str(arduino)+'/'+file,'r')
print file
#Mirem el mes i dia que estem analitzant
Mes=file.split('-')[1].split('.')[0][:-2]
dia=int(file.split('-')[1].split('.')[0][4:6])
#Llegim la primera linea i ens quedem amb el temps
line=f.readline()
camps=line.split(' ')
#comprovem si es una linea plena
while (len(camps)<3) and (line!=""):
line=f.readline()
camps=line.split(' ')
if line=="":
consumL1=0
consumL2=0
consumL3=0
else:
potencia=[float(camps[CampL1]),float(camps[CampL2]),float(camps[CampL3])]
hora=camps[1].split(':')
temps_aux=float(hora[2])+float(hora[1])*60+float(hora[0])*3600
tempstotal=temps_aux
data=camps[0]
#Ens movem linea per linea per llegir les dades
for line in f:
camps=line.split(' ')
if 'str' in line:
break
if len(camps)>2:
lecturapotencia=[float(camps[CampL1]),float(camps[CampL2]),float(camps[CampL3])]
Error=0
for x in lecturapotencia:
Error=Error+abs(x)
#Comprovem si no hi ha valors molt grans que considerem erronis
if Error < 10000:
hora=camps[1].split(':')
temps=float(hora[2])+float(hora[1])*60+float(hora[0])*3600
temps_aux=temps-temps_aux
#Si temps_aux es negatiu vol dir que les dades ja son del dia seguent
if temps_aux < 0:
temps_aux=temps_aux+86400
tempstotal=tempstotal+temps_aux
#Multipliquem el valor de la potencia pel temps transcorregut des del ultim cop.
potencia[:]=[x*temps_aux for x in potencia]
#Afegim el valor de la potencia als consums
for i in range(3):
consum[i].append(potencia[i])
temps_aux=temps
potencia=lecturapotencia
#Afegim aixo per si no estem calculant exactament les 24h dun dia
if tempstotal < 86400:
potencia[:]=[x*(86400-tempstotal) for x in potencia]
for i in range(3):
consum[i].append(potencia[i])
#valorx.append(data)
consumL1=sum(consum[0])/3600000
consumL2=sum(consum[1])/3600000
consumL3=sum(consum[2])/3600000
print consumL1, consumL2, consumL3
#Nomes fem els calculs si estem en lultim mes
if Mes == mesactual:
valorDiariL1[arduino][dia]=consumL1
valorDiariL2[arduino][dia]=consumL2
valorDiariL3[arduino][dia]=consumL3
if Mes == Mes_aux:
consum_mensualL1[arduino][aux_mes]+=consumL1
consum_mensualL2[arduino][aux_mes]+=consumL2
consum_mensualL3[arduino][aux_mes]+=consumL3
else :
aux_mes+=1
consum_mensualL1[arduino][aux_mes]=consumL1
consum_mensualL2[arduino][aux_mes]=consumL2
consum_mensualL3[arduino][aux_mes]=consumL3
Mes_aux=file.split('-')[1].split('.')[0][:-2]
DadesMesos[arduino]=aux_mes+1
print u"Creant gràfics ..."
aux_mes+=1
Colectiu= [0]*3
#for arduino in range(2,8):
for arduino in 5,7:
left = np.arange(1,DIES)
width = 0.95
if arduino == 2:
for i in range(2,4):
for index in range(1,int(DIES)):
valorDiariL1[2][index]+=eval('valorDiariL'+str(i))[2][index]
for index in range(1,aux_mes):
consum_mensualL1[2][index]+=consum_mensualL2[2][index]+consum_mensualL3[2][index]
RANG=2
Colectiu=[u"Taller Col·lectiu","",""]
elif arduino == 3:
Colectiu=["Coworking","Fugart",u"Comú 2a Planta"]
RANG=4
elif arduino == 4:
Colectiu=[u"Comú Taller",u"Comú Sotan",u"Cooperatives"]
RANG=4
elif arduino == 5:
for i in range(2,4):
for index in range(1,int(DIES)):
valorDiariL1[5][index]+=eval('valorDiariL'+str(i))[5][index]
for index in range(1,aux_mes):
consum_mensualL1[5][index]+=consum_mensualL2[5][index]+consum_mensualL3[5][index]
Colectiu=["General de Can Fugarolas","",""]
RANG=2
elif arduino == 6:
for i in range(2,4):
for index in range(1,int(DIES)):
valorDiariL1[6][index]+=eval('valorDiariL'+str(i))[6][index]
for index in range(1,aux_mes):
consum_mensualL1[6][index]+=consum_mensualL2[6][index]+consum_mensualL3[6][index]
Colectiu=["Cronopis","",""]
RANG=2
elif arduino == 7:
Colectiu=["Estudi i Mig","Taller de Cosir","LAC"]
RANG=4
for i in range(1,RANG):
valors=eval('valorDiariL'+str(i))[arduino][1:int(DIES)]
plt.clf()
plt.axis([1,DIES,0,max(valors)*1.1])
plt.bar(left, valors, width, bottom=0.01, color='r',alpha=0.8)
plt.xlabel('Dies')
plt.ylabel('Consum en KWh')
plt.title('Consum Diari del mes '+ Mesos[int(mesactual[2:])-1]+' de l\'any 20'+ mesactual[0:2]+' ('+Colectiu[i-1]+')')
plt.xticks(left+width/2,left)
plt.tight_layout()
plt.savefig('Grafics/'+str(arduino)+'/ConsumL'+str(i)+'Diari_Ard'+str(arduino)+'.png')
left = np.arange(1,len(consum_mensualL1)+1)
for i in range(1,RANG):
valors=eval('consum_mensualL'+str(i))[arduino][0:12]
plt.clf()
plt.axis([0,MESOS+1,0,max(valors)*1.1])
plt.bar(left, valors, width, bottom=0.01, color='y',alpha=0.8)
plt.ylabel('Consum en KWh')
plt.title('Consum Mensual')
plt.title('Consum Mensual de tot un any fins '+ Mesos[int(mesactual[2:])-1]+' del 20'+mesactual[0:2]+' ('+Colectiu[i-1]+')')
plt.xticks(left+width/2,np.roll(Mesos,DadesMesos[arduino]-int(mesactual[2:])))
plt.tight_layout()
plt.savefig('Grafics/'+str(arduino)+'/ConsumL'+str(i)+'Mensual_Ard'+str(arduino)+'.png')
Colectiu=[0]*3
print u"Creant factures ..."
#for arduino in range(2,8):
for arduino in 5,7:
if arduino == 2:
RANG=2
Colectiu=[u"Taller Col·lectiu","",""]
elif arduino == 3:
Colectiu=["Coworking","Fugart",u"Comú 2a Planta"]
RANG=4
elif arduino == 4:
Colectiu=[u"Comú Taller",u"Comú Sotan",u"Cooperatives"]
RANG=4
elif arduino == 5:
Colectiu=["General de Can Fugarolas","",""]
RANG=2
elif arduino == 6:
Colectiu=["Cronopis","",""]
RANG=2
else:
Colectiu=["Estudi i Mig","Taller de Cosir","LAC"]
RANG=4
for i in range(1,RANG):
file='FacturaArduino'+str(arduino)+'_'+str(i)
doc = Document(default_filepath=file,documentclass='report',fontenc='T1',inputenc='utf8',lmodern='True')
doc.packages.append(Package('eurosym'))
doc.packages.append(Package('babel',options='catalan'))
args = Arguments(NoEscape(r'\textheight'),'6cm')
doc.preamble.append(Command('addtolength',args))
args = Arguments(NoEscape(r'\textwidth'),'6.5cm')
doc.preamble.append(Command('addtolength',args))
args = Arguments(NoEscape(r'\voffset'),NoEscape(r'-2.5cm'))
doc.preamble.append(Command('addtolength',args))
args = Arguments(NoEscape(r'\hoffset'),NoEscape(r'-3cm'))
doc.preamble.append(Command('addtolength',args))
doc.append(Command('pagestyle','empty'))
doc.append(Command('begin','center'))
doc.append(Command('begin','flushright'))
doc.append(bold(u"Can Fugarolas\n"))
doc.append(bold('1 de '+Mesos[int(mesactual[2:])]+' de 20'+mesactual[0:2]))
doc.append(Command('end','flushright'))
Text=u"Factura del Col·lectiu "+Colectiu[i-1]+"\n"
doc.append(Command('textbf',Command('Large',Text)))
doc.append(Command('textbf',Command('Large','L'+str(i)+u" de l'Arduino "+str(arduino)+"\n")))
doc.append(Command('end','center'))
doc.append(u'Aquesta factura esta generada automàticament a partir de les dades que hem anat obtenint dels Arduinos que tenim instal·lats a la Fàbrica de Can Fugarolas.\n Tota la informació de com han estat generades aquestes factures les pots trobar a')
doc.append(Command('begin','center'))
doc.append(Command('texttt','http://energia.canfugarolas.org'))
doc.append(Command('end','center'))
doc.append(u'Si vols veure en temps real el consum de tots els col·lectius de la fàbrica, pots consultar-ho a')
doc.append(Command('begin','center'))
doc.append(Command('texttt','http://emon.canfugarolas.org'))
doc.append(Command('end','center'))
doc.append(u'Aquí et deixem el teu consum per dies i per mesos en KWh:')
image_filename = os.path.join('Grafics/'+str(arduino)+'/', 'ConsumL'+str(i)+'Diari_Ard'+str(arduino)+'.png')
image_filename2 = os.path.join('Grafics/5/','ConsumL1Diari_Ard5.png')
with doc.create(Figure(position='h!')) as kittens:
with doc.create(SubFigure(
position='b',
width=NoEscape(r'0.5\linewidth'))) as left_kitten:
left_kitten.add_image(image_filename,
width=NoEscape(r'\linewidth'))
left_kitten.add_caption(u'Últim mes')
with doc.create(SubFigure(
position='b',
width=NoEscape(r'0.5\linewidth'))) as right_kitten:
right_kitten.add_image(image_filename2,
width=NoEscape(r'\linewidth'))
right_kitten.add_caption(u'Últim mes')
kittens.add_caption(u'Consum del Col·lectiu')
image_filename = os.path.join('Grafics/'+str(arduino)+'/', 'ConsumL'+str(i)+'Mensual_Ard'+str(arduino)+'.png')
image_filename2 = os.path.join('Grafics/5/','ConsumL1Mensual_Ard5.png')
with doc.create(Figure(position='h!')) as kittens:
with doc.create(SubFigure(
position='b',
width=NoEscape(r'0.5\linewidth'))) as left_kitten:
left_kitten.add_image(image_filename,
width=NoEscape(r'\linewidth'))
left_kitten.add_caption(u'Últim any')
with doc.create(SubFigure(
position='b',
width=NoEscape(r'0.5\linewidth'))) as right_kitten:
right_kitten.add_image(image_filename2,
width=NoEscape(r'\linewidth'))
right_kitten.add_caption(u'Últim any')
kittens.add_caption(u'Consum de Tot Can Fugarolas')
doc.append(Command('begin','flushright'))
doc.append(Command('textsc',u"Comissió d'Energia de Can Fugarolas"))
doc.append(Command('end','flushright'))
doc.generate_pdf()