Changeset 30931 for branches/eam_branches/ipp-20110213/psModules/src/objects/pmSourceOutputs.c

Timestamp:

Mar 16, 2011, 1:32:43 PM (15 years ago)

Author:

eugene

Message:

create a single test function (pmSourcePositionUseMoments) to check if we should use the moments or peaks for source positions

File:

• branches/eam_branches/ipp-20110213/psModules/src/objects/pmSourceOutputs.c (modified) (3 diffs)

branches/eam_branches/ipp-20110213/psModules/src/objects/pmSourceOutputs.c

@@ -77 +77 @@
 }
 
+// what is the correct postion?
+// * if we have a PSF model fit, use PM_PAR_XPOS,YPOS for X_PSF,Y_PSF 
+// * if we do not have a model:
+// ** if we have moments:
+// *** if the star is saturated, use the moments
+// *** if the moments and peak agree to < DR, use the moments
+// *** otherwise, use the peak
+
 bool pmSourceOutputsSetValues (pmSourceOutputs *outputs, pmSource *source, pmChip *chip, float fwhmMajor, float fwhmMinor, float magOffset) {
 
@@ -112 +120 @@
         outputs->apRadius = source->apRadius;
     } else {
-        bool useMoments = true;
-        useMoments = (useMoments && source->moments);          // can't if there are no moments
-        useMoments = (useMoments && source->moments->nPixels); // can't if the moments were not measured
-        useMoments = (useMoments && !(source->mode && PM_SOURCE_MODE_MOMENTS_FAILURE)); // can't if the moments failed...
+        bool useMoments = pmSourcePositionUseMoments(source);
 
         if (useMoments) {
@@ -175 +180 @@
     return true;
 }
+
+// this function decides if the source position should be based on the peak or the moments.
+// this is only used if we know we should not use a model fit position (eg, no model, or no
+// model yet)
+bool pmSourcePositionUseMoments(pmSource *source) {
+
+    if (!source->moments) return false;          // can't if there are no moments
+    if (!source->moments->nPixels) return false; // can't if the moments were not measured
+    if (source->mode && PM_SOURCE_MODE_MOMENTS_FAILURE) return false; // can't if the moments failed...
+
+    if (source->mode & PM_SOURCE_MODE_SATSTAR) return true; // moments are best for SATSTARs
+
+    float dX = source->moments->Mx - source->peak->xf;
+    float dY = source->moments->My - source->peak->yf;
+    float dR = hypot(dX, dY);
+    
+    // only use the moments position if the moment-peak offset is small or the star is saturated
+    if (dR > 1.5) return false
+
+    return true;
+}