Le fichier fit est un format de fichier image utilisé pour les données d'exploration et les données astronomiques.Il se compose d'une partie appelée en-tête qui contient diverses informations telles que l'heure et la posture, et une partie tableau de la taille de l'image (par exemple, 1024x1024). Il existe également un fichier ajustements en forme de cube, qui se compose de n en-têtes et d'un tableau d'images. https://ja.wikipedia.org/wiki/FITS À propos de la lecture et de l'affichage des données d'ajustement à l'aide de la bibliothèque Python astropy.io.fits. DS9 est l'un des téléspectateurs des ajustements.
PDS Planetary Data System Archive des données d'exploration de la NASA. Tout le monde peut l'utiliser gratuitement. Si vous téléchargez les données d'ici, vous pouvez voir l'image de votre planète préférée! Les images sont publiées au format de fichier ajusté.
Par exemple, les données pour Hayabusa sont ici. Toutes les données d'image prises par Hayabusa sont ouvertes au public.
Installez avec anaconda.
conda install astropy
Lisez les 0èmes données. L'image se compose d'un en-tête + un tableau d'images.
import astropy.io.fits as fits
import matplotlib.pyplot as plt
hdulist=pyfits.open('file.fits')
hdu=hdulist[0]
data=hdu.data #data=fits.getdata('file.fits',0)Mais tu peux.
header=hdu.header
plt.imshow(data)
plt.show()
Lors de la récupération d'un élément d'en-tête, par exemple TI_0
print header["TI_0"]
Spécification des pixels X et Yth
data[Y][X]
hdu = fits.PrimaryHDU(im)
hdulist = fits.HDUList([hdu])
hdulist.writeto('new.fits',overwrite=True)
correspond est un ensemble (hdu) avec en-tête et tableau d'images. Il se compose de hdu = PrimaryHDU (data, header) ''. Lorsque vous utilisez un grand nombre d'extensions, créez plusieurs hdu à partir de la seconde avec
fit.ImageHDU '', et créez-en une avec le sentiment de HDUList ([hdu1, hdu2, hdu3])
. Faites-en un fichier. Autrement dit, il est possible d'enregistrer une pluralité d'images dans un fichier.
Il peut être écrasé en définissant overwrite = True.
hdu.header['NEWKEYWORD']='something'
Recommended Posts