Changeset 6056 for trunk/psphot/src/psphot.c

Timestamp:

Jan 19, 2006, 2:55:33 PM (21 years ago)

Author:

eugene

Message:

adapted everything to work with pmReadout *readout instead of eamReadout *imdata

File:

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

trunk/psphot/src/psphot.c

@@ -5 +5 @@
 
     psMetadata  *config  = NULL;
-    eamReadout  *imdata  = NULL;
+    psMetadata  *header  = NULL;
+    pmReadout   *readout  = NULL;
     psArray     *sources = NULL;
     psArray     *peaks   = NULL;
@@ -18 +19 @@
 
     // load input data (config and images (signal, noise, mask)
-    imdata = psphotSetup (config);
+    readout = psphotSetup (config, &header);
 
     // run a single-model test if desired
-    psphotModelTest (imdata, config);
+    psphotModelTest (readout, config);
 
     // measure image stats for initial guess 
-    sky = psphotImageStats (imdata, config);
+    sky = psphotImageStats (readout, config);
 
     // generate a background model (currently, 2D polynomial)
     // XXX this should be available to be re-added to the original image
-    psphotImageBackground (imdata, config, sky);
+    psphotImageBackground (readout, config, sky);
 
     // find the peaks in the image
-    peaks = psphotFindPeaks (imdata, config, sky);
+    peaks = psphotFindPeaks (readout, config, sky);
 
     // construct sources and measure basic stats
-    sources = psphotSourceStats (imdata, config, peaks);
+    sources = psphotSourceStats (readout, config, peaks);
 
     // classify sources based on moments, brightness
@@ -49 +50 @@
     switch (FITMODE) {
       case 0:
-        psphotEnsemblePSF (imdata, config, sources, psf, sky);
+        psphotEnsemblePSF (readout, config, sources, psf, sky);
         break;
 
       case 1:
-        psphotEnsemblePSF (imdata, config, sources, psf, sky);
-        psphotFullFit (imdata, config, sources, psf, sky);
+        psphotEnsemblePSF (readout, config, sources, psf, sky);
+        psphotFullFit (readout, config, sources, psf, sky);
         psphotReplaceUnfit (sources);
-        psphotApResid (imdata, sources, config, psf);
+        psphotApResid (readout, sources, config, psf);
         break;
 
       case 2:
-        psphotEnsemblePSF (imdata, config, sources, psf, sky);
-        psphotBlendFit (imdata, config, sources, psf, sky);
+        psphotEnsemblePSF (readout, config, sources, psf, sky);
+        psphotBlendFit (readout, config, sources, psf, sky);
         psphotReplaceUnfit (sources);
-        psphotApResid (imdata, sources, config, psf);
+        psphotApResid (readout, sources, config, psf);
         break;
 
       case 3:
-        psphotApplyPSF (imdata, config, sources, psf, sky);
+        psphotApplyPSF (readout, config, sources, psf, sky);
         break;
 
       case 4:
-        psphotApplyPSF (imdata, config, sources, psf, sky);
-        psphotFitExtended (imdata, config, sources, sky);
+        psphotApplyPSF (readout, config, sources, psf, sky);
+        psphotFitExtended (readout, config, sources, sky);
         break;
     }
 
     // write out data in appropriate format
-    psphotOutput (imdata, config, sources, psf, sky);
+    psphotOutput (readout, header, config, sources, psf, sky);
     psLogMsg ("psphot", 3, "wrote output: %f sec\n", psTimerMark ("psphot"));
     psLogMsg ("psphot", 3, "complete psphot run: %f sec\n", psTimerMark ("complete"));