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
|
/*******************************************************************************
* Johns Hopkins University
* Center For Astrophysical Sciences
* FUSE
*******************************************************************************
*
* Synopsis: cf_apply_fes_mask(fitsfile *fesfits, float **image,
* float *mask, int axis1, int axis2)
*
* Description: Apply FES mask in *mask to the raw FES image in **image
* To apply the mask I simply loop through the
* mask image pixel by pixel and
* If the mask pixel value != 0.0
* then fes image pixel value = mask pixel value.
*
* Arguments: fitsfile *fesfits fits file for image, so that we can add
* history/comment lines to the HDU
* float **image FES image
* float *mask FES mask
* int axis1, axis2 size of each axis of the images.
* This presumes that both images are
* the same size when they get here.
*
* Returns: 0 upon successful completion
*
* History: 07/08/98 gak calfes_design.070898 design documented
* 07/22/98 mlr started work
* 04/06/99 mlr Broke subroutine out from cf_fes_mask.c
* 08/10/99 mlr added fesfits to the argument list
*
*
******************************************************************************/
#include "calfuse.h"
#include "cf_calfes.h"
#define CF_PRGM_ID "cf_fes_apply_mask"
#define CF_VER_NUM "1.4"
int cf_fes_apply_mask(fitsfile *fesfits, float **image,
float *mask, int naxis1, int naxis2)
{
int i, j; /*loop counters */
cf_timestamp(CF_PRGM_ID, CF_VER_NUM, "Applying FES pixel mask");
for (j=0; j < naxis2; j++)
for (i=0; i < naxis1; i++)
if (mask[j*naxis1+i] != FES_GOOD_PIX )
(*image)[j*naxis1+i] = mask[j*naxis1+i];
cf_timestamp(CF_PRGM_ID, CF_VER_NUM, "Done applying FES mask");
return(0);
}
|
