Changeset 12933

Timestamp:

Apr 19, 2007, 4:42:44 PM (19 years ago)

Author:

eugene

Message:

adding residuals

Location:

branches/eam_02_branch/psphot/src

Files:

• Makefile.am (modified) (1 diff)
• psphot.h (modified) (1 diff)
• psphotChoosePSF.c (modified) (2 diffs)
• psphotMakeResiduals.c (added)
• psphotRadialPlot.c (modified) (2 diffs)

branches/eam_02_branch/psphot/src/Makefile.am

@@ -56 +56 @@
         psphotDeblendSatstars.c  \
         psphotMosaicSubimage.c   \
+        psphotMakeResiduals.c    \
         psphotAddNoise.c
 

branches/eam_02_branch/psphot/src/psphot.h

@@ -103 +103 @@
 bool            psphotDeblendSatstars (psArray *sources, psMetadata *recipe);
 
+bool            psphotMakeResiduals (psArray *sources, psMetadata *recipe, pmPSF *psf);
+
 #endif

branches/eam_02_branch/psphot/src/psphotChoosePSF.c

@@ -155 +155 @@
     try = models->data[bestN];
 
+    // XXX calculate the PSF residual table here (and add to the pmPSF structure)
+
     // measure and save the median value of dparams[PM_PAR_SXX],dparams[PM_PAR_SYY]
     // these are used by psphotEvalPSF to identify extended sources
@@ -221 +223 @@
     }
 
+    // build a PSF residual image
+    psphotMakeResiduals (try->sources, recipe, try->psf);
+
     // XXX test dump of psf star data and psf-subtracted image
     if (psTraceGetLevel("psphot.psfstars") > 5) { 

branches/eam_02_branch/psphot/src/psphotRadialPlot.c

@@ -8 +8 @@
 static int nCount = 0;
 
-bool psphotRadialPlot (int *kapa, char *filename, pmSource *source) {
+bool psphotRadialPlot (int *kapa, const char *filename, pmSource *source) {
 
     Graphdata graphdata;
@@ -95 +95 @@
     sprintf (pagename, "%02d", nCount);
     if (nCount == 0) {
-        KiiPS (*kapa, false, KAPA_PS_NEWPLOT, filename, pagename);
+        KiiPS (*kapa, filename, false, KAPA_PS_NEWPLOT, pagename);
     } else {
-        KiiPS (*kapa, false, KAPA_PS_NEWPAGE, filename, pagename);
+        KiiPS (*kapa, filename, false, KAPA_PS_NEWPAGE, pagename);
     }