์๋ ํ์ธ์,
๋๋ astropy 2.0.1, Python 2.7์ ์ฌ์ฉํ๊ณ ์์ต๋๋ค. FITS ๋ฐ์ดํฐ๋ฅผ ๋ฐฐ์ด๋ก ๊ฐ์ ธ์ค๋ ค๊ณ ํ ๋( ํ์ผ ๋งํฌ์ ๋ง์ )
from astropy.io import fits
image_data = fits.getdata('NGC_628_NA_MOM2_THINGS.FITS')
๋๋ NaN์ผ๋ก ๊ฐ๋ ์ฐฌ numpy ๋ฐฐ์ด์ ์ป์ต๋๋ค.
array([[[[ nan, nan, nan, ..., nan, nan, nan],
[ nan, nan, nan, ..., nan, nan, nan],
[ nan, nan, nan, ..., nan, nan, nan],
...,
[ nan, nan, nan, ..., nan, nan, nan],
[ nan, nan, nan, ..., nan, nan, nan],
[ nan, nan, nan, ..., nan, nan, nan]]]], dtype=float32)
HDU ๋ชฉ๋ก์ผ๋ก ์ด๊ธฐ
hdu_list = fits.open('NGC_628_NA_MOM2_THINGS.FITS')
image_data = hdu_list[0].data
๊ฐ์ ๊ฒฐ๊ณผ๋ฅผ ๋ฐํํฉ๋๋ค. ๊ทธ๋ฌ๋ SAOImage ds9๋ ์ฑ๊ณต์ ์ผ๋ก ๋ฐ์ดํฐ๋ฅผ ์ด๊ณ ํ์ํ ์ ์์ต๋๋ค.
๊ทธ ์ด๋ฏธ์ง์๋ ๊ฐ์ฅ์๋ฆฌ์ ๋ง์ nan
์์ต๋๋ค. ๊ทธ๋๋ ๊ด์ฐฎ์ ๋ณด์
๋๋ค.
from astropy.io import fits
import matplotlib.pyplot as plt
data = fits.getdata('NGC_628_NA_MOM2_THINGS.FITS')
plt.imshow(np.squeeze(data), origin='lower', cmap='gray')
plt.show()
์, ๋ด ์๋ชป์ด์ผ.
@MSeifert04 ์ด ์ด๊ณ ์ ๋ต๋ณ์ ๊ฐ์ฌ๋๋ฆฝ๋๋ค!
๊ด์ฐฎ์์. ๋ฐ์ดํฐ ์ ํ์ฉํ์ธ์ :)
๊ฐ์ฅ ์ ์ฉํ ๋๊ธ
๊ทธ ์ด๋ฏธ์ง์๋ ๊ฐ์ฅ์๋ฆฌ์ ๋ง์
nan
์์ต๋๋ค. ๊ทธ๋๋ ๊ด์ฐฎ์ ๋ณด์ ๋๋ค.