Changeset 32704 for trunk/psphot/src/psphotStackImageLoop.c

Timestamp:

Nov 18, 2011, 11:14:36 AM (14 years ago)

Author:

eugene

Message:

make sure the psphot output headers are copied to the output container

File:

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

trunk/psphot/src/psphotStackImageLoop.c

@@ -10 +10 @@
 bool UpdateHeadersForFPA (pmConfig *config, pmFPAview *view);
 bool UpdateHeadersForChip (pmConfig *config, pmFPAview *view);
+bool UpdateHeadersForReadout (pmConfig *config, pmFPAview *view);
 
 bool psphotStackImageLoop (pmConfig *config) {
@@ -66 +67 @@
                     return false;
                 }
+
+                UpdateHeadersForReadout(config, view);
+
                 psMemDump("psphot");
             }
@@ -101 +105 @@
 }
 
+bool UpdateHeadersForReadout (pmConfig *config, pmFPAview *view) {
+
+    psMetadata *recipe = psMetadataLookupPtr (NULL, config->recipes, PSPHOT_RECIPE);
+    bool useRaw = psMetadataLookupBool (NULL, recipe, "PSPHOT.STACK.USE.RAW");
+
+    int num = psphotFileruleCount(config, "PSPHOT.INPUT");
+
+    // loop over the available readouts
+    for (int i = 0; i < num; i++) {
+
+        // find the currently selected readout
+        pmFPAfile *output = pmFPAfileSelectSingle(config->files, "PSPHOT.STACK.OUTPUT.IMAGE", i); // File of interest
+        psAssert (output, "missing file?");
+
+        pmFPAfile *inputRaw = pmFPAfileSelectSingle(config->files, "PSPHOT.STACK.INPUT.RAW", i); // File of interest
+        pmFPAfile *inputCnv = pmFPAfileSelectSingle(config->files, "PSPHOT.STACK.INPUT.CNV", i); // File of interest
+        pmFPAfile *input = useRaw ? inputRaw : inputCnv;
+        psAssert (input, "missing input file");
+
+        // just copy the input headers to the output headers, then update version info
+        pmReadout *inReadout = pmFPAviewThisReadout(view, input->fpa); ///< Chip in the input
+        pmReadout *outReadout = pmFPAviewThisReadout(view, output->fpa); ///< Chip in the output
+
+        psMetadata *header = psMetadataLookupPtr (NULL, inReadout->analysis, "PSPHOT.HEADER");
+        psMetadataAdd (outReadout->analysis, PS_LIST_TAIL, "PSPHOT.HEADER",  PS_DATA_METADATA, "header stats", header);
+    }
+    return true;
+}
+
 bool UpdateHeadersForChip (pmConfig *config, pmFPAview *view) {