Changeset 13512 for trunk/ppStack/src/ppStackArguments.c

Timestamp:

May 24, 2007, 11:48:31 AM (19 years ago)

Author:

Paul Price

Message:

Use seeing information provided to perform additional rejection

File:

• trunk/ppStack/src/ppStackArguments.c (modified) (2 diffs)

trunk/ppStack/src/ppStackArguments.c

@@ -85 +85 @@
     psMetadataAddU8(arguments,  PS_LIST_TAIL, "-mask-bad", 0, "Mask value for bad pixels", 0);
     psMetadataAddU8(arguments,  PS_LIST_TAIL, "-mask-blank", 0, "Mask value for blank region", 0);
+    psMetadataAddStr(arguments, PS_LIST_TAIL, "-seeing", 0, "Name of file with seeing FWHMs", NULL);
 
     if (argc == 1 || !psArgumentParse(arguments, &argc, argv) || argc != 4) {
@@ -112 +113 @@
     VALUE_ARG_RECIPE_INT("-mask-blank",     "MASK.BLANK",   U8, 0);
 
+    const char *seeingName = psMetadataLookupStr(NULL, arguments, "-seeing"); // Name of seeing file
+    if (seeingName && strlen(seeingName) > 0) {
+        psArray *seeing = psVectorsReadFromFile(seeingName, "%f"); // Array with seeing vector only
+        if (!seeing || seeing->n != 1) {
+            psError(PS_ERR_UNKNOWN, false, "Unable to read file with seeing FWHMs.");
+            psFree(seeing);
+            goto ERROR;
+        }
+        psMetadataAddVector(config->arguments, PS_LIST_TAIL, "SEEING", 0, "Seeing in each image",
+                            seeing->data[0]);
+        psFree(seeing);
+    }
+
     psTrace("ppStack", 1, "Done reading command-line arguments\n");
     psFree(arguments);