Changeset 15808 for branches/eam_branch_20071212/psModules/src/detrend/pmMaskBadPixels.c

Timestamp:

Dec 13, 2007, 7:54:08 AM (19 years ago)

Author:

eugene

Message:

working on new masking api

File:

• branches/eam_branch_20071212/psModules/src/detrend/pmMaskBadPixels.c (modified) (2 diffs)

branches/eam_branch_20071212/psModules/src/detrend/pmMaskBadPixels.c

@@ -81 +81 @@
 
 psImage *pmMaskFlagSuspectPixels(psImage *out, const pmReadout *readout, float rej,
-                                 psMaskType maskVal, float frac, psRandom *rng)
+                                 psMaskType maskVal, float median, float stdev)
 {
     PS_ASSERT_PTR_NON_NULL(readout, NULL);
@@ -103 +103 @@
     psImage *mask = readout->mask;      // Corresponding mask
 
-    if (rng) {
-        psMemIncrRefCounter(rng);
-    } else {
-        rng = psRandomAlloc(PS_RANDOM_TAUS, 0);
-    }
-
-    // XXX note that this now will accept any of several stats options
-    psStats *stats = psStatsAlloc(PS_STAT_ROBUST_MEDIAN | PS_STAT_ROBUST_STDEV);
-    stats->nSubsample = frac * image->numCols * image->numRows;
-    if (!psImageBackground(stats, NULL, image, mask, maskVal, rng) ||
-            !isfinite(stats->robustMedian) || !isfinite(stats->robustUQ) || !isfinite(stats->robustLQ)) {
-        psError(PS_ERR_UNKNOWN, false, "Unable to measure image statistics.\n");
-        psFree(stats);
-        psFree(rng);
-        return NULL;
-    }
-    psFree(rng);
-
-    float median = stats->robustMedian; // Median value
-    float stdev = stats->robustStdev; // Estimate of the standard deviation
-    psFree(stats);
-
     if (!out) {
         out = psImageAlloc(image->numCols, image->numRows, PS_TYPE_S32);