Changeset 16413

Timestamp:

Feb 12, 2008, 3:41:06 PM (18 years ago)

Author:

Paul Price

Message:

When statistics are not available, it means that the statistics failed --- mask entire image.

File:

• trunk/psModules/src/detrend/pmMaskBadPixels.c (modified) (5 diffs)

trunk/psModules/src/detrend/pmMaskBadPixels.c

@@ -104 +104 @@
     psImage *mask = readout->mask;      // Corresponding mask
 
-    float median = psMetadataLookupF32 (&status, readout->analysis, "READOUT.MEDIAN");
-    if (!status) {
-        psError(PS_ERR_PROGRAMMING, true, "Missing READOUT.MEDIAN from readout analysis data");
-        return NULL;
-    }
-
-    float stdev  = psMetadataLookupF32 (&status, readout->analysis, "READOUT.STDEV");
-    if (!status) {
-        psError(PS_ERR_PROGRAMMING, true, "Missing READOUT.STDEV from readout analysis data");
-        return NULL;
-    }
-
-    psTrace ("psModules.detrend", 3, "suspect: %f +/- %f\n", median, stdev);
-
     if (!out) {
         out = psImageAlloc(image->numCols, image->numRows, PS_TYPE_S32);
         psImageInit(out, 0);
+    }
+
+    bool whole = false;                 // Mask whole image?
+    float median = psMetadataLookupF32 (&status, readout->analysis, "READOUT.MEDIAN");
+    if (!status || !isfinite(median)) {
+        whole = true;
+    }
+    float stdev  = psMetadataLookupF32 (&status, readout->analysis, "READOUT.STDEV");
+    if (!status || !isfinite(stdev)) {
+        whole = true;
+    }
+
+
+    psTrace ("psModules.detrend", 3, "suspect: %f +/- %f\n", median, stdev);
+
+    if (whole) {
+        // If we get down here and the statistics are missing, then we should go and mask the entire image
+        psWarning("Missing statistics --- flagging entire image as suspect.");
+        return (psImage*)psBinaryOp(out, out, "+", psScalarAlloc(1.0, PS_TYPE_S32));
     }
 
@@ -146 +151 @@
       case PM_MASK_ID_VALUE:
         limit = thresh;
-        break;
+        break;
 
       case PM_MASK_ID_FRACTION:
         limit = thresh * nTotal;
-        break;
+        break;
 
       case PM_MASK_ID_SIGMA: {
@@ -162 +167 @@
         }
         limit = thresh * stats->clippedStdev;
-        psTrace ("psModules.detrend", 3, "bad: %f -> %f\n", stats->clippedStdev, limit);
+        psTrace ("psModules.detrend", 3, "bad: %f -> %f\n", stats->clippedStdev, limit);
         psFree(stats);
-        break;
+        break;
       }
 
@@ -195 +200 @@
         limit = max + 1.0 - thresh * sqrtf((float)max + 1.0);
 
-        psTrace ("psModules.detrend", 3, "bad: mode: %d, stdev: %f, limit: %f\n", max, sqrtf((float)max + 1.0), limit);
-        break;
+        psTrace ("psModules.detrend", 3, "bad: mode: %d, stdev: %f, limit: %f\n", max, sqrtf((float)max + 1.0), limit);
+        break;
       }
       default:
@@ -237 +242 @@
     if (!strcasecmp (string, "VALUE")) {
       return PM_MASK_ID_VALUE;
-    } 
+    }
     if (!strcasecmp (string, "FRACTION")) {
       return PM_MASK_ID_FRACTION;
-    } 
+    }
     if (!strcasecmp (string, "SIGMA")) {
       return PM_MASK_ID_SIGMA;
-    } 
+    }
     if (!strcasecmp (string, "POISSON")) {
       return PM_MASK_ID_POISSON;
-    } 
+    }
     return PM_MASK_ID_NONE;
 }