Changeset 13900 for trunk/psphot/src/psphotImageMedian.c

Timestamp:

Jun 19, 2007, 4:40:44 PM (19 years ago)

Author:

Paul Price

Message:

Extensive changes to APIs to allow use of a nominated value to mask
against (the maskVal). Previously, the mask values were either
hard-coded (e.g., PM_MASK_SAT) or taken as anything non-zero. The
code is tested with psModules (which has similar changes) and does not
crash, but neither is it successful in marking all bad pixels. For
this reason, I have left the "gutter" pixels (cell gaps) set to 0
instead of NAN in pmFPAMosaic.

File:

• trunk/psphot/src/psphotImageMedian.c (modified) (3 diffs)

trunk/psphot/src/psphotImageMedian.c

@@ -40 +40 @@
 // generate the median in NxN boxes, clipping heavily
 // linear interpolation to generate full-scale model
-bool psphotImageMedian (pmConfig *config, const pmFPAview *view)
+bool psphotImageMedian (pmConfig *config, const pmFPAview *view, psMaskType maskVal)
 {
     bool status = true;
@@ -195 +195 @@
             // XXX don't bother trying if there are no valid pixels...
 
-            if (psImageBackground(stats, subset, submask, 0xff, rng)) {
+            if (psImageBackground(stats, subset, submask, maskVal, rng)) {
                 if (stats->options & PS_STAT_ROBUST_QUARTILE) {
                     modelData[iy][ix] = stats->robustMedian;
@@ -205 +205 @@
                 psStatsOptions currentOptions = stats->options;
                 stats->options = PS_STAT_ROBUST_MEDIAN | PS_STAT_ROBUST_STDEV;
-                if (!psImageBackground(stats, subset, submask, 0xff, rng)) {
+                if (!psImageBackground(stats, subset, submask, maskVal, rng)) {
                     psLogMsg ("psphot", PS_LOG_WARN, "Failed to estimate background using ROBUST_MEDIAN for "
                                "(%dx%d, (row0,col0) = (%d,%d)",