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.
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()
Opa, minha culpa.
@ MSeifert04 Muito obrigado por esta resposta super-rápida!
Sem problemas. Boa sorte com os dados :)
Comentários muito úteis
Essa imagem apenas tem muitos
nan
s nas bordas. Parece bom, porém: