@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 fecha se usa una función lambda que usa el formato %H%M
y luego eliminamos la información de la fecha usando el método time
. Puedes hacer que esta columna sea el índice.
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]
)
S no te resulta cómodo usar pandas.DataFrame.plot, puedes usar matplotlib directamente como tu haces:
plt.semilogy(df["HHMM"], df["Short"], color='blue')