Changeset 42960

Timestamp:

Feb 2, 2026, 2:09:19 PM (3 months ago)

Author:

tdeboer

Message:

merging ppSub changes into ops tag

Location:

tags/ipp-ops-20220906/ppSub/src

Files:

• ppSubArguments.c (modified) (3 diffs, 1 prop)
• ppSubDefineOutput.c (modified) (2 diffs, 1 prop)

tags/ipp-ops-20220906/ppSub/src/ppSubArguments.c

@@ -104 +104 @@
     psMetadataAddBool(arguments, PS_LIST_TAIL, "-visual", 0, "Show diagnostic plots", NULL);
     psMetadataAddBool(arguments, PS_LIST_TAIL, "-updatemode", 0, "update mode?", false);
+    psMetadataAddStr(arguments, PS_LIST_TAIL, "-skycell_id", 0, "skycell id", NULL);
 
     if (argc == 1 || !psArgumentParse(arguments, &argc, argv) || argc != 2) {
@@ -110 +111 @@
 
     psMetadataAddStr(config->arguments, PS_LIST_TAIL, "OUTPUT", 0, "Name of the output image", argv[1]);
-
 
     const char *inImage = psMetadataLookupStr(NULL, arguments, "-inimage"); // Name of input image
@@ -151 +151 @@
     }
 
+    char *skycell_id = psMetadataLookupStr(NULL, arguments, "-skycell_id"); // Name of skycell_id
+    psMetadataAddStr(config->arguments, PS_LIST_TAIL, "SKYCELL_ID", 0, "skycell id", skycell_id);
+
     data->stamps = psMemIncrRefCounter(psMetadataLookupStr(NULL, arguments, "-stamps"));
 

tags/ipp-ops-20220906/ppSub/src/ppSubDefineOutput.c

@@ -30 +30 @@
 
     bool status = false;
+    float zptObs_in=0.;
+    float zptErr_in=0.;
+    float zptObs_ref=0.;
+    float zptErr_ref=0.;
 
     pmFPAview *view = ppSubViewReadout(); // View to readout
@@ -139 +143 @@
     psFree(inBase);
 
+    //read the zeropoint keywords
+    pmCell *cell_in = pmFPAfileThisCell(config->files, view, "PPSUB.INPUT");
+    pmHDU *hdu_in = pmHDUFromCell(cell_in);
+    zptObs_in = psMetadataLookupF32 (NULL, hdu_in->header, "ZPT_OBS");
+    zptErr_in = psMetadataLookupF32 (NULL, hdu_in->header, "ZPT_ERR");
+
+    pmCell *cell_ref = pmFPAfileThisCell(config->files, view, "PPSUB.REF");
+    pmHDU *hdu_ref = pmHDUFromCell(cell_ref);
+    zptObs_ref = psMetadataLookupF32 (NULL, hdu_ref->header, "ZPT_OBS");
+    zptErr_ref = psMetadataLookupF32 (NULL, hdu_ref->header, "ZPT_ERR");
+
+    if (zptObs_in != 0.) { 
+        psMetadataAddF32 (outHDU->header, PS_LIST_TAIL, "ZPT_OBS_IN", PS_META_REPLACE, "measured zero point",  zptObs_in);
+        psMetadataAddF32 (outHDU->header, PS_LIST_TAIL, "ZPT_ERR_IN", PS_META_REPLACE, "error on zero point",  zptErr_in);
+    }
+
+    if (zptObs_ref != 0.) { 
+        psMetadataAddF32 (outHDU->header, PS_LIST_TAIL, "ZPT_OBS_REF", PS_META_REPLACE, "measured zero point for template image",  zptObs_ref);
+        psMetadataAddF32 (outHDU->header, PS_LIST_TAIL, "ZPT_ERR_REF", PS_META_REPLACE, "error on zero point for template image",  zptErr_ref);
+    }
+
+    //also output the skycell id
+    char *skycell_id = psMetadataLookupStr(NULL, config->arguments, "SKYCELL_ID");
+
+    psMetadataAddStr (outHDU->header, PS_LIST_TAIL, "SKYCELL_ID", PS_META_REPLACE, "skycell_id",  skycell_id);
+
     ppSubVersionHeader(outHDU->header);