Changeset 42183 for branches/eam_branches/ipp-20220316/fpcamera/src/fpcameraAnalysis.c

Timestamp:

Apr 28, 2022, 4:13:55 PM (4 years ago)

Author:

eugene

Message:

fpcamera now reads in the chip image, mask, variance; generates pmSources and loads them on the analysis MD; saves the residual image

File:

• branches/eam_branches/ipp-20220316/fpcamera/src/fpcameraAnalysis.c (modified) (3 diffs)

branches/eam_branches/ipp-20220316/fpcamera/src/fpcameraAnalysis.c

@@ -1 +1 @@
 # include "fpcamera.h"
 
-# define ESCAPE(ERROR, MSG) { psErrorStackPrint(stderr, MSG); return false; }
+# define ESCAPE(ERROR, MSG) { psErrorStackPrint(stderr, MSG); psFree (view); return false; }
 
 /* \brief this function loops over chips and performs forced photometry for the references */
@@ -7 +7 @@
 
     bool status;
+    pmChip *chip;
+    pmCell *cell;
+    pmReadout *readout;
+    pmFPAview *view = NULL;
 
     // measure the total elapsed time in fpcameraAnalysis.
@@ -15 +19 @@
     if (!recipe) ESCAPE (FPCAMERA_ERR_CONFIG, "Can't find FPCAMERA recipe");
 
-    psLogMsg("fpcamera", 3, "TIMEMARK: fpcameraMosaicAstrom: %f sec\n", psTimerMark ("complete"));
+    // loop over the input images
+    pmFPAfile *input = psMetadataLookupPtr (NULL, config->files, "FPCAMERA.INPUT");
+    if (!input) ESCAPE(FPCAMERA_ERR_CONFIG, "Can't find or interpret output file rule FPCAMERA.INPUT!");
+
+    // astrometry reference (smf)
+    pmFPAfile *astrom = psMetadataLookupPtr (NULL, config->files, "FPCAMERA.INPUT.ASTROM");
+    if (!astrom) ESCAPE(FPCAMERA_ERR_CONFIG, "Can't find or interpret output file rule FPCAMERA.INPUT.ASTROM!");
+
+    // only activate input image-type files
+    pmFPAfileActivate (config->files, false, NULL);
+    pmFPAfileActivate (config->files, true, "FPCAMERA.INPUT");
+    pmFPAfileActivate (config->files, true, "FPCAMERA.INPUT.MASK");
+    pmFPAfileActivate (config->files, true, "FPCAMERA.INPUT.VARIANCE");
+    pmFPAfileActivate (config->files, true, "FPCAMERA.RESID");
+    // Note FPCAMERA.RESID is tied to FPCAMERA.INPUT so they must be active in the same block.
+
+    view = pmFPAviewAlloc (0);
+
+    // load images at FPA level (if appropriate)
+    if (!pmFPAfileIOChecks (config, view, PM_FPA_BEFORE)) ESCAPE(FPCAMERA_ERR_DATA, "failed to load images at FPA level");
+
+    // load images at chip level (if appropriate) 
+    while ((chip = pmFPAviewNextChip (view, input->fpa, 1)) != NULL) {
+        psTrace ("fpcamera", 4, "Chip %d: %x %x\n", view->chip, chip->file_exists, chip->process);
+        if (!chip->process || !chip->file_exists) { continue; }
+        if (!pmFPAfileIOChecks (config, view, PM_FPA_BEFORE)) ESCAPE(FPCAMERA_ERR_DATA, "failed to load images at Chip level");
+
+        // loop over all cells in chip
+        while ((cell = pmFPAviewNextCell (view, input->fpa, 1)) != NULL) {
+            psLogMsg ("fpcamera", 5, "Cell %d: %x %x\n", view->cell, cell->file_exists, cell->process);
+            if (!cell->process || !cell->file_exists) { continue; }
+
+            // loop over all readouts in cell
+            while ((readout = pmFPAviewNextReadout (view, input->fpa, 1)) != NULL) {
+                psLogMsg ("fpcamera", 6, "Cell %d: %x %x\n", view->cell, cell->file_exists, cell->process);
+                if (!readout->data_exists) { continue; }
+
+                fpcameraChooseRefstars (input, astrom, view);
+            }
+        }
+        if (!pmFPAfileIOChecks (config, view, PM_FPA_AFTER)) ESCAPE(FPCAMERA_ERR_IO, "failure in IOChecks(AFTER) at Chip");
+    }
+    if (!pmFPAfileIOChecks (config, view, PM_FPA_AFTER)) ESCAPE(FPCAMERA_ERR_IO, "failure in IOChecks(AFTER) at FPA");
+
+    psFree (view);
+
+    psLogMsg("fpcamera", 3, "TIMEMARK: fpcameraAnalysis: %f sec\n", psTimerMark ("analysis"));
     return true;
 }
+
+// NOTES
+// chip->process is set (unset) based on command-line -chip selections (in fpcameraArguments)
+// chip->file_exists is set (in pmFPAFlags.c:pmChipSetFileStatus) by pmFPAfileDefineFromArgs (in pmFPAAddSource...)