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.
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()
Ups, meine Schuld.
@MSeifert04 Vielen Dank für diese superschnelle Antwort!
Kein Problem. Viel Glück mit den Daten :)
Hilfreichster Kommentar
Dieses Bild hat nur viele
nan
s an den Rändern. Sieht aber ok aus: