Astropy: fit.getdata gibt NaNs zurück

Erstellt am 28. Aug. 2017  ·  3Kommentare  ·  Quelle: astropy/astropy

Hi,

Ich verwende Astropy 2.0.1, Python 2.7. Wenn ich versuche, FITS-Daten in das Array zu bekommen ( passt auf den Dateilink ) von

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

Ich bekomme ein numpy Array voller 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)

Als HDU-Liste öffnen

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

das gleiche Ergebnis zurückgeben. Aber SAOImage ds9 kann erfolgreich Daten öffnen und anzeigen.

io.fits

Hilfreichster Kommentar

Dieses Bild hat nur viele nan s an den Rändern. Sieht aber ok aus:

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

Alle 3 Kommentare

Dieses Bild hat nur viele nan s an den Rändern. Sieht aber ok aus:

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

Ups, meine Schuld.
@MSeifert04 Vielen Dank für diese superschnelle Antwort!

Kein Problem. Viel Glück mit den Daten :)

War diese Seite hilfreich?
0 / 5 - 0 Bewertungen