Flexible Image Transport System (FITS)

from astropy.io import fits
import matplotlib.pyplot as plt
from matplotlib.colors import LogNorm
import glob
import pandas as pd

hdul = {}
df = []
for file in sorted(glob.glob('fits/*.fits')):
    tis = fits.open(file)
    hdul[file] = tis
    for nn in [1, 2, 3]:
        assert tis[nn].shape==(3400, 3400)
    label = file.split('coadd_')[1].split('.')[0]
    band = label.split('_')[0]
    field = label.split('_')[1]
    df.append([field, band, len(tis)])
df = pd.DataFrame(df, columns=['field', 'band', 'length'])
df

for kk, vv in hdul.items():
    for mm in [1, 2, 3]:
        label = kk.split('coadd_')[1].split('.')[0]
        band = label.split('_')[0]
        field = label.split('_')[1]
        plt.imshow(vv[mm].data, norm=LogNorm(), cmap='gray')
        plt.savefig(f'png/butler_deep_coadd_{field}_{band}_{mm}.png', bbox_inches='tight')
        plt.axis('off'); plt.close()