Changeset 26523 for branches/eam_branches/20091201/psphot/src/psphotModelBackground.c

Timestamp:

Jan 6, 2010, 11:41:09 AM (17 years ago)

Author:

eugene

Message:

starting work on psphotStack

File:

• branches/eam_branches/20091201/psphot/src/psphotModelBackground.c (modified) (2 diffs)

branches/eam_branches/20091201/psphot/src/psphotModelBackground.c

@@ -328 +328 @@
 }
 
-
+// XXX these two functions are absurdly-named and very similar -- fix!
 
 psImage *psphotBackgroundModel(pmReadout *ro, const pmConfig *config)
@@ -376 +376 @@
     return true;
 }
+
+// XXX supply filename or keep PSPHOT.INPUT fixed?
+bool psphotStackModelBackground (pmConfig *config, const pmFPAview *view)
+{
+    bool status = true;
+
+    int num = psMetadataAddS32 (&status, config->arguments, "INPUTS.NUM");
+    psAbort (!status, "programming error: must define INPUTS.NUM");
+
+    // loop over the available readouts
+    for (int i = 0; i < num; i++) {
+
+        pmFPAfile *file = pmFPAfileSelectSingle(config->files, "PSPHOT.INPUT", i); // File of interest
+
+        // find the currently selected readout
+        pmReadout  *readout = pmFPAviewThisReadout (view, file->fpa);
+        PS_ASSERT_PTR_NON_NULL (readout, false);
+
+        psImageBinning *binning = psphotBackgroundBinning(readout->image, config); // Image binning parameters
+        pmReadout *model = pmFPAGenerateReadout(config, view, "PSPHOT.BACKMDL", inFPA, binning); // this needs to be a MULTI like PSPHOT.INPUT
+        pmReadout *modelStdev = pmFPAGenerateReadout(config, view, "PSPHOT.BACKMDL.STDEV", inFPA, binning); // this needs to be a MULTI like PSPHOT.INPUT
+
+        // XXX save results on model->analysis or readout->analysis??
+        if (!backgroundModel(model->image, modelStdev->image, model->analysis, readout, binning, config)) {
+            psError(PS_ERR_UNKNOWN, false, "Unable to generate background model");
+            return false;
+        }
+        npass ++;
+    }
+    return true;
+}