Changeset 30811 for branches/eam_branches/ipp-20110213/ppStack/src/ppStackUpdateHeader.c

Timestamp:

Mar 6, 2011, 12:39:13 PM (15 years ago)

Author:

eugene

Message:

update output headers to include more info on the input images (airmasses, exptimes, zpts)

File:

• branches/eam_branches/ipp-20110213/ppStack/src/ppStackUpdateHeader.c (modified) (2 diffs)

branches/eam_branches/ipp-20110213/ppStack/src/ppStackUpdateHeader.c

@@ -70 +70 @@
     psMetadataAddStr(hdu->header, PS_LIST_TAIL, "TESS_ID", PS_META_REPLACE, "type of stack", tessID);
 
-    // add data about the inputs
+    psMetadataAddF32(hdu->header, PS_LIST_TAIL, "AIRM_SLP", PS_META_REPLACE, "airmass slope", options->airmassSlope);
+    
+
+    // write the following fields to the output headers:
+    // INP_NNNN : input file names
+    // SCL_NNNN : scale factor for each input file
+    // ZPT_NNNN : zero point of input file 
+    // EXT_NNNN : exptime of input file
+    // AIR_NNNN : airmass of input file
+
+    // use separate loops so they are blocked together in the headers (more readable)
+
+    // save the input filenames
     for (int i = 0; i < options->num; i++) {
         char field[64];
@@ -78 +90 @@
         psMetadataAddStr(hdu->header, PS_LIST_TAIL, field, PS_META_REPLACE, "input image name", basename);
         psFree(basename);
-
-        float inScale = pow(10.0, -0.4*options->norm->data.F32[i]);
-        snprintf (field, 64, "SCL_%04d", i);
-        psMetadataAddF32(hdu->header, PS_LIST_TAIL, field, PS_META_REPLACE, "input image scale factor", inScale);
     }   
 
+    // save the input normalizations
+    for (int i = 0; i < options->num; i++) {
+        char field[64];
+
+        float value = options->norm ? pow(10.0, -0.4*options->norm->data.F32[i]) : NAN;
+        snprintf (field, 64, "SCL_%04d", i);
+        psMetadataAddF32(hdu->header, PS_LIST_TAIL, field, PS_META_REPLACE, "input image scale factor", value);
+    }   
+    // save the input exptimes
+    for (int i = 0; i < options->num; i++) {
+        char field[64];
+
+        float value = options->zpInput ? options->zpInput->data.F32[i] : NAN;
+        snprintf (field, 64, "ZPT_%04d", i);
+        psMetadataAddF32(hdu->header, PS_LIST_TAIL, field, PS_META_REPLACE, "input image zero point", value);
+    }   
+    // save the input normalizations
+    for (int i = 0; i < options->num; i++) {
+        char field[64];
+
+        float value = options->expTimeInput ? options->expTimeInput->data.F32[i] : NAN;
+        snprintf (field, 64, "EXP_%04d", i);
+        psMetadataAddF32(hdu->header, PS_LIST_TAIL, field, PS_META_REPLACE, "input image exptime", value);
+    }   
+    // save the input normalizations
+    for (int i = 0; i < options->num; i++) {
+        char field[64];
+
+        float value = options->airmassInput ? options->airmassInput->data.F32[i] : NAN;
+        snprintf (field, 64, "AIR_%04d", i);
+        psMetadataAddF32(hdu->header, PS_LIST_TAIL, field, PS_META_REPLACE, "input image airmass", value);
+    }   
     return true;
 }