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
|
/*******************************************************************************
* Johns Hopkins University
* Center For Astrophysical Sciences
* FUSE
*******************************************************************************
*
* Synopsis:int cf_fes_get_cal_image(char cal_file_name[], float **cal_image,
* int naxis1, int naxis2)
*
* Description: open the specified calibration file for an FES image
* loop through each extension in the calibration file
* until you find the one that has the same size naxis1 and
* naxis2 as the input image.
* read the image in that extention into *cal_image
* close the calibration file
*
* Arguments: char cal_file_name[],
* float **cal_image,
* int naxis1, int naxis2
*
* Returns: none
*
* History: 04/14/99 mlr started work
* 04/16/99 mlr modified to use libcf and FITSIO.h
*
******************************************************************************/
#include "calfuse.h"
#define CF_PRGM_ID "cf_fes_get_cal_image"
#define CF_VER_NUM "1.4"
int cf_fes_get_cal_image(char cal_file_name[], float **cal_image,
int naxis1, int naxis2)
{
fitsfile *fes_cal_file;
int status;
int i, hdu, hdutype, num_exts;
int cal_axis1, cal_axis2;
char buffer[FLEN_CARD];
cf_timestamp(CF_PRGM_ID, CF_VER_NUM, "Begin getting Cal image");
status = 0;
/* open the input FES CALIBRATION file */
FITS_open_file(&fes_cal_file,cf_cal_file(cal_file_name),READONLY,&status);
FITS_read_key(fes_cal_file, TINT, "NEXTEND", &num_exts, buffer, &status);
for (i = 1; i<= num_exts; i++)
{
hdu = i + 1;
FITS_movabs_hdu(fes_cal_file, hdu, &hdutype, &status);
cf_fes_read(fes_cal_file, cal_image, &cal_axis1, &cal_axis2);
if ((cal_axis1 == naxis1) && (cal_axis2 == naxis2))
break;
}
FITS_close_file(fes_cal_file, &status);
if (i>num_exts)
{
*cal_image=NULL;
cf_timestamp(CF_PRGM_ID, CF_VER_NUM, "No Cal image of requested size");
return(-1);
}
cf_timestamp(CF_PRGM_ID, CF_VER_NUM, "Done getting Cal image");
return(status);
}
|
