Estoy queriendo graficar la amplitud versus el tiempo de todo un día. El archivo donde se encuentra la data es de formato fits, donde el tiempo se muestra de la siguiente manera:
[ 0 1 2 ... 86397 86398 86399]
Este tiempo esta en segundos y, por otro lado, la amplitud esta en dB y se representa como:
[ -28.3 -27.84 -29.88 ... 0.18 -0.92 -0.28]
Mi código en python para poder leer la data y graficarla para todo un día es:
# -*- coding: utf-8 -*-
#!/usr/bin/env python
from astropy.io import fits
import matplotlib.pyplot as plt
import pylab
import os, sys
import numpy as np
import datetime
from datetime import timedelta
import time
import matplotlib
from datetime import datetime
#Dimensiones de la figura a graficar
fig_size = [18, 20]
plt.rcParams['figure.figsize'] = fig_size
data_vlf = fits.open('PLO-20171018.fits')
t = pylab.transpose(data_vlf[0].data)[0]
time_vlf = t.astype(int)
time_format = [str(timedelta(seconds=s)) for s in time_vlf]
t_format = '%H:%M:%S'
date_vlf = [datetime.strptime(i, t_format).time() for i in time_format]
ch = pylab.transpose(data_vlf[0].data)[6]
pylab.plot(date_vlf, ch, color='black')
plt.legend(['Channel 6'])
plt.xlabel('Time UT')
plt.ylabel('Amplitude')
plt.show()
Pero también deseo ver como evoluciona la amplitud cada hora del día, para ello, escribí el siguiente código:
h_i = 0
h_i_2 = 3600
h_f = 21600
for i in range(1, 7):
for num1 in range(h_i, h_f, 3600):
for num2 in range(h_i_2, h_f+1, 3600):
t_1 = date_vlf[num1:num2]
ch_1 = ch[num1:num2]
plt.subplot(3, 2, i)
pylab.plot(t_1, ch_1)
plt.xlabel('1 hour (UT)')
plt.ylabel('Amplitude')
plt.show()
Sin embargo, solo grafica una hora específica y no de cada hora en el día, en que se debe modificar mi código. Sigo aprendiendo el lenguaje de python, si mi código se puede mejorar agradecería cualquier ayuda.
Gracias de antemano