Use astropy not pyfits