Astropy: fit.getdata retorna NaNs

Criado em 28 ago. 2017  ·  3Comentários  ·  Fonte: astropy/astropy

Oi,

Estou usando astropy 2.0.1, Python 2.7. Quando tento obter dados FITS no array ( se encaixa no link do arquivo ) por

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

Eu obtenho um array numpy cheio de NaNs

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)

Abrir como lista de HDU

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

retorna o mesmo resultado. Mas SAOImage DS9 pode abrir e mostrar dados com sucesso.

io.fits

Comentários muito úteis

Essa imagem apenas tem muitos nan s nas bordas. Parece bom, porém:

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

Todos 3 comentários

Essa imagem apenas tem muitos nan s nas bordas. Parece bom, porém:

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

Opa, minha culpa.
@ MSeifert04 Muito obrigado por esta resposta super-rápida!

Sem problemas. Boa sorte com os dados :)

Esta página foi útil?
0 / 5 - 0 avaliações