Changeset 23463 for trunk/ppMerge/src/ppMergeLoop.c

Timestamp:

Mar 20, 2009, 5:56:03 PM (17 years ago)

Author:

eugene

Message:

add ability to handle CTEMASK mode (combine like a flat; mask pixels below a minimum value; output a mask image)

File:

• trunk/ppMerge/src/ppMergeLoop.c (modified) (4 diffs)

trunk/ppMerge/src/ppMergeLoop.c

@@ -118 +118 @@
       case PPMERGE_TYPE_BIAS:
       case PPMERGE_TYPE_DARK:
+      case PPMERGE_TYPE_CTEMASK:
         break;
       default:
@@ -208 +209 @@
             // call the init functions
             switch (type) {
+              case PPMERGE_TYPE_CTEMASK:
               case PPMERGE_TYPE_BIAS:
               case PPMERGE_TYPE_FLAT:
@@ -252 +254 @@
                 // Start a job
                 switch (type) {
+                  case PPMERGE_TYPE_CTEMASK:
                   case PPMERGE_TYPE_BIAS:
                   case PPMERGE_TYPE_FLAT:
@@ -408 +411 @@
             }
 
+            // calculate CTEMASK after stats so stats reflect median image
+            if (type == PPMERGE_TYPE_CTEMASK && outRO) {
+                // need to apply range cuts on the output image
+                psAssert (outRO->mask, "mask is not defined");
+                psAssert (outRO->mask->numCols == outRO->image->numCols, "mismatch between image and mask");
+                psAssert (outRO->mask->numRows == outRO->image->numRows, "mismatch between image and mask");
+
+                // CTEMASK parameters
+                float cteMin = psMetadataLookupF32(NULL, arguments, "CTE.MIN"); // Number of fringe points
+
+                char *cteMaskName = psMetadataLookupStr (&mdok, config->arguments, "MASK.SET.VALUE");
+                psImageMaskType cteMaskValue = pmConfigMaskGet(cteMaskName, config);
+
+                psF32 **outputImage = outRO->image->data.F32; 
+                psImageMaskType **outputMask = outRO->mask->data.PS_TYPE_IMAGE_MASK_DATA; 
+                for (int iy = 0; iy < outRO->image->numRows; iy++) {
+                    for (int ix = 0; ix < outRO->image->numCols; ix++) {
+                        if (outputImage[iy][ix] < cteMin) {
+                            outputMask[iy][ix] |= cteMaskValue;
+                        }
+                    }
+                }
+            }
+
             psFree(outRO);
             cellNum++;