만약 저작권을 포기하는 것이 법률상 불가능하다면 저작권자는 이 저작물을 누구나 자유로이 어떤 목적으로든 제한 없이 사용할 수 있도록 허용합니다.
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)