Changeset 32850 for branches/eam_branches/ipp-20111122/psphot/src/psphotKronIterate.c

Timestamp:

Dec 2, 2011, 11:42:55 AM (15 years ago)

Author:

eugene

Message:

interpolate the sky-limit radius to 0.0 flux; allow faint sources to have a somewhat larger max size; ignore DEFECT and SATURATION (garbage) in masking the neighbors for kron mags

Location:

branches/eam_branches/ipp-20111122/psphot

Files:

• . (modified) (1 prop)
• src (modified) (1 prop)
• src/psphotKronIterate.c (modified) (3 diffs)

branches/eam_branches/ipp-20111122/psphot/src/psphotKronIterate.c

@@ -111 +111 @@
 
         pmSource *source = sources->data[i];
-        if (!source->peak) continue; // XXX how can we have a peak-less source?
-
-        // allocate space for moments
-        if (!source->moments) continue;
 
         // set a window function for each source based on the moments
+        // (this skips really bad sources (no peak, no moments, DEFECT)
         psphotKronWindowSetSource (source, kronWindow, false, true);
     }
@@ -221 +218 @@
 
             // XXX TEST : use a window based on the radial profile numbers: max is skyRadius, min is RADIUS
-            float maxWindow = source->skyRadius;
+            // if we lack the skyRadius (eg MATCHED sources), go to the default value
+            float maxWindow = isfinite(source->skyRadius) ? source->skyRadius : RADIUS;
             float windowRadius = PS_MAX(RADIUS, maxWindow);
 
@@ -400 +398 @@
 
     if (!source) return false;
+    if (!source->peak) return false; // XXX how can we have a peak-less source?
     if (!source->moments) return false;
+    if (source->type == PM_SOURCE_TYPE_DEFECT) return false;
+    if (source->type == PM_SOURCE_TYPE_SATURATED) return false;
     psAssert(kronWindow, "need a window");