can now use either pyfits or astropy for fits handling