def plot_renov(self, data=None): # Plot renovables (sin eólica) datos_plot = data if data is not None else self.data rdp.plot_renov(datos_plot)
errores[errores.anyo != 2009].hist() errores_anyo = errores[dt.date.today().strftime('%Y')] print 'Hay %lu errores en la base de datos, producidos en %lu días (%lu de este año).\nEl último fue en %s, con un delta_Ts de %lu' \ % (len(errores), len(set(errores.ordinal)), len(errores_anyo), errores.index[-1], errores.delta_T[-1]) dias_errores = sorted(set(errores.str_dia)) print data.nuc[data.nuc < 3000] data.nuc.plot() print data.car[data.car < 10] print data['2007-11-24 10:00':'2007-11-24 14:00'] # Plot renovables (sin eólica) rdp.plot_renov(data) # Plot hid, car, gf (buscando problemas en datos) fig, hejes = rdp.get_lienzo() # data.dem.plot(ax=hejes,kind='Area'), plt.show() data.gf.plot(ax=hejes, label='Fuel', alpha=.7) data.hid.plot(ax=hejes, label=u'Hidráulica', alpha=.7) data.car.plot(ax=hejes, label=u'Carbón', alpha=.7) rdp.plt.show() # Días raros (errores en los datos de REE) for d in DIAS_RAROS: rdp.plot_dia_raro(data, d) # print df0.ix['2007-10-26 21:00'] # fig, hejes = rdp.get_lienzo()
errores[errores.anyo != 2009].hist() errores_anyo = errores[dt.date.today().strftime('%Y')] print 'Hay %lu errores en la base de datos, producidos en %lu días (%lu de este año).\nEl último fue en %s, con un delta_Ts de %lu' \ % (len(errores), len(set(errores.ordinal)), len(errores_anyo), errores.index[-1], errores.delta_T[-1]) dias_errores = sorted(set(errores.str_dia)) print data.nuc[data.nuc < 3000] data.nuc.plot() print data.car[data.car < 10] print data['2007-11-24 10:00':'2007-11-24 14:00'] # Plot renovables (sin eólica) rdp.plot_renov(data) # Plot hid, car, gf (buscando problemas en datos) fig, hejes = rdp.get_lienzo() # data.dem.plot(ax=hejes,kind='Area'), plt.show() data.gf.plot(ax=hejes, label='Fuel', alpha=.7) data.hid.plot(ax=hejes, label=u'Hidráulica', alpha=.7) data.car.plot(ax=hejes, label=u'Carbón', alpha=.7) rdp.plt.show() # Días raros (errores en los datos de REE) for d in DIAS_RAROS: rdp.plot_dia_raro(data, d) # print df0.ix['2007-10-26 21:00']