Changeset 31654 for branches/eam_branches/ipp-20110505/psModules/src/objects/pmSourcePhotometry.c

Timestamp:

Jun 21, 2011, 6:18:28 PM (15 years ago)

Author:

eugene

Message:

add source to flag peaks with brighter neighbors in the footprint

File:

• branches/eam_branches/ipp-20110505/psModules/src/objects/pmSourcePhotometry.c (modified) (2 diffs)

branches/eam_branches/ipp-20110505/psModules/src/objects/pmSourcePhotometry.c

@@ -191 +191 @@
     }
 
+    pmSourceNeighborFlags (source);
+
     // measure the aperture magnitude, if (SN > AP_MIN_SN)
     if (!isfinite(SN)) {
@@ -264 +266 @@
 
 */
+
+bool pmSourceNeighborFlags (pmSource *source) {
+
+    return false;
+
+    // source must have a peak to have a footprint
+    if (!source) return false;
+    if (!source->peak) return false;
+    if (!source->peak->footprint) return false;
+    if (!source->peak->footprint->peaks) return false;
+    if (!source->peak->footprint->peaks->n) return false;
+
+    // find the brightest peak (first peak)
+    pmPeak *brightPeak = source->peak->footprint->peaks->data[0];
+
+    // are we the brightest peak?
+    if (source->peak == brightPeak) return true;
+
+    // if not, raise a flag:
+    source->mode2 |= PM_SOURCE_MODE2_HAS_BRIGHTER_NEIGHBOR;
+
+    // but, this is a common situation.  more interesting is if the ratio flux_n / (r^2 flux_o) is large
+
+    float radius2 = PS_SQR(source->peak->xf - brightPeak->xf) + PS_SQR(source->peak->yf - brightPeak->yf);
+
+    float ratio = brightPeak->rawFlux / (source->peak->rawFlux * radius2);
+
+    if (ratio > 1) {
+        source->mode2 |= PM_SOURCE_MODE2_BRIGHT_NEIGHBOR_1;
+    }
+
+    if (ratio > 10) {
+        source->mode2 |= PM_SOURCE_MODE2_BRIGHT_NEIGHBOR_10;
+    }
+
+    return true;
+}
 
 // return source model magnitude