Skip to main content
3 de 3
se eliminaron 6 caracteres en el cuerpo
FJSevilla
  • 62.1k
  • 7
  • 43
  • 68

@abulafia indica correctamente en su respuesta los fallos que tienes en tu código y como parsear la columna con la hora, no voy a incidir en ello. Dado que importas Pandas, aunque no lo uses, puedes usar pandas.read_csv para parsear tu archivo, lo que simplifica en mi opinión mucho la tarea y te evitas trabajar con listas. En esencia tienes un csv con espacios como separador y con algunas filas que son comentarios:

import matplotlib.pyplot as plt
import pandas as pd

df = pd.read_csv('20071018_G10xr_1m.txt', 
                 sep="\s+",
                 comment="#",
                 names=("HHMM", "Short", "Long"), 
                 date_parser=lambda date: pd.datetime.strptime(date, '%H%M').time(),
                 parse_dates=[0]
                 )


fig, ax = plt.subplots()
df.plot(x="HHMM", y="Short", color='blue', ax=ax)
ax.set_yscale('log')
ax.legend(['GOES Xray Flux'])
ax.set_xlabel('Tiempo UTC')
ax.set_ylabel(r'X-ray Flux ($Watts/m^{2}$)')
plt.show()

Para parsear la hora se usa una función lambda que usa el formato %H%M para obtener un objeto datetime y luego eliminamos la información de la fecha usando el método time().

El dataframe que genera es algo así:

       HHMM         Short          Long
0  00:00:00  3.450000e-09  3.730000e-09
1  00:01:00  3.450000e-09  3.730000e-09
2  00:02:00  3.450000e-09  3.730000e-09
3  00:03:00  3.450000e-09  3.730000e-09
4  00:04:00  3.440000e-09  3.730000e-09
5  00:05:00  3.490000e-09  3.730000e-09
6  00:06:00  3.440000e-09  3.730000e-09
7  00:07:00  3.440000e-09  3.730000e-09
8  00:08:00  3.440000e-09  3.730000e-09
9  00:09:00  3.440000e-09  3.730000e-09

Si no te interesa la columna Long, ignorala al leer el archivo usando el argumento usecols:

df = pd.read_csv('20071018_G10xr_1m.txt', 
                 sep="\s+",
                 comment="#",
                 names=("HHMM", "Short", "Long"), 
                 date_parser=lambda date: pd.datetime.strptime(date, '%H%M').time(),
                 parse_dates=[0], 
                 usecols=[0, 1]
                 )

Si no te resulta cómodo usar pandas.DataFrame.plot, puedes usar matplotlib directamente como tu haces:

plt.semilogy(df["HHMM"], df["Short"], color='blue')
FJSevilla
  • 62.1k
  • 7
  • 43
  • 68