Needs pyfits!