Changeset 10331

Timestamp:

Nov 30, 2006, 3:09:15 PM (19 years ago)

Author:

Paul Price

Message:

Add processing data to header.

File:

• trunk/ppImage/src/ppImageDetrendFringe.c (modified) (4 diffs)

trunk/ppImage/src/ppImageDetrendFringe.c

@@ -141 +141 @@
 
 
-psImage *ppImageDetrendFringeGenerate(pmCell *science, pmCell *fringes)
+bool ppImageDetrendFringeGenerate(pmCell *science, pmCell *fringes)
 {
     PS_ASSERT_PTR_NON_NULL(science, false);
@@ -159 +159 @@
         expTime = 1.0;
     }
+
+    pmHDU *hdu = pmHDUFromCell(science);// HDU  of interest
 
     // Construct the fringe image from the scale
@@ -174 +176 @@
             psBinaryOp(sumFringe, sumFringe, "+", fringe);
         }
+
+        psVector *md5 = psImageMD5(fringe); // md5 hash
+        psString md5string = psMD5toString(md5); // String
+        psFree(md5);
+        psStringPrepend(&md5string, "Fringe image %d (scale %.3f) MD5: ",
+                        i, solution->coeff->data.F32[i + 1]);
+        psMetadataAddStr(hdu->header, PS_LIST_TAIL, "HISTORY", PS_META_DUPLICATE_OK,
+                         md5string, "");
+        psFree(md5string);
     }
     if (expTime != 1.0) {
@@ -179 +190 @@
     }
 
+    // Apply the correction to all readouts
+    psArray *readouts = science->readouts; // Array of science readouts
+    for (int i = 0; i < readouts->n; i++) {
+        pmReadout *readout = readouts->data[i]; // Readout of interest
+        if (!readout || !readout->data_exists) {
+            continue;
+        }
+
+        // XXX: Make generic, so subregions may be subtracted as well
+        if (!psBinaryOp(readout->image, readout->image, "-", sumFringe)) {
+            psError(PS_ERR_UNEXPECTED_NULL, false, "Unable to subtract fringe.\n");
+            return false;
+        }
+    }
+    psFree(sumFringe);
+
     return sumFringe;
 }