Astropy: fit.getdata devuelve NaN

Creado en 28 ago. 2017  ·  3Comentarios  ·  Fuente: astropy/astropy

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.

io.fits

Comentario más útil

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()

png image 573 x 429 pixels

Todos 3 comentarios

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()

png image 573 x 429 pixels

Vaya, culpa mía.
@ MSeifert04 ¡ Muchas gracias por esta respuesta súper rápida!

No hay problema. Buena suerte con los datos :)

¿Fue útil esta página
0 / 5 - 0 calificaciones