Changeset 13166

Timestamp:

May 2, 2007, 3:00:56 PM (19 years ago)

Author:

eugene

Message:

adding optional peak limit (useful for dense fields)

File:

• trunk/psphot/src/psphotFindPeaks.c (modified) (2 diffs)

trunk/psphot/src/psphotFindPeaks.c

@@ -112 +112 @@
     }
 
+    // limit the total number of returned peak as specified
+    if (pass == 1) {
+        psArraySort (peaks, pmPeakSortBySN);
+        int NMAX = psMetadataLookupS32 (&status, recipe, "PEAKS_NMAX");
+        if (NMAX && (peaks->n > NMAX)) {
+            psArray *tmpPeaks = psArrayAllocEmpty (NMAX);
+            for (int i = 0; i < NMAX; i++) {
+                psArrayAdd (tmpPeaks, 100, peaks->data[i]);
+            }
+            psFree (peaks);
+            peaks = tmpPeaks;
+        }
+    }    
+
     // optional dump of all peak data
     char *output = psMetadataLookupStr (&status, recipe, "PEAKS_OUTPUT_FILE");
@@ -118 +132 @@
     }
     psLogMsg ("psphot", PS_LOG_INFO, "%ld peaks: %f sec\n", peaks->n, psTimerMark ("psphot"));
-    //
+
     // If they asked us to return a set of pmFootprints, not a raw pmPeak list,
     // go ahead and find the footprints, and assign them their peaks.