파일:Light_Curve_of_RT_Aurigae_Observed_by

1. 기본 정보2. 이미지 설명

1. 기본 정보

출처	TESS 관측 데이터를 이용해 게시자가 직접 작성 관측 ID: tess2023289093419-s0071-0000000171855358-0266-s
날짜	2024/7/11
저작자	게시자 본인
저작권	퍼블릭 도메인
기타	기타 정보가 있으면 삽입해 주세요.

2. 이미지 설명

고전 세페이드 변광성인 마차부자리 RT의 광도곡선 및 파워스펙트럼, TESS 관측 데이터를 내려받아 matplotlib를 활용하여 작성됨.

다운로드 링크: https://mast.stsci.edu/portal/Mashup/Clients/Mast/Portal.html
위 사이트로 접속하여 원하는 천체의 이름을 검색하면, TESS를 통해 관측한 측광 데이터를 포함한 각종 천문관측 데이터를 자유롭게 내려받는 것이 가능함.

작성에 활용된 파이썬 소스코드는 아래와 같음. 이 소스코드를 활용해 광도곡선을 그리고자 할 경우, 내려받은 압축 폴더의 압축을 푼 다음, 폴더 안의 광도곡선 데이터가 저장된 fits 파일(확장자를 제외한 파일명이 lc로 끝나는 fits 파일, 이 데이터의 경우 tess2023289093419-s0071-0000000171855358-0266-s_lc.fits)이 저장된 경로에 해당 소스코드를 파일로 작성해 저장한 후 코드를 실행하면 그래프가 png 이미지로 저장됨.

#!syntax python
import numpy as np
import matplotlib.pyplot as plt
from astropy.io import fits
from astropy.timeseries import LombScargle as ls

hdul=fits.open('tess2023289093419-s0071-0000000171855358-0266-s_lc.fits')

t=hdul[1].data['TIME']
fop=hdul[1].data['SAP_FLUX']
feop=hdul[1].data['SAP_FLUX_ERR']

xunit=hdul[1].header['TUNIT1']
yunit=hdul[1].header['TUNIT8']

dnan=np.isfinite(fop)

t1=t[dnan]
fop1=fop[dnan]
feop1=feop[dnan]

fig,ax=plt.subplots(2,1,figsize=(18,12))
plt.subplots_adjust(hspace=0.35)

ax[0].errorbar(t,fop,feop,label='RT Aur',linewidth=0.7,elinewidth=0.7)

f,p=ls(t1,fop1).autopower()
f1=1/f
ax[1].plot(f1,p,label='Power Spectrum')
ax[1].set_xlim(0,7.5)

ax[0].set_xlabel(xunit,fontsize=15)
ax[0].set_ylabel('Instrumental Flux (%s)'%str(yunit),fontsize=15)
ax[0].legend(loc='upper left',fontsize=17)
ax[0].tick_params(labelsize=15)
ax[0].set_title('Light Curve of RT Aurigae',fontsize=20)
ax[0].grid(color='gray',linestyle=':',linewidth=0.9)
ax[1].set_xlabel('period (days)',fontsize=15)
ax[1].set_ylabel('power',fontsize=15)
ax[1].text(0.2,0.5,'period = %s days' %str(round(f1[p.argmax()],4)),fontsize=19)
ax[1].legend(loc='upper left',fontsize=17)
ax[1].tick_params(labelsize=15)
ax[1].set_title('Power Spectrum of RT Aurigae',fontsize=20)
ax[1].grid(color='gray',linestyle=':',linewidth=0.9)
plt.savefig('Light_Curve_of_RT_Aurigae_Observed_by_TESS.png',facecolor='#ffffff',bbox_inches='tight',dpi=700)

파일:Light_Curve_of_RT_Aurigae_Observed_by_TESS.png

1. 기본 정보

2. 이미지 설명

분류