Astropy: fit.getdata๋Š” NaN์„ ๋ฐ˜ํ™˜ํ•ฉ๋‹ˆ๋‹ค.

์— ๋งŒ๋“  2017๋…„ 08์›” 28์ผ  ยท  3์ฝ”๋ฉ˜ํŠธ  ยท  ์ถœ์ฒ˜: astropy/astropy

์•ˆ๋…•ํ•˜์„ธ์š”,

๋‚˜๋Š” 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()

png image 573 x 429 pixels

๋ชจ๋“  3 ๋Œ“๊ธ€

๊ทธ ์ด๋ฏธ์ง€์—๋Š” ๊ฐ€์žฅ์ž๋ฆฌ์— ๋งŽ์€ 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()

png image 573 x 429 pixels

์•—, ๋‚ด ์ž˜๋ชป์ด์•ผ.
@MSeifert04 ์ด ์ดˆ๊ณ ์† ๋‹ต๋ณ€์— ๊ฐ์‚ฌ๋“œ๋ฆฝ๋‹ˆ๋‹ค!

๊ดœ์ฐฎ์•„์š”. ๋ฐ์ดํ„ฐ ์ž˜ ํ™œ์šฉํ•˜์„ธ์š” :)

์ด ํŽ˜์ด์ง€๊ฐ€ ๋„์›€์ด ๋˜์—ˆ๋‚˜์š”?
0 / 5 - 0 ๋“ฑ๊ธ‰