Changeset 14787

Timestamp:

Sep 7, 2007, 11:00:34 AM (19 years ago)

Author:

eugene

Message:

convert to new psImageBackground API, add option to plot histograms of background cells

File:

• branches/eam_branch_20070830/psphot/src/psphotImageMedian.c (modified) (6 diffs)

branches/eam_branch_20070830/psphot/src/psphotImageMedian.c

@@ -80 +80 @@
                                  PS_STAT_FITTED_MEAN |
                                  PS_STAT_FITTED_MEAN_V2 |
-                                 PS_STAT_FITTED_MEAN_V3))) {
+                                 PS_STAT_FITTED_MEAN_V3 |
+                                 PS_STAT_FITTED_MEAN_V4))) {
         statsOptionLocation = PS_STAT_FITTED_MEAN;
     }
@@ -101 +102 @@
     } else if (statsOptionLocation & PS_STAT_FITTED_MEAN_V3) {
         statsOptionWidth = PS_STAT_FITTED_STDEV_V3;
+    } else if (statsOptionLocation & PS_STAT_FITTED_MEAN_V4) {
+        statsOptionWidth = PS_STAT_FITTED_STDEV_V4;
     } else {
         psAbort("Unable to estimate variance of selected statsOptionLocations 0x%x", statsOptionLocation);
@@ -196 +199 @@
             // XXX don't bother trying if there are no valid pixels...
 
-            if (psImageBackground(stats, subset, submask, maskVal, rng)) {
+            psVector *sample = NULL;
+
+            if (psImageBackground(stats, &sample, subset, submask, maskVal, rng)) {
                 if (stats->options & PS_STAT_ROBUST_QUARTILE) {
                     modelData[iy][ix] = stats->robustMedian;
@@ -203 +208 @@
                 }
                 modelStdevData[iy][ix] = psStatsGetValue(stats, statsOptionWidth);
+
+                // supply sample to plotting routing
+                psphotDiagnosticPlots (config, "IMAGE.BACKGROUND.CELL.HISTOGRAM", ix, iy, modelData[iy][ix], modelStdevData[iy][ix], sample);
+                psFree (sample);
             } else {
                 psStatsOptions currentOptions = stats->options;
                 stats->options = PS_STAT_ROBUST_MEDIAN | PS_STAT_ROBUST_STDEV;
-                if (!psImageBackground(stats, subset, submask, maskVal, rng)) {
+                if (!psImageBackground(stats, &sample, subset, submask, maskVal, rng)) {
                     psLogMsg ("psphot", PS_LOG_WARN, "Failed to estimate background using ROBUST_MEDIAN for "
                                "(%dx%d, (row0,col0) = (%d,%d)",
@@ -214 +223 @@
                     modelData[iy][ix] = psStatsGetValue (stats, PS_STAT_ROBUST_MEDIAN);
                     modelStdevData[iy][ix] = psStatsGetValue(stats, PS_STAT_ROBUST_STDEV);
+
+                    // supply sample to plotting routing
+                    psphotDiagnosticPlots (config, "IMAGE.BACKGROUND.CELL.HISTOGRAM", ix, iy, modelData[iy][ix], modelStdevData[iy][ix], sample);
+                    psFree (sample);
                 }
                 // drop errors caused by psImageBackground failures
@@ -224 +237 @@
             psFree (submask);
         }
+    }
+
+    if (psTraceGetLevel("psphot") > 5) {
+        char name[256];
+        sprintf (name, "backraw.%02d.fits", npass);
+        psphotSaveImage (NULL, model->image, name);
     }