Astropy: fit.getdata renvoie des NaN

Créé le 28 août 2017 · 3Commentaires · Source: astropy/astropy

Salut,

J'utilise Astropy 2.0.1, Python 2.7. Lorsque j'essaie d'obtenir des données FITS dans un tableau ( convient au lien de fichier ) en

from astropy.io import fits
image_data = fits.getdata('NGC_628_NA_MOM2_THINGS.FITS')

Je reçois un tableau numpy plein de NaN

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)

Ouvrir en tant que liste HDU

hdu_list = fits.open('NGC_628_NA_MOM2_THINGS.FITS')
image_data = hdu_list[0].data

renvoie le même résultat. Mais SAOImage ds9 peut ouvrir et afficher des données avec succès.

io.fits

Source

Amarchuk

Commentaire le plus utile

Cette image a juste beaucoup de nan sur les bords. ça a l'air pas mal quand même :

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 le 28 août 2017

👍3

Tous les 3 commentaires

Cette image a juste beaucoup de nan sur les bords. ça a l'air pas mal quand même :

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 le 28 août 2017

👍3

Oups, ma faute.
@MSeifert04 Merci beaucoup pour cette réponse ultra rapide !

Amarchuk le 28 août 2017

Aucun problème. Bonne chance avec les données :)

MSeifert04 le 28 août 2017

Cette page vous a été utile?

0 / 5 - 0 notes