Astropy: fits.getdataはNaNを返します

作成日 2017年08月28日 · 3コメント · ソース: astropy/astropy

やあ、

astropy 2.0.1、Python2.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

同じ結果を返します。ただし、SAOImageds9はデータを正常に開いて表示できます。

io.fits

ソース

Amarchuk

最も参考になるコメント

その画像の端には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 2017年08月28日

👍3

全てのコメント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 2017年08月28日

👍3

おっと、私のせい。
@ MSeifert04この超高速の回答を

Amarchuk 2017年08月28日

問題ない。データで頑張ってください:)

MSeifert04 2017年08月28日

このページは役に立ちましたか？

0 / 5 - 0 評価