Changeset 6862 for trunk/psphot/src/psphotParseCamera.c

Timestamp:

Apr 14, 2006, 11:58:44 AM (20 years ago)

Author:

eugene

Message:

renamed config to recipe in lower-level functions, added pixel weights to pmSource, other cleanups

File:

• trunk/psphot/src/psphotParseCamera.c (modified) (1 diff)

trunk/psphot/src/psphotParseCamera.c

@@ -1 +1 @@
 # include "psphot.h"
 
+// define the needed / desired I/O files
 bool psphotParseCamera (pmConfig *config) {
 
-    bool status;
+    bool status = false;
 
-    status = false;
     pmFPAfile *input = pmFPAfileFromArgs (&status, config, "PSPHOT.INPUT", "INPUT");
-    if (!status) {
-        psAbort (__func__, "missing INPUT entry");
-    }
+    if (!status) { psAbort (__func__, "missing INPUT entry"); }
 
-    // recipe override values (command-line options):
+    // select recipe options supplied on command line
+    psMetadata *recipe  = psMetadataLookupPtr (&status, config->recipes,   "PSPHOT");
     psMetadata *options = psMetadataLookupPtr (&status, config->arguments, "PSPHOT.OPTIONS");
-    psMetadata *recipe = psMetadataLookupPtr (&status, config->recipes, "PSPHOT");
-    psMetadataIterator *iter = psMetadataIteratorAlloc (options, PS_LIST_HEAD, NULL);
-    psMetadataItem *item = NULL;
-    while ((item = psMetadataGetAndIncrement (iter)) != NULL) {
-        psMetadataAddItem (recipe, item, PS_LIST_TAIL, PS_META_REPLACE);
-    }
-    psFree (iter);
+    psMetadataCopy (recipe, options);
+
+    // set default recipe values here
+    psMetadataAddStr (recipe, PS_LIST_TAIL, "FITMODE",         PS_META_NO_REPLACE, "default fitting mode", "NONE");
+    psMetadataAddStr (recipe, PS_LIST_TAIL, "PHOTCODE",        PS_META_NO_REPLACE, "default photcode",     "NONE");
+    psMetadataAddStr (recipe, PS_LIST_TAIL, "BREAK_POINT",     PS_META_NO_REPLACE, "default break point",  "NONE");
+    psMetadataAddF32 (recipe, PS_LIST_TAIL, "ZERO_PT",         PS_META_NO_REPLACE, "default zero point",    25.00);
+    psMetadataAddS32 (recipe, PS_LIST_TAIL, "BACKGROUND.XBIN", PS_META_NO_REPLACE, "default binning",          64);
+    psMetadataAddS32 (recipe, PS_LIST_TAIL, "BACKGROUND.YBIN", PS_META_NO_REPLACE, "default binning",          64);
 
     // these calls bind the I/O handle to the specified fpa
     pmFPAfileDefine (config->files, config->camera, input->fpa, "PSPHOT.OUTPUT");
-    pmFPAfileDefine (config->files, config->camera, input->fpa, "PSPHOT.RESID");
 
-    // pmFPAfileDefine (config->files, format, input, "PSPHOT.PSF_INPUT");
-    // pmFPAfileDefine (config->files, format, input, "PSPHOT.PSF_OUTPUT");
+    // optionally save the residual image 
+    if (psMetadataLookupBool(NULL, recipe, "SAVE.RESID")) {
+        pmFPAfileDefine (config->files, config->camera, input->fpa, "PSPHOT.RESID");
+    }
+
+    int DX = psMetadataLookupS32 (&status, recipe, "BACKGROUND.XBIN");
+    int DY = psMetadataLookupS32 (&status, recipe, "BACKGROUND.YBIN");
 
     // these calls construct a new fpa for the I/O handle 
-    int DX = psMetadataLookupS32 (&status, recipe, "BACKGROUND_XBIN");
-    if (!status) {DX = 64;}
-    int DY = psMetadataLookupS32 (&status, recipe, "BACKGROUND_YBIN");
-    if (!status) {DY = 64;}
-    pmFPAfileFromFPA (config, input->fpa, DX, DY, "PSPHOT.BACKMDL");
-    pmFPAfileFromFPA (config, input->fpa,  1,  1, "PSPHOT.BACKGND");
-    pmFPAfileFromFPA (config, input->fpa,  1,  1, "PSPHOT.BACKSUB");
+
+    // optionally save the background model (small FITS image)
+    if (psMetadataLookupBool(NULL, recipe, "SAVE.BACKMDL")) {
+        pmFPAfileFromFPA (config, input->fpa, DX, DY, "PSPHOT.BACKMDL");
+    }
+    // optionally save the full background image
+    if (psMetadataLookupBool(NULL, recipe, "SAVE.BACKGND")) {
+        pmFPAfileFromFPA (config, input->fpa,  1,  1, "PSPHOT.BACKGND");
+    }
+    // optionally save the background-subtracted image
+    if (psMetadataLookupBool(NULL, recipe, "SAVE.BACKSUB")) {
+        pmFPAfileFromFPA (config, input->fpa,  1,  1, "PSPHOT.BACKSUB");
+    }
+    // optionally save the PSF Model
+    if (psMetadataLookupBool(NULL, recipe, "SAVE.PSF")) {
+        pmFPAfileDefine (config->files, config->camera, input->fpa, "PSPHOT.PSF.SAVE");
+    }
+    // optionally load the PSF Model
+    if (psMetadataLookupBool(NULL, recipe, "LOAD.PSF")) {
+        pmFPAfileDefine (config->files, config->camera, input->fpa, "PSPHOT.PSF.LOAD");
+    }
 
     // pmFPAfileConstruct (config->files, format, config->camera, "PSPHOT.PSF_SAMPLE");
 
-    // psphot is supplied with the output name
-    // this needs to be better: supply to all output (WRITE) files?
+    // supply the output name (from cmd-line) to all output (WRITE) files
     char *output = psMetadataLookupPtr(&status, config->arguments, "OUTPUT");
     pmFPAfileAddFileNames (config->files, "OUTPUT", output, PM_FPA_MODE_WRITE);
 
-    // set default recipe values here
-    psMetadataAddStr (recipe, PS_LIST_TAIL, "FITMODE", PS_META_NO_REPLACE, "default fitting mode", "NONE");
-    psMetadataAddStr (recipe, PS_LIST_TAIL, "PHOTCODE", PS_META_NO_REPLACE, "default photcode", "NONE");
-    psMetadataAddStr (recipe, PS_LIST_TAIL, "BREAK_POINT", PS_META_NO_REPLACE, "default break point", "NONE");
-    psMetadataAddF32 (recipe, PS_LIST_TAIL, "ZERO_PT", PS_META_NO_REPLACE, "default zero point", 25.00);
-
-    // Chip selection: if we are using a single chip, select it for each FPA
-    // Chip numbers to work on: -chip 1,2,3,5,8,10
-    # if (0)
-    char *chips = psMetadataLookupStr(NULL, config->arguments, "CHIP_SELECTIONS"); 
-    if (chips != NULL) {
-        char *p = chips;
-        while (strlen(p)) {
-            char *q = strchr (p, ',');
-            if (q == NULL) { 
-                q = p + strlen(p);
-            } else {
-                *q = 0;
-                q = q + 1;
-            }
-            int chipNum = atoi(p);
-            // XXX EAM : extend this to allow an array of selected chips (by name)
-            if (! pmFPASelectChip(input->fpa, chipNum)) {
+    // Chip selection: turn on only the chips specified
+    char *chipLine = psMetadataLookupStr(NULL, config->arguments, "CHIP_SELECTIONS"); 
+    psArray *chips = psStringSplitArray (chipLine, ",");
+    if (chips->n > 0) {
+        pmFPASelectChip (input->fpa, -1, true); // deselect all chips
+        for (int i = 0; i < chips->n; i++) {
+            int chipNum = atoi(chips->data[i]);
+            if (! pmFPASelectChip(input->fpa, chipNum, false)) {
                 psErrorStackPrint(stderr, "Chip number %d doesn't exist in camera.\n", chipNum);
                 exit(EXIT_FAILURE);
             }
-            psLogMsg("psphot", PS_LOG_INFO, "Operating only on chip %d\n", chipNum);
-            p = q;
         }
     }
-    # endif
 
     psTrace(__func__, 1, "Done with psphotParseCamera...\n");