Changeset 24905

Timestamp:

Jul 22, 2009, 9:52:11 PM (17 years ago)

Author:

Paul Price

Message:

If statistics can't be calculated, kill entire readout.

File:

• trunk/psModules/src/detrend/pmPattern.c (modified) (1 diff)

trunk/psModules/src/detrend/pmPattern.c

@@ -48 +48 @@
     psRandom *rng = psRandomAlloc(PS_RANDOM_TAUS); // Random number generator
     if (!psImageBackground(stats, NULL, ro->image, ro->mask, maskVal, rng)) {
-        psWarning("Unable to calculate statistics on readout.");
+        psWarning("Unable to calculate statistics on readout; masking entire readout.");
+        psErrorClear();
         psFree(stats);
         psFree(rng);
-        return false;
+        psImageInit(image, NAN);
+        if (mask) {
+            psBinaryOp(mask, mask, "|", psScalarAlloc(maskBad, PS_TYPE_IMAGE_MASK));
+        }
+        if (ro->variance) {
+            psImageInit(image, NAN);
+        }
+        return true;
     }
     float lower = stats->robustMedian - thresh * stats->robustStdev; // Lower bound for data