Changeset 12436

Timestamp:

Mar 13, 2007, 4:42:49 PM (19 years ago)

Author:

eugene

Message:

adding SKY_CLIP_SIGMA option, using statsDefaults to reset stats options

File:

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

trunk/psphot/src/psphotImageMedian.c

@@ -62 +62 @@
         stats->options |= PS_STAT_USE_BINSIZE;
     }
+
+    psMetadataConfigPrint (stderr, recipe);
+
+    // optionally set the sigma clipping
+    stats->clipSigma = psMetadataLookupF32 (&status, recipe, "SKY_CLIP_SIGMA");
+    if (!status) {
+        stats->clipSigma = 1.0;
+    }
+
+    // stats is not initialized by psStats???  use this to save the input options
+    psStats *statsDefaults = psStatsAlloc (statsOption);
+    *statsDefaults = *stats;
 
     // find the currently selected readout
@@ -142 +154 @@
             int sx = ix*DX - xOffset;
             int sy = iy*DY - yOffset;
+
+            // XXX test override
+            // sx = 1955;
+            // sy = 4;
+
             region = psRegionSet (sx, sx + 2*DX, sy, sy + 2*DY);
             region = psRegionForImage (image, region);
@@ -152 +169 @@
             psImage *submask = psImageSubset (mask, region);
 
+            // reset the default values
+            *stats = *statsDefaults;
+
             // Use the selected background statistic for the first pass
             // If it fails, fall back on the "ROBUST_MEDIAN" version
             // If both fail, set the pixel to NAN and (below) interpolate
             // XXX psImageBackground will probably be renamed psImageStats
+            // XXX don't bother trying if there are no valid pixels...
             if (psImageBackground(stats, subset, submask, 0xff, rng)) {
                 if (stats->options & PS_STAT_ROBUST_QUARTILE) {
@@ -164 +185 @@
             } else {
                 psStatsOptions currentOptions = stats->options;
-                stats->options |= PS_STAT_ROBUST_MEDIAN;
+                stats->options = PS_STAT_ROBUST_MEDIAN;
                 if (!psImageBackground(stats, subset, submask, 0xff, rng)) {
                     psLogMsg ("psphot", PS_LOG_WARN, "Failed to estimate background using ROBUST_MEDIAN for "
@@ -299 +320 @@
 # ifdef TESTSAVE
     psphotSaveImage (NULL, image, "backsub.fits");
+    exit (1);
 # endif
 
     psLogMsg ("psphot", PS_LOG_INFO, "subtracted background model: %f sec\n", psTimerMark ("psphot"));
     psFree(stats);
+    psFree(statsDefaults);
     psFree(rng);