1 files changed, 56 insertions, 0 deletions
diff --git a/src/fes/cf_fes_apply_mask.c b/src/fes/cf_fes_apply_mask.c
new file mode 100644
index 0000000..c083121
--- /dev/null
+++ b/src/fes/cf_fes_apply_mask.c
@@ -0,0 +1,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);
+}