use astropy.io.fits in preference to pyfits