Changeset 7738 for trunk/ppImage/src/ppFocus.c

Timestamp:

Jun 28, 2006, 2:26:49 PM (20 years ago)

Author:

eugene

Message:

adding ppFocus elements

File:

• trunk/ppImage/src/ppFocus.c (modified) (2 diffs)

trunk/ppImage/src/ppFocus.c

@@ -25 +25 @@
     }
     
-    
-    
-
-    // Image Arithmetic Loop
-    if (!ppImageLoop(config, options)) {
-        psErrorStackPrint(stderr, "");
+    // we search the argument data for the named fileset (argname)
+    psArray *infiles = psMetadataLookupPtr(NULL, config->arguments, "INPUT");
+    if (!infiles) {
+        psTrace("pmFPAfile", 5, "Failed to find INPUT in argument list");
         exit(1);
     }
 
-    psLogMsg ("ppImage", 3, "complete ppImage run: %f sec\n", psTimerMark (TIMERNAME));
+    // allocate vectors for analysis
+    psVector *focus = psVectorAlloc (infiles->n, PS_TYPE_F32);
+    psVector *fwhm = psVectorAlloc (infiles->n, PS_TYPE_F32);
+
+    for (int i = 0; i < infiles->n; i++) {
+
+        // Image Arithmetic Loop
+        if (!ppImageLoop(config, options)) {
+            psErrorStackPrint(stderr, "");
+            exit(1);
+        }
+        
+        // determine FWHM at reference location in image
+        // (also removes PPIMAGE.INPUT from config->files)
+        ppFocusGetFWHM (config, focus, fwhm);
+
+        // silently ignore failure for (i == infiles->n)
+        pmFPAfileDefineSingleFromArgs (NULL, config, "PPIMAGE.INPUT", "INPUT", i+1);
+    }
+
+    ppFocusFitFWHM (config, focus, fwhm);
+
+    psLogMsg ("ppFocus", 3, "complete ppFocus run: %f sec\n", psTimerMark (TIMERNAME));
 
     // Cleaning up
@@ -44 +64 @@
 // - the input list is a set of independent images (not multiple files for a single image)
 // - each pass to ppImageLoop performs the analysis on a different pmFPAfile
-// - 
+// - after each ppImageLoop, grap the input pmFPAfile and extract the FWHM stats