Changeset 30975 for branches/eam_branches/ipp-20110213/psphot/src/psphotBasicDeblend.c

Timestamp:

Mar 18, 2011, 2:10:35 PM (15 years ago)

Author:

eugene

Message:

we had 3 items conceptually equivalent to the brightness of a peak: "value", "flux", and "SN"; Ive modified pmPeak to carry explicitly named elements "rawFlux", "smoothFlux", and "detValue". Ive also added the corresponding variance values for rawFlux and smoothFlux. This lets us choose independently of the peak detection process whether flux comparisons are done relative to the smoothed or unsmoothed image. It also simplifies tests of the peak flux compared to something else

File:

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

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

@@ -63 +63 @@
     for (int i = 0; i < SN->n; i++) {
         source = sources->data[i];
-        SN->data.F32[i] = source->peak->SN;
+        SN->data.F32[i] = source->peak->rawFlux;
     }
     psVector *index = psVectorSortIndex (NULL, SN);
@@ -115 +115 @@
         // threshold is fraction of the source peak flux
         // image is background subtracted; source->moments->Sky should always be 0.0
-        threshold = FRACTION * source->peak->SN;
+        threshold = FRACTION * sqrt(source->peak->detValue);
         // threshold is no less than NSIGMA
         threshold = PS_MAX (threshold, NSIGMA);
@@ -138 +138 @@
             testSource = overlap->data[k];
             if (testSource->mode & PM_SOURCE_MODE_BLEND) continue;
-            if (testSource->peak->value > source->peak->value) continue;
+            if (testSource->peak->rawFlux > source->peak->rawFlux) continue;
             for (int j = 0; j < xv->n; j+=2) {
                 if (fabs(yv->data.F32[j] - testSource->peak->y) > 0.5) continue;