Use astropy.io.fits with fallback to pyfits