Changeset 15807 for branches/eam_branch_20071212/ppMerge/src/ppMergeMaskByImageStats.c

Timestamp:

Dec 13, 2007, 7:53:28 AM (19 years ago)

Author:

eugene

Message:

updates to get new mask api working

File:

• branches/eam_branch_20071212/ppMerge/src/ppMergeMaskByImageStats.c (modified) (2 diffs)

branches/eam_branch_20071212/ppMerge/src/ppMergeMaskByImageStats.c

@@ -15 +15 @@
 #include "ppMergeMask.h"
 
+// this function measures the statistics of the pixel values for a single readout (median and stdev)
+// and identifies pixels which deviate from the median by more than MASK.SUSPECT * stdev.
 
 // Generate a mask image consisting of pixels which are outliers from the image mean
@@ -65 +67 @@
                             view->chip, view->cell, view->readout);
                     float frac = options->sample / (float)(roIn->image->numCols * roIn->image->numRows);
+
+                    // measure the global image stats: median and stdev
+                    // 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);
+
                     suspect = pmMaskFlagSuspectPixels(suspect, roIn, options->maskSuspect,
-                                                      options->combine->maskVal, frac, rng);
+                                                      options->combine->maskVal, median, stdev);
+                    
+                    psMetadataAddF32(roIn->analysis, PS_LIST_TAIL, "MEDIAN", 0, "median", median);
+                    psMetadataAddF32(roIn->analysis, PS_LIST_TAIL, "STDEV", 0, "stdev", stdev);
                 }