In [14]:

    

%matplotlib inline


    

In [12]:

    

import numpy as np
import matplotlib.pyplot as plt
from astropy.io import ascii


    

In [10]:

    

data = ascii.read('/home/glauffer/Dropbox/FURG/final_project/data/kplr001026032-2009166043257_llc.dat')
data
# Time (days), Raw Flux, Raw Flux error, Corrected Flux, Corrected Flux error


    

    
Out[10]:




col1
col2
col3
col4
col5
54964.5127135
12319.002
4.01829672
17661.9863
5.71493769
54964.5331482
12321.8613
4.01851368
17651.8203
5.7225852
54964.5535829
12317.9189
4.01818419
17646.8379
5.72262812
54964.5740174
12307.4531
4.01739311
17634.7559
5.71586132
54964.594452
12308.417
4.01751661
17634.8652
5.72388887
54964.6148867
12315.8311
4.01795244
17649.2188
5.719244
54964.6353212
12316.9014
4.01794434
17653.3203
5.71940374
54964.6557558
12316.9521
4.01808023
17655.0664
5.72260284
54964.6761905
12307.1855
4.01737928
17641.3965
5.72496271
54964.6966251
12309.168
4.01743841
17644.1133
5.72099781
54964.7170596
12318.3076
4.01818132
17658.6465
5.71628428
...
...
...
...
...
54997.7800308
12541.0615
4.0330925
17682.334
5.7287097
54997.8004652
12551.2129
4.03390551
17700.7715
5.73312283
54997.8208995
12539.7217
4.03310013
17682.5957
5.73356152
54997.8413337
12548.2871
4.03369236
17693.873
5.72497129
54997.8617681
12542.6279
4.03333378
17680.5625
5.72649193
54997.8822024
12553.9414
4.03419304
17690.6582
5.72741556
54997.9026366
12541.3213
4.03326321
17676.8809
5.7296176
54997.9230709
12549.7539
4.03379345
17684.1777
5.73302174
54997.9435053
12551.208
4.03387833
17688.9473
5.72724915
54997.9639395
12547.0566
4.0336175
17679.7754
5.73365974
54997.9843738
12547.208
4.03364801
17688.3789
5.72657919

In [15]:

    

plt.plot(data['col1'], data['col4'])


    

    
Out[15]:





[<matplotlib.lines.Line2D at 0x7f0c1a473cf8>]






    

In [19]:

    

#convertendo fluxo em magnitude
mag = -2.5 * np.log(data['col4'])
mag


    

    
Out[19]:





<Column name='col4' unit=None format=None description=None>
array([-24.44792486, -24.44648548, -24.44577973, ..., -24.45173819,
       -24.45044158, -24.45165786])

In [22]:

    

plt.plot(data['col1'], mag)
plt.gca().invert_yaxis()
plt.xlabel('Dias [JD]')
plt.ylabel('Magnitude')


    

    
Out[22]:





<matplotlib.text.Text at 0x7f0c198cc898>






    

In [28]:

    

kpl6b = ascii.read('/home/glauffer/Dropbox/FURG/final_project/data/kplr010874614-2009131105131_llc.dat')
#kpl6b
mag6b = -2.5 * np.log(kpl6b['col4'])

plt.clf()
plt.plot(kpl6b['col1'], mag6b)
#plt.plot(kpl6b['col1'], kpl6b['col4'])
plt.gca().invert_yaxis()
plt.title('Kepler 6B - KIC 10874614')
plt.xlabel('Dias [JD]')
plt.ylabel('Magnitude')


    

    
Out[28]:





<matplotlib.text.Text at 0x7f0c19dd9e10>






    

In [29]:

    

#supernova KSN 2011a
sn = ascii.read('/home/glauffer/Dropbox/FURG/final_project/data/kplr008480662-2011271113734_llc.dat')
magSN = -2.5 * np.log(sn['col4'])

plt.clf()
plt.plot(sn['col1'], magSN)
plt.gca().invert_yaxis()
plt.title('KSN 2011a - KIC 8480662')
plt.xlabel('Dias [JD]')
plt.ylabel('Magnitude')


    

    
Out[29]:





<matplotlib.text.Text at 0x7f0c1a63aac8>






    

In [ ]: