Changeset 32205 for branches/eam_branches/ipp-20110710/psphot/src/psphotModelTestReadout.c

Timestamp:

Aug 29, 2011, 2:59:09 AM (15 years ago)

Author:

eugene

Message:

find the footprint saddles; iterative search for kron (in masked); add iterative search for kron (downweighting); add footprint based mask

File:

• branches/eam_branches/ipp-20110710/psphot/src/psphotModelTestReadout.c (modified) (2 diffs)

branches/eam_branches/ipp-20110710/psphot/src/psphotModelTestReadout.c

@@ -84 +84 @@
 
     // define the source of interest
-    float xObj     = psMetadataLookupF32 (&status, recipe, "TEST_FIT_X");
+   float xObj     = psMetadataLookupF32 (&status, recipe, "TEST_FIT_X");
     float yObj     = psMetadataLookupF32 (&status, recipe, "TEST_FIT_Y");
     if (!isfinite(xObj) || !isfinite(yObj)) psAbort ("object position is not defined");
@@ -102 +102 @@
     status = pmSourceLocalSky (source, PS_STAT_SAMPLE_MEDIAN, INNER, maskVal, markVal);
     if (!status) psAbort("pmSourceLocalSky error");
+
+    { 
+        // XXX I want to test an iterative aperture for brighter sources
+        float radius = mRADIUS;
+        for (int i = 0; i < 10; i++) {
+
+            // get the source moments
+            status = pmSourceMoments (source, radius, 0.0, 0.0, MIN_KRON_RADIUS, maskVal);
+            if (!status) psAbort("psSourceMoments error");
+
+            float oldRadius = radius;
+            radius = source->moments->Mrf * RADIUS;
+            
+            fprintf (stderr, "%d %f  %f  %f\n", i, oldRadius, radius, source->moments->Mrf);
+        }
+        exit (0);
+    }
 
     // get the source moments