어린이망원경

적외선으로 만나는 우주 4 - 허셜 우주망원경 관측 데이터 불러오기

페이지 정보

작성자 깡지쌤 / GETA 강원지구… 작성일21-10-13 00:00 조회5회 댓글0건

본문



허셜 우주망원경이 관측한 데이터를 불러오는 파이썬 코드를 배워볼 것입니다.
전체 코드는 서지자료와 여기에 올려두겠으니 필요하신 분은 드래그하여 가져가시면 됩니다.

import matplotlib.pyplot as plt
from astropy.wcs import WCS
from astropy.io import fits
from astropy.utils.data import get_pkg_data_filename
import numpy as np

filename = get_pkg_data_filename('C:\\\\111\\\\IC434.fits')
hdu = fits.open(filename)[1]
wcs = WCS(hdu.header)
fig = plt.figure(figsize=(18, 12))

plt.subplot(projection=wcs)
plt.grid(color='white', linestyle='--')
bar=plt.imshow(hdu.data, origin='lower', cmap=plt.cm.gist_heat, vmin=0, vmax=1000)
plt.colorbar(bar, label='Herschel 250'r'$\\mu$m (MJy/sr)')
plt.xlabel('RA')
plt.ylabel('Dec')
plt.title('IC 434 dust continuum map(250'r'$\\mu$m)', size=15)
plt.savefig('C:\\\\Users\\\\kalch\\\\Desktop\\\\IC434250.png',
bbox_inches='tight', pad_inches=0.1, dpi=300)
levelinter=np.linspace(0,850,3)
plt.contour(hdu.data, linewidths=2, colors='white', linestyles="--", levels=levelinter)
print(hdu.header)
plt.show()

댓글목록

등록된 댓글이 없습니다.

Total 3,381건 1 페이지
게시물 검색
Copyright © www.estore365.kr. All rights reserved.  Contact : help@oxmail.xyz