Changeset 9907 for trunk/ppMerge/src/ppMergeCombine.c

Timestamp:

Nov 7, 2006, 5:34:39 PM (20 years ago)

Author:

Paul Price

Message:

Pulling out function that got implemented in psLib

File:

• trunk/ppMerge/src/ppMergeCombine.c (modified) (6 diffs)

trunk/ppMerge/src/ppMergeCombine.c

@@ -78 +78 @@
 #endif
 
+// Write fringe extensions
+static void writeFringes(psMetadata *fringes, // Table of extensions with fringe statistics to write
+                         psFits *fits   // FITS file to which to write
+                         )
+{
+    assert(fringes);
+
+    // Write the fringe extensions
+    psMetadataIterator *fringesIter = psMetadataIteratorAlloc(fringes, PS_LIST_HEAD, false); // Iterator
+    psMetadataItem *fringeItem;         // Item from iteration
+    while ((fringeItem = psMetadataGetAndIncrement(fringesIter))) {
+        const char *extname = fringeItem->name; // Extension name
+        psArray *fringe = fringeItem->data.V; // Fringe data (array of fringe components)
+        if (!pmFringesWriteFits(fits, NULL, fringe, extname)) {
+            psWarning("Unable to write fringe data to extension %s\n", extname);
+        }
+    }
+
+    // Clear the list
+    while (psListLength(fringes->list) > 0) {
+        psMetadataRemoveIndex(fringes, PS_LIST_TAIL);
+    }
+
+    return;
+}
+
+
 // Combine the inputs
 bool ppMergeCombine(psImage *scales,    // Scales for each cell of each integration, or NULL
@@ -101 +128 @@
     // Iterate over the FPA
     pmFPA *fpa = data->out;             // Output FPA
+    psMetadata *fringeExtns = psMetadataAlloc(); // List of fringe extensions to write
     pmFPAview *view = pmFPAviewAlloc(0);// View of FPA, for iteration
     int cellNum = -1;                   // Cell number in the whole FPA
@@ -296 +324 @@
                 psStringAppend(&extname, "FRINGE_%s_%s", chipName, cellName);
 
-                if (!pmFringesWriteFits(data->outFile, NULL, fringes, extname)) {
-                    psWarning("Unable to write fringe data to extension %s\n", extname);
-                }
-
+                psMetadataAdd(fringeExtns, PS_LIST_TAIL, extname, PS_DATA_UNKNOWN,
+                              "Fringe data to write", fringes);
                 psFree(fringes);
                 psFree(extname);
@@ -329 +355 @@
                 pmCellWrite(cell, data->outFile, config->database, false);
                 pmCellFreeData(cell);
+
+                if (options->fringe) {
+                    writeFringes(fringeExtns, data->outFile);
+                }
             }
         }
@@ -344 +374 @@
             pmChipWrite(chip, data->outFile, config->database, false, false);
             pmChipFreeData(chip);
+
+            if (options->fringe) {
+                writeFringes(fringeExtns, data->outFile);
+            }
         }
     }
@@ -357 +391 @@
         pmFPA *fpaIn = data->in->data[i]; // Input FPA
         pmFPAFreeData(fpaIn);
+
+        if (options->fringe) {
+            writeFringes(fringeExtns, data->outFile);
+        }
     }
     pmFPAFreeData(data->out);