fits routines with error handling

[martinkilbinger]

We should implement some functions to handle FITS files with better error messaging.
E.g. a wrapper for astropy.io.fits.getdata that checks whether the HDU number exists in the FITS file.
(Standard astropy error is just an IndexError).

Use this in e.g. shear_psf_leakage.