def to_fits(self): primary_hdu = fits.PrimaryHDU() image_hdu = fits.ImageHDU(self.data, self.header) hdulist = fits.HDUList([primary_hdu, image_hdu]) return hdulist