ASCENT-zouho2022
25/32

データ読み込み解析結果のグラフ表示 観測データが MAXIJ1820.csv というファイルに保存されているとき、次のプログラムで読み込むことができます。 matplotlibを用いて解析結果をグラフ表示できます。 import pandas as pd は、pandasを使うためのおまじないです。pd.read_csvは、データを読み込む命令です。読み込んだデータはdfという名前で参照できます。 読み込んだデータの各列はdf[‘4-10keV’]のように各列のラベルを用いて取り出せます。y=df[‘4-10keV’]/df[‘2-4keV’]とすれば、各エネルギー範囲の光子数の比が求まります。このプログラムを実行すると、右図のように、低エネルギー(2-4keV)と高エネルギー(4-10keV)の光子数の時間変化が表示されます。既に公開されているデータの解析からでも、新しい発見があるかもしれないよ!import pandas as pddf = pd.read_csv('./MAXIJ1820.csv')df.head()import matplotlib.pyplot as pltdf.plot(x='MJD',y=['4-10keV','2-4keV'])plt.ylabel('photons cm$^{-2}$ s$^{-1}$',fontsize=12)plt.show()25

元のページ  ../index.html#25

このブックを見る