1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
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);
}
|
