Hola,
Estoy usando astropy 2.0.1, Python 2.7. Cuando trato de obtener datos FITS en la matriz (se ajusta al enlace del archivo ) por
from astropy.io import fits
image_data = fits.getdata('NGC_628_NA_MOM2_THINGS.FITS')
Obtengo una matriz llena 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 HDU
hdu_list = fits.open('NGC_628_NA_MOM2_THINGS.FITS')
image_data = hdu_list[0].data
devuelve el mismo resultado. Pero SAOImage ds9 puede abrir y mostrar datos con éxito.
Esa imagen solo tiene mucho nan
s en los bordes. Aunque se ve bien:
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()
Vaya, culpa mía.
@ MSeifert04 ¡ Muchas gracias por esta respuesta súper rápida!
No hay problema. Buena suerte con los datos :)
Comentario más útil
Esa imagen solo tiene mucho
nan
s en los bordes. Aunque se ve bien: