Changeset 30903 for branches/eam_branches/ipp-20110213/psphot/src/psphotFindFootprints.c

Timestamp:

Mar 14, 2011, 3:07:22 PM (15 years ago)

Author:

eugene

Message:

psphotSignificanceImage now returns a pmReadout so we can keep the smoothed signal image as well as the S/N image; psphotCullPeaks (and pmFootprintCullPeaks) now uses the smooted image to cull the peaks -- this deals better with low-significance detections are large, low-surface brightness areas

File:

• branches/eam_branches/ipp-20110213/psphot/src/psphotFindFootprints.c (modified) (3 diffs)

branches/eam_branches/ipp-20110213/psphot/src/psphotFindFootprints.c

@@ -1 +1 @@
 # include "psphotInternal.h"
 
-bool psphotFindFootprints (pmDetections *detections, psImage *significance, pmReadout *readout, psMetadata *recipe, const float threshold, const int pass, psImageMaskType maskVal) {
+bool psphotFindFootprints (pmDetections *detections, pmReadout *significance, pmReadout *readout, psMetadata *recipe, const float threshold, const int pass, psImageMaskType maskVal) {
 
     bool status;
@@ -18 +18 @@
     PS_ASSERT (status, false);
 
-    // find the raw footprints & assign the peaks to those footprints
-    psArray *footprints = pmFootprintsFind (significance, threshold, npixMin);
+    // find the raw footprints in the smoothed significance image & assign the peaks to those footprints
+    psArray *footprints = pmFootprintsFind (significance->variance, threshold, npixMin);
 
     if (pmFootprintsAssignPeaks(footprints, detections->peaks) != PS_ERR_NONE) {
@@ -56 +56 @@
     }
 
-    psphotCullPeaks(readout, recipe, detections->footprints);
+    psphotCullPeaks(readout, significance, recipe, detections->footprints);
     detections->peaks = pmFootprintArrayToPeaks(detections->footprints);
     psLogMsg ("psphot", PS_LOG_INFO, "%ld peaks, %ld total footprints: %f sec\n", detections->peaks->n, detections->footprints->n, psTimerMark ("psphot.footprints"));