Changeset 28782 for branches/eam_branches/ipp-20100621/psphot/src/psphotExtendedSourceFits.c

Timestamp:

Jul 29, 2010, 2:36:33 PM (16 years ago)

Author:

eugene

Message:

attempting to get a better starting guess; attempting to downweight neighbor detections

File:

• branches/eam_branches/ipp-20100621/psphot/src/psphotExtendedSourceFits.c (modified) (3 diffs)

branches/eam_branches/ipp-20100621/psphot/src/psphotExtendedSourceFits.c

@@ -39 +39 @@
     int NplainPass = 0;
     bool savePics = false;
+    float radius;
 
     // find the currently selected readout
@@ -164 +165 @@
         Next ++;
 
+        // set the radius based on the footprint (also sets the mask pixels)
+        if (!psphotSetRadiusFootprint(&radius, readout, source, markVal, 1.0)) return false;
+
+        // XXX note that this changes the source moments that are published...
+        // recalculate the source moments using the larger extended-source moments radius
+        // at this stage, skip Gaussian windowing, and do not clip pixels by S/N
+        // this uses the footprint to judge both radius and aperture?
+        // XXX save the psf-based moments for output
+        if (!pmSourceMoments (source, radius, 0.0, 0.0, maskVal)) {
+            // subtract the best fit from the object, leave local sky
+            pmSourceSub (source, PM_MODEL_OP_FULL, maskVal);
+            // XXX raise an error of some kind
+            continue;
+        }
+
         // save the modelFlux here in case we need to subtract it (for failure)
         psImage *modelFluxStart = psMemIncrRefCounter (source->modelFlux);
@@ -242 +258 @@
 
           // test for fit quality / result
+          modelFit->fitRadius = radius;
           psArrayAdd (source->modelFits, 4, modelFit);