change from pyfits to astropy