Changeset 34492 for trunk/psphot/src/psphotFindDetections.c

Timestamp:

Oct 2, 2012, 11:18:35 AM (14 years ago)

Author:

eugene

Message:

attempt to reinsert sources on the second pass (actually on the "not-first-pass"). However, it is not an error for no sources to be present if no sources are added and/or subtracted, the significance image will not be regenerated

File:

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

trunk/psphot/src/psphotFindDetections.c

@@ -71 +71 @@
     }
 
+    // first pass vs other: if this is the first pass, the code will use PEAKS_NSIGMA_LIMIT and
+    // only attempt to detect PEAKS_NMAX entries.  If 'firstPass' is false, the code will
+    // attempt to replace the subtracted sources in order to measure the footprints.  After
+    // replacement, it is necessary to regenerate the significance image.  If no sources are
+    // available, the code will skip the significance image regeneration step.
+
     bool replaceSourcesForFootprints = false;
     if (firstPass) {
@@ -126 +132 @@
     if (useFootprints) {
         if (replaceSourcesForFootprints) {
+            bool modified = false;
             // subtract the noise for all sources including satstars
-            psphotAddOrSubNoiseReadout(config, view, filerule, index, recipe, false);
-            psphotReplaceAllSourcesReadout (config, view, filerule, index, recipe, false);
+            modified |= psphotAddOrSubNoiseReadout(config, view, filerule, index, recipe, false);
+            modified | psphotReplaceAllSourcesReadout (config, view, filerule, index, recipe, false);
 
             // add in the satstars
-            psphotAddOrSubSatstarsReadout (config, view, filerule, index, recipe, true);
+            modified |= psphotAddOrSubSatstarsReadout (config, view, filerule, index, recipe, true);
 
-            psFree (significance);
-            significance = psphotSignificanceImage (readout, recipe, maskVal);
+            if (modified) {
+                psFree (significance);
+                significance = psphotSignificanceImage (readout, recipe, maskVal);
+            }
 
             // display the significance image