from astropy.table import Table
import pandas as pd
import matplotlib.pyplot as plt
from datetime import datetime,timedelta,date
pd.set_option('display.max_columns',None)
"""
费米卫星每列数据
START:从2001年1月1日00时00分00秒开始到此刻的总秒数,不算闰秒
STOP:
"""
class FermiPlot():
def __init__(self,fitsName):
self.fitsName=fitsName
self.data=Table.read(self.fitsName).to_pandas()
print(self.data.columns)
print(self.data)
self.data["START"]=self.data["START"].apply(self.modifiedJulianDay)
print(self.data)
def plotFermi(self):
data=self.data
fig,ax=plt.subplots(figsize=(10,5))
yInterval=plt.MultipleLocator(5*10**(-7))
ax.yaxis.set_major_locator(yInterval)
ax.set_ylim((0,3*10**(-6)))
ax.scatter(data['START'],data['FLUX_100_300000'])
ax.set_xlim((54000,60000))
画费米卫星光变曲线
最新推荐文章于 2025-05-28 20:48:46 发布