Changeset 14755

Timestamp:

Sep 5, 2007, 11:38:36 AM (19 years ago)

Author:

eugene

Message:

load externally-supplied PSF before finding peaks, only try for PSF Clump if PSF not supplied

File:

• trunk/psphot/src/psphotReadout.c (modified) (4 diffs)

trunk/psphot/src/psphotReadout.c

@@ -79 +79 @@
     psphotModelTest (config, view, recipe, maskVal, maskMark);
 
-    // find the peaks in the image
+    // load the psf model, if suppled
+    // XXX include code in here to save FWHM_X,FWHM_Y in the recipe
+    pmPSF *psf = psphotLoadPSF (config, view, recipe);
+    
+    // find the peaks in the image.  
+
+    // XXX clean this up into a single function.  if psf is defined, we should treat this as
+    // pass "2", not "1".  re-define this boolean to be "have PSF"
     psArray *peaks;
     psArray *footprints = NULL;
@@ -102 +109 @@
     }
 
-    // psArray *bounds = psphotGetBounds ();
-    // psArray *sources = psphotBoundsToSources ();
-
     // construct sources and measure basic stats
     psArray *sources = psphotSourceStats (readout, recipe, peaks, maskVal, maskMark);
@@ -122 +126 @@
     }
 
+    // use source classes to select PSF candidate stars (from PSF clump)?
+    bool findPSFClump = (psf == NULL);
+
     // classify sources based on moments, brightness
-    if (!psphotRoughClass (sources, recipe, true, maskSat)) {
+    if (!psphotRoughClass (sources, recipe, findPSFClump, maskSat)) {
         psLogMsg ("psphot", 3, "failed to find a valid PSF clump for image");
         return psphotReadoutCleanup (config, readout, recipe, NULL, sources);
@@ -131 +138 @@
     }
 
-    // load the psf model, if suppled
-    pmPSF *psf = psphotLoadPSF (config, view, recipe);
-    if (!psf) {
+    // if we were not supplied a PSF, choose one here
+    if (psf == NULL) {
         // use bright stellar objects to measure PSF
         psf = psphotChoosePSF (readout, sources, recipe, maskVal, maskMark);