Changeset 15865 for trunk/psModules/src/detrend/pmMaskBadPixels.c

Timestamp:

Dec 16, 2007, 12:22:06 PM (18 years ago)

Author:

eugene

Message:

added pmMaskIdentifyMode functions

File:

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

trunk/psModules/src/detrend/pmMaskBadPixels.c

@@ -5 +5 @@
 #include <stdio.h>
 #include <math.h>
+#include <strings.h>
 #include <pslib.h>
 
@@ -125 +126 @@
     psFree(stats);
 
+    psTrace ("psModules.detrend", 3, "suspect: %f +/- %f\n", median, stdev);
+
     if (!out) {
         out = psImageAlloc(image->numCols, image->numRows, PS_TYPE_S32);
@@ -142 +145 @@
 }
 
-
-psImage *pmMaskIdentifyBadPixels(const psImage *suspects, float thresh, psMaskType maskVal)
+psImage *pmMaskIdentifyBadPixels(const psImage *suspects, psMaskType maskVal, int nTotal, float thresh, pmMaskIdentifyMode mode)
 {
     PS_ASSERT_IMAGE_NON_NULL(suspects, NULL);
@@ -150 +152 @@
 
     float limit = NAN;                  // Limit for masking
-    if (thresh > 0) {
+
+    switch (mode) {
+      case PM_MASK_ID_VALUE:
+        limit = thresh;
+        break;
+
+      case PM_MASK_ID_FRACTION:
+        limit = thresh * nTotal;
+        break;
+
+      case PM_MASK_ID_SIGMA: {
         psStats *stats = psStatsAlloc(PS_STAT_CLIPPED_STDEV); // Statistics
         stats->clipSigma = 5.0;
@@ -160 +172 @@
         }
         limit = thresh * stats->clippedStdev;
-        psFree(stats);
-    } else if (thresh < 0) {
+        psTrace ("psModules.detrend", 3, "bad: %f -> %f\n", stats->clippedStdev, limit);
+        psFree(stats);
+        break;
+      }
+
+      case PM_MASK_ID_POISSON: {
         psStats *stats = psStatsAlloc(PS_STAT_MAX); // Statistics
         if (!psImageStats(stats, suspects, NULL, 0)) {
@@ -189 +205 @@
         limit = max + 1.0 - thresh * sqrtf((float)max + 1.0);
 
-    } else {
-        psError(PS_ERR_BAD_PARAMETER_VALUE, true,
-                "Threshold (%f) must be a positive or negative real number.\n", thresh);
+        psTrace ("psModules.detrend", 3, "bad: mode: %d, stdev: %f, limit: %f\n", max, sqrtf((float)max + 1.0), limit);
+        break;
+      }
+      default:
+        psError(PS_ERR_BAD_PARAMETER_VALUE, true, "Invalid mask identify mode");
         return NULL;
     }
@@ -212 +230 @@
     }
 
+    psTrace ("psModules.detrend", 3, "bad pixel threshold: %f", limit);
+
     for (int y = 0; y < suspects->numRows; y++) {
         for (int x = 0; x < suspects->numCols; x++) {
@@ -222 +242 @@
     return badpix;
 }
+
+pmMaskIdentifyMode pmMaskIdentifyModeFromString (const char *string) {
+
+    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;
+}