Changeset 5772 for trunk/psphot/src/psphotEnsemblePSF.c

Timestamp:

Dec 13, 2005, 6:43:44 AM (21 years ago)

Author:

eugene

Message:

substantial work on ensemble PSF fitting, better functions for evaluation, overall code cleanups

File:

• trunk/psphot/src/psphotEnsemblePSF.c (modified) (3 diffs)

trunk/psphot/src/psphotEnsemblePSF.c

@@ -12 +12 @@
     // source analysis is done in spatial order
     sources = psArraySort (sources, psphotSortByY);
-
-    // this should be added to the PM_SOURCE flags:
-    int PM_SOURCE_BLEND = PM_SOURCE_OTHER + 1;
 
     float OUTER_RADIUS     = psMetadataLookupF32 (&status, config, "SKY_OUTER_RADIUS");
@@ -35 +32 @@
         // skip non-astronomical objects (very likely defects)
         // XXX EAM : should we try these anyway?
-        if (inSource->type == PM_SOURCE_BLEND) continue;
+        if (inSource->mode &  PM_SOURCE_BLEND) continue;
         if (inSource->type == PM_SOURCE_DEFECT) continue; 
         if (inSource->type == PM_SOURCE_SATURATED) continue;
@@ -73 +70 @@
         psImage *mask = otSource->mask;
 
-        // XXX EAM : set model to unit peak, zero sky, maybe use peak (x,y)
+        // XXX EAM : use these lines to fit to the peak
         // model->params->data.F32[2] = inSource->peak->x;
         // model->params->data.F32[3] = inSource->peak->y;
+        // XXX EAM : better option: improve the peak with 2D poly fit 3x3
+
+        // set model to unit peak, zero sky (we assume sky is constant)
         model->params->data.F32[0] = 0.0;
         model->params->data.F32[1] = 1.0;