Changeset 27848

Timestamp:

May 3, 2010, 3:03:40 PM (16 years ago)

Author:

eugene

Message:

working towards convolved stack photometry

Location:

trunk/psphot

Files:

• doc/stack.txt (modified) (2 diffs)
• src/psphotStackImageLoop.c (modified) (1 diff)
• src/psphotStackMatchPSFs.c (added)

trunk/psphot/doc/stack.txt

@@ +1 @@
+
+20100503:
+
+  psphotStack : outline
+
+  * load N images
+  * make chisq detection image
+  * detect sources on all (N+1) images
+  * convolve N input images (to match, to target)
+  * measure source properties on convolved images
 
 20100126:
@@ -4 +14 @@
   * watch out for psphotSetMomentsWindow & MOMENTS_SX_MAX,etc
   * watch out for psphotSignificanceImage.c:,psphotEfficiency using the FWHM_MAJ from psphotChoosePSF
-    * ppSimDetections.c : SIGMA_SMOOTH
-ppSmooth/src/ppSmoothReadout.c:    psMetadataAddF32(recipe, PS_LIST_TAIL, "EFFECTIVE_AREA", PS_META_REPLACE, "Effective Area", effArea);
-ppSmooth/src/ppSmoothReadout.c:    psMetadataAddF32(recipe, PS_LIST_TAIL, "SIGNIFICANCE_SCALE_FACTOR", PS_META_REPLACE, "Signicance scale factor", factor);
-
 
 20100120 : more stack processing mods:

trunk/psphot/src/psphotStackImageLoop.c

@@ -41 +41 @@
                 if (! readout->data_exists) { continue; }
 
+# if (0)                
+                // uncomment to generate matched psfs
+                if (!psphotStackMatchPSFs (config, view)) {
+                    psError(psErrorCodeLast(), false, "failure in psphotStackMatchPSFs for chip %d, cell %d, readout %d\n", view->chip, view->cell, view->readout);
+                    psFree (view);
+                    return false;
+                }
+# endif
+
                 // XXX for now, we assume there is only a single chip in the PHU:
                 if (!psphotStackReadout (config, view)) {
-                    psError(psErrorCodeLast(), false, "failure in psphotReadout for chip %d, cell %d, readout %d\n", view->chip, view->cell, view->readout);
+                    psError(psErrorCodeLast(), false, "failure in psphotStackReadout for chip %d, cell %d, readout %d\n", view->chip, view->cell, view->readout);
                     psFree (view);
                     return false;