diff options
Diffstat (limited to 'src/fes/cf_fes_cal.c')
| -rw-r--r-- | src/fes/cf_fes_cal.c | 115 |
1 files changed, 115 insertions, 0 deletions
diff --git a/src/fes/cf_fes_cal.c b/src/fes/cf_fes_cal.c new file mode 100644 index 0000000..0f08bbd --- /dev/null +++ b/src/fes/cf_fes_cal.c @@ -0,0 +1,115 @@ +/******************************************************************************* + * Johns Hopkins University + * Center For Astrophysical Sciences + * FUSE + ******************************************************************************* + * + * Synopsis: int cf_fes_cal(int cal_type, fitsfile *fesfits, float **image, + * int inaxis1, int inaxis2) + * + * Description: Get bias calibration filename for the raw FES image in + * the current hdu pointed at by *infits. + * Call cf_fes_get_cal_image() to open the cal file and get + * the calibration image. + * Call cf_fes_apply_flat(); + * return the modified image to the calling routine. + * 4/10/99 MLR removed int inhdu as an argument and removed + * the call to fits_movabs_hdu. This requires any + * user of this subroutine to do these steps before + * calling cf_fes_flat. + * 4/10/99 MLR removed (fitsfile *outfits, int outhdu) from + * the argument list. We will leave writing the image + * out to the calling routine. + * + * + * Arguments: fitsfile *infits Input FITS file pointer + * float **image already read in FES image + * int inaxis1 image size of the input image + * int inaxis2 + * + * Returns: none + * + * History: 07/08/98 gak calfes_design.070898 design documented + * 07/22/98 mlr started work + * 04/10/99 mlt pulled the actual work out of this + * routine and into cf_fes_apply_flat + * 04/12/99 mlr modified the input arguments(see above) + * use the now generic cf_fes_read to + * get the calibration image. + * 04/14/99 mlr pulled the reading out and put it + * into cf_fes_get_cal_image.c + * 04/16/99 mlr finished modifications to utilize + * libcf and FITSIO.h(error handling). + * + * ToDo: 4/12/99 I have added the naxis1 and naxis2 parameters to the + * argument list. I need to use them to compare binning + * factors in the two images. + * I should combine cf_fes_mask, cf_fes_bias, cf_fes_flat + * into one generic cf_fes_do_cal("cf_fes_mask", *fesfits,...etc) + * + ******************************************************************************/ + +#include <string.h> +#include <stdio.h> +#include <stdlib.h> +#include "calfuse.h" +#include "cf_calfes.h" + +#define CF_VER_NUM "1.4" + +int cf_fes_cal(int cal_type, fitsfile *fesfits, float **image, + int inaxis1, int inaxis2) +{ + char buffer[FLEN_CARD]; + char calfilename[30]; + fitsfile *calfile; + int status; + float *calimage; + char *prog_id, *msg, *keyword, *msg2; + + int (*func)(fitsfile *, float **, float *, int, int); + + status = 0; + + switch (cal_type) { + case FES_FLAT: + prog_id = "cf_fes_flat"; + msg = "Begin FES flatfield"; + msg2 = "Done FES flatfield"; + keyword = "FFLT1FCL"; + func = cf_fes_apply_flat; + break; + case FES_MASK: + prog_id = "cf_fes_mask"; + msg = "Begin FES mask"; + msg2 = "Done FES mask"; + keyword = "MASK_FCL"; + func = cf_fes_apply_mask; + break; + case FES_BIAS: + prog_id = "cf_fes_bias"; + msg = "Begin FES bias"; + msg2 = "Done FES bias"; + keyword = "BIAS1FCL"; + func = cf_fes_apply_bias; + break; + } + + + cf_timestamp(prog_id, CF_VER_NUM, msg); + + FITS_read_key(fesfits, TSTRING, keyword, calfilename, buffer, &status); + + if (cf_fes_get_cal_image(calfilename, &calimage, inaxis1, inaxis2) == -1) + return (-1); + + (*func)(fesfits, image, calimage, inaxis1, inaxis2); + + cf_fes_proc_update(fesfits, prog_id, "COMPLETE"); + + free(calimage); + + cf_timestamp(prog_id, CF_VER_NUM, msg2); + + return(status); +} |