pyfits -> astropy.io.fits