def create_mock_fits(): x = np.ones((5, 5)) prihdu = fits.PrimaryHDU(x) # Single extension FITS img = fits.ImageHDU(data=x) singlehdu = fits.HDUList([prihdu, img]) singlehdu.writeto('image.fits', clobber=True)