from astropy.io import fits
import numpy as np

# Configura numpy per visualizzare tutti i dati senza troncamenti
np.set_printoptions(threshold=np.inf)

# Percorso del file FITS
file_path = '/home/alfonso/Scrivania/THESEUS/xgis_m7-main/fits_matrix/matrix.fits'

# Apri il file FITS
with fits.open(file_path) as hdul:
    # Stampa un sommario delle HDU nel file
    hdul.info()
    
    # Itera su tutte le HDU e stampa l'header e i dati
    for idx, hdu in enumerate(hdul):
        print(f"\nHeader HDU {idx}:")
        print(hdu.header)
        print("\nDati:")
        print(hdu.data)  # Stampa i dati completi senza troncamenti

